平稳相近似

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

\int_{a}^{b} f(t,x)*e^{ivh(t)}dt

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

原理
F(v,x)=\int_{0}^{\pi} cos(v*h(t) dt

位相：vh(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 := \int_{0}^{1}cos(5t-5xsin(t))dt +\int_{1}^{2}cos(5t-5xsin(t))dt +\int_{2}^{\pi}cos(5t-5x*sin(t))dt

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

一阶近似展开
设有下列积分

F(v)=\int_{a}^{b}e^{ivh(t)}f(t)dt ,

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

将h(t)作泰勒展开：

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

于是

F_{k}(v) \approx e^{ivh(t)}f(t_k)\int_{-\infty}^{\infty}e^{ivs^{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(vh*(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(vh*(t_k)))}+\frac{i}{v}(\frac{e^{ivh(a)}f(a)}{h'(a)}-\frac{ivh(b)f(b)}{h'(b)})

参考文献
*D.Richards Advanced Mathematical Methods with Maple,Cambridge 2002.
*Frank J. Oliver, NIST Handbook of Mathematical Functions,Cambridge 2010