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+