уравняться

(Подлежит удалению) Уравнять сигнал с помощью объекта эквалайзера

уравнивание будет удалено в следующем выпуске. Использовать comm.LinearEqualizer или comm.DecisionFeedback вместо этого.

Синтаксис

y = equalize(eqobj,x) y = equalize(eqobj,x,trainsig) [y,yd] = equalize(...) [y,yd,e] = equalize(...)

Описание

y = equalize(eqobj,x) обрабатывает вектор сигнала основной полосы частот x с объектом эквалайзера eqobj и возвращает выровненный вектор сигнала y. В конце процесса, eqobj содержит обновленную информацию о состоянии, такую как значения веса эквалайзера и значения входного буфера. Построить eqobj, используйте lineareq или dfe функция. equalize функция предполагает, что сигнал x отбирается при nsamp выборки на символ, где nsamp - значение nSampPerSym имущество eqobj. Для адаптивных алгоритмов, отличных от CMA, эквалайзер адаптируется в режиме принятия решения с использованием детектора, заданного SigConst имущество eqobj. Задержка эквалайзера (eqobj.RefTap-1)/eqobj.nSampPerSym.

Обратите внимание, что (eqobj.RefTap-1) должно быть целым кратным nSampPerSym. Для компенсатора с дробным интервалом отводы разнесены на доли периода символа. Опорный отвод относится к обучающим символам и, таким образом, должен совпадать с целым числом символов (т.е. целым числом выборок на символ). eqobj.RefTap=1 соответствует первому символу, eqobj.RefTap=nSampPerSym+1 ко второму и так далее. Поэтому (eqobj.RefTap-1) должно быть целым кратным nSampPerSym.

Если eqobj.ResetBeforeFiltering является 0, equalize использует существующую информацию о состоянии в eqobj при запуске операции выравнивания. В результате, equalize(eqobj,[x1 x2]) эквивалентно [equalize(eqobj,x1) equalize(eqobj,x2)]. Перезагружать eqobj вручную, примените reset функция для eqobj.

Если eqobj.ResetBeforeFiltering является 1, equalize сброс eqobj перед началом операции выравнивания перезаписать любую предыдущую информацию о состоянии в eqobj.

y = equalize(eqobj,x,trainsig) первоначально использует обучающую последовательность для адаптации эквалайзера. После обработки обучающей последовательности корректор адаптируется в режиме принятия решения. Длина вектора trainsig должно быть меньше или равно length(x)-(eqobj.RefTap-1)/eqobj.nSampPerSym.

[y,yd] = equalize(...) возвращает вектор yd обнаруженных символов данных.

[y,yd,e] = equalize(...) возвращает результат вычисления ошибки. Для адаптивных алгоритмов, отличных от CMA, e - вектор ошибок между y и опорный сигнал, где опорный сигнал состоит из обучающей последовательности или обнаруженных символов.

Примеры

свернуть все

Настройка линейных выравнивателей

В этом примере настраиваются рекомендуемые comm.LinearEqualizer Системные object™ и устаревшие lineareq с сопоставимыми настройками.

Инициализация переменных и поддерживающих объектов

d = randi([0 3],1000,1);
x = pskmod(d,4,pi/4);
r = awgn(x,25);
sps = 2; %samples per symbol for oversampled cases
nTaps = 6;
txFilter = comm.RaisedCosineTransmitFilter('FilterSpanInSymbols',nTaps, ...
    'OutputSamplesPerSymbol',4);
rxFilter = comm.RaisedCosineReceiveFilter('FilterSpanInSymbols',nTaps, ...
    'InputSamplesPerSymbol',4,'DecimationFactor',2);
x2 = txFilter(x);
r2 = rxFilter(awgn(x2,25,0.5));
filterDelay = txFilter.FilterSpanInSymbols/2 + ...
    rxFilter.FilterSpanInSymbols/2; % Total filter delay in symbols

Чтобы сравнить выровненный выходной сигнал, постройте график совокупностей с использованием кода, например:

% plot(yNew,'*')
% hold on
% plot(yOld,'o')
% hold off; legend('New Eq','Old Eq'); grid on

Использовать алгоритм LMS с линейным эквалайзером

Настроить lineareq и comm.LinearEqualizer объекты с сопоставимыми параметрами. LeakageFactor удалено из алгоритма LMS. comm.LinearEqualizer Система object™ предполагает, что коэффициент утечки всегда равен 1.

eqOld = lineareq(5,lms(0.05),pskmod(0:3,4,pi/4))

eqOld =
  EqType: 'Linear Equalizer'
  AlgType: 'LMS'
  nWeights: 5
  nSampPerSym: 1
  RefTap: 1
  SigConst: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  StepSize: 0.0500
  LeakageFactor: 1
  Weights: [0 0 0 0 0]
  WeightInputs: [0 0 0 0 0]
  ResetBeforeFiltering: 1
  NumSamplesProcessed: 0

