aaexample741.m

自适应信号处理中LMS

function[ems]=aaexample741(mu,M)
%function[ems]=aaexample741(mu,M);
%this function produces figures like in example 7.4.1;
eav=zeros(1,1000);
dn(1)=0;dn(2)=0;x(1)=0;x(2)=0;
for k=1:120
    for n=3:1000
        dn(n)=.96*dn(n-1)-.2*dn(n-2)+2*(rand-.5);
        x(n)=dn(n-1);
    end;
    [w,y,e]=aalms(x,dn,mu,M);
    eav=eav+e.^2;
end;
ems=eav/120;