Day 22 · 2026.07.14

微分方程

Differential Equations — 用变化率书写世界运行的法则
"自牛顿以来,人类逐渐领悟:物理定律总是用微分方程的语言写就。" — Steven Strogatz

常微分方程

Ordinary Differential Equations · 用「此刻的变化率」推演整段未来
常微分方程
直觉版

微分方程不直接告诉你「你在哪」,而是告诉你「在这个位置上,你正以多快的速度改变」。它是一条局部规则;把它一小步一小步顺下去,就长出整条命运轨迹。

看人口。最朴素的假设:新增人数正比于现有人数——人越多生得越多,于是 $\frac{dP}{dt}=rP$,解出来是指数爆炸。但地球装不下无穷多人,资源一紧增长就该刹车。Verhulst 加了个「刹车项」:当 $P$ 逼近环境容量 $K$,增长率被自己按到零。方程成了 $\frac{dP}{dt}=rP\left(1-\frac{P}{K}\right)$,轨迹从爆炸的指数曲线弯成一条优雅的 S 形(logistic)曲线:先缓、再猛、后饱和。没人「设计」这条 S 曲线——它是那条局部规则自己生长出来的形状。

正式定义
$$\frac{dP}{dt}=rP\left(1-\frac{P}{K}\right)$$

$P(t)$ 是随时间变的种群数,$r$ 是内禀增长率,$K$ 是环境容量(可承载上限)。括号里的 $1-P/K$ 是那个刹车:$P$ 很小时它约等于 $1$(近乎指数增长),$P$ 接近 $K$ 时它趋于 $0$(增长停住)。只需给一个初值 $P(0)$,整条曲线就被唯一确定。

K 容量上限 指数 rP logistic S 形 t P
为什么美

同一个 logistic 方程,几乎原封不动地描述了:细菌繁殖、谣言与病毒传播、新技术采纳曲线、肿瘤生长。这些表面毫不相干的现象,共享同一条「自我限制的增长」骨架。这正是微分方程最惊人之处——它不描述某个具体故事,而是捕捉一整类故事背后共同的生长语法。一条方程,写尽万物由盛转饱的宿命。

应用

流行病学的 SIR 模型是一组耦合 ODE,各国封控决策靠它推演;化学动力学、电路、控制系统全是 ODE。在 AI 里,Neural ODE 把神经网络的「层」看成连续时间——残差网络逐层前向传播,本质是数值求解一个 ODE,深度成了「积分时长」。

一句话精华 + 思考题
微分方程给的是局部规则「此刻怎么变」,解出的却是全局命运「整条轨迹」——初值加规则,未来即被锁定。
思考题:logistic 方程完全确定(无随机项),为何仍能描述看似「不可预测」的增长?「确定的规则」与「可预测的未来」是一回事吗?(Day 16 动力系统会给出否定的答案。)

热方程

The Heat Equation · 万物趋于平滑,时间不可倒流
偏微分方程
直觉版

把一根铁棒中间烧红,热会往两边冷处流。流多快?取决于「这点比邻居热多少」。数学把「与邻居的落差」精确化为二阶空间导数(曲率):一个点若像凸起的尖峰(高于邻居)就降温,若像凹陷的谷(低于邻居)就升温。于是所有尖角被磨圆、起伏被抹平,最终归于均匀。

这就是扩散——不只是热,墨水在水中散开、气味在房间弥漫,遵循同一条方程。而它藏着一个秘密:时间有方向。你能看着墨水散开,却从没见过散开的墨水自己聚回一点。热方程正向温和地抹平一切,倒着运行却会疯狂放大最细微的噪声——它给了「时间箭头」一个纯数学的解释。

正式定义
$$\frac{\partial u}{\partial t}=\alpha\,\frac{\partial^2 u}{\partial x^2}$$

$u(x,t)$ 是位置 $x$、时刻 $t$ 的温度,$\alpha>0$ 是热扩散系数。左边「温度随时间怎么变」,右边正是那个曲率:$\frac{\partial^2 u}{\partial x^2}>0$(凹谷)则升温,$<0$(凸峰)则降温。方程说的就是一句话——越不平的地方,被抹平得越急

t=0 t 增大 → 抹平 x
为什么美

