%% 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(