利用哈密顿凯莱定理求解微分方程组两种方法.docVIP

利用哈密顿凯莱定理求解微分方程组两种方法 　　摘要微分方程组在很多实际问题中（尤其是工程技术方面）都有广泛的应用，一些比较复杂的数学模型往往会导出多于一个微分方程的微分方程组，而且通过某些简化的假设和适当的变换，这种方程组又可化为一阶线性微分方程组。 　　关键词 Hamilton-Cayclay 定理基解 矩阵特征 多项式 　　中图分类号：O17文献标志码：A 　　 　　Two Ways for Hamilton-Cayclay to Solve Differential Equations 　　JU Li 　　(Henan Finance and Taxation Vocational College, Zhengzhou, Henan 450000) 　　AbstractDifferential equations are widely used in many practical problems (especially engineering techniques), some of the complex mathematical models tend to export differential equations more than one differential equation, but also through some of the simplifying assumptions and the appropriate transformation, this equation can be turned into a first-order linear differential equations. 　　Key wordsHamilton-Cayclay; theorem-based solution; matrix characteristics; polynomial 　　 　　本文考虑形如 　　（1） 　　的一阶线性常微分方程组，其中已知函数(t) ( , j = 1,2,…,n )在区间a≤t≤b上是连续的，方程组（1）关于x1,x2,…xn,及xi1,xi2,…xin是线性的。所以引进下面的记号 　　 　　这里A（t）是n??n矩阵，它的元是n2个函数(t) ( , j = 1,2,…,n )。 　　 　　这里f (t), x1,xi是n??矩阵或n维列向量。方程组（1）可以写成xi = A(t)x + f (t)，而这里主要研究一阶齐次线性微分方程组的解法，即f (t) = 0的时候xi = A(t)x。 　　而对一阶齐次线性微分方程组一般的解法相对来说比较复杂，在计算过程中容易出错，所以引入以下两种解法。这两种解法的原理是哈密顿?凯莱定理。 　　1 哈密顿?凯莱定理 　　定理1设A是数域P上一个n??n矩阵，f () = |E - A|是A的特征多项式，则 　　f (A) = An - (a11 + a22 + … + ann) An-2 + … + (-1)n |A| = 0 　　2 二种求解方法 　　由于A的特征多项式为 　　f () ≡ det (E - A) ≡ n + cn-1 n-1 + c1 + c0 　　令A的特征值为1,2,…n,不必彼此互异。 　　方法1exp At = j+1 (t) Pj 　　其中P0 = E, Pj =(A - kE), j = 1,2,3,…,n，并且r1(t),r2(t),…rn(t)为初值问题 　　，的解。 　　方法2构造纯量函数Z(t)满足方程Z(n) + Cn-1Z(n-1) +… + C1Z1 + C0Z = 0以及初值条件Z(0) = Z(0) = … = Zn-2(0) = 0, Z(n-1) (0) = 1。 　　定义 　　则exp At = qj (t)Aj。 　　其中q0(t),q1(t),…qn-1(t)为列向量q0(t) = CZ(t)的n个元素。 　　3 两种方法的比较及应用举例 　　求解线性微分方程组的方法有很多。通常利用线性代数的理论求解，Jordan标准形的方法就是其中的一种，其理论较为简洁，按实际操作起来可能比较麻烦。本文利用Hamilton-Cayclay定理以及文献的结论，将基解矩阵的计算问题进行简化，降低了求解方程组的难度，尤其是当对应的特征方程含有重根时，其效果更为突出。其中，方法一将基解矩阵的计算问题归结为求解带下三角形矩阵的齐次线性微分方程组的初值问题： 　　例1求解微分方程组的基解矩阵expAt 　　解方程组对应的特征方程为 　　det (A - E) == 2 - 2 + 2 = 0，解这个代数方程，可知矩阵A有特征