%rls
randn('seed',0);
rand('seed',0);
NoOfDate=8000;
Order=32;
Lambda=0.98;
Delta=0.001;
x=randn(NoOfDate,1);
h=rand(Order,1);
d=filter(h,1,x);
P=Delta*eye(Order,Order);
w=zeros(Order,1);
for n=Order:NoOfDate;
    u=x(n:-1:n-Order+1);
    pi_=u'*P;
    k=Lambda+pi_*u;
    K=pi_'/k;
    e(n)=d(n)-w'*u;
    PPrime=K*pi_;
    P=(P-PPrime)/Lambda;
    w_err(n)=norm(h-w);
end;
figure;
plot(20*log10(abs(e)));
title('學(xué)習(xí)曲線');
xlabel('迭代次數(shù)');
ylabel('輸出誤差估計(jì)');
figure;
semilogy(w_err);
title('矢量誤差估計(jì)');
xlabel('迭代次數(shù)');
ylabel('誤差權(quán)矢量');