本文是《微分方程完全入门》系列的第三篇，覆盖拉普拉斯变换的定义、性质、逆变换及其在解ODE中的应用。阅读本文前需掌握一阶和二阶ODE的基础知识。掌握拉普拉斯变换后，你会发现：解微分方程变成了做代数运算——这是控制理论中传递函数概念的数学起点。

第三篇：拉普拉斯变换——从时域到频域的桥梁

一、概念引入——为什么需要一个”变换”

1.1 问题

解一个二阶ODE需要：

  1. 求特征方程的根
  2. 确定齐次解形式
  3. 猜测特解（待定系数法，常常还要处理”共振”的例外）
  4. 代入初始条件确定任意常数

如果是高阶系统（四阶、六阶），这个流程变得极其繁琐。更糟的是，待定系数法的每一步都需要针对 f(t)f(t) 的形式特殊处理。

有没有一种方法，可以让求解 ODE 变得像解代数方程一样简单？

有：拉普拉斯变换。它的核心思想是：

把时域中的微分方程，通过一个积分变换，变成频域中的一个代数方程。解完代数方程后，再变换回时域。

1.2 历史

时间 人物 贡献
1744 欧拉（Euler） 首次使用了类似拉普拉斯变换的积分形式求解ODE
1785 拉普拉斯（Pierre-Simon Laplace） 在《概率论的分析理论》中系统定义了拉普拉斯变换
1899 赫维赛德（Oliver Heaviside） 独立发明”运算微积分”——本质上就是拉普拉斯变换的工程版本
1930s 多奇（Gustav Doetsch） 严格建立了拉普拉斯变换的数学理论基础

有趣的是，赫维赛德（英国电气工程师）在 1899 年使用的”运算微积分”比数学家们的严格化早了30年——工程界”先用起来”的典型例子。他的方法在电路分析中极其有效，但当时被一些数学家批评为”没有理论基础的直觉”，直到拉普拉斯变换的严格理论建立后才得到证明。

二、定义与基本性质

2.1 定义

一个时域函数 f(t)f(t) 的拉普拉斯变换定义如下：

F(s)=L{f(t)}=∫0∞e−stf(t) dt \boxed{F(s) = \mathcal{L}\{f(t)\} = \int_0^\infty e^{-st} f(t)\,dt}

其中 ss 是一个复变量：s=σ+iωs = \sigma + i\omega。

物理含义：

  • 时域函数 f(t)f(t) 是”随时间变化的量”
  • 拉普拉斯变换将其分解为 不同衰减率（σ\sigma）和不同频率（ω\omega）的指数函数的叠加
  • 变换后的 F(s)F(s) 是复频域中的表示

2.2 基本性质

