clear;
ts=cputime;
clc;
N=2000;
% close all;
sum=zeros(1,N); 
Ci=1000;
h = waitbar(0,'Please wait...');
for i=1:Ci 
 waitbar(i/Ci,h) 
N=2000;
N1=1000;
Nd1=8;
Nd2=5;
 Noise1=randn(1,N);
       Noise3=randn(1,N);
        x2=randn(1,N1);
         x=Noise1+0.01*Noise3;
d=[zeros(1,Nd2) 0.8*x(1:end-Nd2)]+[zeros(1,Nd1) 0.4*x(1:end-Nd1)]+[0.01*x2(1:1000) zeros(1,1000)];%主回路信號(hào)   
  lenth=10;  %濾波器階次
%   w=[0 0 0 0 0 0.7 0 0 0.3 0];
       w=zeros(1,lenth);     %系數(shù)初始化
        u=0.01;
%        u=1/(x*x'/N*lenth)*0.1;             %步長(zhǎng)
       y=zeros(1,N);
       e=ones(1,N); 
%開(kāi)始算法；
 for n=lenth:N;
        xx=x(n:-1:n-lenth+1);
        y(n)=xx*w';
        e(n)=d(n)-y(n);
        w=w+2*u*xx*e(n);
        ee(n)=e(n).^2; 
  end
   sum=sum+10*log10(ee);
end
 close(h)
  h1=figure(1);
 plot(sum/Ci,'r')
  title('學(xué)習(xí)曲線(xiàn)')
%  title('濾波器階數(shù)對(duì)濾波效果影響')
 xlabel('迭代次數(shù)Ci');
 ylabel('誤差/db');
 axis([0,1000,-50,0]);
  hold on