%LSeg1 
u=[-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1]; %系統(tǒng)辨識的輸入信號為一個周期的M序列(模擬白噪聲)
z=zeros(1,16); %定義輸出觀測值的長度
for k=3:16 
    z(k)=1.5*z(k-1)-0.7*z(k-2)+u(k-1)+0.5*u(k-2); %用理想輸出值作為觀測值
end
subplot(3,1,1) %畫三行一列圖形窗口中的第一個圖形
stem(u) %畫出輸入信號u的經(jīng)線圖形
subplot(3,1,2) %畫三行一列圖形窗口中的第二個圖形
i=1:1:16; %橫坐標(biāo)范圍是1到16，步長為1
plot(i,z) %圖形的橫坐標(biāo)是采樣時刻i, 縱坐標(biāo)是輸出觀測值z, 圖形格式為連續(xù)曲線
subplot(3,1,3) %畫三行一列圖形窗口中的第三個圖形
stem(z),grid on%畫出輸出觀測值z的經(jīng)線圖形，并顯示坐標(biāo)網(wǎng)格
u,z%顯示輸入信號和輸出觀測信號
%L=14%數(shù)據(jù)長度
HL=[-z(2) -z(1) u(2) u(1);-z(3) -z(2) u(3) u(2);-z(4) -z(3) u(4) u(3);-z(5) -z(4) u(5) u(4);-z(6) -z(5) u(6) u(5);-z(7) -z(6) u(7) u(6);-z(8) -z(7) u(8) u(7);-z(9) -z(8) u(9) u(8);-z(10) -z(9) u(10) u(9);-z(11) -z(10) u(11) u(10);-z(12) -z(11) u(12) u(11);-z(13) -z(12) u(13) u(12);-z(14) -z(13) u(14) u(13);-z(15) -z(14) u(15) u(14)] %給樣本矩陣HL賦值
ZL=[z(3);z(4);z(5);z(6);z(7);z(8);z(9);z(10);z(11);z(12);z(13);z(14);z(15);z(16)]% 給樣本矩陣zL賦值
%calculating parameters%計算參數(shù)
c1=HL'*HL; c2=inv(c1); c3=HL'*ZL; c=c2*c3 %計算并顯示 
%DISPLAY PARAMETERS
a1=c(1), a2=c(2), b1=c(3), b2=c(4) %從 中分離出并顯示a1 、a2、 b1、 b2
%End