Вот, нашёл Формулу Адаптивного Шумоподавителя - [Ссылки могут видеть только зарегистрированные пользователи. ]

Может кто-нибудь подкинуть Целый КОД этого Фильтра из МатЛаба или МатКада?
Там ещё есть Фильтр Виннера - то же здорово было бы найти готовый КОД из МатЛаба/МатКада...


Фильтр состоит из блока цифрового фильтра с регулируемыми коэффициентами и адаптивного алгоритма для настройки и изменения коэффициентов фильтра.
На фильтр одновременно подаются входные сигналы y(k) и x(k). Сигнал y(k) содержит полезный сигнал s(k) и некоррелированный с ним загрязняющий сигнал g(k). Сигнал x(k) какого-либо источника шума, коррелированный с g(k), который используется для формирования оценки сигнала ^g(k)

. Полезный сигнал оценивается по разности:
^s(k)=y(k)−^g(k)=s(k)+g(k)−^g(k)(1. 1)

Возводим уравнение в квадрат и получаем:
^s2(k)=s2(k)+(g(k)−^g(k))2+2s(k)(g( k)−^g(k))(1.2)

Вычислим математическое ожидание левой и правой части этого уравнения:
M[^s2(k)]=M[s2(k)]+M[g(k)−^g(k)2]+2M[s(k)(g(k)−^g(k))](1.3)

Последнее слагаемое в выражении равно нулю, поскольку сигнал s(k)
не коррелирует с сигналами g(k) и ^g(k)

M[^s2(k)]=M[s2(k)]+M[(g(k)−^g(k)2)]+2M[s(k)(g(k)−^g(k))](1.4)

В этом выражении M[s2(k)]=W(s(k))
– мощность сигнала s(k), M[^s2(k)]=W(^s(k)) – оценка мощности сигнала s(k) и общая выходная мощность, M[g(k)−^g(k)2]=W(ξg)

- остаточная мощность шума, который может содержаться в выходном сигнале. При настройке адаптивного фильтра к оптимальному положению минимизируется мощность остаточного шума, а, следовательно, и мощность выходного сигнала:
minW(^s(k))=W(s(k))+minW(ξg)(1.5)

На мощность полезного сигнала настройка не влияет, поскольку сигнал не коррелирован с шумом. Эффект минимизации общей выходной мощности будет выражаться в максимизации выходного отношения сигнал/шум. Если настройка фильтра обеспечивает равенство ^g(k)=g(k)
, то при этом ^s(k)=s(k). Если сигнал не содержит шума, адаптивный алгоритм должен устанавливать нулевые значения всем коэффициентам цифрового фильтра.