(常微分方程)第4章常系数线性方程.ppt

4.1　常系数齐次线性方程的解法　　在这一节里，我们考虑具有常系数a1,a2,…,an的齐次线性方程: 　　如果能看出它的n个线性无关的解，则可立即得到它的通解.但是怎么看出来呢？欧拉提出了一种方法，即待定指数函数法：设想方程(4.1)有指数形式的解： 　　例4.1　解方程 　　方程(4.1)的特征多项式P(λ)是λ的n次多项式.如果它可以分解为一次因式，即 的解都是方程(4.1)的解.为了解方程(4.4),注意利用归纳法容易证明公式 令i=1,2,…,r，我们便得到方程(4.1)的如下n个解： 　　定理4.1　如果方程(4.1)的特征多项式P(λ)可以分解为 　　例4.2　解方程 　 　解　写出特征多项式，并作因子分解： 　　例4.3　解方程 　　由于特征多项式可能有不是实数的根，因此定理4.1给出的基本解中有的可能不是实值函数.能否进一步得到全部都由实值函数构成的基本解组呢？回答是肯定的.　　事实上，由于方程(4.1)的系数ai都是实的，它的特征多项式是实系数多项式，因此若复数λ0=α+iβ是n0重特征根，则其共轭复数 也是n0重特征根.定理4.1给出的基本解组中将同时出现: 把式(4.9)与(4.10)中各对应函数分别相加后除以2和分别相减后除以2i，便得到2n0个实值解： 　　例4.4　求例4.3中方程的实基本解组.　　解　注意 4.2　常系数齐次线性方程组的解法　　现在转到常系数齐次线性方程组 　　1.矩阵指数函数eAt　　设x0为任意n维常向量.由定理3.1知，方程组(4.11)初值条件x(0)=x0的解φ(t)在-∞t∞上存在且唯一.把它代入方程组(4.11)，并对它从0到t积分，可得 其中,E是单位矩阵.如此迭代m次，便得 为此，对每一给定的t，在以0与t为端点的区间上，|φ(s)|有界.设M是它的一个上界,则 这样定义的矩阵函数具有下列性质：　　(1)eAt是方程组(4.11)的基本解矩阵，或者确切地说，X=eAt是下述初值问题的解： 　　应用矩阵指数函数eA(t-t0)，常系数齐次线性方程组(4.11)满足初值条件 　　2.基本解矩阵的结构　　下面我们分析基本解矩阵eAt的具体结构.　　首先考虑简单情形： 因此 即见 　　最后，考虑一般情形.由线性代数知道，若A的特征多项式det(λE-A)的分解式为 其中 　　证明　取Φ(t)=eAtT，则 　　3.待定系数法　　比较定理4.1和4.2可知，为了计算一个基本解组，对于n阶齐次线性方程(4.1)，只需知道相应的特征根及其重数.而对齐次线性方程组(4.11)，仅仅知道系数矩阵A的特征根及其重数是不够的，还要知道把A化成若尔当标准型的n阶非奇异常矩阵T，以及对应于每个特征根λi的全部若尔当小块的数目mi和这些若尔当小块的阶数kij.然而这样的矩阵T至今还没有比较简捷有效的算法.在具体求方程组(4.11)的基本解组时，下述定理比定理4.2实用. 　　定理4.3　如果方程组(4.11)的系数矩阵A的特征多项式det(λE-A)的分解式为式(4.20)，则对A的每一特征根λi(i=1,…,r)，方程组(4.11)有形如 　　证明　注意，方程组(4.11)的基本解矩阵(4.26)可以表示为 　　假如方程组(4.11)的系数矩阵A的特征多项式det(λE-A)的分解式为式(4.20).由定理4.3知，对应于特征根λi，方程组(4.11)有ni个形如式(4.26)，即形如 即 　　注4.1　无论是定理4.2、4.3,还是待定系数法，都与定理4.1一样，有一个特征多项式的分解问题，特征根中有的可能是复数.此外，用以化矩阵A为若尔当形J的矩阵T，以及从式(4.28)解出的hk，都有可能是复数阵，因此所得到的可能是复数值的基本解组.由于方程组(4.11)的系数矩阵A是实数阵，因此我们很容易进一步求得方程组(4.11)的实数值基本解组. 　　例4.5　求下述方程组 知特征根为2，其重数为3，故该方程组有如下形式的三个线性无关的解： 注意 从式(4.29)的第一式可得到三个线性无关的h0，譬如取 从而得到基本解组： 　　例4.6　求方程组 知特征根为λ1=1,λ２＝１＋2i，λ3=1-2i,它们都是单根.因此该方程组有形如 从而得到复基本解组： 4.3　算子解法与拉氏变换法　　先考虑常系数非齐次线性方程 　　我们知道，求f(t)的原函数，等价于求解微分方程： 　　(2)若P(D)=P1(D)·P2(D)，则 　　(5)设fk(t)是关于t的k次多项式，λ、ω是实数，则 　　性质(1)、(2)和(5)是显然的，我们只证明性质(3)和(4).为证明性质(3