1480 字
7 分钟
常微分方程求解(6)

nn阶常数矩阵A\mathbf{A}中的每一元素aij (i,j=1,,n)a_{ij}\ (i,j=1,\cdots,n)都是常数,则称

dxdt=Ax+f(t)\frac{\mathrm{d}\mathbf{x}}{\mathrm{d}t} = \mathbf{Ax} + \mathbf{f}(t)

为常系数线性微分方程组。

我们先介绍常系数齐次方程组

dxdt=Ax\frac{\mathrm{d}\mathbf{x}}{\mathrm{d}t} = \mathbf{Ax}

的解法,再讨论常系数非齐次方程组的情形。

一、常系数齐次线性方程组的解法#

与解常系数线性方程类似,根据方程组是齐次、线性,A\mathbf{A} 是常数矩阵的特点,我们可设

x=veλt\mathbf{x} = \mathbf{v} e^{\lambda t}

来试解,其中 vv 是常向量,λ\lambda 是常数,二者待定。代入得:

λveλt=Aveλt\lambda \mathbf{v} e^{\lambda t} = \mathbf{Av} e^{\lambda t}

注意到 v=Ev\mathbf{v}=\mathbf{Ev}, E\mathbf{E} 是单位矩阵:

E=(100010001)\mathbf{E} = \begin{pmatrix} 1 & 0 & \cdots & 0\\ 0 & 1 & \cdots & 0\\ \vdots & \vdots & \ddots & \vdots\\ 0 & 0 & \cdots & 1\\ \end{pmatrix}

移项,并约去非零因子 eλte^{\lambda t},有

(AλE)v=0(\mathbf{A} - \lambda \mathbf{E}) \mathbf{v} = 0

这是一个齐次线性代数方程组,v\mathbf{v} 的各分量是未知数。由线性代数可知,有非零解 v\mathbf{v}(即v\mathbf{v}的各分量不全为零)的充要条件是上式的系数行列式等于零,即:

det(AλE)=0\det(\mathbf{A} - \lambda \mathbf{E}) = 0

或:

D(λ)=a11λa12a1na21a22λa2nan1an2annλ=0D(\lambda) = \begin{vmatrix} a_{11}-\lambda & a_{12} & \cdots & a_{1n}\\ a_{21} & a_{22}-\lambda & \cdots & a_{2n}\\ \vdots & \vdots & \ddots & \vdots\\ a_{n1} & a_{n2} & \cdots & a_{nn}-\lambda \end{vmatrix} = 0

在线性代数中,称上面方程为矩阵A\mathbf{A}的特征方程,其根为A\mathbf{A}的特征根(或称特征值)。与此类似,我们称其和其根为常系数齐次线性微分方程组的特征方程和特征根。

如果λ=λk\lambda = \lambda_k是的一个特征根,则将它代入可求得相应的非零解 v=vk\mathbf{v} = \mathbf{v}_k。在线性代数中称这种非零向量vkv_k为矩阵A\mathbf{A}属于λk\lambda_k的特征向量,从而:

x(t)=vkeλkt\mathbf{x}(t) = \mathbf{v}_k e^{\lambda_k t}

是齐次方程的一个解。

1、特征根为单根#

设矩阵A\mathbf{A}的特征根都是单根,即有nn个不同的特征根λ1,,λn\lambda_1,\cdots,\lambda_n。又设vi\mathbf{v}_i是属于特征根λi\lambda_i的特征向量(i=1,,n)(i=1,\cdots,n),则方程组有nn个不同的解:

v1eλ1t, v2eλ2t, , vneλnt\mathbf{v}_1 e^{\lambda_1 t},\ \mathbf{v}_2 e^{\lambda_2 t},\ \cdots,\ \mathbf{v}_n e^{\lambda_n t}

由于它们在线性无关区间(,+)(-\infty, +\infty),因此它们构成的一个基本解组。于是

v=i=1ncivieλit\mathbf{v} = \sum_{i=1}^n c_i \mathbf{v}_i e^{\lambda_i t}

是的通解,其中ci(i=1,,n)c_i(i=1,\cdots,n)是任意常数。

2、特征根有复根#

若特征方程有复数根,设A\mathbf{A}是实数矩阵,复特征根必共轭地成对出现。设λ1=α+iβ\lambda_1 = \alpha + i\beta是一个特征根,v1=p+iq\mathbf{v}_1 = \mathbf{p} + i\mathbf{q}是属于λ1\lambda_1的一个特征向量,则λ2=αiβ\lambda_2 = \alpha - i\beta也是的一个特征根,且v2=v1=piq\mathbf{v}_2 = \mathbf{v}_1^* = \mathbf{p} - i\mathbf{q}是属于λ2\lambda_2的一个特征向量。因此有两个复值解

