%誤差性能曲面和等值曲線圖
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
????%產(chǎn)生白噪聲
????bzs=wgn（1N+1