平稳相近似

平稳相近似是一种处理在给定区间内被积分函数快速振荡的定积分的近似法。下列定积分

$\underset{a}{\overset{b}{\int }}f(t,x)*e^{ivh(t)}dt$

被积分函数$f(t,x)*e^{ivh(t)}$在{a,b}区间快速振荡，而位相在此区间内的变化相对缓慢。这类积分用拉普拉斯近似无效，必须用平稳相近似法。

原理

绿线代表振荡的被积分函数，红线代表位相

主要贡献来自平稳位相区

$F(v,x)=\underset{0}{\overset{\pi }{\int }}cos(v*h(t)dt$

位相：$v*h(t)=v*(t-xsin(t))$

被积分函数：$cos(v(t-xsin(t))$

位相的平稳点由下列微分方程给出：

$\frac{d}{dt}(v*(t-xsin(t))=0$,令$v=5,x=35$得平稳相点：

$t_s=t=arccos(1/35)=1.54$

将积分F(v,x)分为三段：

$sf:=\underset{0}{\overset{1}{\int }}cos(5*t-5*x*sin(t))dt+\underset{1}{\overset{2}{\int }}cos(5*t-5*x*sin(t))dt+\underset{2}{\overset{\pi }{\int }}cos(5*t-5*x*sin(t))dt$

由图可见，在平稳相点左右区间[1,2]，位相平稳，这个平稳相区间贡献了定积分的主要部分，在平稳相区间之外[0,1]和$[2,\pi ]$，由于被积分函数振荡激烈，正负相消，因此贡献可以忽略不计。这就是平稳相近似的原理

一阶近似展开

设有下列积分

$F(v)=\underset{a}{\overset{b}{\int }}{e}^{ivh(t)}f(t)dt$,

并设$h(t)$的平稳相点为$t_k$

将$h(t)$作泰勒展开：

$h(t)\approx h(t_k)+\frac{1}{2}*(t-t_k{)}^2{h}^{″}(t_k)+\cdots$

于是

$F_k(v)\approx e^{ivh(t)}f(t_k)$$\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{iv{s}^2{h}^{″}(t_k)/2}ds$

将所有平稳相点$k=1\cdots N$的贡献相加得

$F(v)\approx \sqrt{\frac{2*\pi }{abs(v)}}$$\sum _{k=1}^N$$\frac{f(t_k)}{\sqrt{abs(h^{″}(t_k))}}exp(i(vh(t)+\frac{\pi }{4}sgn(v{h}^{″}(t_k)))$

一般情形被积分函数是不是积分区间的周期函数，因此积分区间两头的贡献也必须补入，

$F(v)\approx \sqrt{\frac{2*\pi }{abs(v)}}$$\sum _{k=1}^N$$\frac{f(t_k)}{\sqrt{abs(h^{″}(t_k))}}exp(i(vh(t)+\frac{\pi }{4}sgn(v{h}^{″}(t_k)))$$+\frac{i}{v}(\frac{e^{ivh(a)}f(a)}{h^{\prime }(a)}-\frac{ivh(b)f(b)}{h^{\prime }(b)})$

参考文献

• D.Richards Advanced Mathematical Methods with Maple,Cambridge 2002.
• Frank J. Oliver, NIST Handbook of Mathematical Functions,Cambridge 2010