热方程的解有一个惊艳的显式形式:一个初始点热源,随时间散开成一条正态分布(高斯钟形),且越散越宽,标准差正比于 $\sqrt{t}$。这条「热核」把两个看似无关的世界焊在一起:确定性的热扩散随机的布朗运动。无数分子的随机乱撞,宏观上竟精确遵守一条光滑的偏微分方程。概率与分析,在此握手。

应用

图像处理的高斯模糊就是让像素亮度按热方程扩散几步;金融的 Black-Scholes 期权定价方程,换元后正是一个热方程。最惊人的是扩散生成模型(Stable Diffusion、Sora 的底层):训练时把图像扩散加噪「抹平」成纯噪声,生成时再学着反转这个不可逆过程——AI 画画,本质是在数学上倒放一部热方程的电影。

一句话精华 + 思考题
热方程把「变化率正比于与邻居的落差」写成一行公式,抹平一切起伏,并给了时间一个不可逆的方向。
思考题:既然反向热方程会放大噪声、数学上「病态」,扩散模型凭什么还能稳定地反转它、生成清晰图像?靠的是什么额外信息填补了这份不可逆?

波方程

The Wave Equation · 一字之差,命运翻转
偏微分方程
直觉版

波方程和热方程长得几乎一样——只把左边的一阶时间导换成二阶。就这一个字之差,行为天翻地覆。

一阶时间导(热方程)说「曲率决定速度」,系统单调滑向平衡、再不回头。二阶时间导(波方程)说「曲率决定加速度」——这就像弹簧:拉离平衡越远,回拉的力越大,于是它冲过平衡点、荡到另一侧、再荡回来,永不停歇地振荡。琴弦、水波、传播的声与光,都是这份「过冲—回弹」在空间里一站站接力。热方程让扰动消失,波方程让扰动旅行

正式定义
$$\frac{\partial^2 u}{\partial t^2}=c^2\,\frac{\partial^2 u}{\partial x^2}$$

$u(x,t)$ 是位移(弦的高度、气压、场强),$c$ 是波速。左边是加速度,右边仍是曲率。d'Alembert 给出了它惊人干净的通解:$u=f(x-ct)+g(x+ct)$——任意一个右行波加一个左行波,保持形状匀速前进。整个解空间,就是「向左走的形状」与「向右走的形状」的自由叠加。

为什么美

把一阶时间导升成二阶,系统整个性格就翻了个面:耗散变守恒、不可逆变可逆。波方程能量不随时间衰减,且时间可逆——一段波动录像倒放,仍是合法解(热扩散倒放则荒谬)。更震撼的是,Maxwell 把电与磁的方程组合并,得到的正是一个波方程,波速 $c$ 恰好等于光速——光,就此被「预言」成一种电磁波。一个方程结构,统一了声、水、光、地震与引力波。

应用

声学与乐器设计、地震波成像(给地球做 CT)、超声与雷达、光纤通信,全都在解波方程。2015 年 LIGO 探测到的引力波,是爱因斯坦场方程在弱场下退化出的波方程之解,迟到百年的验证。量子力学的薛定谔方程,更是波方程与热方程间用复数搭起的近亲。

一句话精华 + 思考题
把时间导数从一阶提到二阶,抹平的扩散就变成不衰减的振荡——同一副数学骨架,因一个数字而分出「消失」与「旅行」两种世界。
思考题:为什么正是「二阶时间导」带来了时间可逆与能量守恒?物理定律的基本方程大多是二阶时间导($F=ma$ 也是),这暗示了宇宙怎样的一种深层对称?

数值方法

Numerical Methods · 当公式解不存在时,就顺着斜率走
数值分析
直觉版

残酷的事实:绝大多数微分方程没有初等函数写成的解。三体问题、天气、湍流——你永远得不到一个干净的 $u(t)=\dots$ 公式。但方程给了你一样宝贵的东西:每一点的斜率。那就够了。

最朴素的 Euler 法像蒙眼走山路:站在当前点,摸一下脚下的坡度(方程算出的斜率),沿这方向迈一小步 $h$,落到新点,再摸坡度、再迈步……一步步把连续曲线近似成一条折线。步子越小越贴合真解,但每步的微小误差会累积,折线慢慢偏离真曲线。改进的思路是 Runge-Kutta(RK4):迈步前先在步内多处试探斜率、取加权平均,让这一步「预判」曲线的弯度,精度陡增数个数量级。

