常微分方程基础

常微分方程基础与一阶方程

物理问题引入：自然定律为什么总是写成微分方程？

在物理学中，绝大多数基本物理定律（牛顿定律、麦克斯韦方程、热力学传导、薛定谔方程）并不是直接给出物理量本身（如位置 $x$ x 或电荷 $q$ q），而是给出物理量的 变化率：

牛顿第二定律指出：加速度（位置随时间的二阶导数 $\frac{d^{2} x}{d t^{2}}$ d^2 xdt^2）正比于合外力；
放射性原子核衰变指出：单位时间衰变的原子数（数量的一阶导数 $\frac{d N}{d t}$ dNdt）正比于当前存在的原子核总数；
欧姆与法拉第定律指出：感应电动势正比于磁通量的变化率．

将物理量及其各阶导数联系起来的方程，就是 微分方程（Differential Equations）．当未知函数仅依赖于 单一自变量（如仅依赖时间 $t$ t 或空间单坐标 $x$ x）时，称为 常微分方程（ODE）．

1. 常微分方程的基本概念与分类

1.1 阶数与线性

一般 $n$ n 阶常微分方程形式为：

F (x, y, y^{'}, y^{″}, \dots, y^{(n)}) = 0

Fleft(x, y, y', y'', dots, y^(n)right) = 0

阶数（Order）：方程中出现的最高阶导数次数 $n$ n；
线性（Linear）：未知函数 $y$ y 及其所有导数 $y^{'}, y^{″}, \dots$ y', y'', dots 均以一次幂出现，且彼此没有交叉相乘项．线性方程的一般形式为：
$a_{n} (x) y^{(n)} + \dots + a_{1} (x) y^{'} + a_{0} (x) y = f (x)$ a_n(x) y^(n) + dots + a_1(x) y' + a_0(x) y = f(x)
若非齐次项 $f (x) = 0$ f(x) = 0，称为 齐次方程；反之称为 非齐次方程．

1.2 通解与特解

通解：包含与方程阶数相同数目的独立任意常数（ $C_{1}, \dots, C_{n}$ C_1, dots, C_n）的解族；
特解：通过实验给定的 初始条件（如 $t = 0$ t=0 时的位置与速度）唯一确定常数后的解．

2. 常见一阶常微分方程的求解范式

一阶方程的标准形式为 $\frac{d y}{d x} = f (x, y)$ dydx = f(x, y)．

2.1 分离变量法（Separation of Variables）

若方程右端可分解为 $x$ x 与 $y$ y 函数的乘积：

\frac{d y}{d x} = g (x) h (y) ⟹ \frac{d y}{h (y)} = g (x) d x

dydx = g(x) h(y) implies dyh(y) = g(x) dx

两边直接积分即得通解： $\int \frac{d y}{h (y)} = \int g (x) d x + C$ int dyh(y) = int g(x)dx + C．

2.2 一阶线性微分方程与积分因子法

标准一阶线性方程形式为：

\frac{d y}{d x} + P (x) y = Q (x)

dydx + P(x)y = Q(x)

引入积分因子 $μ (x) = \exp (\int P (x) d x)$ mu(x) = expleft(int P(x)dxright)，方程左边可重写为全导数：

\frac{d}{d x} [y e^{\int P (x) d x}] = Q (x) e^{\int P (x) d x}

ddxleft[ y e^int P(x)dx right] = Q(x) e^int P(x)dx

两边积分直接得到通解公式（常数变易法结果）：

y (x) = \exp (- \int P (x) d x) [\int Q (x) \exp (\int P (x) d x) d x + C]

y(x) = expleft(-int P(x)dxright) left[ int Q(x) expleft(int P(x)dxright) dx + C right]

3. 物理建模例题：带黏性阻力的雨滴下落与收尾速度

物理问题：质量为 $m$ m 的雨滴在重力场中自静止下落，空气阻力与瞬时速度成正比： $f_{drag} = - b v$ f_drag = -b v．求雨滴速度随时间演化的规律及收尾速度．

建模与求解

选取向下为正方向，根据牛顿第二定律：

m \frac{d v}{d t} = m g - b v ⟹ \frac{d v}{d t} + \frac{b}{m} v = g

m dvdt = mg - bv implies dvdt + bmv = g

这是一阶线性常系数微分方程，满足 $P (t) = \frac{b}{m}, Q (t) = g$ P(t) = b/m, Q(t) = g．分离变量或利用通解公式：

\frac{d v}{g - \frac{b}{m} v} = d t ⟹ - \frac{m}{b} \ln | g - \frac{b}{m} v | = t + C_{0}

dvg - bmv = dt implies -mb lnleft|g - bmvright| = t + C_0

取指数并利用初值条件 $v (0) = 0$ v(0) = 0：

v (t) = \frac{m g}{b} (1 - \exp (- \frac{b}{m} t))

v(t) = mgbleft( 1 - expleft(-bmtright) right)

结果物理分析

特征弛豫时间： $τ = \frac{m}{b}$ tau = mb，系统在该时间内完成主要加速；
收尾速度（Terminal Velocity）：当 $t ≫ τ$ t gg tau 时，指数衰减项趋近于零，重力与空气阻力达成动力学平衡，速度趋向稳定极限：
$v_{terminal} = lim_{t \to \infty} v (t) = \frac{m g}{b}$ v_terminal = lim_t -> inf v(t) = mgb

4. 学习衔接

下一节：研究力学谐振子、阻尼受迫振动与共振，进入二阶常系数线性微分方程；
更优求解工具：对于带有初始条件的初值问题，亦可直接使用拉普拉斯变换快速求解．

本页面最近更新：2026/10/1 10:49:01，更新历史
发现错误？想一起完善？在 GitHub 上编辑此页！
本页面贡献者：Leafuke, Physics Learning Wiki Team
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供，附加条款亦可能应用

常微分方程基础