%遺忘因子遞推最小二乘參數估計（FFRLS）
clear?all;?close?all;

a=[1?-1.5?0.7]‘;?b=[1?0.5]‘;?d=3;?%對象參數
na=length（a）-1;?nb=length（b）-1;?%na、nb為A、B階次

L=1000;?%仿真長度
uk=zeros（d+nb1）;?%輸入初值：uk（i）表示u（k-i）
yk=zeros（na1）;?%輸出初值
u=randn（L1）;?%輸入采用白噪聲序列
xi=sqrt（0.1）*randn（L1）;?%白噪聲序列

thetae_1=zeros（na+nb+11）;?%thetae初值
P=10^6*eye（na+nb+1）;
lambda=0.98;?%遺忘因子范圍[0.9?1]
for?k=1:L
????if?k==501
????????a=[1?-1?0.4]‘;b=[1.5?0.2]‘;?%對象參數突變
????end
????theta（:k）=[a（2:na+1）;b];?%對象參數真值
????
????phi=[-yk;uk（d:d+nb）];
????y（k）=phi‘*theta（:k）+xi（k）;?%采集輸出數據
???
????%遺忘因子遞推最小二乘法
????K=P*phi/（lambda+phi‘*P*phi）;
??