1.5 两百年演进里程碑

1.5 两百年演进里程碑

前四节把这门语言的语法、出生现场与三大原型讲完了，本节做第一章的收尾：把从 1747 年到今天两百多年的学科演化压进一张时间线。读历史不是为了背年份，而是为了看清一条主线——每一次范式跃迁，都是"解"的概念与"求解"的手段同时升级；这条主线会一路延伸到第五章的前沿与第七章的工程交付。

一条时间线上的四次跃迁

把两百年切成四段，每一段的核心问题都不同：解是什么（经典奠基期）、解在什么空间里存在（现代理论期）、解怎么算出来（计算革命期）、解能不能学会（智能融合期）。

偏微分方程两百年演进四阶段

四段分界年份只取整十，因为四个时期并不接力棒式交接：1928 年 CFL 条件论文出现在"现代理论期"的年头里，却是计算革命的开山作；施瓦茨 1950 年前后的分布理论为弱解正名时，有限元的专利已经备案。时期的划分是记忆的脚手架，不是历史的封条。

经典奠基期：从物理里认出方程

1747 到十九世纪末，学科的主要工作是把物理现象翻译成方程，再为少数幸运情形找到显式解。弦振动之争（本章第二节已复盘）催生了波动方程与三角级数；十九世纪初热传导问题把傅里叶级数送上舞台；拉普拉斯与泊松把势理论系统化，静电、引力、不可压缩无旋流共用同一个椭圆型方程；格林 1828 年的自助论文提出格林函数，把"边界影响"变成可计算的积分算子，这条线索直通第三章的点源响应。

这一期的解法是"手艺"：分离变量、积分变换、特征线法，各管一段，换一个方程常常从头再来。手艺的局限在十九世纪末暴露得越来越明显——大量有物理背景的问题（湍流、激波、自由边界）根本找不到显式解。但手艺沉淀了两件 durable 的资产：三大原型方程的定性行为（本章第三节），以及叠加原理与级数语言。前者到今天仍是给陌生方程画像的参照系，后者是第四章谱方法的直接祖先。

现代理论期：解住进函数空间

二十世纪上半叶，主题从"找解"变成"说清解在哪里"。阿达马 1902 年提出适定性三要求——存在、唯一、连续依赖——把"一个问题算不算良性"变成可检验的判据（第二章第三节展开）。勒雷 1934 年绕开光滑性引进弱解，证明不可压缩流体方程弱解整体存在；施瓦茨 1945 年后的分布理论把"求导"扩展到不可导对象上，弱解从此有了严格语法（第二章第四节）。索伯列夫空间为"解在什么空间里找"提供了标准住址。

同一时期，数值化的前奏已经吹响。1928 年库朗、弗里德里希斯与列维联合发表的论文给出了显式格式的稳定性条件（后来以三人首字母缩写命名），这是"数值解可信吗"第一次有了定量答案。冯·诺依曼在曼哈顿计划期间发展的稳定性分析方法，让格式设计从玄学变成可验算的程序。这期的遗产也最容易被初学者低估：没有适定性与弱解理论，数值解算出来也不知道自己在算什么——差分格式收敛到的那个人为构造的极限对象，恰恰要用函数空间的语言来描述。

计算革命期：解可以算出来

1950 到 1990 年代，计算机把"算出解"从奢望变成日常。有限差分方法在流体与武器物理中成熟；有限元法从波音的三角机翼设计（Turner、Clough 等人 1956 年前后）走出，凭借对复杂几何的适应力成为结构工程的行业标准；有限体积法在可压缩流动与守恒律问题上立住脚跟；1965 年快速傅里叶变换问世，谱方法的每步代价从平方级降到对数级，"用全局基函数拿指数精度"从纸面优势变成实用武器。第四章会把这些方法逐一请上场，每个都配可复算的实验。

理论侧的总结性成果是拉克斯 1956 年的等价定理：对一个线性适定问题，相容加稳定等价于收敛。它把散落各处的经验（格式要像方程、算起来要稳）收紧成一条可操作的定理，也为第七章的验证与确认提供了数学起点——所谓"验证"，很大程度就是在检查这条定理的两个条件在代码里是否还成立。

智能融合期：解可以被学会

2010 年代以后，机器学习进入这个领域。物理信息网络把方程残差写进损失函数，让网络在"满足数据"与"满足方程"之间取平衡；神经算子不再学单个解，而是学"从输入函数到解函数"的映射本身，换边界条件、换系数不必重新训练。第五章最后一节会拆解这两条路线的原理与局限。需要保持的清醒是：智能方法目前不能替代传统求解器，它在正问题上的效率优势多见于高维参数化场景与快速多查询场景，其误差界与泛化边界仍是开放问题。把它当作第四章方法谱系的新成员，而不是终结者，是本册的立场。

用一次手算与一次上机致敬 1928

CFL 条件值得在第一章就亲手碰一下，因为它是"推导与数值互证"最短的一个案例：对波动方程 u_tt = c²u_xx 取显式中心差分格式，理论推导给出稳定性条件 c·dt/dx 不超过 1（第四章第二节会完整推出）。推导只花两行，但亲眼看看违反它会发生什么，记忆才牢固。

import numpy as np

def wave_explode(r, steps=60):
    """r = c*dt/dx。r<=1 稳定，r>1 爆炸。返回最后一步的峰值。"""
    x = np.linspace(0, 1, 101); dx = x[1] - x[0]; dt = r * dx / 1.0   # c=1
    u = np.exp(-300 * (x - 0.5)**2)          # 初始位移：中间一个鼓包
    u_prev = u.copy()                         # 初始速度为零
    for _ in range(steps):
        u_next = np.zeros_like(u)
        u_next[1:-1] = (2*u[1:-1] - u_prev[1:-1]
                        + r*r * (u[2:] - 2*u[1:-1] + u[:-2]))
        u_prev, u = u, u_next
    return np.abs(u).max()

for r in [0.9, 1.0, 1.05, 1.2]:
    print(f"r = {r:4.2f}   最后一步峰值 = {wave_explode(r):.3e}")
​

典型输出形状：r 取 0.9 与 1.0 时峰值维持在 1 附近（解在正常传播）；r 取 1.05 与 1.2 时峰值涨到数十甚至更高，鼓包被锯齿噪声吞没。理论条件与数值现象严丝合缝——1928 年那篇论文预言的，正是你在循环里看到的爆炸。第四章会把这条条件从泰勒展开推起，并把同样的实验升级为带误差阶数对照的完整复盘。

本节要点

• 四次跃迁：经典奠基（认出方程）、现代理论（解住进函数空间）、计算革命（解可以算出）、智能融合（解可以被学会），四期叠加而非替代；
• 经典遗产：三大原型的定性行为与级数语言，分别是第二章分类学与第四章谱方法的祖先；
• 现代遗产：适定性判据、弱解与索伯列夫空间，为"数值解收敛到了谁"提供了描述语言；
• 计算里程碑：CFL 条件、冯·诺依曼稳定性分析、有限元、FFT 与拉克斯等价定理，构成第七章工程规范的理论底座；
• 智能定位：物理信息网络与神经算子是方法谱系的新成员，误差与泛化问题未关卷；
• 一分钟的互证：CFL 条件的手算与四行代码实验，是全册方法论的缩影。

第一章到此收束：语言已经发明，骨架已经立起。下一章给方程画像——分类判别式与适定性三问，是把本章的直觉变成可操作的判断流程的地方。