目录

第七章 线性微分方程组

7.1 引言

在实际问题中,常常需要同时考虑多个相互关联的未知函数,这就导致了微分方程组的研究。线性微分方程组是微分方程理论中最重要的部分之一,它不仅在数学理论上完善,在控制理论、电路分析、生态系统等领域也有广泛应用。

7.2 线性微分方程组的一般理论

7.2.1 方程组的标准形式

一阶线性微分方程组的标准形式为: $\begin{cases} x_1' = a_{11}(t)x_1 + a_{12}(t)x_2 + \cdots + a_{1n}(t)x_n + f_1(t) \\ x_2' = a_{21}(t)x_1 + a_{22}(t)x_2 + \cdots + a_{2n}(t)x_n + f_2(t) \\ \vdots \\ x_n' = a_{n1}(t)x_1 + a_{n2}(t)x_2 + \cdots + a_{nn}(t)x_n + f_n(t) \end{cases}$

其中 $x_i = x_i(t)$ 是未知函数,$a_{ij}(t)$ 和 $f_i(t)$ 是已知连续函数。

7.2.2 矩阵表示

引入向量记号: $\mathbf{x} = \begin{pmatrix} x_1 \\ x_2 \\ \vdots \\ x_n \end{pmatrix}, \quad A(t) = \begin{pmatrix} a_{11}(t) & \cdots & a_{1n}(t) \\ \vdots & \ddots & \vdots \\ a_{n1}(t) & \cdots & a_{nn}(t) \end{pmatrix}, \quad \mathbf{f}(t) = \begin{pmatrix} f_1(t) \\ f_2(t) \\ \vdots \\ f_n(t) \end{pmatrix}$

方程组可简写为: $\mathbf{x}' = A(t)\mathbf{x} + \mathbf{f}(t)$

当 $\mathbf{f}(t) \equiv \mathbf{0}$ 时,称为齐次线性方程组;否则称为非齐次线性方程组

初值条件为 $\mathbf{x}(t_0) = \mathbf{x}_0$。

例 7.1 将方程组写成矩阵形式: $\begin{cases} x_1' = 2x_1 + 3x_2 + e^t \\ x_2' = -x_1 + 4x_2 \end{cases}$

解: $\mathbf{x}' = \begin{pmatrix} 2 & 3 \\ -1 & 4 \end{pmatrix}\mathbf{x} + \begin{pmatrix} e^t \\ 0 \end{pmatrix}$

7.2.3 高阶方程与方程组的等价性

任何高阶线性微分方程都可以化为一阶线性方程组。

例 7.2 将 $y'' + p(t)y' + q(t)y = f(t)$ 化为方程组。

解: 令 $x_1 = y, x_2 = y'$,则: $\begin{cases} x_1' = x_2 \\ x_2' = -q(t)x_1 - p(t)x_2 + f(t) \end{cases}$

即 $\mathbf{x}' = \begin{pmatrix} 0 & 1 \\ -q(t) & -p(t) \end{pmatrix}\mathbf{x} + \begin{pmatrix} 0 \\ f(t) \end{pmatrix}$

7.3 齐次线性方程组的通解结构

7.3.1 叠加原理

定理 7.1 若 $\mathbf{\varphi}_1(t), \ldots, \mathbf{\varphi}_k(t)$ 是齐次方程组 $\mathbf{x}' = A(t)\mathbf{x}$ 的解,则它们的线性组合 $c_1\mathbf{\varphi}_1 + \cdots + c_k\mathbf{\varphi}_k$ 也是解。

7.3.2 基本解组与基本解矩阵

定义 7.1 齐次方程组的 $n$ 个线性无关的解称为基本解组

定义 7.2 以基本解组为列向量构成的矩阵 $\Phi(t) = \begin{pmatrix} \mathbf{\varphi}_1(t) & \mathbf{\varphi}_2(t) & \cdots & \mathbf{\varphi}_n(t) \end{pmatrix}$ 称为基本解矩阵

定理 7.2 齐次方程组的通解为: $\mathbf{x}(t) = \Phi(t)\mathbf{c} = c_1\mathbf{\varphi}_1 + c_2\mathbf{\varphi}_2 + \cdots + c_n\mathbf{\varphi}_n$ 其中 $\mathbf{c} = (c_1, \ldots, c_n)^T$ 是任意常向量。

例 7.3 验证 $\mathbf{\varphi}_1 = \begin{pmatrix} e^t \\ e^t \end{pmatrix}, \mathbf{\varphi}_2 = \begin{pmatrix} e^{2t} \\ 2e^{2t} \end{pmatrix}$ 是方程组 $\mathbf{x}' = \begin{pmatrix} 0 & 1 \\ -2 & 3 \end{pmatrix}\mathbf{x}$ 的基本解组。

解: 验证 $\mathbf{\varphi}_1' = \begin{pmatrix} e^t \\ e^t \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -2 & 3 \end{pmatrix}\begin{pmatrix} e^t \\ e^t \end{pmatrix} = \begin{pmatrix} e^t \\ e^t \end{pmatrix}$ ✓

类似验证 $\mathbf{\varphi}_2$。

计算Wronski行列式: $\det\begin{pmatrix} e^t & e^{2t} \\ e^t & 2e^{2t} \end{pmatrix} = 2e^{3t} - e^{3t} = e^{3t} \neq 0$

故线性无关,是基本解组。

