Предсказанная среднеквадратическая ошибка для адаптивного фильтра LMS
[mmseemse] = msepred(lmsFilt,x,d)[mmseemse,meanw,mse,tracek] = msepred(lmsFilt,x,d)[mmseemse,meanw,mse,tracek] = msepred(lmsFilt,x,d,m)[ предсказывает установившиеся значения в сходимости минимальной среднеквадратической ошибки, mmseemse] = msepred(lmsFilt,x,d)mmse, и избыточной среднеквадратической ошибки, emse, учитывая вход и желаемые последовательности сигнала ответа в x и d и количествах в Системе dsp.LMSFilter object™, lmsFilt.
[ задает дополнительный фактор десятикратного уменьшения для вычисления mmseemse,meanw,mse,tracek] = msepred(lmsFilt,x,d,m)meanw, mse и tracek. Если m> 1, каждый m th ожидаемое значение каждой из этих последовательностей сохранен. Если не использовано, значение m является значением по умолчанию, которое является тем.
Среднеквадратическая ошибка (MSE) измеряет среднее значение квадратов ошибок между желаемым сигналом и первичным входом сигнала к адаптивному фильтру. Сокращение этой ошибки сходится первичный входной параметр к желаемому сигналу. Определите ожидаемое значение MSE и моделируемое значение MSE в каждый раз мгновенное использование функции msesim и msepred. Сравните эти значения MSE друг с другом и относительно минимального MSE и установившихся значений MSE. Кроме того, вычислите сумму квадратов содействующих ошибок, данных трассировкой содействующей ковариационной матрицы.
Примечание: Если вы используете R2016a или более ранний релиз, заменяете каждый вызов объекта с эквивалентным синтаксисом шага. Например, obj(x) становится step(obj,x).
Инициализация
Создайте Систему dsp.FIRFilter object™, который представляет неизвестную систему. Передайте сигнал, x, к КИХ-фильтру. Вывод неизвестной системы является желаемым сигналом, d, который является суммой вывода неизвестной системы (КИХ-фильтр) и аддитивный шумовой сигнал, n.
num = fir1(31,0.5); fir = dsp.FIRFilter('Numerator',num); iir = dsp.IIRFilter('Numerator',sqrt(0.75),... 'Denominator',[1 -0.5]); x = iir(sign(randn(2000,25))); n = 0.1*randn(size(x)); d = fir(x) + n;
Фильтр LMS
Создайте Системный объект dsp.LMSFilter, чтобы создать фильтр, который адаптируется, чтобы вывести желаемый сигнал. Установите длину адаптивного фильтра к 32 касаниям, размера шага к 0,008, и фактор десятикратного уменьшения для анализа и симуляции к 5. Переменная simmse представляет моделируемый MSE между выводом неизвестной системы, d, и выводом адаптивного фильтра. Переменная mse дает соответствующее ожидаемое значение.
l = 32; mu = 0.008; m = 5; lms = dsp.LMSFilter('Length',l,'StepSize',mu); [mmse,emse,meanW,mse,traceK] = msepred(lms,x,d,m); [simmse,meanWsim,Wsim,traceKsim] = msesim(lms,x,d,m);
Постройте результаты MSE
Сравните значения моделируемого MSE, предсказал MSE, минимальный MSE и итоговый MSE. Итоговое значение MSE дано суммой минимального MSE и избыточного MSE.
nn = m:m:size(x,1); semilogy(nn,simmse,[0 size(x,1)],[(emse+mmse)... (emse+mmse)],nn,mse,[0 size(x,1)],[mmse mmse]) title('Mean Squared Error Performance') axis([0 size(x,1) 0.001 10]) legend('MSE (Sim.)','Final MSE','MSE','Min. MSE') xlabel('Time Index') ylabel('Squared Error Value')
Предсказанный MSE следует за той же траекторией как моделируемый MSE. Обе этих траектории сходятся с установившимся (итоговым) MSE.
Постройте содействующие траектории
meanWsim является средним значением моделируемых коэффициентов, данных msesim. meanW является средним значением предсказанных коэффициентов, данных msepred.
Сравните моделируемые и предсказанные средние значения коэффициентов фильтра LMS 12,13,14, и 15.
plot(nn,meanWsim(:,12),'b',nn,meanW(:,12),'r',nn,... meanWsim(:,13:15),'b',nn,meanW(:,13:15),'r') PlotTitle ={'Average Coefficient Trajectories for';... 'W(12), W(13), W(14), and W(15)'}
PlotTitle = 2x1 cell array
{'Average Coefficient Trajectories for'}
{'W(12), W(13), W(14), and W(15)' }
title(PlotTitle) legend('Simulation','Theory') xlabel('Time Index') ylabel('Coefficient Value')
В устойчивом состоянии оба сходятся траектории.
Сумма содействующих ошибок в квадрате
Сравните сумму содействующих ошибок в квадрате, данных msepred и msesim. Эти значения даны трассировкой содействующей ковариационной матрицы.
semilogy(nn,traceKsim,nn,traceK,'r') title('Sum-of-Squared Coefficient Errors') axis([0 size(x,1) 0.0001 1]) legend('Simulation','Theory') xlabel('Time Index') ylabel('Squared Error Value')
Функция msepred вычисляет каждый элемент последовательности среднеквадратических ошибок с помощью следующих уравнений:
где,
mmse – Минимальная среднеквадратическая ошибка.
dk – Диагональные элементы содействующей ковариационной матрицы. Начальным значением dk дают
.
meanw – Средние значения вектора коэффициентов даны lmsFilt.InitialConditions.
Wopt – Оптимальные Винеровские коэффициенты фильтра.
Q L-by-L матрица, столбцы которой являются собственными векторами входной матрицы автокорреляции, R, так, чтобы R Q = Q Lam. Lam является диагональной матрицей, содержащей соответствующие собственные значения.
diagA – Фрагмент аналитической матрицы перехода MSE.
dk ◦diagA – Адамар или entrywise продукт dk и diagA.
λ – Вектор-столбец, содержащий собственные значения входной матрицы автокорреляции, R.
μ – Размер шага дан lmsFilt.StepSize.
kurt – Среднее значение эксцесса собственного вектора отфильтровало сигналы.
λ ◦2 – Адамар или entrywise степень вектор-столбца, содержащего собственные значения.
D Управление термином для анализа MSE.
[1] Hayes, M.H. Статистическая цифровая обработка сигналов и моделирование. Нью-Йорк: John Wiley & Sons, 1996.