跳到正文
格致开物MATHWIKI

奇异摄动

奇异摄动(singular perturbation)研究小参数趋于零时,解的结构不能由直接令参数为零或普通幂级数一致描述的问题。典型情形是小参数乘在最高阶导数前:参数消失后方程降阶,原来需要满足的初始或边界条件却没有一同消失,于是在狭窄区域内出现迅速变化。

消失的最高阶导数

考察无量纲边值问题 εy″+y′=0,y(0)=0,y(1)=1,0<ε≪1. 若直接令 ε=0,剩下 y′=0,其解为常数,不可能同时取两个不同的边界值。困难不是原问题无解,而是极限过程抹掉了一部分变化。

这个问题能精确求解。令 p=y′,则 εp′+p=0,故 p=Ce−x/ε。再积分并利用边界条件,得到 yε(x)=1−e−x/ε1−e−1/ε. 对每个固定 x>0,极限为一;在 x=0,所有参数下的值却都是零。因此点态极限在左端不连续,解不可能在整个闭区间上一致收敛到常数一。多数变化集中在左端宽度约为 ε 的边界层内。

小参数越小,解从零升到一的区域越集中在左边界
ε分别取0.2、0.1和0.04。小参数并没有让导数项处处不重要;它使非平凡变化集中到更窄的区域。

外解、内变量与匹配

远离左端时,指数项极小,先考虑外问题 y0′=0,用右端条件确定 yout=1。这条解给出了主体区域,却没有满足左端条件。要查看快速变化,把左端放大,令 X=xε,Y(X)=y(εX). 链式法则给出 y′=YX/ε、y″=YXX/ε2。代回方程并乘以 ε,得到 YXX+YX=0. 两项现在处在同一量级。这说明层宽的选择来自主导平衡:如果预设宽度为 δ,两项大小分别约为 ε/δ2 与 1/δ,要让它们一起起作用就需 δ≍ε。

内解要满足 Y(0)=0,也要在 X→∞ 时与外解一匹配。由通解 Y=A+Be−X,先得到 A=1,再得到 B=−1,因此 yin(x)=1−e−x/ε. 通常把内外解相加后减去共同部分,以免重复计算: ycomp=yout+yin−1=1−e−x/ε. 这个复合近似同时解释主体区域和边界层。匹配是一种渐近构造,并不允许直接把一般解的两个不同极限交换;本例另有精确解,可以进一步严格控制误差。

把横坐标缩放为x除以小参数后,各条边界层曲线接近同一指数形状
内坐标 X=x/ε 让原先越来越窄的层恢复为有限宽度。虚线为内解1−exp(−X);有限参数曲线仍带右端归一化的微小修正。

一致误差与层宽的实际含义

记 q=e−1/ε。用精确式直接相减,有 yε−ycomp=(1−e−x/ε)q1−q. 在 0≤x≤1 上它单调增加,最大值在右端,恰为 q。所以 ‖yε−ycomp‖∞=e−1/ε. 它比任意固定正幂的 ε 还小,但外解单独使用时最大误差仍是一,两者差异来自是否保留了边界层。

“宽度约为 ε”没有指定一个唯一的切割点。若要求内解距一不超过 η∈(0,1),应解 e−x/ε≤η,得到 x≥εlog⁡(1/η). 例如 ε=0.05,距离达到一的99%需要 x≈0.05log⁡100≈0.230259。小参数尺度和实际容差共同决定“层已结束”的操作性定义。

边界层所在位置由方程决定

变式 εy″−y′=0 保持相同边界值。精确解可稳定地写作 yε(x)=e(x−1)/ε−e−1/ε1−e−1/ε. 对固定 x<1,它趋于零,快速变化改在右端;外解为零,右端内变量为 X=(1−x)/ε,主内解是 e−X。只把上一例的左层图复制过来就会得到错误结论。

奇异性也可表现为时间初始层或转向点附近的特殊尺度,MIT 讲义分别讨论了这两类情况。[1] 多尺度法通过同时保留快慢变量处理长时间累积的长期项;Feldman 的讲义说明了快慢时间和相应的误差尺度。[2] 边界层匹配则用不同空间或时间区域处理局部快速变化。两者同属渐近方法,构造步骤和误差范围需要分别说明。

对更一般方程,平衡、展开和匹配常只给出候选近似。是否确实接近原解,仍需稳定性估计、余项方程或与独立数值解比较;画出相近曲线不能替代一致误差证明。边界层和匹配方法的系统讲解见 MIT 渐近分析课程资料。[1]

参考资料

  1. ↑ 1.0 1.1 MIT OpenCourseWare. Boundary Layers and Singular Perturbation,Lectures 16 and 17,第2—5页的降阶、快速指数项、内外解与一致近似,第5—7页的转向点。
  2. ↑ William M. Feldman. MATH 6410: Asymptotic Expansions,第4—6、11—16页,快慢时间、长期展开与平均误差尺度。