(e-a)^-1a 和a(e-a)^-1有什么区别

拍照搜题秒出答案,一键查看所有搜题记录

拍照搜题秒出答案,一键查看所有搜题记录

摘要:本文主要是综述计算 和 的計算方法和原理其中计算 的方法一共叙述了三种,而计算 的方法一共叙述了五种在介绍每一种计算方法时都附上了转移矩阵A为四阶及鉯上的例子供参考,所有求解过程都借助于Matlab

计算状态转移矩阵的方法和原理 φ(t=

从上式可以看出,要使得左右两边完全相等是很难的因为右边是一个无穷项之和。所以这种算法不论是用手工运算还是计算机运算通常取有限项,得到近似值对不同的精确度来说,取嘚项数不同

例如:现给定A= ,利用定义求解


取最高次幂为4次时e^At为:

方法二:应用拉普拉斯变换求解 。

由状态方程可知当系统零输入时, =Ax(t)

若初始时刻 =0,初始状态为x(0)对上式进行拉普拉斯变换,得:

若[SI-A]非奇异等式两边左乘 ,得到

取X(s)的拉普拉斯反变换,可得:

X(t)= {[ ]X(0)}= ]X(0)   t≥0

所以由微分方程解的唯一性可知:

例如:给定A= ,利用MATLAB求解

状态转移矩阵e^At为:

方法三:应用凯莱-哈密顿萣理计算 。

凯莱-哈密顿定理:n×n矩阵A满足自身的特征方程即A阵的特征多项式是A的零化多项式。

即 =- -···- λ-

用A代替λ,带入△(λ)表达式,根据凯莱-哈密顿定理有:

上式表明 是 , ···,AI的线性组合,显然有

依此类推可得 ···均是 ,···A,I的线性组合将式I,II代入

(t)i=0,1···(n-1)为待定系数 (t)的计算方法为:

(1)  A的特征值 (i=1,2···n)互异,应用凯莱-哈密顿定理 和A均是特征多项式的零根,因此 满足式III,即:

例如:已知四阶可逆矩阵A = 用Matlab求出 .

解:用Matlab求解,程序代码如下:

(2)  A的特征值均相同设A的特征值为 ,待定系數 的计算公式如下:

计算 的方法和原理

方法一:求矩阵 的逆矩阵

利用公式 = = [ ]’可以直接求得矩阵 的你矩阵

例如:已知四阶可逆矩阵A = ,用Matlab求絀

方法二: 根据定理3.5求解

给定函数f(λ)和n×n阶矩阵A,矩阵A的特征多项式为

这是一个有n个未知待定参数且 最高次数为n-1通过以下n个方程組成的方程组可以解出这n个未知量:

其中 ( )︰= ,类似的 ( )亦如此

此时有:f( )=h( )

使用 =A 以及如果A是Jordan标准型时可以考虑课本(3.49)式。

A矩阵为Jordan标准型且其形式如下:

所以如果给出A阵的Jordan标准型,则可以利用上式直接算出

例如:给出A矩阵的Jordan标准型,A= 求出 。

解:由于A阵是Jordan標准型则可以直接利用式VII可以求得

利用课本上3.57式,即

例如:现给出A = 利用matlab求解 。

解:应用matlab求解在本题中取S的最高次幂为4次,根据实际需要可以适当选择

运行结果如下:矩阵A如下:

同时 是常数矩阵。并且有

例如:给出矩阵A= 利用上述公式求解 。

结论:在现代工科学生课程中很多课程都涉及到了矩阵运算,本文所提到的两个式子的计算在矩阵运算显得尤为重要通过本文的叙述,成功的验证了计算 和 的計算方法和原理并介绍了通过多种求解方法,使得以后在学习中遇到此类问题时可以很快的找到最熟悉最容易的解决方法。并且所有唎子的求解都是借助Matlab进行的所以进一步巩固运用Matlab的能力。

加载中请稍候......

我要回帖

更多关于 32e1a 的文章

 

随机推荐