4.2 高阶线性微分方程

上一节讨论了几种一阶方程以及三种可降阶的高阶方程的求解问题.本节要讨论的是另一种在科学技术中经常遇到的高阶方程——高阶线性微分方程.主要内容有:高阶线性微分方程的有关概念以及解的性质和结构,常系数高阶线性微分方程的求解方法.

2.1 高阶线性微分方程举例

例 2.1

弹簧的机械振动 图 4.4 表示一个简单的减振装置. 质量为 $m$ 的物体安置在弹簧上, 当物体稳定在位置 $O$ 时, 作用在物体上的重力和弹性力大小相等、方向相反, 这个位置是物体的平衡位置. 设在垂直方向有一随时间周期变化的外界强迫力 $f_1(t)=H\sin pt$ 作用在物体上, 物体将受外力驱使而上下振动, 求振动过程中位移与时间所满足的关系方程.

图4.4

解

取物体的平衡位置为坐标原点, $x$ 轴的方向垂直向下. 设振动开始时刻 t=0, 时刻 $t$ 物体离开平衡位置的位移为 $x(t)$ , 先来建立 $x(t)$ 所应满足的方程.

物体在运动过程中受到外界强迫力、弹性力和介质阻力三个力的作用.

由 Hooke 定律知弹性力

$$f=-kx,$$

其中 $k$ 为弹簧的劲度系数. 设物体在振动过程中所受介质(如空气)的阻力 $f_0$ 与运动速度 $v$ 成正比, 即

$$f_0=-\mu v=-\mu \frac{\mathrm{d}x}{\mathrm{d}t},$$

其中 $\mu$ 为介质的阻尼系数, 负号表示阻力方向与速度方向相反. 于是根据 Newton 第二定律可得公式

$$ma=-kx-\mu \frac{\mathrm{d}x}{\mathrm{d}t}+f_1(t).$$

由于加速度 $a=\frac{d^{2}x}{d{t}^2}$ ，故 $x(t)$ 应满足的微分方程为

$$m\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{t}^2}+\mu \frac{\mathrm{d}x}{\mathrm{d}t}+kx=H\sin pt,$$

将微分方程改写为

$$\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{t}^2}+2\delta \frac{\mathrm{d}x}{\mathrm{d}t}+{\omega }^{2}x=h\sin pt,\text{\hspace{1em}(2.1)}$$

其中 $\delta =\frac{\mu }{2m},\omega =\sqrt{\frac{k}{m}},h=\frac{H}{m}.$

由于已取物体由平衡位置开始运动的时刻为 t=0 且平衡位置为原点, 因此运动还应满足条件

$$x{\mid }_{t=0}=0,v{\mid }_{t=0}=\frac{\mathrm{d}x}{\mathrm{d}t}{|}_{t=0}=0.\text{\hspace{1em}(2.2)}$$

于是,振动过程中位移随时间的变化规律应满足微分方程(2.1)与初值条件(2.2),方程(2.1)是一个二阶线性微分方程.

例 2.2

RLC 电路中的电压变化规律

图 4.5 是 RLC 电路图, 其中 $R$ 为电阻, $L$ 为电感, $C$ 为电容. 设电容器已经充电, 它的两极板间电压为 $E$. 当开关 $S$ 闭合后, 电容器放电, 此时电路中将有电流 $i$ 通过, 产生电磁振荡. 求电容器两极板间电压 $u_C$ 随时间变化规

图4.5

律所满足的方程.

解

根据回路电压定律可知,电容、电感、电阻上的电压 $u_C,u_L,u_R$ 应有如下关系:

$$u_L+u_R+u_C=0.\text{\hspace{1em}(2.3)}$$

由于 $i=C\frac{d{u}_c}{dt}$ ，故

$$u_R=Ri=RC\frac{\mathrm{d}{u}_C}{\mathrm{d}t},u_L=L\frac{\mathrm{d}i}{\mathrm{d}t}=LC\frac{{\mathrm{d}}^2{u}_C}{\mathrm{d}{t}^2},$$

代入(2.3)式,得

$$LC\frac{{\mathrm{d}}^2{u}_c}{\mathrm{d}{t}^2}+RC\frac{\mathrm{d}{u}_c}{\mathrm{d}t}+u_c=0,$$

或

$$\frac{{\mathrm{d}}^2{u}_c}{\mathrm{d}{t}^2}+\frac{R}{L}\frac{\mathrm{d}{u}_c}{\mathrm{d}t}+\frac{1}{LC}{u}_c=0.\text{\hspace{1em}(2.4)}$$

这也是一个二阶线性微分方程. 设开关闭合时刻 t=0, 按所给条件 $u_c$ 尚应满足下列初值条件

$$u_C{\mid }_{t=0}=E,\frac{\mathrm{d}{u}_C}{\mathrm{d}t}{|}_{t=0}=\frac{1}{C}i{|}_{t=0}=0.$$

二阶线性微分方程的一般形式为

$$\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{t}^2}+P_1(t)\frac{\mathrm{d}x}{\mathrm{d}t}+P_2(t)x=F(t),\text{\hspace{1em}(2.5)}$$

其中未知函数 $x(t)$ 及其一阶、二阶导数均为一次的. 当 $F(t)\equiv 0$ 时, 方程(2.5)称为二阶线性齐次微分方程; 当 $F(t)$ 不恒为零时, 称为二阶线性非齐次微分方程. 例如, 方程(2.1)就是一个二阶线性非齐次微分方程, 而方程(2.4)是一个二阶线性齐次微分方程.

一般来说,一个 $n$ 阶微分方程,如果其中的未知函数及其各阶导数都是一次的,则称其为 $n$ 阶线性微分方程,简称为 $n$ 阶线性方程.它的一般形式为

$$p_0(t)x^{(n)}(t)+p_1(t)x^{(n-1)}(t)+\cdots +p_{n-1}(t)\dot{x}(t)+p_n(t)x(t)=f(t).$$

若 $p_0(t)\ne 0$ ，则上式可写成

$$x^{(n)}+P_1(t)x^{(n-1)}+\cdots +P_{n-1}(t)\dot{x}+P_n(t)x=F(t),\text{\hspace{1em}(2.6)}$$

其中

$$P_i(t)=\frac{p_i(t)}{p_0(t)},i=1,2,\dots ,n;F(t)=\frac{f(t)}{p_0(t)}.$$

当 $F(t)$ 不恒为零时, 方程(2.6)称为 $n$ 阶线性非齐次方程; 当 $F(t)\equiv 0$ 时, 方程(2.6)变为

$$x^{(n)}+P_1(t)x^{(n-1)}+\cdots +P_{n-1}(t)\dot{x}+P_n(t)x=0,\text{\hspace{1em}(2.7)}$$

称为与方程(2.6)对应的线性齐次方程.

在本章第一节中我们已经指出， $n$ 阶微分方程的初值条件由未知函数及其从一阶到 $n-1$ 阶的导数在初始点的值所构成。对于方程(2.6)与(2.7)而言，其初值条件的一般形式可写成

$$x(t_0)=x_0,\dot{x}(t_0)={\dot{x}}_0,\dots ,x^{(n-1)}(t_0)=x_0^{(n-1)},\text{\hspace{1em}(2.8)}$$

其中 $t_0$ 为初始时刻, $x_0,{\dot{x}}_0,\cdots ,x_0^{(n-1)}$ 为一组确定的常数.

可以证明(从略):若方程(2.6)中的系数 $P_1(t),P_2(t),\cdots ,P_n(t)$ 以及 $F(t)$ 均在区间 $(a,b)$ 内连续,则方程(2.6)存在唯一的满足初值条件(2.8)的解 $x(t)$ , $t\in (a,b)$ ,其中 $t_0\in (a,b),(x_0,{\dot{x}}_0,\cdots ,x_0^{(n-1)})\in {\mathbf{R}}^n$ .此结论称为$n$阶线性微分方程解的存在唯一性定理.

2.2 线性微分方程解的结构

1. 线性齐次微分方程解的结构

为书写简便起见,我们引入记号

$$L(x)=\frac{{\mathrm{d}}^{n}x}{\mathrm{d}{t}^n}+P_1(t)\frac{{\mathrm{d}}^{n-1}x}{\mathrm{d}{t}^{n-1}}+\cdots +P_{n-1}(t)\frac{\mathrm{d}x}{\mathrm{d}t}+P_n(t)x,\text{\hspace{1em}(2.9)}$$

从而 $n$ 阶线性齐次方程(2.7)可简洁地写为

$$L(x)=0,\text{\hspace{1em}(2.10)}$$

这里

$$L()=\frac{{\mathrm{d}}^n}{\mathrm{d}{t}^n}+P_1(t)\frac{{\mathrm{d}}^{n-1}}{\mathrm{d}{t}^{n-1}}+\cdots +P_{n-1}(t)\frac{\mathrm{d}}{\mathrm{d}t}+P_n(t)$$

称为线性微分算子,它的含义是如果将它作用于某一函数 $x(t)$ , 就意味着对 $x(t)$ 求出直到 $n$ 阶的各阶导数 $x^{(0)},x^{(1)},\cdots ,x^{(n)}$ , 再分别乘 $P_n(t),P_{n-1}(t),\cdots ,P_1(t)$ , 1 后相加, 即得 (2.9) 式 (这里零阶导数 $x^{(0)}$ 即为 $x$ 自身).

容易看出,线性微分算子具有以下性质:

(1) $L(0)=0;$

$$L\left(C_1{x}_1+C_2{x}_2+\cdots +C_n{x}_n\right)=C_{1}L\left(x_1\right)+C_{2}L\left(x_2\right)+\cdots +C_{n}L\left(x_n\right),$$

其中 $C_1,C_2,\cdots ,C_n$ 为任意常数.

定理 2.1 解的叠合性

若 $x_1,x_2,\cdots ,x_n$ 均是线性齐次方程(2.10)(或方程(2.7))的解，则

$$x=C_1{x}_1+C_2{x}_2+\cdots +C_n{x}_n\text{\hspace{1em}(2.11)}$$

