用户工具

站点工具


常微分方程:第七章_线性微分方程组

第七章 线性微分方程组

7.1 引言

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

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

7.2.1 方程组的标准形式

一阶线性微分方程组的标准形式为: <math>\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}</math>

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

7.2.2 矩阵表示

引入向量记号: <math>\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}</math>

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

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

初值条件为 <math>\mathbf{x}(t_0) = \mathbf{x}_0</math>。

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

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

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

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

例 7.2 将 <math>y + p(t)y' + q(t)y = f(t)</math> 化为方程组。 解: 令 <math>x_1 = y, x_2 = y'</math>,则: <math>\begin{cases} x_1' = x_2
x_2' = -q(t)x_1 - p(t)x_2 + f(t) \end{cases}</math> 即 <math>\mathbf{x}' = \begin{pmatrix} 0 & 1
-q(t) & -p(t) \end{pmatrix}\mathbf{x} + \begin{pmatrix} 0
f(t) \end{pmatrix}</math> ===== 7.3 齐次线性方程组的通解结构 ===== ==== 7.3.1 叠加原理 ==== 定理 7.1 若 <math>\mathbf{\varphi}_1(t), \ldots, \mathbf{\varphi}_k(t)</math> 是齐次方程组 <math>\mathbf{x}' = A(t)\mathbf{x}</math> 的解,则它们的线性组合 <math>c_1\mathbf{\varphi}_1 + \cdots + c_k\mathbf{\varphi}_k</math> 也是解。 ==== 7.3.2 基本解组与基本解矩阵 ==== 定义 7.1 齐次方程组的 <math>n</math> 个线性无关的解称为基本解组定义 7.2 以基本解组为列向量构成的矩阵 <math>\Phi(t) = \begin{pmatrix} \mathbf{\varphi}_1(t) & \mathbf{\varphi}_2(t) & \cdots & \mathbf{\varphi}_n(t) \end{pmatrix}</math> 称为基本解矩阵定理 7.2 齐次方程组的通解为: <math>\mathbf{x}(t) = \Phi(t)\mathbf{c} = c_1\mathbf{\varphi}_1 + c_2\mathbf{\varphi}_2 + \cdots + c_n\mathbf{\varphi}_n</math> 其中 <math>\mathbf{c} = (c_1, \ldots, c_n)^T</math> 是任意常向量。 例 7.3 验证 <math>\mathbf{\varphi}_1 = \begin{pmatrix} e^t
e^t \end{pmatrix}, \mathbf{\varphi}_2 = \begin{pmatrix} e^{2t}
2e^{2t} \end{pmatrix}</math> 是方程组 <math>\mathbf{x}' = \begin{pmatrix} 0 & 1
-2 & 3 \end{pmatrix}\mathbf{x}</math> 的基本解组。 解: 验证 <math>\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}</math> ✓ 类似验证 <math>\mathbf{\varphi}_2</math>。 计算Wronski行列式: <math>\det\begin{pmatrix} e^t & e^{2t}
e^t & 2e^{2t} \end{pmatrix} = 2e^{3t} - e^{3t} = e^{3t} \neq 0</math> 故线性无关,是基本解组。 ==== 7.3.3 Wronski行列式与Liouville公式 ==== 定义 7.3 设 <math>\mathbf{\varphi}_1, \ldots, \mathbf{\varphi}_n</math> 是 <math>n</math> 个解,称 <math>W(t) = \det(\mathbf{\varphi}_1, \ldots, \mathbf{\varphi}_n)</math> 为它们的Wronski行列式定理 7.3 <math>n</math> 个解线性无关 <math>\Leftrightarrow W(t) \neq 0</math> 对某点成立。 定理 7.4 (Liouville公式) <math>W(t) = W(t_0)e^{\int_{t_0}^{t} \text{tr}A(s)ds}</math> 其中 <math>\text{tr}A = a_{11} + a_{22} + \cdots + a_{nn}</math> 是矩阵 <math>A</math> 的迹。 ===== 7.4 非齐次线性方程组的通解 ===== ==== 7.4.1 通解结构 ==== 定理 7.5 设 <math>\Phi(t)</math> 是齐次方程组的基本解矩阵,<math>\mathbf{\psi}(t)</math> 是非齐次方程组的特解,则非齐次方程组的通解为: <math>\mathbf{x}(t) = \Phi(t)\mathbf{c} + \mathbf{\psi}(t)</math> ==== 7.4.2 常数变易法 ==== 设特解形式: <math>\mathbf{\psi}(t) = \Phi(t)\mathbf{c}(t)</math> 代入非齐次方程: <math>\Phi'\mathbf{c} + \Phi\mathbf{c}' = A\Phi\mathbf{c} + \mathbf{f}</math> 由于 <math>\Phi' = A\Phi</math>,得: <math>\Phi\mathbf{c}' = \mathbf{f}</math> 故: <math>\mathbf{c}(t) = \int \Phi^{-1}(t)\mathbf{f}(t)dt</math> 特解: <math>\mathbf{\psi}(t) = \Phi(t)\int \Phi^{-1}(t)\mathbf{f}(t)dt</math> 例 7.4 求解 <math>\mathbf{x}' = \begin{pmatrix} 0 & 1
-1 & 0 \end{pmatrix}\mathbf{x} + \begin{pmatrix} 0
\sin t \end{pmatrix}</math>。 解: 齐次方程的基本解矩阵: <math>\Phi(t) = \begin{pmatrix} \cos t & \sin t
-\sin t & \cos t \end{pmatrix}</math> <math>\Phi^{-1}(t) = \begin{pmatrix} \cos t & -\sin t
\sin t & \cos t \end{pmatrix}</math> <math>\Phi^{-1}\mathbf{f} = \begin{pmatrix} -\sin^2 t
\sin t \cos t \end{pmatrix}</math> 积分得: <math>\mathbf{c}(t) = \begin{pmatrix} \frac{t}{2} - \frac{\sin 2t}{4}
-\frac{\cos 2t}{4} \end{pmatrix}</math> 特解: <math>\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}</math> 通解:<math>\mathbf{x} = \Phi(t)\mathbf{c} + \mathbf{\psi}(t)</math>。 ===== 7.5 习题 ===== 习题 7.1 将下列高阶方程化为方程组: a) <math>y
' - 2y + y' - 2y = 0</math> b) <math>y + \omega^2 y = f(t)</math>

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

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

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

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

7.6 参考答案

习题 7.1

a) <math>\mathbf{x}' = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 2 & -1 & 2 \end{pmatrix}\mathbf{x}</math>
b) <math>\mathbf{x}' = \begin{pmatrix} 0 & 1 \\ -\omega^2 & 0 \end{pmatrix}\mathbf{x} + \begin{pmatrix} 0 \\ f(t) \end{pmatrix}</math>

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

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

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

7.7 本章小结

本章主要内容:

  • 矩阵表示:线性方程组可用矩阵形式 <math>\mathbf{x}' = A\mathbf{x} + \mathbf{f}</math> 简洁表示
  • 基本解组与基本解矩阵:齐次方程组的 <math>n</math> 个线性无关解构成基本解组
  • Wronski行列式:判断解的线性无关性
  • 通解结构:非齐次通解 = 齐次通解 + 非齐次特解
  • 常数变易法:求非齐次特解的一般方法
  • Liouville公式:Wronski行列式的演化规律

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

常微分方程/第七章_线性微分方程组.txt · 最后更改: 127.0.0.1