梯度校正参数辨识程序

%随机性问题梯度校正参数辨识
%z（k）=1.75*z（k-1）-0.8*z（k-2）+1.4*u（k-1）+0.2*u（k-2）+s（k）;
%======================================== 
clc
clear
% 产生4位m序列作为输入数据
L=800; % L输入数据长度
y1=1;y2=0;y3=1;y4=0; %4个移位寄存器的输出初始值
for i=1:L;
    x1=xor（y3y4）; %第一个移位寄存器的输入信号
    x2=y1; %第二个移位寄存器的输入信号
    x3=y2; %第三个移位寄存器的输入信号
    x4=y3; %第四个移位寄存器的输入信号
    y（i）=y4; %第四个移位寄存器的输出信号，M序列 幅值“0“和“1“， 
    if y（i）>0.5u（i）=-1; %M序列的值为“1“时辨识的输入信号取“-1”
    else u（i）=1; %M序列的值为“0“时辨识的输入信号取“1”
    end
    y1=x1;y2=x2;y3=x3;y4=x4; %为下一次的输入信号准备
end
s=normrnd（00.31800）; %生成输入噪声，均值为0，方差为0.1，500个

% 产生输出数据
z（2）=0; z（1）=0; %系统最初几个输出的初值 
for k=3:L; %开始求K 
    z（k）=1.75*z（k-1）-0.8*z（k-2）+1.4*u（k-1）+0.2*u（