也是方程(2.10)的解,其中 $C_1,C_2,\cdots ,C_n$ 为任意常数(可以是复数).

证明

将由(2.11)式所表示的函数 $x(t)$ 代入线性齐次方程(2.10)，得

由线性微分算子 $L()$ 的性质(1)可知, $x(t)\equiv 0$ 必定是线性齐次微分方程的解,称为平凡解.线性齐次方程必有平凡解以及解的叠合性是此类方程解的两个重要特性.

$$\begin{array}{rl}L(x)& =L\left(C_1{x}_1+C_2{x}_2+\cdots +C_n{x}_n\right)\\ & =C_{1}L\left(x_1\right)+C_{2}L\left(x_2\right)+\cdots +C_{n}L\left(x_n\right),\end{array}$$

由于 $x_i$ （ $i=1,2,\cdots ,n$ ）是方程（2.10）的解，故有 $L(x_i)=0$ （ $i=1,2,\cdots ,n$ ），从而 $L(x)=0$ ，即由（2.11）式所表示的函数 $x(t)$ 也是方程（2.10）的解。

对于 $n$ 阶线性齐次方程(2.10)而言,既然解(2.11)中含有 $n$ 个任意常数 $C_1$ , $C_2$ , $\cdots$ , $C_n$ , 是否它就是方程(2.10)的通解呢?一般说来,这是不一定的.例如,容易验证 $x_1=\sin t$ 与 $x_2=2\sin t$ 都是二阶线性齐次方程 $\ddot{x}+x=0$ 的解,由定理2.1 知

$$x=C_1\sin t+C_{2}2\sin t=(C_1+2C_2)\sin t\text{\hspace{1em}(2.12)}$$

也是其解. 表达式(2.12)表面上含有两个任意常数 $C_1,C_2$ ，但是它们可以合并成为一个任意常数 $C=C_1+2C_2$ ，因而它并非 $\ddot{x}+x=0$ 的通解. 那么具有 $n$ 个任意常数的解在什么情况下才是通解呢？要回答这一问题，我们还得引入下列关于函数组线性无关的新概念.

定义 2.1 线性相关与线性无关

设 $f_i(t)$ ( $i=1,2,\cdots ,n$ ) 是定义在区间 $I$ 上的 $n$ 个函数，若存在 $n$ 个不全为零的常数 $C_i$ ( $i=1,2,\cdots ,n$ ) 使关系式

$$C_1{f}_1(t)+C_2{f}_2(t)+\cdots +C_n{f}_n(t)=0\text{\hspace{1em}(2.13)}$$

对区间 $I$ 上任何 $t$ 值均成立，则称函数组 $f_i(t)$ (i=1,2, $\cdots$ ,$n$) 在区间 $I$ 上线性相关；若(2.13)式仅当 $C_i$ ( $i=1,2,\cdots ,n$ ) 全为零时成立，则称函数组 $f_i(t)$ ( $i=1,2,\cdots ,n$ ) 在 $I$ 上线性无关或线性独立.

当 $f_i$ 在$I$上线性相关时,由于(2.13)式中 $C_i$ 不全为零,不妨假设 $C_n\ne 0$ ,从而(2.13)式可改写为

$$f_n(t)=-\frac{C_1}{C_n}{f}_1(t)-\frac{C_2}{C_n}{f}_2(t)-\cdots -\frac{C_{n-1}}{C_n}{f}_{n-1}(t),\text{\hspace{1em}(2.14)}$$

函数组线性相关(无关)的定义与线性代数中向量组线性相关(无关)的定义在形式上是一样的. 但是向量组 ${\alpha }_i,i=1,\cdots ,n$ 线性相关是要求存在不全为零的 $C_i,i=1,2,\cdots ,n$ , 使

也就是 $f_i(t)$ 中的某一个函数可由其他n-1个函数的线性组合来表示.反之若(2.14)成立,则(2.13)式成

$$C_1{\alpha }_1+C_2{\alpha }_2+\cdots +C_n{\alpha }_n=0$$

对向量的每一分量均成立；而函数组 $f_i(t),t\in I,i=1,\cdots ,n$ 线性相关是要求(2.13)式在区间$I$上恒成立，线性无关也具有这样的差异.

立.因此, $n$ 个函数是否线性相关与是否能将其中一个函数表示为其他 $n-1$ 个函数的线性组合是等价的.

例 2.3

证明函数组 $1,x,x^2,\cdots ,x^{n-1}$ 在任何区间 $I$ 上线性无关.

证明

使用反证法,若它们线性相关,则必存在 $n$ 个不全为零的常数 $C_i$ (i=0, 1, $\cdots$ ,n-1),使

$$C_0+C_{1}x+C_2{x}_2+\cdots +C_{n-1}{x}^{n-1}=0\text{\hspace{1em}(2.15)}$$

对区间 $I$ 上的一切 $x$ 值成立.但由于方程(2.15)是 $x$ 的 n-1 次代数方程,由代数学基本定理可知,它至多只有 n-1 个实根,也就是说,至多只有 $I$ 上的 n-1 个点使(2.15)式成立.这一矛盾说明,要使(2.15)式在区间 $I$ 上成立,只能是所有 $C_i$ ( $i=0,1,\cdots ,n-1$ ) 均为零.故所给函数组在任何区间 $I$ 线性无关.

由(2.14)式可见,对于两个函数 $f_1(t)$ 与 $f_2(t),t\in I$ 而言,它们线性相关还是线性无关是很容易判别的.若它们线性相关,则其中至少有一个函数可由另一函数线性表示,比如 $f_2(t)=C{f}_1(t)$ ,其中$C$为常数, $\forall t\in I$ ,即

$$\frac{f_2(t)}{f_1(t)}=C,\forall t\in I.\text{\hspace{1em}(2.16)}$$

因此,若两函数之比在某一区间 $I$ 上为一常数,则此两函数在区间 $I$ 上线性相关;若 $\frac{f_2}{f_1}$ 在 $I$ 上不等于常数,则在区间 $I$ 上线性无关.

例如,函数 $e^{2t}$ 与 $e^{3t}$ 在任何区间线性无关,因为

$$\frac{{\mathrm{e}}^{2t}}{{\mathrm{e}}^{3t}}={\mathrm{e}}^{-t}\ne C.$$

当函数组 $f_1$ 与 $f_2$ 在$I$上线性相关时,由(2.16)式可见

$$\begin{array}{rl}{C}_1{f}_1+C_2{f}_2& =C_1{f}_1+C_{2}C{f}_1=(C_1+C_{2}C)f_1,\\ & \forall t\in I.\end{array}$$

即表达式 $C_1{f}_1+C_2{f}_2$ 中的两个任意常数 $C_1,C_2$ 可以合并成为一个. 若 $f_1$ 与 $f_2$ 在 $I$ 上线性无关, 则由于其中任何一个函数都不是另一个函数的倍数, 故表达式 $C_1{f}_1+C_2{f}_2$ 中的两个任意常数不能合并, 或者说 $C_1$ 与 $C_2$ 是两个独立的任意常数. 因此, 对于二阶线性齐次方程

由此可见,对于任一线性表达式 $C_1{\phi }_1(t)+C_2{\phi }_2(t)$ 来说,常数 $C_1$ 与 $C_2$ 能否合并成一个,即 $C_1{\phi }_1(t)+C_2{\phi }_2(t)\equiv \overline{C}{\phi }_1(t)$ 并不决定于这两个常数本身,而取决于函数 ${\phi }_1$ 与 ${\phi }_2$ 线性相关与否,即 $C_1$ 与 $C_2$ 能合并为一个常数的充要条件是 ${\phi }_1$ 与 ${\phi }_2$ 线性相关,此结论可以推广到有限个常数.

$$\ddot{x}+P_1(t)\dot{x}+P_2(t)x=0,$$

如果我们能求得它的任意两个线性无关的特解 $x_1(t),x_2(t)$ ，那么它的通解就可表示为

$$x(t)=C_1{x}_1(t)+C_2{x}_2(t),$$

其中 $C_1$ 与 $C_2$ 为任意常数. 所以求线性齐次方程通解的关键在于寻找它的线性无关的特解.

证明这一结论对于 $n$ 阶线性齐次方程也是正确的. 易得

定理 2.2

对于 $n$ 阶线性齐次微分方程(2.7)，若求得它 $n$ 个线性无关的特解 $x_1(t)$ , $x_2(t)$ , $\cdots$ , $x_n(t)$ ，则它的通解就是

$$x(t)=C_1{x}_1(t)+C_2{x}_2(t)+\cdots +C_n{x}_n(t),$$

其中 $C_1,\cdots ,C_n$ 为任意常数.

证明

由解的叠合性(定理2.1)知, $x(t)$ 也是方程(2.7)的解.要说明它是通解,只需证明常数 $C_1,\cdots ,C_n$ 中任意多个都不能合并.用反证法,不妨假设 $C_k\ne 0$ 且 $C_1$ , $C_2,\cdots ,C_k$ (1<k<n)可以合并为 ${\overline{C}}_k$ (否则可调整解的顺序),这意味着

$$C_1{x}_1(t)+C_2{x}_2(t)+\cdots +C_k{x}_k(t)\equiv {\overline{C}}_k{x}_k(t).$$

于是有

$$\left({\overline{C}}_k-C_k\right)x_k(t)\equiv C_1{x}_1(t)+C_2{x}_2(t)+\cdots +C_{k-1}{x}_{k-1}(t).$$

若 ${\overline{C}}_k-C_k=0$ ，则 $x_1(t),\cdots ,x_{k-1}(t)$ 线性相关，从而 $x_1(t),\cdots ,x_n(t)$ 也必线性相关，矛盾。若 ${\overline{C}}_k-C_k\ne 0$ ，则 $x_k(t)\equiv \frac{1}{{\overline{C}}_k-C_k}\left[C_1{x}_1(t)+\cdots +C_{k-1}{x}_{k-1}(t)\right]$ ，从而 $x_1(t),x_2(t),\cdots ,x_k(t)$ 线性相关，同样矛盾。