性质 公式 工程含义
线性 L{af+bg}=aF+bG\mathcal{L}\{af+bg\} = aF+bG 线性系统的叠加原理
微分 L{f′}=sF(s)−f(0)\mathcal{L}\{f'\} = sF(s) - f(0) 微分变成乘法！
二阶微分 L{f′′}=s2F(s)−sf(0)−f′(0)\mathcal{L}\{f''\} = s^2F(s) - sf(0) - f'(0) 高阶微分也一样
积分 L{∫0tf(τ)dτ}=F(s)/s\mathcal{L}\{\int_0^t f(\tau)d\tau\} = F(s)/s 积分变成除法
位移 L{eatf(t)}=F(s−a)\mathcal{L}\{e^{at}f(t)\} = F(s-a) 频域平移
延迟 L{f(t−a)u(t−a)}=e−asF(s)\mathcal{L}\{f(t-a)u(t-a)\} = e^{-as}F(s) 时域延迟=频域乘e−ase^{-as}
卷积 L{f∗g}=F(s)G(s)\mathcal{L}\{f * g\} = F(s)G(s) 时域卷积=频域乘法

其中微分性质是最核心的——它把求解ODE变成了代数运算。

2.3 常见函数的拉普拉斯变换对

时域 f(t)f(t) 频域 F(s)F(s) 备注
δ(t)\delta(t)（单位脉冲） 11 脉冲输入
11（单位阶跃 u(t)u(t)） 1/s1/s 阶跃输入
tt 1/s21/s^2 斜坡输入
tnt^n n!/sn+1n!/s^{n+1} 一般公式
e−ate^{-at} 1/(s+a)1/(s+a) 一阶系统响应
sin⁡ωt\sin\omega t ω/(s2+ω2)\omega/(s^2+\omega^2) 振荡
cos⁡ωt\cos\omega t s/(s2+ω2)s/(s^2+\omega^2) 振荡
e−atsin⁡ωte^{-at}\sin\omega t ω/[(s+a)2+ω2]\omega/[(s+a)^2+\omega^2] 衰减振荡（欠阻尼）
e−atcos⁡ωte^{-at}\cos\omega t (s+a)/[(s+a)2+ω2](s+a)/[(s+a)^2+\omega^2] 衰减振荡（欠阻尼）

三、用拉普拉斯变换求解 ODE

3.1 标准流程

1
2
3
时域ODE   ──拉普拉斯变换──→  频域代数方程
    ↑                              ↓
时域解   ←──拉普拉斯逆变换──   解代数方程

3.2 手算例子 1：一阶ODE（RC电路）

求解：

dvdt+2v=1,v(0)=0 \frac{dv}{dt} + 2v = 1, \quad v(0) = 0

步骤 1：拉普拉斯变换

[sV(s)−v(0)]+2V(s)=1s [sV(s) - v(0)] + 2V(s) = \frac{1}{s} (s+2)V(s)=1s (s+2)V(s) = \frac{1}{s}

步骤 2：解代数方程

V(s)=1s(s+2) V(s) = \frac{1}{s(s+2)}

步骤 3：部分分式展开

1s(s+2)=As+Bs+2 \frac{1}{s(s+2)} = \frac{A}{s} + \frac{B}{s+2}

两边乘以 s(s+2)s(s+2)：

1=A(s+2)+Bs 1 = A(s+2) + Bs

令 s=0s=0：1=2A⇒A=1/21 = 2A \Rightarrow A = 1/2

令 s=−2s=-2：1=−2B⇒B=−1/21 = -2B \Rightarrow B = -1/2

所以：

V(s)=1/2s−1/2s+2 V(s) = \frac{1/2}{s} - \frac{1/2}{s+2}

步骤 4：拉普拉斯逆变换

v(t)=12−12e−2t v(t) = \frac{1}{2} - \frac{1}{2}e^{-2t}

验证：v(0)=0.5−0.5=0v(0) = 0.5 - 0.5 = 0 ✓，v(∞)=0.5v(\infty) = 0.5（稳态值）✓

3.3 手算例子 2：二阶ODE（RLC电路）

求解：

d2ydt2+4y=sin⁡t,y(0)=1,y′(0)=0 \frac{d^2y}{dt^2} + 4y = \sin t, \quad y(0)=1, y'(0)=0

步骤 1：拉普拉斯变换

[s2Y(s)−sy(0)−y′(0)]+4Y(s)=1s2+1 [s^2Y(s) - s y(0) - y'(0)] + 4Y(s) = \frac{1}{s^2+1} s2Y(s)−s+4Y(s)=1s2+1 s^2Y(s) - s + 4Y(s) = \frac{1}{s^2+1} (s2+4)Y(s)=s+1s2+1 (s^2+4)Y(s) = s + \frac{1}{s^2+1}

步骤 2：解代数方程

Y(s)=ss2+4+1(s2+4)(s2+1) Y(s) = \frac{s}{s^2+4} + \frac{1}{(s^2+4)(s^2+1)}

步骤 3：部分分式展开（第二部分）

1(s2+4)(s2+1)=As+Bs2+4+Cs+Ds2+1 \frac{1}{(s^2+4)(s^2+1)} = \frac{As+B}{s^2+4} + \frac{Cs+D}{s^2+1}

展开比较系数得 A=0,B=−1/3,C=0,D=1/3A=0, B=-1/3, C=0, D=1/3：

Y(s)=ss2+4−1/3s2+4+1/3s2+1 Y(s) = \frac{s}{s^2+4} - \frac{1/3}{s^2+4} + \frac{1/3}{s^2+1}

步骤 4：逆变换

y(t)=cos⁡2t−16sin⁡2t+13sin⁡t y(t) = \cos 2t - \frac{1}{6}\sin 2t + \frac{1}{3}\sin t

检验方便性：如果不使用拉普拉斯变换，用待定系数法求解同样的方程需要先求齐次解、猜测特解、处理共振情况…拉普拉斯变换把所有这些变成了有理函数的部分分式展开——一个纯粹代数问题。

四、传递函数——控制理论的核心

4.1 定义

对于一个线性时不变系统，传递函数定义为零初始条件下输出与输入的拉普拉斯变换之比：

G(s)=Y(s)U(s)=输出输入 G(s) = \frac{Y(s)}{U(s)} = \frac{\text{输出}}{\text{输入}}

4.2 从微分方程到传递函数

将系统的ODE两侧做拉普拉斯变换（初始条件为零）：

andnydtn+⋯+a0y=bmdmudtm+⋯+b0u a_n\frac{d^ny}{dt^n} + \cdots + a_0y = b_m\frac{d^mu}{dt^m} + \cdots + b_0u

变成：

(ansn+an−1sn−1+⋯+a0)Y(s)=(bmsm+⋯+b0)U(s) (a_n s^n + a_{n-1}s^{n-1} + \cdots + a_0)Y(s) = (b_m s^m + \cdots + b_0)U(s)

传递函数：

G(s)=Y(s)U(s)=bmsm+⋯+b0ansn+⋯+a0 G(s) = \frac{Y(s)}{U(s)} = \frac{b_m s^m + \cdots + b_0}{a_n s^n + \cdots + a_0}

关键：系统的全部动态特性都编码在 G(s)G(s) 的极点和零点中。

4.3 极点与系统稳定性

  • 分子多项式的根 = 零点（zeros）
  • 分母多项式的根 = 极点（poles）

稳定性判据：所有极点的实部必须为负。

如果有任何极点的实部为正，系统不稳定（响应会发散）。

4.4 手算例子：一阶系统

一阶系统 G(s)=Kτs+1G(s) = \frac{K}{\tau s + 1}，对单位阶跃输入 U(s)=1/sU(s)=1/s：

Y(s)=Kτs+1⋅1s=K(1s−1s+1/τ) Y(s) = \frac{K}{\tau s + 1} \cdot \frac{1}{s} = K\left(\frac{1}{s} - \frac{1}{s + 1/\tau}\right)

逆变换：

y(t)=K(1−e−t/τ) y(t) = K(1 - e^{-t/\tau})

这就是你在PID调参文章中看到的”一阶响应的上升时间 tr≈2.2τt_r \approx 2.2\tau“的数学来源。

五、工程意义——为什么拉普拉斯变换无处不在

拉普拉斯变换在无人机和控制工程中的出现：

场景 用到的拉普拉斯概念 示例文章
PID控制 传递函数 G(s)G(s)、Kp+Ki/s+KdsK_p + K_i/s + K_d s PID调参
滤波 一阶低通 1/(τs+1)1/(\tau s+1) 飞控传感器融合
飞行动力学 纵向模态极点、螺旋/滚转模态 六自由度稳定性
系统辨识 从频率响应估计传递函数 参数辨识

六、核心公式速查卡

公式 含义
F(s)=∫0∞e−stf(t)dtF(s) = \int_0^\infty e^{-st}f(t)dt 拉普拉斯变换定义
L{f′}=sF(s)−f(0)\mathcal{L}\{f'\} = sF(s) - f(0) 微分性质
L{f′′}=s2F(s)−sf(0)−f′(0)\mathcal{L}\{f''\} = s^2F(s) - sf(0) - f'(0) 二阶微分
G(s)=Y(s)/U(s)G(s) = Y(s)/U(s) 传递函数定义
Y(s)=G(s)U(s)Y(s) = G(s)U(s) 输出=传递函数×输入
所有极点 Re(pi)<0\text{Re}(p_i) < 0 稳定性条件

参考文献

  1. Perko, L. (2001). Differential Equations and Dynamical Systems (3rd ed., Chapter 1-2). Springer.
  2. Ogata, K. (2010). Modern Control Engineering (5th ed., Chapter 2-3). Prentice Hall.
  3. Doetsch, G. (1971). Introduction to the Theory and Application of the Laplace Transformation. Springer.
  4. Heaviside, O. (1899). Electromagnetic Theory (Vol. 2). — 运算微积分的原始工作
  5. Laplace, P. S. (1812). Théorie Analytique des Probabilités. — 拉普拉斯变换的首次系统定义

下一节：微分方程组与相平面分析

从一维走向多元：ODE系统、相平面、平衡点分类、李雅普诺夫稳定性。这是理解六自由度无人机方程和控制理论的数学基础。