当把旋钮拧到零，问题反而爆炸
在上一篇关于正则摄动的指南里，你见过那个乐观的情形：一个问题与某个你能解的问题只差一点点，你把这个差别变成一个小旋钮 epsilon，把答案写成幂级数 y = y_0 + epsilon y_1 + epsilon^2 y_2 + ...，然后求解一串整洁的标准问题。首项 y_0 不过是未扰动的答案，而每个修正都是同一类带已知驱动项的方程。当它奏效时，效果极佳，级数在整个定义域上都一致地好。本篇要讲的，是它失灵的那一刻——以及当那个微小参数拒绝被忽略时，它会做出什么。
下面是那个警告信号。奇异摄动是这样一类问题：epsilon 虽小，却不能简单地置零，因为这么做改变的是问题的种类，而非程度——参见奇异摄动。最戏剧性的标志，是一个小参数乘在最高阶导数上，比如区间 0 到 1 上的 epsilon y'' + a(x) y' + b(x) y = 0，两端都有边界条件 y(0) = A 与 y(1) = B。把旋钮彻底关掉，epsilon = 0，那个二阶方程就塌缩成一个一阶方程 a(x) y' + b(x) y = 0。一阶方程只带一个自由常数，所以只能满足两个边界条件中的一个。这个朴素的极限，字面意义上丢掉了一个边界条件。
边界层：丢失的导数回归之处
物理图像是一架机翼切过空气。远离表面处，空气几乎像没有摩擦一样掠过，于是简单的无摩擦模型工作得极漂亮。但空气会黏在机翼上——在金属表面，它必须静止。从「自由来流速度」一路降到「贴壁处完全停止」的全部急剧变化，都被压进一层紧贴表面、微观地薄的薄片里。这层薄片就是边界层——参见边界层。外层那个无摩擦的故事几乎处处正确；它唯独在这层薄皮里出错，而那里正是我们丢掉的那一项（黏性，或我们的 epsilon y''）称王的地方。
这层有多厚？你用主导平衡原理来弄清楚——这是一门艺术：猜出方程里真正要紧、彼此同量级的两项，而其余可忽略；参见主导平衡。在层内，解变化迅速，于是引入拉伸坐标 X = x / delta，其中 delta 是未知的层宽。由链式法则，每个对 x 的求导都带出一个因子 1/delta，于是 y' 变成 (1/delta) dy/dX，而 y'' 变成 (1/delta^2) d^2y/dX^2。这时 epsilon y'' 项是 epsilon/delta^2 乘以 d^2y/dX^2，而 a(x) y' 项是 (1/delta) 乘以 a dy/dX。要求这两项平衡：epsilon/delta^2 必须与 1/delta 相当，这就逼出 delta = epsilon。层宽是 epsilon 量级——正是那个丢失的二阶导数重新发声之处。
在那个拉伸坐标下，方程在主导阶变成常系数的内层方程 d^2y/dX^2 + a(0) dy/dX = 0——一个你能用特征方程 r^2 + a(0) r = 0 求解的问题。它的根是 r = 0 与 r = -a(0)，所以内层解是 C_1 + C_2 e^{-a(0) X}。若 a(0) > 0，这个指数在你离开壁面时衰减，这恰恰是你想要的那种急升后稳定：解在几倍层宽之内从壁值蹿升到外层值，然后变平。（若 a(0) < 0，层会改坐落在另一端 x = 1——系数的符号告诉你层贴在哪一面壁上。）
匹配渐近展开：把两个世界粘起来
现在我们有两幅局部图像：一个在整段区间大部分上有效、缓慢变化的外层解，它只能满足远端的边界条件；以及一个住在层内、快速变化的内层解，它能满足近端的那个。各自都带着待定常数。匹配渐近展开法就是把它们熔合成一个处处有效的近似的机制——参见匹配渐近展开。诀窍在于存在一个重叠区，一片三不管地带：在层的拉伸尺度上它远离壁面，而在外层尺度上它仍贴近壁面，在那里两套描述本应都成立。逼它们在那里彼此一致，就钉死了那些松动的常数。
- 外层展开。令 epsilon = 0（或对它展开），求解降阶后的一阶方程。只施加远离层的那个边界条件。这给出缓变的外层解 y_out(x)，在主体区域有效，但在薄皮里出错。
- 内层展开。用拉伸变量 X = x/epsilon 放大进入层内，使快速变化看起来是量级为一的。重新标度后的方程现在保留了你丢掉的那个最高阶导数项。求解它，并施加层内的边界条件。这给出 y_in(X)，它仍带着一个自由常数。
- 在重叠区匹配。取外层解的内极限（在 y_out 中令 x -> 0）和内层解的外极限（在 y_in 中令 X -> 无穷），要求二者一致。这一个匹配条件钉死剩下的常数。匹配不是边界条件——它是两个区域之间的桥。
- 构造一致近似。把外层解与内层解相加，再减去它们共同的重叠值（被数了两遍的匹配部分）：y_uniform = y_out + y_in - （共同部分）。结果是一个在整段区间（含层在内）都精确的单一公式。
Model problem: epsilon y'' + y' + y = 0 , y(0) = 0 , y(1) = 1 , 0 < epsilon << 1
OUTER (set epsilon = 0): y' + y = 0 => y_out = C e^{-x}
apply the FAR condition y(1)=1: y_out(x) = e^{1-x}
(it CANNOT also satisfy y(0)=0 -- that is the lost boundary condition)
layer is at x=0; thickness from balance epsilon y'' ~ y' => delta = epsilon
INNER (stretch X = x/epsilon): y_XX + y_X = 0 => y_in = A + B e^{-X}
apply the NEAR condition y(0)=0: A + B = 0 => y_in = A(1 - e^{-X})
MATCH: X -> infinity gives y_in -> A ; x -> 0 gives y_out -> e
so A = e
UNIFORM: y ~ e^{1-x} + e(1 - e^{-x/epsilon}) - e
= e^{1-x} - e * e^{-x/epsilon}另一种病：无界增长的长期项
边界层是空间中的奇异摄动。在时间里有一个孪生的麻烦，而且它会在问题看起来完全温顺时伏击你。取一个弱非线性振子，杜芬方程 y'' + y + epsilon y^3 = 0，这是一个微分方程，描述一根越拉伸刚度越强一点的弹簧。未扰动的运动（epsilon = 0）是诚实的简谐振子 y_0 = cos t，永远有界。于是你搬出正则展开 y = y_0 + epsilon y_1 + ...，代入，整理出 epsilon 阶的方程。麻烦出在驱动它的东西上。
y_1 方程的驱动项是 -y_0^3 = -cos^3 t，而恒等式 cos^3 t = (3/4) cos t + (1/4) cos 3t 藏着一颗毒丸：那个 (3/4) cos t 的部分，恰好以振子自身的固有频率振荡。你在共振地驱动这个系统。对共振驱动的响应不是有界的摆动，而是振幅随时间线性增长的那种——一个正比于 t sin t 的项。这个失控的部分就是长期项，得名于 saeculum「一个漫长的世代」，因为在天体力学里这类项要历经数百年才悄悄占据主导；参见长期项。
仔细读懂出了什么错。杜芬方程的真实解是一个完全有界的振荡——能量守恒，弹簧绝不会失控。长期项 epsilon t sin t 是我们展开方式的赝象，而非物理的。对固定的小 epsilon，当 t 不大时它无害；可一旦 t 达到 1/epsilon 量级，那个「小修正」epsilon t sin t 就已涨到量级为一的大小，近似便一文不值。长期项是一面红旗，表明朴素级数被组织错了方向——而在本例中，深层原因是：真实的非线性振子以一个略微偏移的频率振动，而我们却固执地在原始频率上展开。
多重尺度：给缓慢漂移配上它自己的时钟
第一剂解药是为周期运动找到的：庞加莱-林斯泰特方法干脆把频率拉伸。引入伸缩时间 tau = omega t，其中未知的 omega = 1 + epsilon omega_1 + ...，对 tau 展开，并在每一阶恰好选取频率修正 omega_k，使共振驱动相消——参见庞加莱-林斯泰特方法。对杜芬方程，这给出 omega = 1 + (3/8) epsilon（乘以振幅的平方），正是刚化弹簧那个著名的、依赖于振幅的频率。拉伸频率不是表面上的重新贴标签：在错误的频率上展开，本就是长期项的真正根源，所以纠正它是除去病因，而非压制症状。
但庞加莱-林斯泰特只能修正一个偏移的频率；它无法描述振幅本身缓慢漂移——比如一个弱阻尼摆，快速摆动而摆幅在许多周期中渐渐衰减。要处理这个，你需要更强大的多重尺度分析——参见多重尺度分析。想法很大胆：假装快时间 t 与慢时间 T = epsilon t 是独立变量，把解写成两者的函数 y(t, T)。那唯一的物理时钟被劈成两根——一根快指针计数一个个周期，一根慢指针追踪振幅与相位在长程中如何游移。
在操作上，链式法则把时间导数 d/dt 变成偏导组合 d/dt + epsilon d/dT——缓慢漂移正是通过那个额外的 epsilon 块溜进来的，是偏导数的一个干净用法。展开 y = y_0 + epsilon y_1 + ... 并整理各阶。首项解是熟悉的振荡，但它的振幅 A(T) 和相位现在是慢时间的未知函数。在下一阶，同样那个产生长期项的共振驱动再度出现——回报就在这里：要求它消失并没有丢掉信息，反而交出了一个支配 A(T) 如何演化的方程。慢尺度方程来自扼杀长期项，而非来自手工塞进去的任何额外物理。
你该带走的东西
退后一步，本篇的两半其实是同一个想法穿着两套衣裳。正则摄动在小参数的效应一致地小时奏效；当把 epsilon 置零会改变问题的性质时，它就奇异地失败了——要么是降低了一个微分方程的阶（边界层的故事），要么是在长时间里使共振失谐（长期项的故事）。两种情形的解药在精神上是一样的：别假装一个近似能覆盖整个定义域。找到那个被朴素级数压制掉的第二尺度，给它一个诚实的、属于自己的变量，让两个尺度对话。
最后一点诚实。这些方法是建设性的配方，而不是已替你核验过假设的定理。主导平衡给出的是一个候选层宽，你必须验证它自洽——被丢弃的项确实必须算出来更小。匹配解是一个渐近近似，在 epsilon -> 0 时最锐利，而非精确答案；对固定的 epsilon，它能达到一个最佳精度，再无法更进一步。而边界层可以坐落在任一端，甚至在内部，取决于你必须核查、而非臆断的符号。在这般审慎之下使用，这是整个应用数学中最强大的近似机器之一——同一套工具箱处理空气动力学、半导体器件方程、化学反应动力学，以及一个真实摆动钟摆的缓慢衰减。