Matlab系统辨识应用例子.docxVIP

例1、考虑仿真对象 z(k) 1.5z(k 1) 0.7z(k 2) u(k 1) 0.5u(k 2) v(k) 其中，v(k)是服从正态分布的白噪声N (0,1)。输入信号采用4阶M序 列，幅度为1。选择如下形式的辨识模型 z(k) a1z(k 1) a2z(k 2) bu(k 1) bu(k 2) v(k) 设输入信号的取值是从k =1到k =16的M序列，则待辨识参 数？ls 数？ls为L = ( T ) 1 Tz。其中，被辨识参数 和、观测矩阵 的表达 式为: a1 a2 b1 b2z(3)z(4) a1 a2 b1 b2 z(3) z(4) III z(16) z(2) z(1) i(2) u(1) z(3) z(2) u(3) u(2) III III z(15) z(14)u(15) u(14) 程序框图如图 1所示。Matlab仿真程序如下: 图1最小二乘一次完成算法程序框图 %二阶系统的最小二乘一次完成算法辨识程序，文件名： LS.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(14) u(15) u(14)] % 给样本矩阵 赋值 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 c1=HL*HL; c2=inv(c1); c3=HL*ZL; c=c2*c3 % 计算并显示 L %Display Parameters a1=c(1), a2=c(2), b1=c(3),b2=c(4) % 从和中分离出并显示 ai、 a2 、 b1、 b2 %End 程序运行结果： u =[ -1 , 1 , -1 , 1 ,1,1 , 1 , -1 , -1 , -1 , 1 , -1 , -1 , 1 , 1] z =[ 0 , 0, 0.5000 , 0.2500 , 0.5250 , 2.1125 , 4.3012 , 6.4731 , 6.1988 , 3.2670 , -0.9386 , -3.1949 , -4.6352 , 6.2165 , -5.5800 , -2.5185] HL = 0 0 1.0000 -1.0000 1 -0.5000 0 -1.0000 1.0000 -0.2500 -0.5000 1.0000 -1.0000 -0.5250 -0.2500 1.0000 1.0000 -2.1125 -0.5250 1.0000 1.0000 -4.3012 -2.1125 1.0000 1.0000 -6.4731 -4.3012 -1.0000 1.0000 -6.1988 -6.4731 -1.0000 -1.0000 -3.2670 -6.1988 -1.0000 -1.0000 0.9386 -3.2670 1.0000 -