系统辨识 递推最小二乘法

L=3000;



z（2）=0;z（1）=0;%取z的前二个初始值为零
for k=3:L+4; 
    u（k）=randn（1）;
    z（k）=1.5*z（k-1）-0.7*z（k-2）+1.0*u（k-1）+1.5*u（k-2）+ normrnd（0111）;%给出理想的辨识输出采样信号
    zn（k1）=z（k）;
end

c0=[0 0 0 0 ]‘; %a1 a2 b1 b2给出被辨识参数的初始值

p0=10^6*eye（44）;

c=[c0zeros（42999）];%被辨识参数矩阵的初始值及大小


for k=3:3000; %开始求K 
    h1=[-z（k-1）-z（k-2）u（k-1）u（k-2）]‘; x=h1‘*p0*h1+1; x1=inv（x）; %开始求K（k）
    k1=p0*h1*x1;%求出K的值
    d1=z（k）-h1‘