x1(t)=(p+iq)e(α+iβ)t=eαt(pcosβtqsinβt)+ieαt(psinβt+qcosβt)x2(t)=(piq)e(αiβ)t=eαt(pcosβtqsinβt)ieαt(psinβt+qcosβt)\begin{aligned} \mathbf{x}_1(t) &= (\mathbf{p} + i\mathbf{q}) e^{(\alpha + i\beta)t} \\ &= e^{\alpha t} (\mathbf{p}\cos\beta t - \mathbf{q}\sin\beta t) + i e^{\alpha t} (\mathbf{p}\sin\beta t + \mathbf{q}\cos\beta t) \\ \mathbf{x}_2(t) &= (\mathbf{p} - i\mathbf{q}) e^{(\alpha - i\beta)t} \\ &= e^{\alpha t} (\mathbf{p}\cos\beta t - \mathbf{q}\sin\beta t) - i e^{\alpha t} (\mathbf{p}\sin\beta t + \mathbf{q}\cos\beta t) \end{aligned}

它们的实部和虚部分别

eαt(pcosβtqsinβt),eαt(psinβt+qcosβt)e^{\alpha t} (\mathbf{p}\cos\beta t - \mathbf{q}\sin\beta t),\quad e^{\alpha t} (\mathbf{p}\sin\beta t + \mathbf{q}\cos\beta t)

是两个实值解。在基本解组中,λ1\lambda_1λ2\lambda_2所对应的两个复值解用这两个实值解来代替,所得解组依然是一个基本解组。

3、特征根有重根#

如果矩阵A\mathbf{A}的特征根有重根,则不一定能得如单根情况下的nn个线性无关的解,但有如下引理:

设矩阵A\mathbf{A}的特征方程有kk重特征根λ0\lambda_0,则对应于λ0\lambda_0,方程组有下述形式的kk个线性无关的解:

x(t)=(v0+t1!v1+t22!v2++tk1(k1)!vk1)eλ0t\mathbf{x}(t) = (\mathbf{v}_0 + \frac{t}{1!}\mathbf{v}_1 + \frac{t^2}{2!}\mathbf{v}_2 + \cdots + \frac{t^{k-1}}{(k-1)!}\mathbf{v}_{k-1}) e^{\lambda_0 t}

其中vi(i=0,1,,k1)v_i (i=0,1,\cdots,k-1)是某些常向量。

其证明思路是将上式代入,整理得到一系列递推方程:

(Aλ0E)v0=v1(Aλ0E)v1=v2(Aλ0E)vk2=vk1(Aλ0E)vk1=0\begin{aligned} (\mathbf{A}-\lambda_0 \mathbf{E})\mathbf{v}_0 &= \mathbf{v}_1 \\ (\mathbf{A}-\lambda_0 \mathbf{E})\mathbf{v}_1 &= \mathbf{v}_2 \\ &\vdots \\ (\mathbf{A}-\lambda_0 \mathbf{E})\mathbf{v}_{k-2} &= \mathbf{v}_{k-1} \\ (\mathbf{A}-\lambda_0 \mathbf{E})\mathbf{v}_{k-1} &= \mathbf{0} \end{aligned}

由线性代数理论知,这样可求得kk个线性无关的向量v0,v1,,vk1\mathbf{v}_0,\mathbf{v}_1,\cdots,\mathbf{v}_{k-1},从而得kk个线性无关的解。

4、通解形式总结#

设方程组

dxdt=Ax\frac{\mathrm{d}\mathbf{x}}{\mathrm{d}t} = \mathbf{Ax}

的系数矩阵A\mathbf{A}ss个不同的特征根λ1,,λs\lambda_1,\cdots,\lambda_s,其重数分别为n1,,nsn_1,\cdots,n_sn1++ns=nn_1+\cdots+n_s=n。则:

  1. 对于每一个根λi\lambda_i,方程组(3.24)存在形如 pi(j)(t)eλit(j=1,2,,ni)\mathbf{p}_i^{(j)}(t) e^{\lambda_i t}\quad (j=1,2,\cdots,n_i)nin_i个线性无关的解,其中pi(j)(t)\mathbf{p}_i^{(j)}(t)是向量函数,其分量为tt的次数不超过ni1n_i-1的多项式;
  2. 这些解线性无关,构成基本解组;
  3. 方程组的通解为 x(t)=i=1sj=1nicijpi(j)(t)eλit\mathbf{x}(t) = \sum_{i=1}^s \sum_{j=1}^{n_i} c_{ij} \mathbf{p}_i^{(j)}(t) e^{\lambda_i t} 其中cijc_{ij}是任意常数。

