clear;
clf;
w=[2.9];
M=11;
var=0.001;?%噪聲方差
td=7;??%參考信號延時
N=400;?%?訓(xùn)練次數(shù)
ii=500;%迭代次數(shù)
L=ii+td+M-1;?%單個輸入信號長度
u=0.038;?%迭代步長
value=zeros（length（w）L-M+1-td）;

for?ww=1:length（w）
????h=ones（13）;
????h（1）=0.5*（1+cos（2*pi/w（ww）））;
????h（3）=h（1）;
????e2=zeros（NL-M+1-td）;
????for?n=1:N
????????rand（‘seed‘n*N）;??
????????X=sign（2*rand（1L）-1）;?;%產(chǎn)生隨機信號序列
????????D=zeros（size（X））;
????????for?mm=td+1:L?
????????????D（mm）=X（mm-td）;
????????end?
????????U=conv（Xh）;
????????randn（‘seed‘n*N）;
????????V=randn（size（U））.*sqrt（var）;?%產(chǎn)生高斯噪聲
????????R=U+V;?%自適應(yīng)濾波器輸入信號
????????W=zeros（M1）;?%濾波器參數(shù)的初始值為0
????????for?m=1:ii
????????????r=R（td+m:td+m+M-1）‘;
????????????y=r‘*W;
????????????e=D（td+m+M-1）-y;?%?誤差信號
????????????e2（nm）=e.^2;?
????????????W=W+2*u*e*r;?%濾波器參數(shù)迭代
???????
??????
????????end
????end
????value（ww:）=mean（e2）;?%均方誤差
????
end

semilogy（value（1:）‘g‘）;
title（‘CMA算法（u=0.038）‘）;
xlabel（‘迭代次數(shù)‘）;
ylabel（‘均方誤差‘）;