拉普拉斯变换在求解常微分方程初值问题中的操作过程
第 1 阶段：理解拉普拉斯变换的核心逻辑
明确 拉普拉斯变换能将“求导运算”转换成“乘法运算”，把常微分方程转换为代数方程。
记住 定义式：函数 $f(t)$ 的拉普拉斯变换记作 $F(s)$，其定义为：
$$F(s) = \int_0^{+\infty} f(t) e^{-st} dt$$
观察 此积分将自变量从 $t$（通常代表时间）转换为复变量 $s$。原函数 $f(t)$ 称为“时域”函数，变换后的 $F(s)$ 称为“频域”函数。
理解 初值问题中的初始条件 $f(0)$、$f'(0)$ 会作为“初始值”自动嵌入到代数方程中，这是该方法省去复杂常数求解过程的关键原因。
审视 拉普拉斯变换的适用对象是以 $t$ 为自变量的函数，且要求 $t \ge 0$ 时有定义。对于初值问题，默认研究区间为 $[0, +\infty)$。
第 2 阶段：准备基本变换对与性质
掌握 以下最常用的变换对。求解线性常系数微分方程时，这些对应关系几乎覆盖全部需求。
记住 常见函数的变换对：
| 时域函数 $f(t)$ | 拉普拉斯变换 $F(s)$ | 使用场景 |
|---|---|---|
| $1$ | $\dfrac{1}{s}$ | 常数函数 |
| $t^n$（$n$ 为正整数） | $\dfrac{n!}{s^{n+1}}$ | 多项式函数 |
| $e^{at}$ | $\dfrac{1}{s-a}$ | 指数函数 |
| $\sin(\omega t)$ | $\dfrac{\omega}{s^2 + \omega^2}$ | 正弦振动 |
| $\cos(\omega t)$ | $\dfrac{s}{s^2 + \omega^2}$ | 余弦振动 |
| $t^n e^{at}$ | $\dfrac{n!}{(s-a)^{n+1}}$ | 指数与多项式乘积 |
补充 线性性质：拉普拉斯变换是线性运算。设 $c_1$、$c_2$ 为常数，则：
$$\mathcal{L}\{c_1 f_1(t) + c_2 f_2(t)\} = c_1 F_1(s) + c_2 F_2(s)$$
记忆 这个性质允许对方程逐项取变换，然后分别查找变换对。
第 3 阶段：对微分方程两端执行拉普拉斯变换
写出 给定常微分方程的标准形式。以二阶线性常系数方程为例：
$$a y''(t) + b y'(t) + c y(t) = g(t)$$
其中 $y'(0)$、$y(0)$ 由初值条件给出。
执行 对方程逐项取拉普拉斯变换，记 $\mathcal{L}\{y(t)\} = Y(s)$。
使用 以下导数变换公式——这是核心运算规则：
一阶导数：
$$\mathcal{L}\{y'(t)\} = s Y(s) - y(0)$$
二阶导数：
$$\mathcal{L}\{y''(t)\} = s^2 Y(s) - s y(0) - y'(0)$$
代入 已知初始条件 $y(0) = A$、$y'(0) = B$，得到包含 $Y(s)$ 的代数表达式。右端 $g(t)$ 同样取变换，得到 $G(s)$。原微分方程变为：
$$a[s^2 Y(s) - sA - B] + b[sY(s) - A] + cY(s) = G(s)$$
整理 将所有含 $Y(s)$ 的项集中到等号左侧，其余项移至右侧：
$$(a s^2 + b s + c) Y(s) = G(s) + aA s + aB + bA$$
此时微分方程已被完全转化为代数方程，目标是解出 $Y(s)$。
第 4 阶段：在频域中解代数方程
求解 $Y(s)$，用 $s$ 的有理分式来表示：
$$Y(s) = \frac{G(s) + aA s + aB + bA}{a s^2 + b s + c}$$
注意 如果 $g(t)$ 是零，即齐次方程，则分子仅包含初值项的贡献。
检查 分母多项式 $a s^2 + b s + c$ 称为特征多项式。它的根决定了系统本身的动力学特性，也决定了后续反变换中会出现什么类型的函数（指数、正弦或多项式）。
分解 将右端写成若干个简单分式的和。分解方式取决于分母因式的类型：
- 互异实根：分母因式为 $(s - r_1)(s - r_2)$，拆解为：
$$\frac{A_1}{s - r_1} + \frac{A_2}{s - r_2}$$
- 共轭复根：因式形如 $(s - \alpha)^2 + \omega^2$，拆解为：
$$\frac{A(s - \alpha) + B\omega}{(s - \alpha)^2 + \omega^2}$$
- 重根：因式 $(s - r)^k$，拆解为：
$$\frac{A_1}{s - r} + \frac{A_2}{(s - r)^2} + \cdots + \frac{A_k}{(s - r)^k}$$
求解 待定系数 $A_i$、$B$ 等。通分后对比等式两端分子中 $s$ 各次幂的系数，建立线性方程组，解出所有未知系数。
第 5 阶段：执行逆拉普拉斯变换
获得 $Y(s)$ 的部分分式展开式后，进入最后一步：将 $Y(s)$ 转换回时域函数 $y(t)$。这一操作记作：
$$y(t) = \mathcal{L}^{-1}\{Y(s)\}$$
查找 基本逆变换对应关系，是第 2 阶段变换表的反向使用：
| $F(s)$ | 时域函数 $f(t)$ |
|---|---|
| $\dfrac{1}{s-a}$ | $e^{at}$ |
| $\dfrac{n!}{(s-a)^{n+1}}$ | $t^n e^{at}$ |
| $\dfrac{\omega}{(s-\alpha)^2 + \omega^2}$ | $e^{\alpha t} \sin(\omega t)$ |
| $\dfrac{s-\alpha}{(s-\alpha)^2 + \omega^2}$ | $e^{\alpha t} \cos(\omega t)$ |
执行 对部分分式分解结果逐项进行逆变换。利用 逆变换的线性性质：各项分别变换后相加，即可得到最终解。
代入 变换结果：
$$y(t) = \sum_{i} \mathcal{L}^{-1}\{\text{第 } i \text{ 项}\}$$
核实 将 $y(t)$ 代回原方程，并验证初始条件是否被满足。若符合，则求解完成。
第 6 阶段：完整操作示例
求解 以下二阶初值问题：
$$y'' + 3y' + 2y = 0, \quad y(0) = 1, \quad y'(0) = 2$$
执行 对等式两端取拉普拉斯变换：
$$\mathcal{L}\{y''\} + 3\mathcal{L}\{y'\} + 2\mathcal{L}\{y\} = 0$$
代入 初值条件 $y(0) = 1$、$y'(0) = 2$，并使用导数变换公式：
$$[s^2 Y(s) - s \cdot 1 - 2] + 3[sY(s) - 1] + 2Y(s) = 0$$
展开 括号并整理：
$$(s^2 + 3s + 2) Y(s) = s + 2 + 3 = s + 5$$
解得：
$$Y(s) = \frac{s + 5}{s^2 + 3s + 2}$$
分解 分母：$s^2 + 3s + 2 = (s+1)(s+2)$，故设：
$$\frac{s + 5}{(s+1)(s+2)} = \frac{A}{s+1} + \frac{B}{s+2}$$
求解 系数。通分后比较分子：$s + 5 = A(s+2) + B(s+1)$。代入 $s = -1$，得 $A = 4$；代入 $s = -2$，得 $B = -3$。
回代：
$$Y(s) = \frac{4}{s+1} - \frac{3}{s+2}$$
执行 逆变换。查表得 $\mathcal{L}^{-1}\{\frac{1}{s+1}\} = e^{-t}$，$\mathcal{L}^{-1}\{\frac{1}{s+2}\} = e^{-2t}$，因此：
$$y(t) = 4e^{-t} - 3e^{-2t}$$
验证 当 $t = 0$ 时：$y(0) = 4 - 3 = 1$，符合初值。计算 $y'(t) = -4e^{-t} + 6e^{-2t}$，当 $t = 0$ 时：$y'(0) = -4 + 6 = 2$，符合初值条件。求解完成。
暂无评论，快来抢沙发吧！