文档介绍：该【线性多变量系统的运动分析 】是由【毒药 Posion】上传分享，文档一共【11】页，该文档可以免费在线阅读，需要了解更多关于【线性多变量系统的运动分析 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。《现代控制理论基础》第二章(讲义)2II、分析部分第二章线性多变量系统的运动分析在讨论了状态方程的描述、标准形和模型转换后,本章将讨论线性多变量系统的运动分析,即线性状态方程的求解。对于线性定常系统,为保证状态方程解的存在性和唯一性,系统矩阵A和输入矩阵B中各元必须有界。一般来说,在实际工程中,这个条件是一定满足的。 Σ: ()其中,,且初始条件为。将方程()写为在上式两边左乘e-At,可得将上式由O积分到t,得故可求出其解为 ()或 ()式中为系统的状态转移矩阵。对于线性时变系统非齐次状态方程, ()类似可求出其解为 ()一般说来,线性时变系统的状态转移矩阵只能表示成一个无穷项之和,只有在特殊情况下,才能写成矩阵指数函数的形式。《现代控制理论基础》第二章(讲义) ()的解。下面不加证明地给出线性时变系统状态转移矩阵的几个重要性质:1、;2、;3、;4、当A给定后,唯一;5、计算时变系统状态转移矩阵的公式 ()上式一般不能写成封闭形式,可按精度要求,用数值计算的方法取有限项近似。特别地,只有当满足即在矩阵乘法可交换的条件下,才可表示为如下矩阵指数函数形式 ()显然,定常系统的状态转移矩阵不依赖于初始时刻,其性质仅是上述时变系统的特例。------------------------------------------------------------------------------[]试求如下线性定常系统的状态转移矩阵Ф(t)和状态转移矩阵的逆Ф-1(t)。[解]对于该系统,《现代控制理论基础》第二章(讲义)3其状态转移矩阵由下式确定由于其逆矩阵为因此=由于Ф-1(t)=Ф(-t),故可求得状态转移矩阵的逆为------------------------------------------------------------------------------[]求下列系统的时间响应:式中,u(t)为t=0时作用于系统的单位阶跃函数,即u(t)=1(t)。[解],即因此,系统对单位阶跃输入的响应为:《现代控制理论基础》第二章(讲义)4或如果初始状态为零,即X(0)=0,可将X(t)简化为------------------------------------------------------------------------------ ,即着重讨论Caley-Hamilton定理和最小多项式。-哈密尔顿(Caley-Hamilton)定理 在证明有关矩阵方程的定理或解决有关矩阵方程的问题时,凯莱-哈密尔顿定理是非常有用的。 考虑n×n维矩阵A及其特征方程凯莱-哈密尔顿定理指出,矩阵A满足其自身的特征方程,即 ()为了证明此定理,注意到(λI-A)的伴随矩阵adj(λI-A)是λ的n-1次多项式,即式中,。由于可得从上式可看出,A和(i=1,2,…,n)相乘的次序是可交换的。因此,如果(λI-A)及其伴随矩阵adj(λI-A)中有一个为零,则其乘积为零。如果在上式中用A代替λ,显然λI-A为零。这样《现代控制理论基础》第二章(讲义)5即证明了凯莱-哈密尔顿定理。-哈密尔顿定理,任一n×n维矩阵A满足其自身的特征方程,然而特征方程不一定是A满足的最小阶次的纯量方程。我们将矩阵A为其根的最小阶次多项式称为最小多项式,也就是说,定义n×n维矩阵A的最小多项式为最小阶次的多项式φ(λ),即使得φ(A)=0,或者最小多项式在n×n维矩阵多项式的计算中起着重要作用。假设λ的多项式d(λ)是(λI-A)的伴随矩阵adj(λI-A)的所有元素的最高公约式。可以证明,如果将d(λ)的λ最高阶次的系数选为1,则最小多项式φ(λ)由下式给出: ()注意,n×n维矩阵A的最小多项式φ(λ)可按下列步骤求出:1、根据伴随矩阵adj(λI-A),写出作为λ的因式分解多项式的adj(λI-A)的各元素;2、确定作为伴随矩阵adj(λI-A)各元素的最高公约式d(λ)。选取d(λ)的λ最高阶次系数为1。如果不存在公约式,则d(λ)=1;3、最小多项式φ(λ)可由|λI-A|除以d(λ)得到。 前已指出,状态方程的解实质上可归结为计算状态转移矩阵,即矩阵指数函数eAt。如果给定矩阵A中所有元素的值,MATLAB将提供一种计算eAT的简便方法,其中T为常数。除了上述方法外,对eAt的计算还有几种分析方法可供使用。这里我们将介绍其中的四种计算方法。:直接计算法(矩阵指数函数) ()可以证明,对所有常数矩阵A和有限的t值来说,这个无穷级数都是收敛的。:对角线标准形与Jordan标准形法 若可将矩阵A变换为对角线标准形,那么eAt可由下式给出《现代控制理论基础》第二章(讲义)6式中,P是将A对角线化的非奇异线性变换矩阵。类似地,若矩阵A可变换为Jordan标准形,则eAt可由下式确定出eAt=SeJtS–1 ()------------------------------------------------------------------------------[]考虑如下矩阵A[解]该矩阵的特征方程为因此,矩阵A有三个相重特征值λ=1。可以证明,矩阵A也将具有三重特征向量(即有两个广义特征向量)。易知,将矩阵A变换为Jordan标准形的变换矩阵为矩阵S的逆为于是注意到《现代控制理论基础》第二章(讲义)7可得eAt=SeJtS–:拉氏变换法 ()为了求出eAt,关键是必须首先求出(sI-A)的逆。一般来说,当系统矩阵A的阶次较高时,可采用递推算法。------------------------------------------------------------------------------[]考虑如下矩阵AA试用前面介绍的两种方法计算eAt。[解]方法一由于A的特征值为0和-2(λ1=0,λ2=-2),故可求得所需的变换矩阵P为P=因此,由式()可得 《现代控制理论基础》第二章(讲义)8方法二由于可得因此------------------------------------------------------------------------------:化eAt为A的有限项法(Caley-Hamilton定理法) 第四种是利用凯莱-哈密尔顿定理,化为A的有限项,然后通过求待定时间函数获得的方法。必须指出,这种方法相当系统,而且计算过程简单。设A的最小多项式阶数为m。可以证明,采用赛尔维斯特内插公式,通过求解行列式 ()即可求出。利用式()求解时,所得是以(k=0,1,2,…,m-1)和(i=1,2,3,…,m)的形式表示的。此外,也可采用如下等价的方法。将式()按最后一行展开,容易得到 ()从而通过求解下列方程组:《现代控制理论基础》第二章(讲义)9·· ()·可确定出(k=0,1,2…,m-1),进而代入式()即可求得。如果A为n×n维矩阵,且具有相异特征值,则所需确定的的个数为m=n,即有 ()如果A含有相重待征值,但其最小多项式有单根,则所需确定的的个数小于n,这里将不再进一步介绍。------------------------------------------------------------------------------[]考虑如下矩阵A试用化为A的有限项法计算。[解]矩阵A的特征方程为可得相异特征值为λ1=0,λ2=-2。由式(),可得即将上述行列式展开,可得或《现代控制理论基础》第二章(讲义)10另一种可选用的方法是采用式()。首先,由确定待定时间函数和。由于λ1=0,λ2=-2,上述两式变为求解此方程组,可得因此,------------------------------------------------------------------------------