matlab仿真ALE算法

%% this program verifies the ALE technology
% x1 is a narrow band signal while x2 is a wide band signal

%%generate data module
N=2000;
delta=0.01;
sum=0;
for i=1:100
en=sqrt(1)*randn(1,N);
for i=3:N
    x1(i)=sin(0.05*pi*i+2*pi*rand(1));
    x2(i)=en(i)+en(i-1)+en(i-2);
end
x=x1+x2;%从x(3)开始
X=x(:,3:N);%X的原始数据
RXX=xcorr(X);
Y=X(:,10:N-2);%取得X的后2，时延取10

%%data process module
 H1=zeros(2,N);
 H1(:,1)=[3 -4]';
 for i=1:N-20
     e(i+1)=Y(i+1)-H1(:,i)'*[X(i+1); X(i)];
     H1(:,i+1)=H1(:,i)+delta*e(i+1)*[X(i+1); X(i)];
     Yest(i)=H1(:,i)'*[X(i+1); X(i)];
 end
 sum=sum+Yest;
end
Yest=sum/100;
%%display module as follows

 figure(1);
 plot(x1)
 axis([0  2000 -2 2]);
 xlabel(