匹配渐近展开法

来自testwiki
imported>InternetArchiveBot2024年1月8日 (一) 13:46的版本 (Add 1 book for verifiability (20240107)) #IABot (v2.0.9.5) (GreenC bot
(差异) ←上一版本 | 最后版本 (差异) | 下一版本→ (差异)
跳转到导航 跳转到搜索

匹配渐近展开法Template:Lang-en)是数学中用于获得方程方程组高精度近似解的一种常用方法,尤其常用于奇异摄动微分方程的求解。

对于许多奇异摄动问题而言,可以将定义域分成两个或多个部分。其中一部分(通常是范围最大的部分)可以通过正则摄动理论获得渐近展开级数解。然而这个解在其他较小的部分则十分不精确。如果这些部分处于定义域边界上被称为边界层,处于定义域中间则称为内层。可以将边界层或内层内的求解问题当作一个独立的摄动问题处理,以获得相应的“内解”(之前通过正则摄动获得的则称为“外解”)。最后再将内解与外解通过“匹配”的办法合并,以得到在整个定义域内都适用的近似解。[1][2][3]

示例

考虑边值问题

ϵy+(1+ϵ)y+y=0,

其中y为时间t的函数,定义域从0到1,边界条件为y(0)=0y(1)=1ϵ是一个小参数,满足0<ϵ1

外解(t = O(1))

ϵ十分小,故可以当作正则摄动问题处理。取ϵ=0,有

y+y=0.

该方程的解为

y=Aet

其中A为常数。使用边界条件y(0)=0,有A=0。而如果使用另一个边界条件y(1)=1,则有A=e。这说明该解不可能满足所有边界条件,意味着ϵ=0的假设不能在整个定义域中都适用(即奇异摄动问题)。于是,我们能够知道定义域中必定存在一个边界层,其中ϵ与自变量t相比不能再忽略不计。这个边界层位于t=0一侧。于是我们使用另一个边界条件y(1)=1得到适用于边界层以外区域的外解yO=e1t

内解(t = O(ε))

在边界层之内,tϵ都很小,但它们大小相若,故可以定义一个新的O(1) 时间变量τ=t/ϵ。于是原先的边值问题可以改写为

1ϵy(τ)+(1+ϵ)1ϵy(τ)+y(τ)=0,

将两边同乘ϵ再取ϵ=0,得到

y+y=0.

该方程的解为

y=BCeτ

其中BC为常数。使用边界层内的边界条件y(0)=0,得到B=C。故内解为

yI=B(1eτ)=B(1et/ϵ).

匹配

由于对于中间大小的tϵt1)需同时满足内解和外解,故可以令内解的外极限与外解的内极限相等,即limτyI=limt0yO。由此得到常数B=e

合并

ϵ取不同值时的近似解

最后,将匹配好的内解与外解合并,以得到适用于整个定义域的近似解。具体而言,即是将内解与外解相加,再减去内、外解重合部分的值yoverlap(即外解的内极限,或内解的外极限)。此问题中,重合部分的值为e。故可以得到原边值问题的最终近似解为

y(t)=yI+yOyoverlap=e(1et/ϵ)+e1te=e(etet/ϵ).

参考文献

Template:Reflist