多级最小二乘法Matlab程序

clear
clc
%==========================================
%最小二乘法辨识
% Z（k+3）=-0.9*Z（k+2）-0.15*Z（k+1）-0.02*z（k）+0.7*u（k+2）-1.5*u（k+1）+e（k）
%e（k+2）+1.0*e（k+1）+0.41*e（k）=r*v（k+2）

%==========产生M序列作为输入===============
x=[0 1 0 1 1 0 1 1 1];  %initial value
n=405; %n为脉冲数目
M=[];  %存放M序列

for i=1:n
    temp=xor（x（4）x（9））;
    M（i）=x（9）;
    for j=9:-1:2
        x（j）=x（j-1）;
    end
    x（1）=temp;
end


%===========产生均值为0，方差为1的高斯白噪声=============
v=randn（1405）;
e=[];
e（1）=0.3;
e（2）=0.7;
r=0.9;    %控制信噪比
for i=3:405
    e（i）=-1.0*e（i-1）-0.41*e（i-2）+r*v（i）;
end

%=================产生观测序列===================
z=[];
z（1）=-1;
z（2）=0;
z（3）=1.5;

for i=4:405
	z（i）=-0.9*z（i-1）-0.15*z（i-2）-0.02*z（i-3）+0.7*M（i-1）-1.5*M（i-2）+e（i）;
end

%================第一级辨识  辅助模型参数辨识==================

H=zeros（4009）;
for i=1:400
    H（i1）=-z（i+4）;
    H（i2）=-z（i+3）;
    H（i3）=-z（i+2）;
    H（i4）=-z（i+