eqNew = comm.LinearEqualizer('NumTaps',5,'Algorithm','LMS','StepSize',0.05, ...
    'Constellation',pskmod(0:3,4,pi/4),'ReferenceTap',1)

eqNew = comm.LinearEqualizer with properties:
  Algorithm: 'LMS'
  NumTaps: 5
  StepSize: 0.0500
  Constellation: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  ReferenceTap: 1
  InputDelay: 0
  InputSamplesPerSymbol: 1
  TrainingFlagInputPort: false
  AdaptAfterTraining: true
  InitialWeightsSource: 'Auto'
  WeightUpdatePeriod: 1

Вызовите эквалайзеры.

yOld = equalize(eqOld,r,x(1:100));
yNew = eqNew(r,x(1:100));

Использование алгоритма RLS с линейным эквалайзером

Настроить lineareq и comm.LinearEqualizer объекты с сопоставимыми параметрами.

eqOld = lineareq(5,rls(0.95),pskmod(0:3,4,pi/4))

eqOld =
  EqType: 'Linear Equalizer'
  AlgType: 'RLS'
  nWeights: 5
  nSampPerSym: 1
  RefTap: 1
  SigConst: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  InvCorrInit: 0.1000
  InvCorrMatrix: [5×5 double]
  Weights: [0 0 0 0 0]
  WeightInputs: [0 0 0 0 0]
  ResetBeforeFiltering: 1
  NumSamplesProcessed: 0

eqNew = comm.LinearEqualizer('NumTaps',5,'Algorithm','RLS', ...
    'ForgettingFactor',0.95,'Constellation',pskmod(0:3,4,pi/4),'ReferenceTap',1)

eqNew = comm.LinearEqualizer with properties:
  Algorithm: 'RLS'
  NumTaps: 5
  ForgettingFactor: 0.9500
  InitialInverseCorrelationMatrix: 0.1000
  Constellation: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  ReferenceTap: 1
  InputDelay: 0
  InputSamplesPerSymbol: 1
  TrainingFlagInputPort: false
  AdaptAfterTraining: true
  InitialWeightsSource: 'Auto'
  WeightUpdatePeriod: 1

Вызовите эквалайзеры. Когда ResetBeforeFiltering имеет значение true, каждый вызов equalize объект сбрасывает эквалайзер. Получение вызова эквивалентного поведения reset после каждого вызова comm.LinearEqualizer объект.

yOld1 = equalize(eqOld,r,x(1:100));
yOld2 = equalize(eqOld,r,x(1:100));

yNew1 = eqNew(r,x(1:100));
reset(eqNew)
yNew2 = eqNew(r,x(1:100));

Настроить lineareq и comm.LinearEqualizer объекты с сопоставимыми параметрами. Для comm.LinearEqualizer объект, установить начальную матрицу обратной корреляции в eye(5)*0.2.

eqOld = lineareq(5,rls(0.95),pskmod(0:3,4,pi/4))

eqOld =
  EqType: 'Linear Equalizer'
  AlgType: 'RLS'
  nWeights: 5
  nSampPerSym: 1
  RefTap: 1
  SigConst: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  InvCorrInit: 0.1000
  InvCorrMatrix: [5×5 double]
  Weights: [0 0 0 0 0]
  WeightInputs: [0 0 0 0 0]
  ResetBeforeFiltering: 1
  NumSamplesProcessed: 0

eqNew = comm.LinearEqualizer('NumTaps',5,'Algorithm','RLS', ...
    'ForgettingFactor',0.95,'Constellation',pskmod(0:3,4,pi/4),'ReferenceTap',1, ...
    'InitialInverseCorrelationMatrix',eye(5)*0.2)

eqNew = comm.LinearEqualizer with properties:
  Algorithm: 'RLS'
  NumTaps: 5
  ForgettingFactor: 0.9500
  InitialInverseCorrelationMatrix: [5×5 double]
  Constellation: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  ReferenceTap: 1
  InputDelay: 0
  InputSamplesPerSymbol: 1
  TrainingFlagInputPort: false
  AdaptAfterTraining: true
  InitialWeightsSource: 'Auto'
  WeightUpdatePeriod: 1

yOld1 = equalize(eqOld,r,x(1:100));
yOld2 = equalize(eqOld,r,x(1:100));

yNew1 = eqNew(r,x(1:100));
reset(eqNew)
yNew2 = eqNew(r,x(1:100));

Использование алгоритма CMA с линейным эквалайзером

eqOld = lineareq(5,cma(0.05),pskmod(0:3,4,pi/4))

eqOld =
  EqType: 'Linear Equalizer'
  AlgType: 'Constant Modulus'
  nWeights: 5
  nSampPerSym: 1
  SigConst: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  StepSize: 0.0500
  LeakageFactor: 1
  Weights: [1 0 0 0 0]
  WeightInputs: [0 0 0 0 0]
  ResetBeforeFiltering: 1
  NumSamplesProcessed: 0

