Matlab系统辨识仿真CH3,CH4,CH6程序注释与剖析VIP免费

第3 章程序及注释 例3.3 考虑仿真对象 )()2(5.0)1()2(7.0)1(5.1)(kvkukukzkzkz (3.41) 其中，)(kv是服从正态分布的白噪声N)1,0(。输入信号采用4 阶M 序列，幅度为1。选择如下形式的辨识模型 )()2()1()2()1()(2121kvkubkubkzakzakz (3.42) 设输入信号的取值是从k =1 到k =16 的M 序列，则待辨识参数LSθˆ为LSθˆ =LτL1LτLzH)HH(。其中，被辨识参数LSθˆ 、观测矩阵 z L、H L 的表达式为 2121ˆbbaaLSθ , )16()4()3(zzzLz， )14()2()1()15()3()2()14()2()1()15()3()2(uuuuuuzzzzzzLH (3.43) 程序框图如图 3.2 所示。Matlab6.0 仿真程序如下： 赋输入信号初值u 定义输出观测值的长度并计算系统的输出值 画出输入和输出观测值的图形 给样本矩阵 HL 和 zL 赋值 根据式(3.27)计算参数LSθˆ 从LSθˆ中分离出并显示出被辨识参数a1, a2, b1, b2 停机 图 3.2 最小二乘一次完成算法程序框图 %二阶系统的最小二乘一次完成算法辨识程序，在光盘中的文件名：FLch3LSeg1.m 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 subplot(3,1,1) %画三行一列图形窗口中的第一个图形 stem(u) %画输入信号 u 的径线图形 subplot(3,1,2) %画三行一列图形窗口中的第二个图形 i=1:1:16; %横坐标范围是 1 到 16，步长为 1 plot(i,z) %图形的横坐标是采样时刻 i, 纵坐标是输出观测值 z, 图形格式为连续曲线 subplot(3,1,3) %画三行一列图形窗口中的第三个图形 stem(z),grid on %画出输出观测值 z的径线图形，并显示坐标网格 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(1...