RLS算法及遗忘因子对RLS的影响

%rls algorithm
clear all;
x=randn（5001）;%高斯白噪声
for i=3:500
    d（i）=x（i）-1.6*x（i-1）+0.8*x（i-2）;%输出期望
end
w1=zeros（5013）;%估计参数
w2=zeros（5013）;
w3=zeros（5013）;
%参数初始化
R1=0;R2=0;R3=0;
t1=0;t2=0;t3=0;
%lmta=0.9;
%算法更新
for i=3:500
    u=[x（i）x（i-1）x（i-2）]‘;
    k1=t1*u/（0.9+u‘*t1*u）;
    e1=d（i）-w1（i:）*u;
    w1（i+1:）=w1（i:）+（k1*e1）‘;
    R1=0.9*R1+u‘*u;
    t1=1/R1;
    %t=inv（R）;
    u=[x（i）x（i-1）x（i-2）]‘;
    k2=t2*u/（0.8+u‘*t2*u）;
    e2=d（i）-w2（i:）*u;
    w2（i+1:）=w2（i:）+（k2*e2）‘;