实验二线性连续定常系统的运动分析.docx

《实验二线性连续定常系统的运动分析.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验二线性连续定常系统的运动分析一、实验目的1．掌握线性连续定常系统的状态转移矩阵的求法，学会用MATLAB求解状态转移矩阵。2．掌握线性连续定常系统的状态方程的求解方法，学会用MATLAB求解线性连续定常系统的时间响应，并绘制相应的状态响应曲线和输出响应曲线。二、实验原理1．线性连续定常系统状态转移矩阵的计算设线性连续定常系统的状态空间表达式为xAxBu，则其状态转移矩yCxDu阵为Φ(t)eAt从时间角度看，状态转移矩阵使状态向量随着时间的推移不断地作坐标变换，不断地在状态空间中作转移，故称为状态转移矩阵。对于线性连续定常系统，其状态转移矩阵与其矩阵指数函数相同，可利用直接求解法、

2、拉氏变换法、标准型法和待定系数法等方法对其进行求解。（1）直接求解法Φ(t)eAt1AktkIAt1A2t2L1AktkLk0k!2!k!（2）拉氏变换法Φ(t)eAtL1sIA1（3）标准型法对系统矩阵A进行线性非奇异变换，将其变换为对角线矩阵或约旦矩阵AP1AP，其中P为非奇异变换阵。状态转移矩阵为Φ(t)eAtPeAtP1，其中AP1AP若A的特征值1,2,L,n两两互异，则A为对角线矩阵，此时e1t0Φ(t)eAtPeAtP1POP10ent若A有n重特征值i，则A为约旦矩阵，此时eitteitLΦAtAtQ1Q0OO(t)eQeMOO0L0（4）待定系数法1tn1eit(n

3、1)!MQ1ttei根据凯莱－哈密顿（以下简称C-H）定理，线性连续定常系统的状态转移矩阵为n1Φ(t)eAtaj(t)Aja0(t)Ia1(t)ALan1(t)An1j0其中，a0(t),a1(t),L,an1(t)为t的标量函数，可按A的特征值确定。若A的特征值1,2,L,n两两互异，则a0(t)12Ln11t111e1a1(t)12Ln1e2t222MMMMOMMan1(t)12Ln1entnnn若A有n重特征值i，则a0(t)00L01111)!tn1e1t(na1(t)00L1(n1)1(n12)!tn2e1tMMMNMMMan2(t)0121Lan1(t)12L11n1n

4、21!11!1te1tn1e1t12．线性连续定常系统的时间响应设线性连续定常系统的状态空间表达式为xAxBu，满足初始条件x(t)t0x(0)yCxDu则其状态方程的解（状态响应）为x(t)Φ(t)x(0)t)Bu()dΦ(t（直接求解法）0x(t)L1sI1x(0)sIA1（拉氏变换法）ABU(s)其输出方程的解（输出响应）为y(t)CΦ(t)x(0)t)Bu()dDu(t)或CΦ(t0y(t)L1CsIA1C1BU(s)Du(t)x(0)sIA3．MATLAB相关函数命令简介（1）矩阵指数函数expm()phi=expm(A*t)：返回常矩阵A的矩阵指数函数eAt，其中phi的维

5、数和矩阵A的维数相同。（2）拉氏逆变换ilaplace()g=ilaplace(G)：返回函数G(s)的拉氏逆变换g(t)，即ftL1Gs（3）任意输入下系统响应lsim()[y,x]=lsim(sys,u,t,x0)：直接求取线性系统在任意输入信号下的响应，其中sys为线性系统模型，u为与时间向量t对应的输入向量，y表示输出响应y(t)，x表示状态响应x(t)。（4）单位阶跃响应step()[y,t,x]=step(sys)：返回系统在零状态条件下的单位阶跃响应，其中y表示单位阶跃输出响应y(t)，x表示单位阶跃状态响应y(t)；step(sys)：此时不返回任何变量，而自动地绘制单

6、位阶跃响应输出曲线。相关响应函数：脉冲响应函数impulse()，零输入响应函数initial()（5）二维图像绘制plot()plot(x,y)：绘制y=y(x)函数的二维图形，其中x和y具有相同维数的列向量或行向量。三、实验内容1．状态转移矩阵的计算（1）已知某线性连续定常系统的系统矩阵为A0165a．根据直接求解法（expm函数）和拉氏变换法（ilaplace函数）分别求解该系统的状态转移矩阵，并比较这2种方法获得的状态转移矩阵是否相同；b．计算t0.3s时的状态转移矩阵的值。2．线性连续定常系统的状态响应和输出响应（1）已知某线性连续定常系统的状态空间表达式为&012ux23

7、x1，utett00，其中x0y12x0a．利用ilaplace()函数求解该系统的状态响应x(t)和输出响应y(t)的解析解；b．利用lsim()函数计算0t10s该系统的状态响应x(t)和输出响应y(t)的数值解；c．根据状态响应x(t)和输出响应y(t)的数值解结果，利用plot()函数绘制其状态响应曲线和输出响应曲线。&010ux2x1（2）已知某线性连续定常系统的状态方程为3，其中y12x初始状态x02,利用step()函数求解单位阶