设n阶常数矩阵A中的每一元素aij (i,j=1,⋯,n)都是常数,则称
dtdx=Ax+f(t)
为常系数线性微分方程组。
我们先介绍常系数齐次方程组
dtdx=Ax
的解法,再讨论常系数非齐次方程组的情形。
一、常系数齐次线性方程组的解法#
与解常系数线性方程类似,根据方程组是齐次、线性,A 是常数矩阵的特点,我们可设
x=veλt来试解,其中 v 是常向量,λ 是常数,二者待定。代入得:
λveλt=Aveλt注意到 v=Ev, E 是单位矩阵:
E=10⋮001⋮0⋯⋯⋱⋯00⋮1移项,并约去非零因子 eλt,有
(A−λE)v=0这是一个齐次线性代数方程组,v 的各分量是未知数。由线性代数可知,有非零解 v(即v的各分量不全为零)的充要条件是上式的系数行列式等于零,即:
det(A−λE)=0或:
D(λ)=a11−λa21⋮an1a12a22−λ⋮an2⋯⋯⋱⋯a1na2n⋮ann−λ=0在线性代数中,称上面方程为矩阵A的特征方程,其根为A的特征根(或称特征值)。与此类似,我们称其和其根为常系数齐次线性微分方程组的特征方程和特征根。
如果λ=λk是的一个特征根,则将它代入可求得相应的非零解 v=vk。在线性代数中称这种非零向量vk为矩阵A属于λk的特征向量,从而:
x(t)=vkeλkt是齐次方程的一个解。
1、特征根为单根#
设矩阵A的特征根都是单根,即有n个不同的特征根λ1,⋯,λn。又设vi是属于特征根λi的特征向量(i=1,⋯,n),则方程组有n个不同的解:
v1eλ1t, v2eλ2t, ⋯, vneλnt由于它们在线性无关区间(−∞,+∞),因此它们构成的一个基本解组。于是
v=i=1∑ncivieλit是的通解,其中ci(i=1,⋯,n)是任意常数。
2、特征根有复根#
若特征方程有复数根,设A是实数矩阵,复特征根必共轭地成对出现。设λ1=α+iβ是一个特征根,v1=p+iq是属于λ1的一个特征向量,则λ2=α−iβ也是的一个特征根,且v2=v1∗=p−iq是属于λ2的一个特征向量。因此有两个复值解
x1(t)x2(t)=(p+iq)e(α+iβ)t=eαt(pcosβt−qsinβt)+ieαt(psinβt+qcosβt)=(p−iq)e(α−iβ)t=eαt(pcosβt−qsinβt)−ieαt(psinβt+qcosβt)它们的实部和虚部分别
eαt(pcosβt−qsinβt),eαt(psinβt+qcosβt)是两个实值解。在基本解组中,λ1和λ2所对应的两个复值解用这两个实值解来代替,所得解组依然是一个基本解组。
3、特征根有重根#
如果矩阵A的特征根有重根,则不一定能得如单根情况下的n个线性无关的解,但有如下引理:
设矩阵A的特征方程有k重特征根λ0,则对应于λ0,方程组有下述形式的k个线性无关的解:
x(t)=(v0+1!tv1+2!t2v2+⋯+(k−1)!tk−1vk−1)eλ0t其中vi(i=0,1,⋯,k−1)是某些常向量。
其证明思路是将上式代入,整理得到一系列递推方程:
(A−λ0E)v0(A−λ0E)v1(A−λ0E)vk−2(A−λ0E)vk−1=v1=v2⋮=vk−1=0由线性代数理论知,这样可求得k个线性无关的向量v0,v1,⋯,vk−1,从而得k个线性无关的解。
4、通解形式总结#
设方程组
dtdx=Ax的系数矩阵A有s个不同的特征根λ1,⋯,λs,其重数分别为n1,⋯,ns,n1+⋯+ns=n。则:
- 对于每一个根λi,方程组(3.24)存在形如
pi(j)(t)eλit(j=1,2,⋯,ni)
的ni个线性无关的解,其中pi(j)(t)是向量函数,其分量为t的次数不超过ni−1的多项式;
- 这些解线性无关,构成基本解组;
- 方程组的通解为
x(t)=i=1∑sj=1∑nicijpi(j)(t)eλit
其中cij是任意常数。
二、常系数非齐次线性方程组的解法#
对于右端特殊的常系数非齐次线性方程组,也可采用与线性方程类似的待定系数法求解。但一般性我们介绍变动任意常数法(变参数法)。
考虑一般非齐次线性方程组
dtdx=Ax+f(t)设A及f(t)在区间(a,b)内连续,且已知对应的齐次线性方程组
dtdx=Ax的一个基本解矩阵为X(t),则其通解为
x(t)=X(t)c其中c是n维任意常向量。设
x=X(t)c(t)是非齐次线性方程组的解,则
dtdx=dtdXc(t)+X(t)dtdc=AX(t)c(t)+X(t)dtdc代入原方程得
X(t)dtdc=f(t)由于X(t)为基本解矩阵,其行列式不为零,故可求逆X−1(t),于是
dtdc=X−1(t)f(t)积分得
c(t)=∫t0tX−1(τ)f(τ)dτ因此非齐次项的一个特解为
x∗(t)=X(t)∫t0tX−1(τ)f(τ)dτ从而通解为
x(t)=X(t)c+X(t)∫t0tX−1(τ)f(τ)dτ其中c为n维任意常向量。
若给定初值条件x(t0)=x0,则可得c=X−1(t0)x0,于是初值问题
dtdx=Ax+f(t),x(t0)=x0的解为
x(t)=X(t)X−1(t0)x0+X(t)∫t0tX−1(τ)f(τ)dτ三、消元法#
对于未知函数个数不多(例如n=2)的常系数齐次或非齐次方程组,我们也可以采取反其道而行之的办法,将其化为高阶线性方程来求解。