这是本文档旧的修订版!
第四章 高阶微分方程
4.1 引言
高阶微分方程是指阶数 $n \geq 2$ 的微分方程。在物理学和工程学中,许多问题需要用高阶微分方程来描述,如力学中的Newton第二定律导出的运动方程通常是二阶的。
本章介绍高阶微分方程的一般理论,重点是线性高阶方程的基本理论,以及降阶法等求解技巧。
4.2 高阶微分方程的一般形式
$n$ 阶微分方程的一般形式为: $F(x, y, y', y'', \ldots, y^{(n)}) = 0$
若能解出最高阶导数,可写成: $y^{(n)} = f(x, y, y', \ldots, y^{(n-1)})$
初值条件为: $y(x_0) = y_0, y'(x_0) = y_1, \ldots, y^{(n-1)}(x_0) = y_{n-1}$
4.3 可降阶的高阶方程
某些特殊类型的高阶方程可以通过变量替换降为低阶方程求解。
4.3.1 不含 $y$ 及其低阶导数的方程
形如 $y^{(n)} = f(x)$ 的方程,直接积分 $n$ 次即可。
例 4.1 求解 $y''' = e^x$。
解: $y'' = e^x + C_1$ $y' = e^x + C_1x + C_2$ $y = e^x + \frac{C_1x^2}{2} + C_2x + C_3$
4.3.2 不含 $y$ 的方程
形如 $F(x, y^{(k)}, y^{(k+1)}, \ldots, y^{(n)}) = 0$,令 $z = y^{(k)}$,降为 $n-k$ 阶方程。
特别地,不含 $y$ 的二阶方程 $F(x, y', y'') = 0$,令 $z = y'$,则 $y'' = z'$,方程变为 $F(x, z, z') = 0$。
例 4.2 求解 $xy'' + y' = x$。
解: 令 $z = y'$,则 $z' = y''$,方程变为: $xz' + z = x$ 即 $z' + \frac{1}{x}z = 1$
这是线性方程,通解为: $z = \frac{C_1}{x} + \frac{x}{2}$
积分得: $y = C_1\ln|x| + \frac{x^2}{4} + C_2$
4.3.3 不含 $x$ 的方程(自治方程)
形如 $F(y, y', y'') = 0$ 的方程,令 $z = y'$,则: $y'' = \frac{dz}{dx} = \frac{dz}{dy}\frac{dy}{dx} = z\frac{dz}{dy}$
方程变为 $F(y, z, z\frac{dz}{dy}) = 0$,这是关于 $z$ 和 $y$ 的一阶方程。
例 4.3 求解 $yy'' + (y')^2 = 0$。
解: 令 $z = y'$,则 $y'' = z\frac{dz}{dy}$。
代入得: $y \cdot z\frac{dz}{dy} + z^2 = 0$ $z(y\frac{dz}{dy} + z) = 0$
情况一: $z = 0$,即 $y' = 0$,得 $y = C$。
情况二: $y\frac{dz}{dy} + z = 0$,分离变量: $\frac{dz}{z} = -\frac{dy}{y}$ $\ln|z| = -\ln|y| + C_1$ $z = \frac{C_1}{y}$
即 $y' = \frac{C_1}{y}$,分离变量: $ydy = C_1dx$ $\frac{y^2}{2} = C_1x + C_2$
通解:$y^2 = C_1x + C_2$(包含了 $y = C$ 的情况)
4.4 高阶线性微分方程
4.4.1 定义与一般理论
定义 4.1 形如 $y^{(n)} + p_1(x)y^{(n-1)} + \cdots + p_{n-1}(x)y' + p_n(x)y = f(x)$ 的方程称为 $n$ 阶线性微分方程,其中 $p_1(x), \ldots, p_n(x), f(x)$ 是已知连续函数。
当 $f(x) \equiv 0$ 时,称为齐次线性方程;否则称为非齐次线性方程。
4.4.2 线性微分算子
引入微分算子 $D = \frac{d}{dx}$,定义 $L[y] = D^n y + p_1(x)D^{n-1}y + \cdots + p_n(x)y$ 或 $L = D^n + p_1(x)D^{n-1} + \cdots + p_n(x)$
则方程可写为 $L[y] = f(x)$。
定理 4.1 $L$ 是线性算子,即:
- $L[y_1 + y_2] = L[y_1] + L[y_2]$ (可加性)
- $L[cy] = cL[y]$ (齐次性)
4.4.3 齐次线性方程的通解结构
定义 4.2 设 $y_1(x), y_2(x), \ldots, y_n(x)$ 是 $n$ 个函数,如果存在不全为零的常数 $c_1, c_2, \ldots, c_n$ 使得 $c_1y_1 + c_2y_2 + \cdots + c_ny_n = 0$ 则称这 $n$ 个函数线性相关;否则称为线性无关。
定理 4.2 (叠加原理) 若 $y_1, y_2, \ldots, y_k$ 是齐次线性方程 $L[y] = 0$ 的解,则它们的线性组合 $c_1y_1 + c_2y_2 + \cdots + c_ky_k$ 也是解。
定理 4.3 $n$ 阶齐次线性方程有 $n$ 个线性无关的解(基本解组),通解为 $y = c_1y_1 + c_2y_2 + \cdots + c_ny_n$
4.4.4 Wronski行列式
定义 4.3 设 $y_1, y_2, \ldots, y_n$ 是 $n$ 个 $n-1$ 次可微函数,称 $W(x) = \begin{vmatrix} y_1 & y_2 & \cdots & y_n \\ y_1' & y_2' & \cdots & y_n' \\ \vdots & \vdots & \ddots & \vdots \\ y_1^{(n-1)} & y_2^{(n-1)} & \cdots & y_n^{(n-1)} \end{vmatrix}$ 为这 $n$ 个函数的Wronski行列式。
定理 4.4 若 $y_1, y_2, \ldots, y_n$ 是齐次线性方程的解,则:
- 它们线性无关 $\Leftrightarrow W(x) \neq 0$ 对某点成立
- 它们线性相关 $\Leftrightarrow W(x) \equiv 0$
定理 4.5 (Liouville公式) $W(x) = W(x_0)e^{-\int_{x_0}^{x} p_1(t)dt}$
4.4.5 非齐次线性方程的通解结构
定理 4.6 设 $y^*$ 是非齐次方程 $L[y] = f(x)$ 的一个特解,$y_1, y_2, \ldots, y_n$ 是对应齐次方程的基本解组,则非齐次方程的通解为: $y = y^* + c_1y_1 + c_2y_2 + \cdots + c_ny_n$
定理 4.7 (叠加原理) 若 $y_1^*$ 是 $L[y] = f_1(x)$ 的特解,$y_2^*$ 是 $L[y] = f_2(x)$ 的特解,则 $y_1^* + y_2^*$ 是 $L[y] = f_1(x) + f_2(x)$ 的特解。
4.5 二阶线性微分方程
4.5.1 一般形式
$y'' + p(x)y' + q(x)y = f(x)$
齐次方程: $y'' + p(x)y' + q(x)y = 0$
4.5.2 已知一个特解求通解
若已知齐次方程的一个非零解 $y_1$,可用降阶法求另一个线性无关的解。
设 $y_2 = y_1 \cdot v(x)$,代入方程确定 $v(x)$。
例 4.4 已知 $y_1 = x$ 是方程 $x^2y'' - 2xy' + 2y = 0$ 的解,求通解。
解: 设 $y = xv$,则 $y' = v + xv'$,$y'' = 2v' + xv''$。
代入方程: $x^2(2v' + xv'') - 2x(v + xv') + 2xv = 0$ $x^3v'' + 2x^2v' - 2xv - 2x^2v' + 2xv = 0$ $x^3v'' = 0$
故 $v'' = 0$,$v = C_1x + C_2$。
通解:$y = x(C_1x + C_2) = C_1x^2 + C_2x$
另一个线性无关解为 $y_2 = x^2$。
4.6 习题
习题 4.1 求解下列可降阶方程:
a) $y''' = x + \sin x$
b) $y'' = y' + x$
c) $yy'' = (y')^2$
习题 4.2 验证 $y_1 = e^x, y_2 = e^{2x}$ 是方程 $y'' - 3y' + 2y = 0$ 的基本解组,并写出通解。
习题 4.3 已知 $y_1 = e^x$ 是方程 $xy'' - (x+1)y' + y = 0$ 的解,求通解。
习题 4.4 计算下列函数组的Wronski行列式,并判断线性相关性:
a) $e^x, e^{2x}, e^{3x}$
b) $1, x, x^2$
c) $e^x, xe^x, x^2e^x$
习题 4.5 证明:若 $y_1, y_2$ 是二阶齐次线性方程的两个线性无关解,则它们的Wronski行列式 $W(y_1, y_2) = y_1y_2' - y_2y_1'$ 在定义域内恒不为零。
4.7 参考答案
习题 4.1
a) $y = \frac{x^4}{24} - \sin x + C_1x^2 + C_2x + C_3$
b) $y = C_1e^x - \frac{x^2}{2} - x + C_2$
c) $y = C_2e^{C_1x}$ 或 $y = C$
习题 4.2 通解 $y = C_1e^x + C_2e^{2x}$
习题 4.3 通解 $y = C_1e^x + C_2(x+1)$
习题 4.4
a) $W = 2e^{6x} \neq 0$,线性无关
b) $W = 2 \neq 0$,线性无关
c) $W = 2e^{3x} \neq 0$,线性无关
4.8 本章小结
本章主要内容:
- 可降阶方程:通过变量替换降为低阶方程
- 不含 $y$:令 $z = y'$
- 不含 $x$:令 $z = y', y'' = z\frac{dz}{dy}$
- 高阶线性方程理论
- 线性微分算子
- 基本解组与通解结构
- Wronski行列式判断线性无关性
- 非齐次方程通解 = 特解 + 齐次通解
- 降阶法:已知一个解可求出另一个线性无关解
这些理论为后续学习常系数线性方程和特殊函数奠定基础。