%遞推最小二乘參數估計（RLS）
clear?all;?close?all;

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

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

theta=[a（2:na+1）;b];?%對象參數真值

thetae_1=zeros（na+nb+11）;?%thetae初值
P=10^6*eye（na+nb+1）;?
for?k=1:L
????phi=[-yk;uk（d:d+nb）];?%此處phi為列向量
????y（k）=phi‘*theta+xi（k）;?%采集輸出數據
???
????