微分方程完全入门（一）：从变化率说起——一阶常微分方程与工程建模

本文是《微分方程完全入门》系列的第一篇，覆盖一阶常微分方程（ODE）的基本概念和求解方法。阅读本文仅需掌握一元微积分（导数、基本积分）。

---

第一部分：什么是微分方程

一、概念引入——一个数值计算问题的启示

1.1 问题

假设你有一架四旋翼无人机在悬停，关掉油门让它自由下落。忽略空气阻力，只受重力$g=9.8\ \text{m/s}^2$的作用。

它的速度 $v(t)$ 随时间变化吗？当然。每一秒它都加速 $9.8\ \text{m/s}$。用微积分的语言：

$$\frac{dv}{dt} = g$$

这个方程中包含了未知函数 $v(t)$ 及其导数 $\frac{dv}{dt}$。它描述的不是某个特定的速度值，而是速度如何随时间变化。

这就是微分方程。

1.2 定义

微分方程（Differential Equation） 是包含未知函数及其导数的方程。

根据自变量的个数分为两类：

• 常微分方程（ODE）：未知函数只有一个自变量（通常为时间 $t$ 或空间 $x$）
• 偏微分方程（PDE）：未知函数有多个自变量（如温度 $T(x,y,z,t)$）

阶数（Order）：方程中最高阶导数的阶数

$$\frac{dv}{dt} = g \quad \text{——一阶 ODE}$$ $$\frac{d^2x}{dt^2} + \omega^2 x = 0 \quad \text{——二阶 ODE}$$ $$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2} \quad \text{——二阶 PDE（热方程）}$$

1.3 性质的表格

性质	含义	举例
阶数	最高导数的阶	$\frac{dv}{dt}=g$ 一阶；$\frac{d^2x}{dt^2}+x=0$ 二阶
线性	未知函数及其导数以一次方形式出现	$\frac{dy}{dt}+2y=0$ 线性；$\frac{dy}{dt}+y^2=0$ 非线性
齐次性	所有项都含未知函数或其导数	$y'+y=0$ 齐次；$y'+y=1$ 非齐次
通解	含有任意常数的解族	$v = gt + C$
特解	给定初始条件后确定的唯一解	$v(0)=0 \implies v=gt$

1.4 历史的5分钟

时间	人物	贡献
1671	牛顿（Newton）	在《流数法》中首次系统研究了含”流数”（即导数）的方程——微分方程的起点
1684	莱布尼茨（Leibniz）	发明了现代的导数符号 $\frac{dy}{dx}$，使得方程更易操作
1693	约翰·伯努利（Johann Bernoulli）	首次使用”微分方程”（differential equation）这个术语
1743	欧拉（Euler）	引入积分因子法，建立线性ODE的一般理论
1824	阿贝尔（Abel）	证明五次方程无根式解——这推动了从”求解析解”到”定性分析”的转变

关键转折：早期数学家认为所有微分方程都应该能用初等函数表达。19世纪中叶李雅普诺夫和李群的工作展示了另一种可能——不求出解的表达式，而是通过分析方程的对称性和稳定性来理解解的行为。这种”定性分析”的思路在当代控制理论和飞行动力学中至关重要。

1.5 在已发表文章中的出现

文章	微分方程类型	使用方式
无人机飞行物理学（6-DOF推导）	12个一阶ODE系统	平动+转动方程、四元数微分方程
数值分析完全入门	ODE数值解法	欧拉法、RK4求解
PID调参原理与方法	PID控制器方程	误差的积分/微分方程
势流理论	PDE（拉普拉斯方程）	$\nabla^2 \phi = 0$

---

二、一阶 ODE 的几何含义——方向场

2.1 概念引入

一个一阶 ODE 可以写成一般形式：

$$\frac{dy}{dx} = f(x, y)$$

它的几何含义非常直观：在平面上的每一点 $(x,y)$，我们都可以计算出一个斜率 $\frac{dy}{dx}$。如果把每个点关联一个短线段（箭头），就形成了方向场（slope field / direction field）。

微分方程的解就是一条处处与该方向场相切的曲线。

2.2 手算例子

考虑最简单的 ODE：

$$\frac{dy}{dx} = -y$$

方向场：

• 在 $y=1$ 处，斜率为 $-1$（下坡）
• 在 $y=2$ 处，斜率为 $-2$（更陡的下坡）
• 在 $y=-1$ 处，斜率为 $1$（上坡）
• 在 $y=0$ 处，斜率为 $0$（水平）

通解：$y = Ce^{-x}$（指数衰减）

无论从哪个点出发，解曲线都趋向于 $y=0$。这就是平衡点的概念——稍后在第四篇中讲相平面时会详细讨论。

---

三、一阶 ODE 的三种基本解法

3.1 可分离变量法——最简单的类型

形式：

$$\frac{dy}{dx} = g(x)h(y)$$

即等号右边是 $x$ 的函数乘以 $y$ 的函数（可因子分解）。

解法：

$$\frac{dy}{h(y)} = g(x)\,dx$$

两边分别积分：

$$\int \frac{dy}{h(y)} = \int g(x)\,dx + C$$