二、常系数非齐次线性方程组的解法#

对于右端特殊的常系数非齐次线性方程组,也可采用与线性方程类似的待定系数法求解。但一般性我们介绍变动任意常数法(变参数法)。

考虑一般非齐次线性方程组

dxdt=Ax+f(t)\frac{\mathrm{d}\mathbf{x}}{\mathrm{d}t} = \mathbf{Ax}+ \mathbf{f}(t)

A\mathbf{A}f(t)f(t)在区间(a,b)(a,b)内连续,且已知对应的齐次线性方程组

dxdt=Ax\frac{\mathrm{d}\mathbf{x}}{\mathrm{d}t} = \mathbf{Ax}

的一个基本解矩阵为X(t)\mathbf{X}(t),则其通解为

x(t)=X(t)c\mathbf{x}(t) = \mathbf{X}(t) \mathbf{c}

其中c\mathbf{c}nn维任意常向量。设

x=X(t)c(t)\mathbf{x} = \mathbf{X}(t) \mathbf{c}(t)

是非齐次线性方程组的解,则

dxdt=dXdtc(t)+X(t)dcdt=AX(t)c(t)+X(t)dcdt\begin{align*} \frac{\mathrm{d}\mathbf{x}}{\mathrm{d}t} &= \frac{\mathrm{d}\mathbf{X}}{\mathrm{d}t} \mathbf{c}(t) + \mathbf{X}(t) \frac{\mathrm{d}\mathbf{c}}{\mathrm{d}t} \\ &= \mathbf{A} \mathbf{X}(t) \mathbf{c}(t) + \mathbf{X}(t) \frac{\mathrm{d}\mathbf{c}}{\mathrm{d}t} \end{align*}

代入原方程得

X(t)dcdt=f(t)\mathbf{X}(t) \frac{\mathrm{d}\mathbf{c}}{\mathrm{d}t} = \mathbf{f}(t)

由于X(t)\mathbf{X}(t)为基本解矩阵,其行列式不为零,故可求逆X1(t)\mathbf{X}^{-1}(t),于是

dcdt=X1(t)f(t)\frac{\mathrm{d}\mathbf{c}}{\mathrm{d}t} = \mathbf{X}^{-1}(t) \mathbf{f}(t)

积分得

c(t)=t0tX1(τ)f(τ)dτ\mathbf{c}(t) = \int_{t_0}^{t} \mathbf{X}^{-1}(\tau) \mathbf{f}(\tau) \,\mathrm{d}\tau

因此非齐次项的一个特解为

x(t)=X(t)t0tX1(τ)f(τ)dτ\mathbf{x}^*(t) = \mathbf{X}(t) \int_{t_0}^{t} \mathbf{X}^{-1}(\tau) \mathbf{f}(\tau) \,\mathrm{d}\tau

从而通解为

x(t)=X(t)c+X(t)t0tX1(τ)f(τ)dτ\mathbf{x}(t) = \mathbf{X}(t) \mathbf{c} + \mathbf{X}(t) \int_{t_0}^{t} \mathbf{X}^{-1}(\tau) \mathbf{f}(\tau) \,\mathrm{d}\tau

其中c\mathbf{c}nn维任意常向量。

若给定初值条件x(t0)=x0\mathbf{x}(t_0) = \mathbf{x}_0,则可得c=X1(t0)x0\mathbf{c} = \mathbf{X}^{-1}(t_0)\mathbf{x}_0,于是初值问题

dxdt=Ax+f(t),x(t0)=x0\frac{\mathrm{d}\mathbf{x}}{\mathrm{d}t} = \mathbf{Ax} + \mathbf{f}(t),\qquad \mathbf{x}(t_0)=\mathbf{x}_0

的解为

x(t)=X(t)X1(t0)x0+X(t)t0tX1(τ)f(τ)dτ\mathbf{x}(t) = \mathbf{X}(t) \mathbf{X}^{-1}(t_0) \mathbf{x}_0 + \mathbf{X}(t) \int_{t_0}^{t} \mathbf{X}^{-1}(\tau) \mathbf{f}(\tau)\,\mathrm{d}\tau

三、消元法#

对于未知函数个数不多(例如n=2n=2)的常系数齐次或非齐次方程组,我们也可以采取反其道而行之的办法,将其化为高阶线性方程来求解。

分享

如果这篇文章对你有帮助,欢迎分享给更多人!

常微分方程求解(6)
https://www.laoguantx.cn/posts/odesolutions6/
作者
老官童鞋gogo
发布于
2025-06-10
许可协议
CC BY-NC-SA 4.0

部分信息可能已经过时

封面
Sample Song
Sample Artist
封面
Sample Song
Sample Artist
0:00 / 0:00