7.3.3 Wronski行列式与Liouville公式

定义 7.3 设 $\mathbf{\varphi}_1, \ldots, \mathbf{\varphi}_n$ 是 $n$ 个解,称 $W(t) = \det(\mathbf{\varphi}_1, \ldots, \mathbf{\varphi}_n)$ 为它们的Wronski行列式

定理 7.3 $n$ 个解线性无关 $\Leftrightarrow W(t) \neq 0$ 对某点成立。

定理 7.4 (Liouville公式) $W(t) = W(t_0)e^{\int_{t_0}^{t} \text{tr}A(s)ds}$

其中 $\text{tr}A = a_{11} + a_{22} + \cdots + a_{nn}$ 是矩阵 $A$ 的迹。

7.4 非齐次线性方程组的通解

7.4.1 通解结构

定理 7.5 设 $\Phi(t)$ 是齐次方程组的基本解矩阵,$\mathbf{\psi}(t)$ 是非齐次方程组的特解,则非齐次方程组的通解为: $\mathbf{x}(t) = \Phi(t)\mathbf{c} + \mathbf{\psi}(t)$

7.4.2 常数变易法

设特解形式: $\mathbf{\psi}(t) = \Phi(t)\mathbf{c}(t)$

代入非齐次方程: $\Phi'\mathbf{c} + \Phi\mathbf{c}' = A\Phi\mathbf{c} + \mathbf{f}$

由于 $\Phi' = A\Phi$,得: $\Phi\mathbf{c}' = \mathbf{f}$

故: $\mathbf{c}(t) = \int \Phi^{-1}(t)\mathbf{f}(t)dt$

特解: $\mathbf{\psi}(t) = \Phi(t)\int \Phi^{-1}(t)\mathbf{f}(t)dt$

例 7.4 求解 $\mathbf{x}' = \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix}\mathbf{x} + \begin{pmatrix} 0 \\ \sin t \end{pmatrix}$。

解: 齐次方程的基本解矩阵: $\Phi(t) = \begin{pmatrix} \cos t & \sin t \\ -\sin t & \cos t \end{pmatrix}$

$\Phi^{-1}(t) = \begin{pmatrix} \cos t & -\sin t \\ \sin t & \cos t \end{pmatrix}$

$\Phi^{-1}\mathbf{f} = \begin{pmatrix} -\sin^2 t \\ \sin t \cos t \end{pmatrix}$

积分得: $\mathbf{c}(t) = \begin{pmatrix} \frac{t}{2} - \frac{\sin 2t}{4} \\ -\frac{\cos 2t}{4} \end{pmatrix}$

特解: $\mathbf{\psi}(t) = \Phi(t)\mathbf{c}(t) = \begin{pmatrix} \frac{t\cos t}{2} - \frac{\sin t}{4} \\ -\frac{t\sin t}{2} \end{pmatrix}$

通解:$\mathbf{x} = \Phi(t)\mathbf{c} + \mathbf{\psi}(t)$。

7.5 习题

习题 7.1 将下列高阶方程化为方程组:

a) $y''' - 2y'' + y' - 2y = 0$

b) $y'' + \omega^2 y = f(t)$

习题 7.2 验证 $\mathbf{\varphi}_1 = \begin{pmatrix} 1 \\ 1 \end{pmatrix}e^{3t}, \mathbf{\varphi}_2 = \begin{pmatrix} 1 \\ -1 \end{pmatrix}e^{-t}$ 是方程组 $\mathbf{x}' = \begin{pmatrix} 1 & 2 \\ 2 & 1 \end{pmatrix}\mathbf{x}$ 的基本解组。

习题 7.3 设 $\Phi(t) = \begin{pmatrix} e^{2t} & te^{2t} \\ 0 & e^{2t} \end{pmatrix}$,验证它是某齐次方程组的基本解矩阵,并求对应的 $A(t)$。

习题 7.4 利用常数变易法求解 $\mathbf{x}' = \begin{pmatrix} 2 & 1 \\ 0 & 2 \end{pmatrix}\mathbf{x} + \begin{pmatrix} e^{2t} \\ 0 \end{pmatrix}$。

习题 7.5 证明Liouville公式:若 $\Phi(t)$ 是 $\mathbf{x}' = A(t)\mathbf{x}$ 的基本解矩阵,则 $\det\Phi(t) = \det\Phi(t_0)e^{\int_{t_0}^t \text{tr}A(s)ds}$。

7.6 参考答案

习题 7.1

a) $\mathbf{x}' = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 2 & -1 & 2 \end{pmatrix}\mathbf{x}$

b) $\mathbf{x}' = \begin{pmatrix} 0 & 1 \\ -\omega^2 & 0 \end{pmatrix}\mathbf{x} + \begin{pmatrix} 0 \\ f(t) \end{pmatrix}$

习题 7.2 验证代入方程和Wronski行列式非零。

习题 7.3 $A(t) = \Phi'(t)\Phi^{-1}(t) = \begin{pmatrix} 2 & 1 \\ 0 & 2 \end{pmatrix}$

习题 7.4 通解 $\mathbf{x} = \begin{pmatrix} (C_1 + C_2t + \frac{t^2}{2})e^{2t} \\ C_2e^{2t} \end{pmatrix}$

7.7 本章小结

本章主要内容:

这些理论是下一章研究常系数线性方程组的基础。