%隨機(jī)性問(wèn)題梯度校正參數(shù)辨識(shí)
%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
%?產(chǎn)生4位m序列作為輸入數(shù)據(jù)
L=800;?%?L輸入數(shù)據(jù)長(zhǎng)度
y1=1;y2=0;y3=1;y4=0;?%4個(gè)移位寄存器的輸出初始值
for?i=1:L;
????x1=xor（y3y4）;?%第一個(gè)移位寄存器的輸入信號(hào)
????x2=y1;?%第二個(gè)移位寄存器的輸入信號(hào)
????x3=y2;?%第三個(gè)移位寄存器的輸入信號(hào)
????x4=y3;?%第四個(gè)移位寄存器的輸入信號(hào)
????y（i）=y4;?%第四個(gè)移位寄存器的輸出信號(hào)，M序列?幅值“0“和“1“，?
????if?y（i）>0.5u（i）=-1;?%M序列的值為“1“時(shí)辨識(shí)的輸入信號(hào)取“-1”
????else?u（i）=1;?%M序列的值為“0“時(shí)辨識(shí)的輸入信號(hào)取“1”
????end
????y1=x1;y2=x2;y3=x3;y4=x4;?%為下一次的輸入信號(hào)準(zhǔn)備
end
s=normrnd（00.31800）;?%生成輸入噪聲，均值為0，方差為0.1，500個(gè)

%?產(chǎn)生輸出數(shù)據(jù)
z（2）=0;?z（1）=0;?%系統(tǒng)最初幾個(gè)輸出的初值?
for?k=3:L;?%開(kāi)始求K?
????z（k）=1.75*z（k-1）-0.8*z（k-2）+1.4*u（k-1）+0.2*u（