nlms.m

来自「基于变步长LMS在干扰抵消中的应用」· M 代码 · 共 14 行

function [y,w]=nlms(s,order,d)
w=zeros(order,1);
s=s(:);
N=size(s,1);
alpha=0.95;
beita=0.001;
M=5;
for n=order:N
    u=s(n:-1:n-order+1);
    y(n)=w'*u;
    e(n)=d(n)-y(n);
    mu=alpha*2/(M*(beita+u'*u));
    w=w+mu*u*e(n);
end