资源描述：
《最小二乘参数辨识的matlab仿真程序注释与剖析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、最小二乘一次完成算法的MATLAB仿真例2・1考虑仿真对象z伙)—1.5?伙-1)+o.7z伙—2)=1心-1)+0.5比伙-2)+v(k)(2-1)其中，W灯是服从正态分布的白噪声N(0,l)o输入信号采用4阶M序列，幅度为1。选择如下形式的辨识模型z(k)+仗-1)+dpZ伙一2)=b{u(k-1)+b2u(k-2)+v(k)(2-2)(2-3)LS其中,设输入信号的取值是从k=1到A=16的M序列，则待辨识参数(2-4)被辨识参数％、观测矩阵"、丹匚的表达式为■込⑶■Zl=2(4)•••hl=$刃6)_9一z(l»一z(l令u(5)w(14)-z(2)-z(l)心)i
2、/(l)_z(3)_z(2)u(3)u(2)例2-1程序框图如图2.1所示:图2.1最小二乘一次完成算法程序框图例2-lMatlab仿真程序如下：%二阶系统的最小二乘一次完成算法辨识程序，文件名：FLch3LSegl.m11=卜1,1,・1,1,1,-1,1,-1,1,1];%系统辨识的输入信号为一个周期的M序列z=zeros(l,16);%定义输出观测值的长度fork=3:16z(k)=1.5*z(k-l)-0.7*z(k-2)+u(k-l)+0.5*u(k-2);%用理想输出值作为观测值subplot(3,l,l)%画三行一列图形窗口中的第一个图形stem(u)%画输入信
3、号u的径线图形subplot(3,l,2)%画三行一列图形窗口中的第二个图形i=l:l:16;%横坐标范是1到16,步长为1plot(i,z)%图形的横坐标是采样时刻纵坐标是输出观测值Z,图形格式为连续曲线subplot⑶1,3)%画三行一列图形窗口中的第三个图形stem(z),gridon%画出输出观测值z的径线图形，并显示坐标网格u,z%显示输入信号和输出观测信号%L=14%数据长度HL=[-z(2)-z(l)u(2)u(l);-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(
4、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(ll)-z(10)u(ll)u(10);-z(12)-z(ll)u(12)u(ll);-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(ll);z(12);z(13);z(14);z(15);z(16)]54给样本
5、矩阵"赋值%CalculatingParameterscl=HLf*HL;c2=inv(cl);c3=HL*ZL;c=c2*c3%计算并显示6拓%DisplayParametersal=c(l),a2=c(2),bl=c(3),b2=c(4)%从血中分离出并显示偽、a2>方1、b2%End例2-1程序运行结果:U=[-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-I,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,・
6、5.5800,-2.5185]00-0.50000-0.2500-0.5000-0.5250-0.250021125-0.5250-4.3012-2.1125-6.4731-4.3012-6.1988•6.4731-3.2670-6.19880.9386-3.26703.19490.93864.63523.19496.21654.63525.58006.21651.0000-1.0000-1.00001.00001.0000-1.00001.00001.00001.00001.00001.00001.0000-1.00001.0000-1.0000-1.0000-1.0000
7、-1.00001.0000-1.0000-1.00001.0000-1.0000•1.00001.0000-1.00001.00001.0000丿14*4)ZL=[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]7c=[-1.500(),0.7000,1.0000,().50()0]ral=-1.5000a2=0.7000bl=1.0000b2=
