自适应滤波器matlab代码

%误差性能曲面和等值曲线图
clear all;
h0=[-5:0.2:5];
h1=[-5:0.2:5];
[mn]=meshgrid（h0h1）;
J=0.55+m.^2+n.^2+2*m.*n.*cos（pi/8）-2^（1/2）*m.*cos（pi/10）-2^（1/2）*n.*cos（9*pi/40）;
figure（1）;
surf（mnJ）;
xlabel（‘h0‘）;ylabel（‘h1‘）;zlabel（‘J‘）;title（‘误差性能曲面‘）
figure（2）;
contour（mnJ50）;hold on
xlabel（‘h0‘）;ylabel（‘h1‘）;title（‘等值曲线‘）


%最陡下降法和LMS算法，HH（:n）是最陡下降法的H（n），HH1（:n）是LMS的H（n）。
b=0.4;N=222;M=100;e=zeros（N+11）;Hs=zeros（2N+1）;Js=zeros（N+11）;J= zeros（N+11）;
for k=0:1
    ryx（k+1）=2^（-1/2）*cos（pi*k/8+pi/10）;
    Rxx（k+1）=cos（pi*k/8）;
end
r_yx=[ryx（1）;ryx（2）];
R_xx=[Rxx（1）Rxx（2）;Rxx（2）Rxx（1）];
for n=0:N
    if n==0
        H=[3;-4];
    else
        H1=H+b*（r_yx-R_xx*H）;
        H=H1;
    end
    HH（:n+1）=H;
end
    
for i=1:M
    %产生白噪声
    bzs=wgn（1N+1