程序：u=[-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1]; %系统辨识的输入信号为一个周期的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
>> su
bp
lot(3 ,1,1) %画三行一列图形窗口中的第一个图形%
>> stem(u) %画输入信号u的径线图形%
>> title('stemplot of input series')
>> subplot(3,1,2) %画三行一列图形窗口中的第二个图形%
>> i=1:1:16; %横坐标范围是1到16，步长为1%
>> plot(i,z) %图形的横坐标是采样时刻i，纵坐标是输出观测值z，图形格式为连续曲线%
>> title('wave of measurements of output')
>> subplot(3,1,3) %画三行一列图形窗口中的第三个图形%
>> stem(z),grid on %画出输出观测值z的径线图形，并显示坐标网格%
>> title('stemplot of measurements of output')
>> u,z %显示输入信号和输出观测信号%
%L=14 数据长度%
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赋值 %
c1=HL'*HL;c2=inv(c1);c3=HL'*ZL;c=c2*c3 %计算并显示参数%
a1=c(1),a2=c(2),b1=c(3),b2=c(4) %从参数中分离出并显示a1,a2,b1,b2%