奇异摄动
奇异摄动(singular perturbation)研究小参数趋于零时,解的结构不能由直接令参数为零或普通幂级数一致描述的问题。典型情形是小参数乘在最高阶导数前:参数消失后方程降阶,原来需要满足的初始或边界条件却没有一同消失,于是在狭窄区域内出现迅速变化。
消失的最高阶导数
考察无量纲边值问题 若直接令 ,剩下 ,其解为常数,不可能同时取两个不同的边界值。困难不是原问题无解,而是极限过程抹掉了一部分变化。
这个问题能精确求解。令 ,则 ,故 。再积分并利用边界条件,得到 对每个固定 ,极限为一;在 ,所有参数下的值却都是零。因此点态极限在左端不连续,解不可能在整个闭区间上一致收敛到常数一。多数变化集中在左端宽度约为 的边界层内。
外解、内变量与匹配
远离左端时,指数项极小,先考虑外问题 ,用右端条件确定 。这条解给出了主体区域,却没有满足左端条件。要查看快速变化,把左端放大,令 链式法则给出 、。代回方程并乘以 ,得到 两项现在处在同一量级。这说明层宽的选择来自主导平衡:如果预设宽度为 ,两项大小分别约为 与 ,要让它们一起起作用就需 。
内解要满足 ,也要在 时与外解一匹配。由通解 ,先得到 ,再得到 ,因此 通常把内外解相加后减去共同部分,以免重复计算: 这个复合近似同时解释主体区域和边界层。匹配是一种渐近构造,并不允许直接把一般解的两个不同极限交换;本例另有精确解,可以进一步严格控制误差。
一致误差与层宽的实际含义
记 。用精确式直接相减,有 在 上它单调增加,最大值在右端,恰为 。所以 它比任意固定正幂的 还小,但外解单独使用时最大误差仍是一,两者差异来自是否保留了边界层。
“宽度约为 ”没有指定一个唯一的切割点。若要求内解距一不超过 ,应解 ,得到 例如 ,距离达到一的99%需要 。小参数尺度和实际容差共同决定“层已结束”的操作性定义。
边界层所在位置由方程决定
变式 保持相同边界值。精确解可稳定地写作 对固定 ,它趋于零,快速变化改在右端;外解为零,右端内变量为 ,主内解是 。只把上一例的左层图复制过来就会得到错误结论。
奇异性也可表现为时间初始层或转向点附近的特殊尺度,MIT 讲义分别讨论了这两类情况。[1] 多尺度法通过同时保留快慢变量处理长时间累积的长期项;Feldman 的讲义说明了快慢时间和相应的误差尺度。[2] 边界层匹配则用不同空间或时间区域处理局部快速变化。两者同属渐近方法,构造步骤和误差范围需要分别说明。
对更一般方程,平衡、展开和匹配常只给出候选近似。是否确实接近原解,仍需稳定性估计、余项方程或与独立数值解比较;画出相近曲线不能替代一致误差证明。边界层和匹配方法的系统讲解见 MIT 渐近分析课程资料。[1]
参考资料
- ↑ 1.0 1.1 MIT OpenCourseWare. Boundary Layers and Singular Perturbation,Lectures 16 and 17,第2—5页的降阶、快速指数项、内外解与一致近似,第5—7页的转向点。
- ↑ William M. Feldman. MATH 6410: Asymptotic Expansions,第4—6、11—16页,快慢时间、长期展开与平均误差尺度。