贡献者: JierPeter
定义 1
设各 $a_i(x)$ 和 $f(t)$ 都是区间 $[a, b]$ 上的连续函数,则称
\begin{equation}
\left(\frac{\mathrm{d}^n y}{ \,\mathrm{d}{x} ^n} \right) +a_1 \left(\frac{\mathrm{d}^{n-1} y}{ \,\mathrm{d}{x} ^{n-1}} \right) +a_2 \left(\frac{\mathrm{d}^{n-2} y}{ \,\mathrm{d}{x} ^{n-2}} \right) +\cdots+a_ny=f(x)~
\end{equation}
为区间 $[a, b]$ 上的
线性微分方程(linear differential equation)。
当 $f(x)=0$ 时,称式 1 为齐次的(homogeneous),否则称为非齐次的(inhomogeneous)。
任取 $x_0\in[a, b]$,如果已知该点处 $y, \left(\frac{\mathrm{d} y}{ \,\mathrm{d}{x} } \right) , \left(\frac{\mathrm{d}^{2} y}{ \,\mathrm{d}{x} ^{2}} \right) , \cdots, \left(\frac{\mathrm{d}^{n-1} y}{ \,\mathrm{d}{x} ^{n-1}} \right) $ 的 $n$ 个值,那么我们可以唯一确定式 1 满足这些初值的特解。这一解的存在与唯一性定理将在线性方程组相关章节讨论。
未完成:引用相关的存在与唯一性定理。
1. 齐次线性微分方程的解的性质与结构
区间 $[a, b]$ 上的齐次线性方程形如
\begin{equation}
\left(\frac{\mathrm{d}^n y}{ \,\mathrm{d}{x} ^n} \right) +a_1 \left(\frac{\mathrm{d}^{n-1} y}{ \,\mathrm{d}{x} ^{n-1}} \right) +a_2 \left(\frac{\mathrm{d}^{n-2} y}{ \,\mathrm{d}{x} ^{n-2}} \right) +\cdots+a_ny=0~.
\end{equation}
定理 1 解的线性叠加原理
如果 $y_1(x), y_2(x), \cdots, y_k(x)$ 是式 2 的 $k$ 个解,而 $c_1, c_2, \cdots, c_k$ 是 $k$ 个常数,那么 $c_1y_1(x)+c_2y_2(x)+\cdots+c_ky_k(x)$ 也是一个解。
证明很简单,利用求导的线性性即可。
线性叠加原理,就是解进行 “数乘” 和 “加法” 后仍为解,意味着式 2 的全体解构成了一个线性空间,称为解空间。如果我们能找到这个线性空间的一组基,也就相当于了解了整个线性空间。这样的基向量,被称为基解,一切解都可以表示为它们的线性组合。基解构成的集合,被称为基本解组。
定义 2 线性相关性
对于区间 $[a, b]$ 上的函数 $y_i(x)$,如果存在一组不全为零的常数 $c_i$,使得 $c_1y_1(x)+\cdots+c_ky_k(x)=0$,那么称 $y_1, \cdots, y_k$ 是线性相关的,否则是线性无关的。
定理 2
如果 $y_1(x), y_2(x), \cdots, y_n(x)$ 是式 2 的 $n$ 个线性无关的解,那么它们构成了一组基解。
定理 2 证明思路提示:我们已经知道,式 2 的每个解由 $n$ 个初值唯一确定。在初值点 $x_0$,给定初值条件 $y(x_0), \left(\frac{\mathrm{d} y}{ \,\mathrm{d}{x} } \right) (x_0), \left(\frac{\mathrm{d}^{2} y}{ \,\mathrm{d}{x} ^{2}} \right) (x_0), \cdots, \left(\frac{\mathrm{d}^{n-1} y}{ \,\mathrm{d}{x} ^{n-1}} \right) (x_0)$,如果 $c_1y_1(x)+\cdots+c_ny_n(x)$ 满足全部 $n$ 个初值条件,那么
\begin{equation}
\left\{\begin{aligned}
&c_1y_1(x_0)+\cdots+c_ny_n(x_0)=y(x_0)\\
&c_1y_1'(x_0)+\cdots+c_ny_n'(x_0)=y'(x_0)\\
&\phantom{123456789012345678}\vdots\\
&c_1 \left(\frac{\mathrm{d}^{n-1} y_1}{ \,\mathrm{d}{x} ^{n-1}} \right) (x_0)+\cdots+c_n \left(\frac{\mathrm{d}^{n-1} y_n}{ \,\mathrm{d}{x} ^{n-1}} \right) (x_0)=y(x_0)~.
\end{aligned}\right.
\end{equation}
而由于 $y_1(x), y_2(x), \cdots, y_n(x)$ 在给定区间上是线性无关的,
式 3 的各方程也就彼此独立
1
,因此能够解出唯一的常数组 $\{c_i\}$。“能够解出” 意味着 “能够把满足初值条件的解表示成 $c_1y_1(x)+\cdots+c_ny_n(x)$”。
如何判断一组 $y_1(x), y_2(x), \cdots, y_n(x)$ 在给定区间 $[a, b]$ 上的线性相关性呢?我们需要借助以下朗斯基行列式。
定义 3 朗斯基行列式
对于区间 $[a, b]$ 上存在 $k-1$ 次导函数的 $k$ 个函数 $f_1, f_2, \cdots, f_k$,定义在 $[a, b]$ 上的函数 $W[f_1, f_2, \cdots, f_k](x)$ 为如下行列式:
\begin{equation}
W[f_1, f_2, \cdots, f_k](x)= \begin{vmatrix}f_1&f_2&\cdots&f_k\\f'_1&f'_2&\cdots&f'_k\\ \vdots&\vdots&\ddots&\vdots\\f_1^{(k-1)}&f_2^{(k-1)}&\cdots&f_k^{(k-1)}\end{vmatrix} (x)~,
\end{equation}
称之为给定函数组的
朗斯基行列式(Wronskian)。
朗斯基行列式,本质上就是式 3 的系数矩阵行列式。式 3 解的存在性和 $\{y_i\}$ 的线性相关性紧密相连,因此我们有以下定理:
定理 3
对于 $[a, b]$ 上的函数 $f_i(x)$,如果它们在给定区间上线性相关,那么 $W[f_1, f_2, \cdots, f_k](x)$ 在 $[a, b]$ 上恒为 $0$。
注意,定理 3 的逆定理(W 恒为 $0$ 则 $\{f_i\}$ 线性相关)不成立。我们试举一例反例来说明:
例 1
考虑 $\mathbb{R}$ 上的函数:
\begin{equation}
f_1(x)= \left\{\begin{aligned}
&0 &\quad &(x<0)\\
& \mathrm{e} ^{-\frac{1}{x}} &\quad &(x\geq 0)~.
\end{aligned}\right.
\end{equation}
和
\begin{equation}
f_2(x)= \left\{\begin{aligned}
& \mathrm{e} ^{-\frac{1}{x}} &\quad&(x<0)\\
&0 &\quad& (x\geq 0)~.
\end{aligned}\right.
\end{equation}
简单计算可得 $W[f_1, f_2]$ 在 $\mathbb{R}$ 上恒为 $0$,但 $f_1, f_2$ 在 $\mathbb{R}$ 上线性无关,而且它们还是无穷次可导的函数,是绝佳的反例。
定理 4
如果 $y_1(x), y_2(x), \cdots, y_n(x)$ 是式 2 的 $n$ 个解且它们线性无关,那么它们的朗斯基行列式在 $[a, b]$ 上处处不为 $0$。
定理 4 的证明思路提示:反设存在某 $x_0\in[a, b]$,使得上述朗斯基行列式为 $0$,那么在 $x_0$ 处的初值条件无法总是给出式 3 的解,这与解的存在性矛盾。
2. 非齐次线性微分方程
容易验证,非齐次方程式 1 的解 ${\phi}_1(x), \phi_2(x)$ 和齐次方程式 2 的解 $\overline{\phi}(x)$ 有如下关系:
- $\phi_1(x)+\overline{\phi}(x)$ 是非齐次方程式 1 的解;
- $\phi_1(x)-\phi_2(x)$ 是齐次方程式 2 的解。
因此,我们容易得到以下定理
定理 5 非齐次方程的解集
设 $\phi_1(x), \phi_2(x), \cdots, \phi_n(x)$ 为式 2 的一组基解,$\varphi(x)$ 为式 1 的一个特解,那么式 1 的通解可以表示为
\begin{equation}
c_0\varphi(x)+c_1\phi_1(x)+c_2\phi_2(x)+\cdots+c_n\phi_n(x)~.
\end{equation}
知道了式 2 的基本解组以后,我们可以用常数变易法来解出式 1 .
假设 $\phi_1(x), \phi_2(x), \cdots, \phi_n(x)$ 为式 2 的一组基解,那么式 2 的任意一个解都可以表示为 $c_1\phi_1(x)+c_2\phi_2(x)+\cdots+c_n\phi_n(x)$ 的形式。现在我们把各常数 $c_i$ 写成待定函数 $c_i(x)$,使得
\begin{equation}
c_1(x)\phi_1(x)+c_2(x)\phi_2(x)+\cdots+c_n(x)\phi_n(x)~
\end{equation}
是
式 1 的解。
例 2 非齐次线性微分方程的常数变易法
考虑方程
\begin{equation}
\frac{\mathrm{d}^2 y}{ \,\mathrm{d}{x} ^2}+y=\tan x~,
\end{equation}
且已知其齐次形式的基本解组为 $\{\cos x, \sin x\}$。
设式 9 的通解为
\begin{equation}
y=c_1(x)\cos x+c_2(x)\sin x~,
\end{equation}
将
式 10 代入
式 9 得到
\begin{equation}
c''_1(x)\cos x+c''_2(x)\sin x-c'_1(x)\sin x+c'_2(x)\cos x=\tan x~.
\end{equation}
式 11 有两个未知函数,所以没有唯一确定的解。我们可以有多种方法指定第二个约束式子,实践中当然是怎么方便计算怎么来。
通常,我们令
\begin{equation}
c'_1(x)\cos x+c'_2(x)\sin x=0~,
\end{equation}
这就意味着 $\frac{ \,\mathrm{d}{y} }{ \,\mathrm{d}{x} }=-c_1(x)\sin x+c_2(x)\cos x$,进而
\begin{equation}
\frac{\mathrm{d}^2 y}{ \,\mathrm{d}{x} ^2}=-c_1(x)\cos x-c_2(x)\sin x-c'_1(x)\sin x+c'_2(x)\cos x~.
\end{equation}
重新将式 13 代回式 9 和式 10 得到
\begin{equation}
-c'_1(x)\sin x+c'_2(x)\cos x=\tan x~.
\end{equation}
联立式 12 和式 14 ,解得
\begin{equation}
\left\{\begin{aligned}
c'_1(x)&=\sin^2x/\cos x\\
c'_2(x)&=\sin x~.
\end{aligned}\right.
\end{equation}
解出 $c_1$ 和 $c_2$:
\begin{equation}
\left\{\begin{aligned}
c_1(x)&=C_1-\sin x+ \ln\left(\frac{\cos\frac{x}{2}+\sin\frac{x}{2}}{\cos\frac{x}{2}-\sin\frac{x}{2}}\right) \\
c_2(x)&=C_2-\cos x~,
\end{aligned}\right.
\end{equation}
其中 $C_1, C_2$ 是积分常数。
代回式 10 ,整理后可得式 9 的通解为
\begin{equation}
y=C_1\cos x+C_2\sin x-\sin 2x+\cos x \ln\left(\frac{\cos\frac{x}{2}+\sin\frac{x}{2}}{\cos\frac{x}{2}-\sin\frac{x}{2}}\right) ~.
\end{equation}
例 2 中关键的简化步骤是规定了 $c'_1(x)\cos x+c'_2(x)\sin x=0$。事实上,同一个思路在用常数变易法解非齐次方程时很常用,因为它能避免出现 $c_i(x)$ 的高阶导函数。
1. ^ 式 3 是关于未知量 $c_i$ 的代数方程组。如果 $\{y_i\}$ 线性相关,那么至少有一个 $y_i$ 可以表示为其它几个 $y_i$ 的线性组合,而 $y'_i, y''_i$ 等也可以表示为相同的线性组合,这就使得式 3 各方程中至少有一个能被其它方程表示出来,从而不独立。由于初值条件的任意性,这就会导致多数初值条件下式 3 无法解出 $\{c_i\}$。反过来,$\{y_i\}$ 线性无关,那么式 3 各方程就彼此独立。
致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者
热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。