eqNew = comm.LinearEqualizer('NumTaps',5,'Algorithm','CMA','StepSize',0.05, ...
    'Constellation',pskmod(0:3,4,pi/4),'ReferenceTap',1)

eqNew = comm.LinearEqualizer with properties:
  Algorithm: 'CMA'
  NumTaps: 5
  StepSize: 0.0500
  Constellation: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  ReferenceTap: 1
  InputSamplesPerSymbol: 1
  AdaptWeightsSource: 'Property'
  AdaptWeights: true
  InitialWeightsSource: 'Auto'
  WeightUpdatePeriod: 1

Вызовите эквалайзеры.

yOld = equalize(eqOld,r);
yNew = eqNew(r);

Использование линейных выравнивателей с учетом задержек сигнала

Настроить lineareq и comm.LinearEqualizer объекты с сопоставимыми параметрами. Фильтры передачи и приема приводят к задержке сигнала между сигналами передачи и приема. Учтите эту задержку, установив RefTap имущества lineareq к значению, близкому к значению задержки в выборках. Дополнительно, nWeights должно быть установлено значение больше, чем RefTap.

eqOld = lineareq(filterDelay*sps+4,lms(0.01),pskmod(0:3,4,pi/4),sps);
eqOld.RefTap = filterDelay*sps+1 % Adjust to synchronize with delayed signal

eqOld =
  EqType: 'Linear Equalizer'
  AlgType: 'LMS'
  nWeights: 16
  nSampPerSym: 2
  RefTap: 13
  SigConst: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  StepSize: 0.0100
  LeakageFactor: 1
  Weights: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
  WeightInputs: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
  ResetBeforeFiltering: 1
  NumSamplesProcessed: 0

eqNew = comm.LinearEqualizer('NumTaps',16,'Algorithm','LMS','StepSize',0.01, ...
    'Constellation',pskmod(0:3,4,pi/4),'InputSamplesPerSymbol',sps, ...
    'ReferenceTap',filterDelay*sps+1,'InputDelay',0)

eqNew = comm.LinearEqualizer with properties:
  Algorithm: 'LMS'
  NumTaps: 16
  StepSize: 0.0100
  Constellation: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  ReferenceTap: 13
  InputDelay: 0
  InputSamplesPerSymbol: 2
  TrainingFlagInputPort: false
  AdaptAfterTraining: true
  InitialWeightsSource: 'Auto'
  WeightUpdatePeriod: 1

yOld1 = equalize(eqOld,r,x(1:100));
yOld2 = equalize(eqOld,r,x(1:100));

yNew1 = eqNew(r,x(1:100));
reset(eqNew)
yNew2 = eqNew(r,x(1:100));

В comm.LinearEqualizer объект, InputDelay используется для синхронизации с задержанным сигналом. NumTaps и ReferenceTap не зависят от значения задержки. Мы можем уменьшить количество отводов, используя InputDelay для синхронизации вместо ссылки. Уменьшение числа отводов также уменьшает собственный шум эквалайзера.

eqNew = comm.LinearEqualizer('NumTaps',11,'Algorithm','LMS','StepSize',0.01, ...
    'Constellation',pskmod(0:3,4,pi/4),'InputSamplesPerSymbol',sps, ...
    'ReferenceTap',6,'InputDelay',filterDelay*sps)

eqNew = comm.LinearEqualizer with properties:
  Algorithm: 'LMS'
  NumTaps: 11
  StepSize: 0.0100
  Constellation: [0.7071 + 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071i 0.7071 - 0.7071i]
  ReferenceTap: 6
  InputDelay: 12
  InputSamplesPerSymbol: 2
  TrainingFlagInputPort: false
  AdaptAfterTraining: true
  InitialWeightsSource: 'Auto'
  WeightUpdatePeriod: 1

yNew1 = eqNew(r2,x(1:100));
reset(eqNew)
yNew2 = eqNew(r2,x(1:100));

Вопросы совместимости

развернуть все

`equalize` будет удален

Предупреждает, начиная с R2020a

уравнивание будет удалено в следующем выпуске. Использовать comm.LinearEqualizer или comm.DecisionFeedback вместо этого. Для примеров, сравнивающих установку comm.LinearEqualizer кому lineareq, см. раздел Настройка линейных выравнивателей.

Документация

уравняться

Синтаксис

Описание

Примеры

Настройка линейных выравнивателей

Вопросы совместимости

`equalize` будет удален

См. также

Объекты

Темы

Документация по инструментам связи

Поддержка

Документация

уравняться

Синтаксис

Описание

Примеры

Настройка линейных выравнивателей

Вопросы совместимости

equalize будет удален

См. также

Объекты

Темы

Документация по инструментам связи

Поддержка

`equalize` будет удален