2.平均迭代轨迹曲线.m

自适应滤波器LMS算法的误差曲线和误差平方曲线

[h0,h1]=meshgrid([-4:0.2:4]);   %确定滤波器系数H(n)的范围
J=0.55+h0.^2+h1.^2+2*h0.*h1*cos(pi/8)-sqrt(2)*h0*cos(pi/10)-sqrt(2)*h1*cos(pi/40*9);
                                %计算各组系数下的误差性能J(H)
axis([-4,4,-4,4]);
hold on;
contour(h0,h1,J,100);           %绘制误差性能曲面的二维等高线
xlabel('h0');
ylabel('h1');

A=[linspace(0,0,102); linspace(0,0,102)];
A=A';
for i=0:99                      %100次LMS迭代实验
    S=sqrt(0.05)*randn(1000,1);
    H=[3; -4];
    T=H';
    e=0;
    for n=0:1:100
        y=sin(pi/8*(n+1)+pi/10)+S(n+1,:);
        e=y-H'*[sqrt(2)*sin(pi/8*(n+1)); sqrt(2)*sin(pi/8*n)];
        H=H+0.4*e*[sqrt(2)*sin(pi/8*(n+1)); sqrt(2)*sin(pi/8*n)];
        T=[T; H'];
    end
    A=A+T/100;                  %所有H(n)轨迹取平均值
end
A=A';
plot(A(1,:),A(2,:),'r');        %绘制H(n)平均值的轨迹曲线