正式定义
$$y_{n+1}=y_n+h\,f(t_n,y_n)$$

这是 Euler 法:$y_n$ 是当前值,$f(t_n,y_n)$ 是方程给出的当前斜率,$h$ 是步长。「当前值 + 斜率 × 步长」= 下一点,它是所有数值积分器的原子。RK4 只是把这步的斜率 $f$ 换成四次探测的加权平均,误差从与 $h^2$ 成正比压到与 $h^5$ 成正比。

真解 Euler 折线(误差累积) h
为什么美

数值方法是微积分的逆向工程:微积分把离散求和推向连续极限,数值方法再把连续切回离散,让沉默的方程说出具体数字。它逼你直面一个三角权衡——精度、步长、稳定性:步子太大会失真,太小则误差累积且算得慢;而某些方法步长一旦超过阈值,误差还会指数爆炸、彻底发散。连续与离散的这场对话,是全部科学计算的底色。

应用

天气预报、流体与热仿真、游戏物理引擎、火箭轨道(阿波罗与 SpaceX 的制导),全靠数值求解微分方程。在 AI 里,扩散模型的每一次采样,就是在数值求解一个反向 ODE/SDE——DDIM、DPM-Solver 这些「加速采样器」,本质是更聪明的 Runge-Kutta。

一句话精华 + 思考题
解不出公式,就顺着方程给的斜率一步步走——数值方法把连续切回离散,代价是精度、步长与稳定性之间永恒的权衡。
思考题:既然步长越小越准,为何不取无穷小的步长?(提示:浮点舍入误差与计算成本会从另一头反扑——见 Day 62。)

深入思考

开放问题 · 推向边界与跨学科
为什么热方程不可逆、波方程却可逆?差别只在一个时间导数的阶数?
关键在时间反演 $t\to -t$。二阶导 $\partial^2/\partial t^2$ 在此替换下不变(负负得正),波方程原封不动——倒放的波仍是合法解,能量守恒。一阶导 $\partial/\partial t$ 却变号,正向抹平的热方程一倒放就成了放大噪声的「反扩散」,数学上病态、物理上从不发生。方程的时间对称性,直接决定了它描述的世界有没有「时间箭头」——一个符号之差,划开可逆与不可逆两个宇宙。
扩散模型如何「反转」本该不可逆的热方程?
诀窍是它并不真的反解那个病态方程,而是额外学了一个东西:一个神经网络在海量图像上训练,学会在每个噪声水平下估计「往哪个方向去噪最像真实图像」(即数据分布对数密度的梯度,score)。这份从数据偷来的先验,恰好填补了不可逆过程丢失的信息,把病态反演变成一个有唯一稳定解的反向随机微分方程。不可逆性并未被违背——是数据的知识替时间倒流付了账。
什么是「适定问题」(well-posed)?为什么它是应用数学的生死线?
Hadamard 提出三条:解存在唯一、且连续依赖初值(输入微扰只引起输出微扰)。第三条最关键:反向热方程正栽在这里——初始数据一点点测量噪声,会被放大成解里的滔天巨浪,任何数值计算都不可信。现实中的「反问题」(从模糊照片反推清晰图、从地表数据反推地下结构)大多病态,必须靠正则化——人为注入先验假设——才能驯服。适定与否,决定一个模型能不能真的拿来算。
数值误差在混沌系统里为何注定失控,这是否意味着长期预测不可能?
混沌系统对初值指数敏感(蝴蝶效应,见 Day 16)。数值方法每步引入的微小舍入误差,等于不断给系统注入新的初值扰动,而它们以指数速度被放大。于是无论算法多精、步长多小,超过某个「可预测时长」后,数值轨迹与真解必然面目全非——这正是天气预报到两周就失效的根本原因。但这不等于全盘无知:单条轨迹算不准,其统计规律(吸引子形状、长期平均)却往往稳定可算。确定性方程 → 不可长期预测的个体 + 可预测的整体,是微分方程最深的一课。