然而对于线性齐次方程(2.7)，当我们求得它的 $n$ 个特解后，怎样来判别它们是否线性无关呢？下面的定理将给出一个简便的方法。

定理 2.3 解的线性无关判别法

设 $x_1(t),\cdots ,x_n(t)$ 是 $n$ 阶线性齐次方程 (2.7) 的 $n$ 个定义在区间 $I$ 上的解，则它们在 $I$ 上线性无关的充要条件是，在 $I$ 上存在一点 $t_0$ ，使由这 $n$ 个解及其各阶导数在 $t_0$ 处所构成的行列式

$$w(t_0)=\left|\begin{array}{cccc}{x}_1(t_0)& x_2(t_0)& \dots & x_n(t_0)\\ {\dot{x}}_1(t_0)& {\dot{x}}_2(t_0)& \dots & {\dot{x}}_n(t_0)\\ ⋮& ⋮& & ⋮\\ x_1^{(n-1)}(t_0)& x_2^{(n-1)}(t_0)& \dots & x_n^{(n-1)}(t_0)\end{array}\right|\ne 0^{\stackrel{◯}{1}},$$

$w(t_0)$ 称为解组 $x_1(t),\cdots ,x_n(t)$ 在 $t_0$ 处的 Wronski 行列式.

* 证 充分性 设 $w(t_0)\ne 0$ , 欲证 $x_1(t),\cdots ,x_n(t)$ 在 $I$ 上线性无关. 设有 $n$ 个常数 $C_i$ ( $i=1,2,\cdots ,n$ ) 使

$$C_1{x}_1(t)+C_2{x}_2(t)+\cdots +C_n{x}_n(t)\equiv 0,\forall t\in I.$$

意到 $x_i(t)$ 为方程(2.7)的解，它们必 $n$ 阶可导，上式两端对 $t$ 逐次求导，得

$$\begin{array}{l}{C}_1{\dot{x}}_1(t)+C_2{\dot{x}}_2(t)+\cdots +C_n{\dot{x}}_n(t)\equiv 0,\forall t\in I,\\ C_1{\ddot{x}}_1(t)+C_2{\ddot{x}}_2(t)+\cdots +C_n{\ddot{x}}_n(t)\equiv 0,\forall t\in I,\\ \dots \dots \\ C_1{x}_1^{(n-1)}(t)+C_2{x}_2^{(n-1)}(t)+\cdots +C_n{x}_n^{(n-1)}(t)\equiv 0,\forall t\in I.\end{array}$$

特别地,对于 $t_0\in I$ 有

