最小二乘法Matlab代码

%递推最小二乘参数估计（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）; %采集输出数据