手算例子 1：RC 电路放电

一个 RC 电路中的电压满足：

$$\frac{dV}{dt} = -\frac{1}{RC}V$$

其中 $R$ 是电阻，$C$ 是电容，$RC$ 是时间常数。

给定初始电压 $V(0) = V_0$。

解：

这是可分离变量的：

$$\frac{dV}{V} = -\frac{1}{RC}dt$$

两边积分：

$$\int \frac{dV}{V} = \int -\frac{1}{RC}dt$$ $$\ln|V| = -\frac{t}{RC} + C_1$$ $$V = e^{-t/RC + C_1} = Ce^{-t/RC}$$

代入初始条件 $V(0) = V_0$：

$$V_0 = C \cdot e^0 = C$$

特解：

$$\boxed{V(t) = V_0 e^{-t/RC}}$$

物理意义：电压按指数衰减到 0。$RC$ 越大，衰减越慢。在飞控电路中，电容滤波的原理正是这个方程——时间常数决定了滤波器的截止频率。

手算例子 2：自由落体（带空气阻力）

更真实的情况是考虑空气阻力——阻力与速度成正比：

$$m\frac{dv}{dt} = mg - kv$$

即：

$$\frac{dv}{dt} = g - \frac{k}{m}v$$

解：

分离变量：

$$\frac{dv}{g - (k/m)v} = dt$$

令 $a = g$，$b = k/m$：

$$\int \frac{dv}{a - bv} = \int dt$$ $$-\frac{1}{b}\ln|a - bv| = t + C$$

代入初始条件 $v(0) = 0$，得到 $C = -\frac{1}{b}\ln a$，整理：

$$\boxed{v(t) = \frac{mg}{k}\left(1 - e^{-(k/m)t}\right)}$$

物理意义：速度不会无限增大，而是趋向于终端速度（terminal velocity） $v_T = mg/k$。这在四旋翼的下降速率限制、或者截击机的速度边界中都有体现。

3.2 一阶线性 ODE——积分因子法

标准形式：

$$\frac{dy}{dx} + P(x)y = Q(x)$$

如果 $Q(x)=0$，是齐次方程；$Q(x) \neq 0$，是非齐次方程。

解法——积分因子法：

1. 计算积分因子 $\mu(x) = e^{\int P(x)\,dx}$
2. 两边乘以 $\mu(x)$：$$\mu(x)\frac{dy}{dx} + \mu(x)P(x)y = \mu(x)Q(x)$$
3. 左边变成 $(\mu y)'$：$$\frac{d}{dx}(\mu y) = \mu Q$$
4. 两边积分：$$\mu y = \int \mu Q\,dx + C$$
5. 解出 $y$。

手算例子：RL 电路

一个 RL 电路（电阻+电感串联，施加电压 $E$）：

$$L\frac{di}{dt} + Ri = E$$

标准形式：

$$\frac{di}{dt} + \frac{R}{L}i = \frac{E}{L}$$

解：

$P(t) = R/L$，$\mu(t) = e^{\int R/L\,dt} = e^{Rt/L}$

两边乘 $\mu$：

$$e^{Rt/L}\frac{di}{dt} + \frac{R}{L}e^{Rt/L}i = \frac{E}{L}e^{Rt/L}$$

即：

$$\frac{d}{dt}\left(i e^{Rt/L}\right) = \frac{E}{L}e^{Rt/L}$$

积分：

$$i e^{Rt/L} = \frac{E}{L}\int e^{Rt/L}dt = \frac{E}{L}\cdot\frac{L}{R}e^{Rt/L} + C = \frac{E}{R}e^{Rt/L} + C$$

所以：

$$i(t) = \frac{E}{R} + Ce^{-Rt/L}$$

代入初始条件 $i(0) = 0$ 得 $C = -E/R$：

$$\boxed{i(t) = \frac{E}{R}\left(1 - e^{-(R/L)t}\right)}$$

物理意义：电流以时间常数 $\tau = L/R$ 趋近稳态值 $E/R$。在无人机电机驱动器中，电机的电学模型正是 RL 电路——这和 BLDC 电机的电流响应直接相关。

3.3 恰当微分方程

形式：

$$M(x,y)\,dx + N(x,y)\,dy = 0$$

恰条件：

$$\frac{\partial M}{\partial y} = \frac{\partial N}{\partial x}$$

如果条件满足，存在一个函数 $\psi(x,y)$ 使得 $d\psi = Mdx + Ndy$，通解为 $\psi(x,y) = C$。

找 $\psi$ 的方法：从 $\frac{\partial \psi}{\partial x} = M$ 积分得到 $\psi = \int M dx + g(y)$，再用 $\frac{\partial \psi}{\partial y} = N$ 确定 $g(y)$。

手算例子

检验并求解：

$$(2xy + 3)\,dx + (x^2 - 1)\,dy = 0$$

检验恰条件：

$$\frac{\partial M}{\partial y} = 2x, \quad \frac{\partial N}{\partial x} = 2x$$

相等，恰当！

求 $\psi$：