$$\{\begin{array}{l}{C}_1{x}_1(t_0)+C_2{x}_2(t_0)+\cdots +C_n{x}_n(t_0)=0,\\ C_1{\dot{x}}_1(t_0)+C_2{\dot{x}}_2(t_0)+\cdots +C_n{\dot{x}}_n(t_0)=0,\\ \dots \dots \\ C_1{x}_1^{(n-1)}(t_0)+C_2{x}_2^{(n-1)}(t_0)+\cdots +C_n{x}_n^{(n-1)}(t_0)=0.\end{array}$$

将这一方程组看作是以 $C_1,\cdots ,C_n$ 为未知数的 $n$ 元线性齐次代数方程组，由于其系数行列式为 $w(t_0)\ne 0$ ，由 Cramer 法则知仅有零解，即 $C_1=C_2=\cdots =C_n=0$ ，所以 $x_1(t),\cdots ,x_n(t)$ 在 $I$ 线性无关.

必要性 只需证明若 $w(t)\equiv 0(\forall t\in I)$ ，则 $x_1(t),\cdots ,x_n(t)$ 必在 $I$ 上线性相关即可。实际上，我们可以证明只要存在一点 $t_0\in I$ 使 $w(t_0)=0$ ，则 $x_1(t),\cdots ,x_n(t)$ 必在 $I$ 上线性相关。考察以 $C_1,C_2,\cdots ,C_n$ 为未知数的线性代数方程组

$$\{\begin{array}{l}{C}_1{x}_1(t_0)+C_2{x}_2(t_0)+\cdots +C_n{x}_n(t_0)=0,\\ C_1{\dot{x}}_1(t_0)+C_2{\dot{x}}_2(t_0)+\cdots +C_n{\dot{x}}_n(t_0)=0,\\ \dots \dots \\ C_1{x}_1^{(n-1)}(t_0)+C_2{x}_2^{(n-1)}(t_0)+\cdots +C_n{x}_n^{(n-1)}(t_0)=0.\end{array}$$

由于其系数行列式 $w(t_0)=0$ ，故此方程组有非零解，设 $C_1^0,C_2^0,\cdots ,C_n^0$ 为其一组非零解，其中 $C_i^0$ （ $i=1,2,\cdots ,n$ ）不全为零.

(2.17)

构造解 $x_1(t)$ , $x_2(t)$ , $\cdots$ , $x_n(t)$ 的线性组合注意:定理2.3的充要条件,仅当 $x_i(t),i=1,\cdots ,n$ 是$n$阶线性方程的解时才成立.对于一般的n-1阶可导函数组 $x_i(t),t\in I,i=1,\cdots ,n.$ $\exists t_0\in I$ 有 $w(t_0)\ne 0$ 只是 $x_i(t)$ 线性无关的充分条件,而非必要条件.换句话说,即使 $w(t)\equiv 0,t\in I$ ,函数组 $x_i(t)$上线性无关.本章定理2.3之后给出了例子.其原因在于,若 $I$上线性无关.本章定理2.3之后给出了例子.其原因在于,若 $w(t_0)=0$ 由(2.17)式求得的非零解 $C_i$ ( $i=1,\cdots ,n$ )只说明 $x_i(t)$ 在 $t_0$ 处所构成的向量组 $x_i(t_0),i=1,\cdots ,n$ 线性相关.尽管 $w(t)\equiv 0,t\in I$ ,对 $\forall t_0\in I$ 均可求得不全为零的常数 $C_i$ ,但对不同的 $t_0$ ,所求得的 $C_i$ 可能不同,不能保证对同一组 $C_i$ ,使 $C_1{x}_1+\cdots +C_n{x}_n\equiv 0$ 在$I$上成立.当 $x_i(t)$ 是方程的解时,上述结论是由解的唯一性所保证的.

$$x(t)=C_1^0{x}_1(t)+C_2^0{x}_2(t)+\cdots +C_n^0{x}_n(t),t\in I,$$

由定理2.1 可知, 它是方程(2.7)的解, 而且由(2.17)知

$$x\left(t_0\right)=C_1^0{x}_1\left(t_0\right)+C_2^0{x}_2\left(t_0\right)+\cdots +C_n^0{x}_n\left(t_0\right)=0,$$ $$x^{(k)}\left(t_0\right)=C_1^0{x}_1^{(k)}\left(t_0\right)+C_2^0{x}_2^{(k)}\left(t_0\right)+\cdots +C_n^0{x}_n^{(k)}\left(t_0\right)=0,k=1,2,\dots ,n-1.$$

所以, 解 $x(t)$ 满足初值条件

$$x\left(t_0\right)=\dot{x}\left(t_0\right)=\ddot{x}\left(t_0\right)=\cdots =x^{(n-1)}\left(t_0\right)=0.\text{\hspace{1em}(2.18)}$$

但是,显然 $x\equiv 0$ 也是方程(2.7)的解,且满足初值条件(2.18).由解的唯一性可知

$$x(t)=C_1^0{x}_1(t)+C_2^0{x}_2(t)+\cdots +C_n^0{x}_n(t)\equiv 0,\forall t\in I,$$

从而 $x_1(t),\cdots ,x_n(t)$ 在 $I$ 线性相关.

由此定理的证明可见，对于由方程(2.7)的 $n$ 个解 $x_1(t),x_2(t),\cdots ,x_n(t)$ 所构成的Wronski行列式 $w(t)$ 只有两种可能，或者 $w(t)\ne 0,\forall t\in I$ ；或者 $w(t)\equiv 0,\forall t\in I$ . 换句话说，只要存在一点 $t_0\in I$ 使 $w(t_0)=0$ ，则必有 $w(t)\equiv 0,\forall t\in I$ . 因为在必要性的证明中，我们实际上已证得：只要 $w(t_0)=0$ ，则 $x_1(t),\cdots ,x_n(t)$ 在 $I$ 上线性相关，若同时还存在 $t_1\in I$ 使 $w(t_1)\ne 0$ ，由充分性知 $x_1(t),\cdots ,x_n(t)$ 在 $I$ 上线性无关，这将产生矛盾。但要特别指出，这一结论仅当 $x_1(t),\cdots ,x_n(t)$ 为线性齐次方程的解时才正确。例如对于函数组

$$x_1=\{\begin{array}{ll}{t}^2,& t⩾0,\\ 0,& t<0,\end{array}{x}_2=\{\begin{array}{ll}0,& t⩾0,\\ t^2,& t<0,\end{array}$$

尽管

$$w(t)=\{\begin{array}{ll}\left|\begin{array}{cc}{t}^2& 0\\ 2t& 0\end{array}\right|=0,& t⩾0,\\ \left|\begin{array}{cc}0& t^2\\ 0& 2t\end{array}\right|=0,& t<0,\end{array}$$

即 $w(t)\equiv 0,\forall t\in \mathbb{R}$ , 但不难证明 (留作练习) $x_1(t)$ 与 $x_2(t)$ 在 $(-\infty ,+\infty )$ 却是线性无关的.

定理 2.4 线性齐次方程的通解结构

若 $x_1,x_2,\cdots ,x_n$ 是 $n$ 阶线性齐次方程 (2.7) 的 $n$ 个线性无关的特解，则它的任一解 $x$ 均可表示为

$$x=C_1{x}_1+C_2{x}_2+\cdots +C_n{x}_n,\text{\hspace{1em}(2.19)}$$

其中 $C_1,C_2,\cdots ,C_n$ 为任意常数.

证明

首先,由解的叠合性(定理2.1)知表示式(2.19)中的$x$为方程(2.7)的解.

其次,我们证明它包括了方程(2.7)的所有解.设方程(2.7)满足任一初值条件:

$$x(t_0)=x_0,\dot{x}(t_0)={\dot{x}}_0,\dots ,x^{(n-1)}(t_0)=x_0^{(n-1)},\text{\hspace{1em}(2.20)}$$

(其中 $x_0,{\dot{x}}_0,\cdots ,x_0^{(n-1)}$ 为任一组常数) 的解为 $x(t)$ . 我们要证, 可适当选取常数

$C_1,\cdots ,C_n$ 使此解 $x(t)$ 可表示为(2.19)的形式. 对解(2.19)求直至 n-1 阶导数, 并将初值条件代入得:

$$\{\begin{array}{l}{x}_0=C_1{x}_1(t_0)+C_2{x}_2(t_0)+\cdots +C_n{x}_n(t_0),\\ {\dot{x}}_0=C_1{\dot{x}}_1(t_0)+C_2{\dot{x}}_2(t_0)+\cdots +C_n{\dot{x}}_n(t_0),\\ \dots \dots \\ x_0^{(n-1)}=C_1{x}_1^{(n-1)}(t_0)+C_2{x}_2^{(n-1)}(t_0)+\cdots +C_n{x}_n^{(n-1)}(t_0).\end{array}\text{\hspace{1em}(2.21)}$$

这是一个以 $C_1,\cdots ,C_n$ 为未知数的线性代数方程组，它的系数行列式为 Wronski 行列式 $w(t_0)$ . 由于 $x_1,\cdots ,x_n$ 是 (2.7) 的线性无关的解，由定理2.3 及其注意:由定理的结论可见,形如(2.19)的解中包含了$n$个独立的任意常数,所以它就是方程(2.7)的通解.该定理还表明:线性齐次方程的通解包括了它的全部解.所以,我们也可用包括全部解的表达式去定义齐次线性方程的通解.不仅如此,后面我们即将看到,非齐次线性方程的通解也包括它的全部解.这是线性方程所特有的性质.

脚注知 $w(t_0)\ne 0$ ，从而由 Cramer 法则知，方程组(2.21)存在唯一一组解 $C_1^0,C_2^0,\cdots ,C_n^0$ . 于是

$$x=C_1^0{x}_1+C_2^0{x}_2+\cdots +C_n^0{x}_n\text{\hspace{1em}(2.22)}$$

便是方程(2.7)满足初值条件(2.20)的解.

2. 线性非齐次微分方程解的结构

定理 2.5 线性非齐次方程解的结构

设 $\bar{x}$ 是 $n$ 阶线性非齐次方程(2.6)的任一特解， $X=C_1{x}_1+C_2{x}_2+\cdots +C_n{x}_n$ 是(2.6)所对应的齐次方程(2.7)的通解，则方程(2.6)的通解为

$$x=X+\bar{x},\text{\hspace{1em}(2.23)}$$

而且通解(2.23)包括了此方程的全部解.

我们看到,与一阶线性非齐次方程解的结构相同,高阶线性非齐次方程的通解等于它的任一特解与其对应的齐次方程的通解之和.

证明

首先,证明由(2.23)所表示的 $x$ 是方程(2.6)的解.将表达式(2.23)代入方程(2.6),注意到 $X$ 与 $\bar{x}$ 分别为方程(2.7)与(2.6)的解,得

$$L(x)=L(X+\bar{x})=L(X)+L(\bar{x})=0+F(t)=F(t),$$

故 $x$ 为方程(2.6)的解.

由于表达式(2.23)中含有 $n$ 个独立的任意常数,故它是方程(2.6)的通解.

其次,证明非齐次方程(2.6)包括了它的全部解.即方程(2.6)的任一解 $\tilde{x}$ 均可由(2.23)式表示.换句话说,要证明 $\tilde{x}$ 可从(2.23)式中适当选取$X$中的常数 $C_1,\cdots ,C_n$ 获得.由于 $\tilde{x}$ 与 $\bar{x}$ 均为非齐次方程(2.6)的解,故有

$$L(\tilde{x}-\bar{x})=L(\tilde{x})-L(\bar{x})=F-F\equiv 0,$$

从而 $\tilde{x}-\overline{x}$ 是齐次方程(2.7)的一个解.由定理2.4可知,它可以由其通解$X$中适当选取常数 $C_1,\cdots ,C_n$ 得到,即

$$\tilde{x}-\overline{x}=X,\text{或}\tilde{x}=\overline{x}+X.\text{I}$$

定理 2.6

若 $x_1$ 与 $x_2$ 分别为线性非齐次方程

$$L(x)=F_1\text{与}L(x)=F_2$$

的解,则 $x_1+x_2$ 必为方程

$$L(x)=F_1+F_2\text{\hspace{1em}(2.24)}$$

的解.

证明

将 $x_1+x_2$ 代入方程(2.24)，得

$$L(x_1+x_2)=L(x_1)+L(x_2)=F_1+F_2.$$

这就表示 $x_1+x_2$ 满足方程(2.24).

应当指出,对于线性微分方程,尽管解的结构已十分清楚,有关理论也颇为完整,但是要把解具体求出来,一般来说是很困难的,这是因为要求得 $n(n⩾2)$ 阶线性齐次方程的那些线性无关的特解是难以做到的.即使对于二阶线性齐次方程也只能对某些特殊情况,才能求出其通解.然而,可喜的是,对于自然科学和工程技术中经常碰到的常系数线性微分方程,我们找到了一种非常简便的方法,只需要通过代数运算就能求出它的通解.

数字资源 4.2.1·已知二阶齐次线性方程的一个特解求另一线性无关特解的 Liouville 公式

查看资源

2.3 高阶常系数线性齐次微分方程的解法

当方程(2.6)中的系数 $P_1,\cdots ,P_n$ 均为常数时,此方程称为常系数线性微分方程,简称为常系数线性方程,常记为

$$x^{(n)}+a_1{x}^{(n-1)}+a_2{x}^{(n-2)}+\cdots +a_{n}x=f(t),$$

其中 $a_1,a_2,\cdots ,a_n$ 均为常数. 本段我们将根据 2.2 段中所述的线性方程的一般理论来讨论常系数线性齐次方程的求解法. 为简洁起见, 我们仅以二阶方程为例来进行讲解, 所述方法可类似地推广至 $n$ 阶方程.

常系数二阶线性齐次方程的一般形式为

数字资源 4.2.2·配套数字化资源

查看资源

求高阶非齐次线性方程特解的常数变易法.

$$\ddot{x}+a_1\dot{x}+a_{2}x=0,$$

(2.25)

其中 $a_1,a_2$ 为常数. 我们欲用待定系数法来求解, 即选用适当的函数 $x(t)$ 代入方程 (2.25), 确定该函数中的系数使其满足此方程.

由于指数函数求导后仍为指数函数,利用这一性质,可设方程(2.25)的解为 x=

$e^{\lambda t}$ ( $\lambda$ 为待定的实的或复的常数).于是 $\dot{x}=\lambda e^{\lambda t\stackrel{◯}{1}},\ddot{x}={\lambda }^2{e}^{\lambda t}$ ，代入方程(2.25)得

$${\mathrm{e}}^{\lambda t}\left({\lambda }^2+a_1\lambda +a_2\right)=0.$$

由于 ${\mathrm{e}}^{\lambda t}\ne 0$ ，故有

$${\lambda }^2+a_1\lambda +a_2=0.\text{\hspace{1em}(2.26)}$$

显然,对于二次代数方程(2.26)的每一个根 $\lambda$ ,就有微分方程(2.25)的一个解 $e^{\lambda t}$ .代数方程(2.26)称为线性齐次微分方程(2.25)的特征方程.它的根称为特征根或特征值.下面根据(2.26)式根的不同情况分别讨论.

对比(2.26)与(2.25)两式可见,求方程(2.25)的特征方程只需将(2.25)中的$x$换成 $\lambda$ ,将$x$的 $i(i=1,2)$ 阶导数换成 $\lambda$ 的$i$次方即可.

（1）设特征方程(2.26)有两个不相等的实根 ${\lambda }_1$ 与 ${\lambda }_2$ . 这时, 所得到的微分方程 (2.25) 的两个特解为

$$x_1={\mathrm{e}}^{{\lambda }_{1}t},x_2={\mathrm{e}}^{{\lambda }_{2}t}.$$

由于它们之比

$$\frac{{\mathrm{e}}^{{\lambda }_{1}t}}{{\mathrm{e}}^{{\lambda }_{2}t}}={\mathrm{e}}^{({\lambda }_1-{\lambda }_2)t}\ne C(\text{常数}),$$

所以,这两个解线性无关.于是由定理2.4知,方程(2.25)的通解为

$$x=C_1{\mathrm{e}}^{{\lambda }_{1}t}+C_2{\mathrm{e}}^{{\lambda }_{2}t},$$

其中 $C_1,C_2$ 为任意常数.

（2）设特征方程(2.26)有重根.这时有 ${\lambda }_1={\lambda }_2=-\frac{1}{2}{a}_1$ . 从而，只能得到微分方程(2.25)的一个特解 $x_1={\mathrm{e}}^{{\lambda }_{1}t}$ . 我们用待定函数法来求与其线性无关的另一个特解 $x_2(t)$ . 由于 $x_2$ 与 $x_1$ 必须线性无关，故 $x_2$ 与 $x_1$ 之比不为常数，从而必为一个 $t$ 的函数，设其为 $h(t)$ . 于是

$$\frac{x_2}{x_1}=h(t)\text{或}{x}_2=h(t)x_1.$$

由于

$${\dot{x}}_2=\dot{h}{x}_1+h{\dot{x}}_1,{\ddot{x}}_2=\ddot{h}{x}_1+2\dot{h}{\dot{x}}_1+h{\ddot{x}}_1,$$

代入方程(2.25)，注意到 $L(x_1)={\ddot{x}}_1+a_1{\dot{x}}_1+a_2{x}_1=0$ ，可得

$$\left(2{\dot{x}}_1+a_1{x}_1\right)\dot{h}+x_1\ddot{h}=0.$$

将 $x_1={\mathrm{e}}^{{\lambda }_{1}t}$ 代入上式并消去因式 ${\mathrm{e}}^{{\lambda }_{1}t}$ ，得

$$\left(2{\lambda }_1+a_1\right)\dot{h}+\ddot{h}=0.$$

由 ${\lambda }_1$ 是方程的重根可知 $2{\lambda }_1+a_1=0$ ，故有

$$\ddot{h}=0,$$

从而

$$h=C_{1}t+C_2.$$

取 $C_1=1,C_2=0$ 便得到方程(2.25)的另一个与 $e^{{\lambda }_{1}t}$ 线性无关的特解

$$x_2=t{\mathrm{e}}^{{\lambda }_{1}t}.$$

在 $h=C_{1}t+C_2$ 中, 只要 $C_1\ne 0$ , 可以任选 $C_1$ 与 $C_2$ 的值均可得另一线性无关特解, 这里取 $C_2=0,C_1=1$ 可使特解最为简单.

因此,齐次方程(2.25)的通解为

$$x=C_1{x}_1+C_2{x}_2=\left(C_1+C_{2}t\right){\mathrm{e}}^{{\lambda }_{1}t},$$

其中 $C_1,C_2$ 为任意常数.

（3）设特征方程(2.26)有一对共轭复根 ${\lambda }_1=\alpha +i\beta ,{\lambda }_2=\alpha -i\beta .$ 此时，齐次方程(2.25)有两个特解

$$x_1={\mathrm{e}}^{(\alpha +\mathrm{i}\beta )t},x_2={\mathrm{e}}^{(\alpha -\mathrm{i}\beta )t}.$$

它们显然线性无关.但是这种复数形式的解不便使用,为了得到便于使用的实数形式的解,我们利用 Euler 公式(参见附录 5),可知

$$x_1={\mathrm{e}}^{\alpha t}(\cos \beta t+\mathrm{i}\sin \beta t),x_2={\mathrm{e}}^{\alpha t}(\cos \beta t-\mathrm{i}\sin \beta t).$$

由解的叠合性可知

$$\frac{1}{2}(x_1+x_2)={\mathrm{e}}^{\alpha t}\cos \beta t,\frac{1}{2\mathrm{i}}(x_1-x_2)={\mathrm{e}}^{\alpha t}\sin \beta t$$

均仍为方程(2.25)的解.它们显然线性无关,故此时齐次方程(2.25)实值形式的通解为

$$x={\mathrm{e}}^{\alpha t}(C_1\cos \beta t+C_2\sin \beta t),$$

其中 $C_1,C_2$ 为任意常数.

综上所述,求二阶常系数线性齐次方程

$$\ddot{x}+a_1\dot{x}+a_{2}x=0$$

通解的步骤可归纳如下：

第一步,根据微分方程写出它的特征方程 ${\lambda }^2+a_1\lambda +a_2=0$ ;

第二步,求解特征方程得两个特征根 ${\lambda }_1$ 与 ${\lambda }_2$ ;

第三步,根据特征根的不同情况参照下表写出微分方程的通解:

特征根 ${\lambda }_1,{\lambda }_2$	二阶常系数齐次方程的通解
两个不同实根 ${\lambda }_1,{\lambda }_2$	$x=C_1{e}^{{\lambda }_{1}t}+C_2{e}^{{\lambda }_{2}t}$
两个相等的实根 ${\lambda }_1={\lambda }_2$	$x=e^{{\lambda }_{1}t}(C_{1}t+C_2)$
一对共轭复根 ${\lambda }_{1,2}=\alpha \pm i\beta$	$x=e^{\alpha t}(C_1\cos \beta t+C_2\sin \beta t)$

例 2.4

求方程 $\ddot{x}+7\dot{x}+12x=0$ 的通解.

解

所给方程的特征方程为

$${\lambda }^2+7\lambda +12=0,$$

其特征根为 ${\lambda }_1=-3,{\lambda }_2=-4$ . 于是原方程的通解为

$$x=C_1{\mathrm{e}}^{-3t}+C_2{\mathrm{e}}^{-4t}.\text{ I }$$

例 2.5

求方程 $\ddot{x}-12\dot{x}+36x=0$ 的通解与满足初值条件 $x(0)=1,\dot{x}(0)=0$ 的特解.

解

其特征方程为

$${\lambda }^2-12\lambda +36=0,$$

解之得

$${\lambda }_1={\lambda }_2=6.$$

所以,通解为

$$x={\mathrm{e}}^{6t}(C_1+C_{2}t).$$

将初值条件代入,得

$$\{\begin{array}{l}1={\mathrm{e}}^0(C_1+C_2\cdot 0)=C_1,\\ 0=6{\mathrm{e}}^0(C_1+C_2\cdot 0)+C_2{\mathrm{e}}^0=6C_1+C_2,\end{array}$$

解得 $C_1=1,C_2=-6$ . 从而所求特解为

$$x={\mathrm{e}}^{6t}(1-6t).\text{ I }$$

例 2.6

求方程 $\ddot{x}+2\dot{x}+5x=0$ 的通解.

解

其特征方程为

$${\lambda }^2+2\lambda +5=0,$$

解之得, ${\lambda }_1=-1+2i,{\lambda }_2=-1-2i.$ 于是通解为

$$x={\mathrm{e}}^{-t}(C_1\cos 2t+C_2\sin 2t).$$

例 2.7

求方程 $x^{(4)}-2\ddot{x}-3\ddot{x}+8\dot{x}-4x=0$ 的通解.

解

其特征方程为

$${\lambda }^4-2{\lambda }^3-3{\lambda }^2+8\lambda -4=0,$$

解之可得特征根为

$${\lambda }_1=2,{\lambda }_2=-2,{\lambda }_3=1,{\lambda }_4=1.$$

对应于 ${\lambda }_1$ 与 ${\lambda }_2$ 可得原方程的两个特解为 $e^{2t},e^{-2t}$ ; 对应于重根 ${\lambda }_3={\lambda }_4=1$ , 利用情形 (2) 中所用的待定函数法, 可得两个特解 $e^t,t{e}^t$ . 解组 $\{e^{2t},e^{-2t},e^t,t{e}^t\}$ 的 Wronski 行列式为

$$w(t)=\left|\begin{array}{cccc}{\mathrm{e}}^{2t}& {\mathrm{e}}^{-2t}& {\mathrm{e}}^t& t{\mathrm{e}}^t\\ 2{\mathrm{e}}^{2t}& -2{\mathrm{e}}^{-2t}& {\mathrm{e}}^t& {\mathrm{e}}^t(1+t)\\ 4{\mathrm{e}}^{2t}& 4{\mathrm{e}}^{-2t}& {\mathrm{e}}^t& {\mathrm{e}}^t(2+t)\\ 8{\mathrm{e}}^{2t}& -8{\mathrm{e}}^{-2t}& {\mathrm{e}}^t& {\mathrm{e}}^t(3+t)\end{array}\right|,$$

易见

$$\begin{array}{rl}w(0)& =\left|\begin{array}{cccc}1& 1& 1& 0\\ 2& -2& 1& 1\\ 4& 4& 1& 2\\ 8& -8& 1& 3\end{array}\right|=\left|\begin{array}{cccc}1& 0& 0& 0\\ 2& -4& -1& 1\\ 4& 0& -3& 2\\ 8& -16& -7& 3\end{array}\right|=\left|\begin{array}{ccc}-4& -1& 1\\ 0& -3& 2\\ -16& -7& 3\end{array}\right|\\ & =\left|\begin{array}{ccc}-4& -1& 1\\ 0& -3& 2\\ 0& -3& -1\end{array}\right|=-4\left|\begin{array}{cc}-3& 2\\ -3& -1\end{array}\right|=-36\ne 0.\end{array}$$

所以这四个解线性无关.从而原方程的通解为

$$x=C_1{\mathrm{e}}^{2t}+C_2{\mathrm{e}}^{-2t}+C_3{\mathrm{e}}^t+C_{4}t{\mathrm{e}}^t.\text{ I }$$

一般来说,为求解 $n$ 阶线性齐次方程

$$x^{(n)}+a_1{x}^{(n-1)}+a_2{x}^{(n-2)}+\cdots +a_{n}x=0,$$

可先写出它的特征方程

$${\lambda }^n+a_1{\lambda }^{n-1}+a_2{\lambda }^{n-2}+\cdots +a_n=0,$$

求出其特征根,然后求出与这些特征根相对应的线性无关的 $n$ 个特解,将所有特征根对应的线性无关的特解与任意常数线性组合便得线性齐次方程的通解.

为方便于我们将各种类型特征根及其在通解中的对应项列表如下：

特征根	通解中的对应项
单实根 $\lambda$	$C{e}^{\lambda t}$
$k$ 重实根 $\lambda$	$e^{\lambda t}(C_1+C_{2}t+\cdots +C_k{t}^{k-1})$
一对共轭单复根 ${\lambda }_{1,2}=\alpha \pm i\beta$	$e^{\alpha t}(C_1\cos \beta t+C_2\sin \beta t)$
一对 $k$ 重共轭复根 $\alpha \pm i\beta$	$e^{\alpha t}\left[(C_{11}+C_{12}t+\cdots +C_{1k}{t}^{k-1})\cos \beta t+(C_{21}+C_{22}t+\cdots +C_{2k}{t}^{k-1})\sin \beta t\right]$

例 2.8

求方程 $y^{(4)}(x)-4y^{\prime \prime \prime }(x)+10y^{\prime \prime }(x)-12y^{\prime }(x)+5y(x)=0$ 的通解.

解

其特征方程为

$${\lambda }^4-4{\lambda }^3+10{\lambda }^2-12\lambda +5=0.$$

全部特征根为 $1,1,1\pm 2i$ . 于是所对应的线性无关的特解为

$${\mathrm{e}}^x,x{\mathrm{e}}^x,{\mathrm{e}}^x\cos 2x,{\mathrm{e}}^x\sin 2x,$$

所以所求通解为

$$\begin{array}{rl}y& =(C_1+C_{2}x){\mathrm{e}}^x+{\mathrm{e}}^x(C_3\cos 2x+C_4\sin 2x)\\ & ={\mathrm{e}}^x(C_1+C_{2}x+C_3\cos 2x+C_4\sin 2x).\end{array}$$

2.4 高阶常系数线性非齐次微分方程的解法

我们由定理2.5 已经知道, 线性非齐次方程的通解等于它的任一特解与对应齐次方程通解之和. 齐次方程通解的求法已由上段解决, 因此求线性非齐次方程通解的关键在于寻求它的一个特解. 对于常系数方程而言, 当非齐次项 $F(t)$ 是某些特殊类型的函数时, 其特解可用待定系数法获得. 为简洁起见, 我们仍以二阶方程为例来讲解.

二阶常系数线性非齐次方程的一般形式为

$$\ddot{x}+a_1\dot{x}+a_{2}x=F(t),\text{\hspace{1em}(2.27)}$$

下面对 $F(t)$ 的几种常见的特殊形式分别进行讨论.

1. $F(t)=\phi (t)e^{\mu t}$

其中 $\mu$ 为常数, $\phi (t)$ 是一个 $m(m⩾0)$ 次多项式, 即

$$\phi (t)=b_m{t}^m+b_{m-1}{t}^{m-1}+\cdots +b_{1}t+b_0.$$

要寻求方程(2.27)的一个特解,就是要找一个函数 $x^{\ast }(t)$ ,使它满足方程(2.27).由于方程(2.27)的右端是一多项式与 $e^{\mu t}$ 的乘积,所以 $x^{\ast }(t)$ 中应含有 $e^{\mu t}$ ;又由于多项式与指数函数乘积的导数仍是多项式与指数函数的乘积,因此,我们想到应令

$$x^{\ast }(t)=Z(t){\mathrm{e}}^{\mu t},\text{\hspace{1em}(2.28)}$$

其中 $Z(t)$ 是一个待定的多项式. 将 (2.28) 式求导后代入方程 (2.27) 消去 $e^{\mu t}$ 后可得

$$\left({\mu }^2+a_1\mu +a_2\right)Z(t)+\left(2\mu +a_1\right)Z^{\prime }(t)+Z^{\prime \prime }(t)=\phi (t).\text{\hspace{1em}(2.29)}$$

如果(2.28)式是方程(2.27)的解,那么(2.29)式应是一个恒等式.因此,应选取 $Z(t)$ 的次数,使(2.29)式左端多项式的次数与右端多项式的次数相等,然后比较两端$t$的同次幂的系数,从而确定待定多项式 $Z(t)$ .以下分三种情形讨论:

(1) $\mu$ 不是特征方程的根. 这时 ${\mu }^2+a_1\mu +a_2\ne 0$ , 于是可令 $Z(t)$ 是一个与 $\phi (t)$ 同次数的多项式

$$Z(t)=B_m{t}^m+B_{m-1}{t}^{m-1}+\cdots +B_{1}t+B_0.\text{\hspace{1em}(2.30)}$$

把(2.30)式代入(2.29)式并比较(2.29)式两端$t$的同次幂系数,定出 $B_i$ (i=0,1,…,$m$)后代入(2.28)式,便得到常系数线性非齐次微分方程(2.27)的一个特解.

(2) $\mu$ 是特征方程的单根.这时有

$${\mu }^2+a_1\mu +a_2=0,2\mu +a_1\ne 0.$$

要使(2.29)式两端多项式的次数相等, $Z(t)$ 应是一个 $m+1$ 次的多项式,因此可令

$$Z(t)=t\left(B_m{t}^m+B_{m-1}{t}^{m-1}+\cdots +B_{1}t+B_0\right).$$

(3) $\mu$ 是特征方程的重根. 这时有

$${\mu }^2+a_1\mu +a_2=0,2\mu +a_1=0,$$

按照上面同样的想法,可令

$$Z(t)=t^2\left(B_m{t}^m+B_{m-1}{t}^{m-1}+\cdots +B_{1}t+B_0\right).$$

综上所述,为求二阶常系数线性非齐次方程

$$\ddot{x}+a_1\dot{x}+a_{2}x=\phi (t){\mathrm{e}}^{\mu t}$$

的一个特解 $x^{\ast }(t)$ ，可令

$$x^{\ast }(t)=t^k{\mathrm{e}}^{\mu t}\Phi (t),$$

其中 $\Phi (t)=B_m{t}^m+B_{m-1}{t}^{m-1}+\cdots +B_{1}t+B_0$ 是与 $\phi (t)$ 同次数的多项式. 当 $\mu$ 不是特征根时, 取 k=0; 当 $\mu$ 是单重特征根时, 取 k=1; 当 $\mu$ 是重特征根时, 取 k=2.

例 2.9

求微分方程 $\ddot{x}-5\dot{x}+6x=t{e}^{2t}$ 的通解.

解

对应于齐次线性方程的特征方程为

$${\lambda }^2-5\lambda +6=0,$$

故特征根为 ${\lambda }_1=2,{\lambda }_2=3$ ，从而齐次方程的通解为

$$x=C_1{\mathrm{e}}^{2t}+C_2{\mathrm{e}}^{3t}.$$

现求非齐次线性方程的一个特解,由于 $\mu =2$ 是特征方程的单根,故应令

$$x^{\ast }=t(B_0+B_{1}t){\mathrm{e}}^{2t},$$

求导后代入原方程,化简得

$$-2B_{1}t+2B_1-B_0\equiv t,$$

比较同次幂的系数,得

$$-2B_1=1,2B_1-B_0=0,$$

从而

$$B_0=-1,B_1=-\frac{1}{2},$$

于是,特解为

$$x^{\ast }=-\left(t+\frac{1}{2}{t}^2\right){\mathrm{e}}^{2t}.$$

所以,原方程的通解为

$$x=C_1{\mathrm{e}}^{2t}+C_2{\mathrm{e}}^{3t}-\left(t+\frac{1}{2}{t}^2\right){\mathrm{e}}^{2t}.i$$

2. $F(t)={\mathrm{e}}^{\mu t}\phi (t)\cos \nu t$ 或 $F(t)={\mathrm{e}}^{\mu t}\phi (t)\sin \nu t$

读者容易证明,若函数 $x=x_{\mathrm{R}}(t)\pm \mathrm{i}{x}_{\mathrm{I}}(t)$ 是方程

$$\ddot{x}+a_1(t)\dot{x}+a_2(t)x=f_1(t)\pm \mathrm{i}{f}_2(t)$$

的解,其中 $i$ 为虚数单位,则其实部 $x_{\mathrm{R}}(t)$ 与虚部 $x_{\mathrm{I}}(t)$ 分别是方程

$$\ddot{x}+a_1(t)\dot{x}+a_2(t)x=f_1(t)\text{与}\ddot{x}+a_1(t)\dot{x}+a_2(t)x=f_2(t)$$

的解.

因此,我们可以先求微分方程

$$\ddot{x}+a_1\dot{x}+a_{2}x={\mathrm{e}}^{\mu t}\phi (t)\cos \nu t\pm {\mathrm{ie}}^{\mu t}\phi (t)\sin \nu t={\mathrm{e}}^{(\mu \pm \mathrm{i}\nu )t}\phi (t)\text{\hspace{1em}(2.31)}$$

的特解,然后分出其实部、虚部便可得出所求微分方程的特解,这样,就把类型2化成了类型1.

例 2.10

求方程 $\ddot{x}+x=t{e}^t\cos t$ 的通解.

解

对应线性齐次方程的特征方程为

$${\lambda }^2+1=0,$$

故特征根为

$$\lambda =\pm \mathrm{i},$$

从而线性齐次方程的通解为

$$x=C_1\cos t+C_2\sin t.$$

为求原方程的一个特解,先求微分方程

$$\ddot{x}+x=t{\mathrm{e}}^{(1+\mathrm{i})t}\text{\hspace{1em}(2.32)}$$

的特解,由于 $1+i$ 不是特征根,故令

$$x^{\ast }=\left(B_0+B_{1}t\right){\mathrm{e}}^{(1+\mathrm{i})t}.$$

求导后代入相应方程(2.32)并化简得

$$B_1(1+2\mathrm{i})t+[2B_1(1+\mathrm{i})+2B_0\mathrm{i}+B_0]\equiv t.$$

比较系数得

$$B_1=\frac{1}{1+2\mathrm{i}}=\frac{1-2\mathrm{i}}{5},B_0=-\frac{2(1+\mathrm{i})}{1+2\mathrm{i}}{B}_1=\frac{-2+14\mathrm{i}}{25},$$

从而

$$x^{\ast }={\mathrm{e}}^t\left(\frac{1-2\mathrm{i}}{5}t+\frac{-2+14\mathrm{i}}{25}\right)(\cos t+\mathrm{i}\sin t).$$

它的实部

$$x_{\mathrm{R}}^{\ast }={\mathrm{e}}^t\left[\left(\frac{1}{5}t-\frac{2}{25}\right)\cos t+\left(\frac{2}{5}t-\frac{14}{25}\right)\sin t\right]=\frac{{\mathrm{e}}^t}{25}[(5t-2)\cos t+(10t-14)\sin t]$$

就是原方程的一个特解,于是原方程的通解为

$$x(t)=C_1\cos t+C_2\sin t+\frac{{\mathrm{e}}^t}{25}[(5t-2)\cos t+(10t-14)\sin t].$$

由此例可见,当非齐次项如类型2所示时,用上述方法求特解时要进行复数运算,颇不方便,下面我们由此导出另一种待定系数法,它可以仅在实数域内作运算.

实数域内的待定系数法.由例2.10我们看到为求方程(2.31)的特解,可令

$$x^{\ast }=t^{k}Z(t){\mathrm{e}}^{(\mu \pm \mathrm{i}\nu )t},$$

其中当 $\mu \pm i\nu$ 不是特征根时, 取 k=0; 当 $\mu +i\nu$ 是特征根时, 取 k=1. 通过代入方程, 比较系数可定出与 $\phi (t)$ 同次数的复值多项式

$$Z(t)=R(t)\pm \mathrm{i}I(t).$$

于是

$$\begin{array}{rl}{x}^{\ast }& =t^k[R(t)+\mathrm{i}I(t)]{\mathrm{e}}^{\mu t}(\cos \nu t+\mathrm{i}\sin \nu t)\\ & =t^k{\mathrm{e}}^{\mu t}\left\{\left[R(t)\cos \nu t-I(t)\sin \nu t\right]+\mathrm{i}[I(t)\cos \nu t+R(t)\sin \nu t]\right\},\end{array}$$

从而其实部与虚部分别为

$$x_{\mathrm{R}}^{\ast }=t^k{\mathrm{e}}^{\mu t}[R(t)\cos \nu t-I(t)\sin \nu t]\text{\hspace{1em}(2.33)}$$

与

$$x_{\mathrm{I}}^{\ast }=t^k{\mathrm{e}}^{\mu t}[I(t)\cos \nu t+R(t)\sin \nu t],\text{\hspace{1em}(2.34)}$$

它们分别是方程

$$\ddot{x}+a_1\dot{x}+a_{2}x={\mathrm{e}}^{\mu t}\phi (t)\cos \nu t\text{与}\ddot{x}+a_1\dot{x}+a_{2}x={\mathrm{e}}^{\mu t}\phi (t)\sin \nu t\text{\hspace{1em}(2.35)}$$

的特解.

(2.33)与(2.34)右端方括号内的表达式都是 $\cos \nu t$ 和 $\sin \nu t$ 分别乘与 $\phi (t)$ 同次数的实值多项式后的代数和.由此可见,对于形如(2.35)的线性非齐次微分方程,我们可以直接令其特解为

$$x^{\ast }=t^k{\mathrm{e}}^{\mu t}\left[Z_1(t)\cos \nu t+Z_2(t)\sin \nu t\right],$$

其中 $Z_1(t)$ 与 $Z_2(t)$ 均为与 $\phi (t)$ 同次数的实系数多项式. 当 $\mu \pm i\nu$ 不是特征根时, 取 k=0; 当 $\mu +i\nu$ 是特征根时, 取 k=1.

例 2.11

求方程 $\ddot{x}-3\dot{x}+2x=t\cos t$ 满足初值条件 $x(0)=\frac{22}{25},\dot{x}(0)=\frac{19}{25}$ 的特解.

解

容易求得对应线性齐次方程的通解为

$$x(t)=C_1{\mathrm{e}}^t+C_2{\mathrm{e}}^{2t}.$$

由于 $\mu +i\nu =i$ 不是特征根, 故令原方程的一个特解为

$$x^{\ast }=(A_0+A_{1}t)\cos t+(B_0+B_{1}t)\sin t.$$

代入原方程后整理得

$$\left(A_0-3A_1-3B_0+2B_1\right)\cos t+\left(3A_0+B_0-2A_1-3B_1\right)\sin t+$$$$\left(A_1-3B_1\right)t\cos t+\left(3A_1+B_1\right)t\sin t\equiv t\cos t,$$

比较两端同类项系数,得

$$A_0-3A_1-3B_0+2B_1=0,3A_0+B_0-2A_1-3B_1=0,$$$$A_1-3B_1=1,3A_1+B_1=0.$$

从而解得

$$A_1=\frac{1}{10},B_1=-\frac{3}{10},A_0=-\frac{3}{25},B_0=-\frac{17}{50},$$

所以得

$$x^{\ast }=\left(-\frac{3}{25}+\frac{1}{10}t\right)\cos t-\left(\frac{17}{50}+\frac{3}{10}t\right)\sin t,$$

故原方程的通解为

$$x=C_1{\mathrm{e}}^t+C_2{\mathrm{e}}^{2t}+\left(-\frac{3}{25}+\frac{1}{10}t\right)\cos t-\left(\frac{17}{50}+\frac{3}{10}t\right)\sin t.$$

为求适合初值条件的特解,先将上式求导,得

$$\dot{x}=C_1{\mathrm{e}}^t+2C_2{\mathrm{e}}^{2t}-\frac{12}{50}\cos t-\frac{9}{50}\sin t-\frac{3}{10}t\cos t-\frac{1}{10}t\sin t.$$

将初值条件代入上面两式, 得

$$\int C_1+C_2-\frac{3}{25}=x(0)=\frac{22}{25},$$$$[C_1+2C_2-\frac{6}{25}=\dot{x}(0)=\frac{19}{25},$$

解得 $C_1=1,C_2=0.$

于是，所求特解为

$$x={\mathrm{e}}^t+\left(-\frac{3}{25}+\frac{1}{10}t\right)\cos t-\left(\frac{17}{50}+\frac{3}{10}t\right)\sin t.$$

例 2.12

求微分方程 $\ddot{x}+3x=\sin 2t$ 的一个特解.

解

我们可以利用例2.11 的方法来求解, 然而, 这个方程不含 $\dot{x}$ 项, 由于正弦函数的二阶导数仍是正弦函数, 容易看出, 可令特解为

$$x^{\ast }=A\sin 2t.$$

代入原方程,得

$$(-4A+3A)\sin 2t\equiv \sin 2t,$$

由此得

$$A=-1,$$

从而原方程的一个特解为

$$x^{\ast }=-\sin 2t.$$

上述求特解的待定系数法,可以推广到一般的 $n$ 阶常系数线性微分方程

$$x^{(n)}+a_1{x}^{(n-1)}+a_2{x}^{(n-2)}+\cdots +a_{n-1}\dot{x}+a_{n}x=F(t),$$

其中 $F(t)$ 是类型 1(或类型 2) 中所示的函数. 若 $\mu (\text{或 }\mu +\mathrm{i}\nu )$ 是对应特征方程的 $k$ 重根, $k=0,1,\cdots ,n$ （或 $k=0,1,\cdots ,\left[\frac{n}{2}\right]$ ），则可令特解为

$$x^{\ast }=t^{k}Z(t){\mathrm{e}}^{\mu t}(\text{或}{x}^{\ast }=t^k{\mathrm{e}}^{\mu t}[Z_1(t)\cos \nu t+Z_2(t)\sin \nu t]),$$

其中 $Z(t)$ ( $Z_1(t)$ 与 $Z_2(t)$ ) 均是与 $\phi (t)$ 同次数的待定多项式.

为方便读者使用,现就 $F(t)$ 的常见类型,针对不同情况将用待定系数法应设置的特解的形式列表如下(其中 $Z(t)$ , $Z_1(t)$ , $Z_2(t)$ 均为与 $\phi$ 同次数的多项式):

$F(t)$ 的类型	应设置特解 $x^{\ast }(t)$ 的形式
$m$ 次多项式 $\phi (t)$	0不是特征值	$x^{\ast }=Z(t)$
	0是 $k$ 重特征值	$x^{\ast }=t^{k}Z(t)$
$\phi (t)e^{\mu t}$	$\mu$ 不是特征值	$x^{\ast }=Z(t)e^{\mu t}$
	$\mu$ 是 $k$ 重特征值	$x^{\ast }=t^{k}Z(t)e^{\mu t}$
$\phi (t)e^{\mu t}\cos \nu t$ 或 $\phi (t)e^{\mu t}\sin \nu t$	$\mu +i\nu$ 不是特征值	$x^{\ast }=e^{\mu t}[Z_1(t)\cos \nu t+Z_2(t)\sin \nu t]$
	$\mu +i\nu$ 是 $k$ 重特征值 $\left(1⩽k⩽\left[\frac{n}{2}\right]\right)$	$x^{\ast }=t^k{e}^{\mu t}[Z_1(t)\cos \nu t+Z_2(t)\sin \nu t]$

例 2.13

求微分方程 $y^{(6)}+y^{(5)}-2y^{(4)}=x-1$ 的通解.

解

对应于线性齐次方程的特征方程为

$${\lambda }^6+{\lambda }^5-2{\lambda }^4={\lambda }^4(\lambda -1)(\lambda +2)=0,$$

故特征根为 ${\lambda }_1=0$ （4 重根）， ${\lambda }_2=1$ （单根）， ${\lambda }_3=-2$ （单根）。于是齐次线性方程的通解为

$$y=C_1{\mathrm{e}}^x+C_2{\mathrm{e}}^{-2x}+C_3{x}^3+C_4{x}^2+C_{5}x+C_6.$$

由于 $\mu =0$ 是 4 重特征根, 故令线性非齐次方程的特解为

$$y^{\ast }=x^4(Ax+B),$$

代入原方程,得

$$120A-240Ax-48B\equiv x-1.$$

比较同次幂系数,得

$$-240A=1,120A-48B=-1,$$

从而

$$A=-\frac{1}{240},B=\frac{1}{96}.$$

于是原方程的通解为

$$y=C_1{\mathrm{e}}^x+C_2{\mathrm{e}}^{-2x}+C_3{x}^3+C_4{x}^2+C_{5}x+C_6+\frac{1}{96}{x}^4-\frac{1}{240}{x}^5.I$$

例 2.14

求解在本章例2.1 中所建立的关于弹簧机械振动的微分方程并加以讨论.

解

在例2.1中我们所建立的微分方程为

$$\ddot{x}+2\delta \dot{x}+{\omega }^{2}x=h\sin pt,\text{\hspace{1em}(2.36)}$$

其中 $\delta =\frac{\mu }{2m},\omega =\sqrt{\frac{k}{m}},h=\frac{H}{m},\mu$ 为介质的阻尼系数，$m$ 为物体的质量，$k$ 为弹簧的劲度系数， $H\sin pt$ 为周期变化的外界强迫力。(2.36) 是一个二阶常系数线性非齐次微分方程，称为强迫振动的微分方程。它所对应的齐次方程是被称为自由振动的微分方程

$$\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{t}^2}+2\delta \frac{\mathrm{d}x}{\mathrm{d}t}+{\omega }^{2}x=0.\text{\hspace{1em}(2.37)}$$

首先,我们来求线性齐次方程(2.37)的通解,它的特征方程为

$${\lambda }^2+2\delta \lambda +{\omega }^2=0,$$

故特征根为

$${\lambda }_1=-\delta +\sqrt{{\delta }^2-{\omega }^2},{\lambda }_2=-\delta -\sqrt{{\delta }^2-{\omega }^2}.$$

下面分三种情况讨论：

(1) ${\delta }^2-{\omega }^2>0.$ 这时, 线性齐次方程的通解为

$$x=C_1{\mathrm{e}}^{-(\delta -\sqrt{{\delta }^2-{\omega }^2})t}+C_2{\mathrm{e}}^{-(\delta +\sqrt{{\delta }^2-{\omega }^2})t}.$$

由此可见，当 $t\to +\infty$ 时， $x(t)\to 0$ 。这说明当外力为零且阻尼 $\mu$ 很大时，有 $\delta >\omega$ ，物体运动按指数函数规律迅速衰减，不会产生振动，如图 4.6所示。

(2) ${\delta }^2-{\omega }^2=0.$ 这时, 特征方程有重根 ${\lambda }_1={\lambda }_2=-\delta .$ 于是线性齐次方程的通解为

$$x=\left(C_1+C_{2}t\right){\mathrm{e}}^{-\delta t},$$

物体仍按指数规律衰减,图形与图 4.6 相似,这时

图4.6

的阻尼 $\mu$ 称为临界阻尼.

(3) ${\delta }^2-{\omega }^2<0.$ 这时, 特征根是共轭复根

$${\lambda }_{1,2}=-\delta \pm \mathrm{i}\sqrt{{\omega }^2-{\delta }^2},$$

线性齐次方程的通解为

$$x={\mathrm{e}}^{-\delta t}\left(C_1\cos \sqrt{{\omega }^2-{\delta }^2}t+C_2\sin \sqrt{{\omega }^2-{\delta }^2}t\right)=A{\mathrm{e}}^{-\delta t}\sin \left(\sqrt{{\omega }^2-{\delta }^2}t+\phi \right).$$

其中 $A=\sqrt{C_1^2+C_2^2}$ 与 $\phi =\arctan \frac{C_1}{C_2}$ 是两个任意常数.

由通解可见,一方面,运动的振幅 $A{e}^{-\delta t}$ 随时间的增大而减少,所以运动是衰减的;另一方面,由于通解中含有正弦函数,所以运动又是振荡的.因此,当阻尼 $\mu$ 很小时, $\delta <\omega$ ,物体做衰减振荡运动,如图 4.7 所示.

为求解线性非齐次方程(2.36)，关键在于求它的任意一个特解.设阻尼 $\mu \ne 0$ ，从而ip不是特征根，于是令特解为

图4.7

$$x^{\ast }=A_1\cos pt+A_2\sin pt,$$

代入方程(2.36)，比较同类项系数，得

$$\{\begin{array}{l}({\omega }^2-p^2)A_1+2\delta p{A}_2=0,\\ ({\omega }^2-p^2)A_2-2\delta p{A}_1=h.\end{array}$$

解之得

$$A_1=-\frac{2\delta ph}{{\left({\omega }^2-p^2\right)}^2+4{\delta }^2{p}^2},A_2=\frac{h\left({\omega }^2-p^2\right)}{{\left({\omega }^2-p^2\right)}^2+4{\delta }^2{p}^2}.$$

于是

$$\begin{array}{rl}{x}^{\ast }& =\frac{h({\omega }^2-p^2)}{({\omega }^2-p^2{)}^2+4{\delta }^2{p}^2}\sin pt-\frac{2\delta ph}{({\omega }^2-p^2{)}^2+4{\delta }^2{p}^2}\cos pt\\ & =B\sin (pt-\psi ),\end{array}$$

其中

$$B=\frac{h}{\sqrt{{\left({\omega }^2-p^2\right)}^2+4{\delta }^2{p}^2}},\psi =\arctan \left(\frac{2p\delta }{{\omega }^2-p^2}\right).$$

所以强迫振动方程(2.36)的通解为

$$x=A{\mathrm{e}}^{-\delta t}\sin (\sqrt{{\omega }^2-{\delta }^2}t+\phi )+B\sin (pt-\psi ),\delta <\omega .$$

容易看出,随时间 $t$ 增大,上式右端第一项很快地衰减为零,故称它为暂态项;第二项

不随时间衰减,故称它为稳态项.由此可见,弹簧的强迫振动规律,主要取决于特解,即

$$x\approx x^{\ast }=B\sin \left(pt-\psi \right),\text{当}t\gg 1\text{时}.$$

换句话说,主要取决于外界强迫力的作用,其中 $p$ 就是强迫力周期运动的频率.

当阻尼很小时,振动的振幅为

$$B=\frac{h}{\sqrt{{\left({\omega }^2-p^2\right)}^2+4{\delta }^2{p}^2}}\approx \frac{h}{\left|{\omega }^2-p^2\right|}.$$

由此可见，当外界强迫力的频率 $p$ 与弹簧的固有频率 $\omega =\sqrt{\frac{k}{m}}$ 相差不大时，振幅 $B$ 将会很大，这时就会产生所谓共振现象.

共振现象在很多问题中有很大的破坏作用,它可能引起机器损坏、桥梁折断、建筑物倒塌等严重事故.例如,1831年一队士兵以整齐的步伐通过英国Lanchester附近的布劳顿吊桥时,由于整齐的步伐产生了周期性的外力,而且这个力的频率非常接近于吊桥振动的固有频率,从而引起了共振,导致了吊桥倒塌.因此,在一些工程问题中,常需算出固有频率,调整有关参数和采取各种措施,以避免共振现象的发生.但是,事物都是一分为二的,有时人们却恰恰需要利用共振.例如,收音机、电视机必须要调节频率使之与所接收电台、电视台的发射频率相同,产生共振,才能收到需要的信息.

2.5 高阶变系数线性微分方程的求解问题

高于一阶的变系数线性微分方程,一般来说是难以用初等方法求解的.但是,对于某些特殊类型的微分方程,我们可以通过适当的变换把它化成常系数线性微分方程,从而使求解问题得到解决.下面将要介绍的后继课中颇有用处的 Euler 微分方程就是其中的一种.

Euler 微分方程 Euler 微分方程的一般形式是

$$t^n\frac{{\mathrm{d}}^{n}x}{\mathrm{d}{t}^n}+a_1{t}^{n-1}\frac{{\mathrm{d}}^{n-1}x}{\mathrm{d}{t}^{n-1}}+\cdots +a_{n-1}t\frac{\mathrm{d}x}{\mathrm{d}t}+a_{n}x=f(t),$$

其中 $a_1,a_2,\cdots ,a_n$ 都是常数. 对于这类微分方程, 通过变换

$$t={\mathrm{e}}^{\tau \stackrel{◯}{1}}\text{或}\tau =\ln t(t>0)$$

就可化为常系数线性微分方程,下面通过例子来具体说明.

例 2.15

求微分方程 $t^2\ddot{x}-t\dot{x}+x=0$ 的通解.

解

这是一个 Euler 微分方程, 令 $t=e^{\tau }$ , 即 $\tau =\ln t$ , 从而有

$$\frac{\mathrm{d}x}{\mathrm{d}t}=\frac{\mathrm{d}x}{\mathrm{d}\tau }\cdot \frac{\mathrm{d}\tau }{\mathrm{d}t}=\frac{1}{t}\frac{\mathrm{d}x}{\mathrm{d}\tau },\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{t}^2}=\frac{1}{t^2}\left(\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{\tau }^2}-\frac{\mathrm{d}x}{\mathrm{d}\tau }\right),$$

代入原方程,化简得

$$\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{\tau }^2}-2\frac{\mathrm{d}x}{\mathrm{d}\tau }+x=0.$$

这是一个二阶常系数线性齐次微分方程,容易求得它的通解为

$$x=\left(C_1\tau +C_2\right){\mathrm{e}}^{\tau },$$

把 $\tau$ 代换成 $\ln t$ ，即得原方程的通解为

$$x=\left(C_1\ln t+C_2\right)t.$$

例 2.16

求微分方程 $(x+2{)}^2\frac{{\mathrm{d}}^{3}y}{\mathrm{d}{x}^3}+(x+2)\frac{{\mathrm{d}}^{2}y}{\mathrm{d}{x}^2}+\frac{\mathrm{d}y}{\mathrm{d}x}=1$ 的通解.

解

令 $x+2=t$ ，原方程变为

$$t^2\frac{{\mathrm{d}}^{3}y}{\mathrm{d}{t}^3}+t\frac{{\mathrm{d}}^{2}y}{\mathrm{d}{t}^2}+\frac{\mathrm{d}y}{\mathrm{d}t}=1,$$

这不是 Euler 方程,但两端乘 $t$, 即得 Euler 微分方程

$$t^3\frac{{\mathrm{d}}^{3}y}{\mathrm{d}{t}^3}+t^2\frac{{\mathrm{d}}^{2}y}{\mathrm{d}{t}^2}+t\frac{\mathrm{d}y}{\mathrm{d}t}=t.$$

再令 $t=e^{\tau }$ ，上方程可化为

$$\frac{{\mathrm{d}}^{3}y}{\mathrm{d}{\tau }^3}-2\frac{{\mathrm{d}}^{2}y}{\mathrm{d}{\tau }^2}+2\frac{\mathrm{d}y}{\mathrm{d}\tau }={\mathrm{e}}^{\tau },$$

容易求得它的通解为

$$y=C_0+{\mathrm{e}}^{\tau }(C_1\cos \tau +C_2\sin \tau )+{\mathrm{e}}^{\tau }.$$

所以,原方程的通解为

$$y=C_0+(x+2)\left[C_1\cos \ln (x+2)+C_2\sin \ln (x+2)+1\right].$$

4.2 高阶线性微分方程 课后真题与自测练习

习题与真题

收录本小节对应期中、期末与精选习题· 提供推导解析、试卷溯源与 LaTeX 导出

进入自测