%隨機性問題梯度校正參數辨識
%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（