$$\psi = \int (2xy + 3)\,dx = x^2y + 3x + g(y)$$ $$\frac{\partial \psi}{\partial y} = x^2 + g'(y) = x^2 - 1$$

所以 $g'(y) = -1$，$g(y) = -y + C_1$

因此：

$$\psi(x,y) = x^2y + 3x - y = C$$

这就是通解。

---

四、三种解法的对比

方法	适用形式	核心步骤	难度	工程出场率
可分离变量	$y' = g(x)h(y)$	分离+积分	⭐	最高（RC充放电、指数增长/衰减、重力+阻力）
积分因子法	$y' + P(x)y = Q(x)$	乘$\mu$+积分	⭐⭐	高（RL电路、混合器模型、一阶控制系统）
恰当方程	$Mdx+Ndy=0$	检验+找$\psi$	⭐⭐⭐	低（热力学中的恰当微分）

---

五、一阶 ODE 在飞控与仿真中的应用

5.1 一阶控制系统

PID 控制器中的 PI 环节本质上是一阶 ODE。一个比例-积分控制器：

$$u(t) = K_p e(t) + K_i \int_0^t e(\tau)\,d\tau$$

写成微分方程形式：

$$\frac{du}{dt} = K_p\frac{de}{dt} + K_i e(t)$$

5.2 气压计高度滤波

无人机气压计测量高度时使用了低通滤波——本质上就是求解一阶ODE：

$$\frac{dh_{\text{filtered}}}{dt} = \frac{1}{\tau}(h_{\text{raw}} - h_{\text{filtered}})$$

这是一个一阶线性ODE，解就是指数趋近真实高度。$\tau$ 越大，滤波越平滑但延迟越大。

5.3 电机响应

四旋翼电机的转速响应近似为一阶系统：

$$\frac{d\omega}{dt} = -\frac{1}{\tau}\omega + \frac{K}{\tau}u$$

其中 $u$ 是 PWM 输入，$\tau$ 是电机时间常数（通常 50-200 ms）。这个模型直接决定了飞控的响应速度。

---

六、完整概念地图

微分方程
├── 常微分方程 (ODE) — 单自变量
│   ├── 一阶
│   │   ├── 可分离变量: dy/dx = g(x)h(y)
│   │   ├── 一阶线性: y'+P(x)y=Q(x) → 积分因子
│   │   └── 恰当方程: Mdx+Ndy=0
│   ├── 二阶 (本篇预告)
│   │   ├── 齐次 ← 特征根法
│   │   └── 非齐次 ← 待定系数法
│   └── 高阶/方程组 (第四篇)
│       └── 相平面分析、稳定性
└── 偏微分方程 (PDE) — 多自变量 (第六篇)
    └── 拉普拉斯方程、热方程、波动方程

工程应用树:
一阶 ODE ─→ RC/RL 电路、电机响应、低通滤波、自由落体
二阶 ODE ─→ 弹簧-阻尼系统、飞行器模态、PID
方程组 ──→ 六自由度无人机方程、卡尔曼滤波
PDE ────→ 势流、热传导、声波传播

---

七、核心公式速查卡

公式	方法	适用条件
$\frac{dy}{dx} = g(x)h(y) \Rightarrow \int\frac{dy}{h(y)} = \int g(x)dx$	分离变量	$y'=g(x)h(y)$
$y' + P(x)y = Q(x) \Rightarrow y = \frac{1}{\mu}\int\mu Q\,dx$，$\mu=e^{\int P dx}$	积分因子	一阶线性
$Mdx+Ndy=0$，$\partial M/\partial y = \partial N/\partial x \Rightarrow \psi(x,y)=C$	恰当方程	恰当条件满足
$\dot{v} = g - (k/m)v \Rightarrow v(t) = \frac{mg}{k}(1-e^{-(k/m)t})$	分离变量	重力+线性阻力
$\dot{i} = (E-Ri)/L \Rightarrow i(t) = \frac{E}{R}(1-e^{-(R/L)t})$	积分因子	RL电路
$\dot{V} = -V/RC \Rightarrow V(t) = V_0e^{-t/RC}$	分离变量	RC电路放电

---

参考文献

1. Perko, L. (2001). Differential Equations and Dynamical Systems (3rd ed., Chapter 1). Springer. ISBN: 978-0-387-95116-4.
2. Boyce, W. E., & DiPrima, R. C. (2017). Elementary Differential Equations and Boundary Value Problems (11th ed., Chapter 2). Wiley.
3. Newton, I. (1671). De Methodis Serierum et Fluxionum (On the Methods of Series and Fluxions). — 微分方程的起点
4. Euler, L. (1743). “De integratione aequationum differentialium”. Opera Omnia.
5. Strogatz, S. H. (2018). Nonlinear Dynamics and Chaos (2nd ed., Chapter 1-2). CRC Press. — 一阶ODE的定性分析入门

---

下一节：二阶线性方程：弹簧、阻尼与共振

从RLC电路和弹簧-质量-阻尼系统切入，系统讲解二阶线性齐次/非齐次ODE。特征根法、阻尼比与自然频率→这正是四旋翼飞行动力学和PID控制器设计的数学基础。