Используйте Систему фильтра RLS object™, чтобы определить значение сигналов.

Примечание: Если вы используете R2016a или более ранний релиз, заменяете каждый вызов объекта с эквивалентным step синтаксис. Например, obj(x) становится step(obj,x).

rls1 = dsp.RLSFilter(11, 'ForgettingFactor', 0.98);
filt = dsp.FIRFilter('Numerator',fir1(10, .25)); % Unknown System
x = randn(1000,1);                       % input signal
d = filt(x) + 0.01*randn(1000,1); % desired signal
[y,e] = rls1(x, d);
w = rls1.Coefficients;
subplot(2,1,1), plot(1:1000, [d,y,e]);
title('System Identification of an FIR filter');
legend('Desired', 'Output', 'Error');
xlabel('time index'); ylabel('signal value');
subplot(2,1,2); stem([filt.Numerator; w].');
legend('Actual','Estimated');
xlabel('coefficient #'); ylabel('coefficient value');

Подавление помех

rls2 = dsp.RLSFilter('Length', 11, 'Method', 'Householder RLS');
filt2 = dsp.FIRFilter('Numerator',fir1(10, [.5, .75]));
x = randn(1000,1);                           % Noise
d = filt2(x) + sin(0:.05:49.95)';     % Noise + Signal
[y, err] = rls2(x, d);
subplot(2,1,1), plot(d), title('Noise + Signal');
subplot(2,1,2), plot(err), title('Signal');

Продемонстрируйте адаптивный алгоритм RLS с помощью обратной модели системы идентификации, показанной здесь.

Расположение каскадом адаптивного фильтра с неизвестным фильтром заставляет адаптивный фильтр сходиться к решению, которое является инверсией неизвестной системы.

Если передаточная функция неизвестной системы и адаптивного фильтра является H (z) и G (z), соответственно, ошибка, измеренная между желаемым сигналом и сигналом от каскадной системы, достигает своего минимума, когда продукт H (z) и G (z) равняется 1, G (z) *H (z) = 1. Для этого отношения, чтобы быть верным, G (z) должен равняться 1/H (z), инверсия передаточной функции неизвестной системы.

Чтобы продемонстрировать, что это верно, создайте сигнал ввести к каскадной паре фильтра, s.

s = randn(3000,1);

В каскадном случае фильтров неизвестный фильтр приводит к задержке сигнала, прибывающего в точку суммирования после обоих фильтров. Чтобы препятствовать тому, чтобы адаптивный фильтр пытался адаптироваться к сигналу, это еще не видело (эквивалентный предсказанию будущего), задержите желаемый сигнал 12 выборками, порядком неизвестной системы.

Обычно вы не знаете порядка системы, которую вы пытаетесь идентифицировать. В этом случае задержите желаемый сигнал количеством выборок, равных половине порядка адаптивного фильтра. Задержка входа требует предварительного ожидания 12 выборок нулевых значений к входу, s.

delay = zeros(12,1);
d = [delay; s(1:2988)]; % Concatenate the delay and the signal.

Необходимо сохранить желаемый сигнальный вектор d та же длина как x, следовательно настройте количество элемента сигнала, чтобы допускать выборки задержки.

Несмотря на то, что не обычно верный, для этого примера вы знаете порядок неизвестного фильтра, таким образом, вы добавляете задержку, равную порядку неизвестного фильтра.

Для неизвестной системы используйте lowpass, КИХ-фильтр 12-го порядка.

filt = dsp.FIRFilter;
filt.Numerator = fir1(12,0.55,'low');

Фильтрация s обеспечивает сигнал входных данных для адаптивной функции алгоритма.

x = filt(s);

Чтобы установить алгоритм RLS, инстанцируйте dsp.RLSFilter объект и набор его Length, ForgettingFactor, и InitialInverseCovariance свойства.

Для получения дополнительной информации о входных условиях подготовить объект алгоритма RLS, обратитесь к dsp.RLSFilter.

p0 = 2 * eye(13);
lambda = 0.99;
rls = dsp.RLSFilter(13,'ForgettingFactor',lambda,...
   'InitialInverseCovariance',p0);

Большая часть процесса к этой точке совпадает с предыдущими примерами. Однако, поскольку этот пример стремится разработать обратное решение, необходимо быть осторожными, о котором сигнал несет данные и который является желаемым сигналом.

Более ранние примеры адаптивных фильтров используют отфильтрованный шум в качестве желаемого сигнала. В этом случае, отфильтрованный шум (x) несет неизвестную информацию о системе. С Распределением Гаусса и отклонением 1, неотфильтрованный шумовой d желаемый сигнал. Код, чтобы запустить этот адаптивный пример фильтра

[y,e] = rls(x,d);

где, y возвращает отфильтрованный выходной параметр и e содержит сигнал ошибки, когда фильтр адаптируется, чтобы найти инверсию неизвестной системы.

Предполагаемые коэффициенты фильтра RLS получены путем ввода rls.Coefficients в командной строке MATLAB.

b = rls.Coefficients;

Просмотрите частотную характеристику адаптированного фильтра RLS (обратная система, G (z)) использование freqz. Обратная система похожа на фильтр highpass с линейной фазой.

freqz(b,1)

Просмотрите частотную характеристику неизвестной системы, H (z). Ответ является ответом фильтра нижних частот с частотой среза 0,55.

freqz(filt.Numerator,1)

Результатом каскада неизвестной системы и адаптированного фильтра является компенсированная система с расширенной частотой среза 0,8.

overallCoeffs = conv(filt.Numerator,b);
freqz(overallCoeffs,1)