%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）‘;
?