一维驻相法 | Codeboy's Blog

快速振荡的积分往往很小，因为相邻位置的复数方向不断旋转、互相抵消。在相位导数为零的地方，旋转会暂时放慢，一小段区域的贡献得以保留。驻相法计算这些区域留下的大小和方向。

幅度按平方根尺度衰减，复相位却不能丢掉。多个驻点的贡献可能增强，也可能相互抵消；只把各个绝对值相加会改变答案。

形式陈述

一个没有端点干扰的定理

设实相位 $\phi\in C^\infty(\mathbb R)$ ，振幅 $A\in C_c^\infty(\mathbb R)$ ，即无限光滑且在一个有界区间外为零。假设 $A$ 的支撑附近只有一个驻点 $t_0$ ，并且相位的导数满足

\phi'(t_0)=0,\qquad d=\phi''(t_0)\ne0.

由于振幅紧支撑，下式的积分可在包含支撑的有限区间按 Riemann 积分理解。当 $\lambda\to+\infty$ ，

I(\lambda)=\int_{\mathbb R}A(t)e^{i\lambda\phi(t)}dt =e^{i\lambda\phi(t_0)} e^{i\pi\operatorname{sgn}(d)/4} A(t_0)\sqrt{\frac{2\pi}{\lambda|d|}} +O(\lambda^{-3/2}).

这里 $\operatorname{sgn}(d)$ 是 $d$ 的正负号，余项按渐近展开的加法形式估计。 $d>0$ 时额外转过 $\pi/4$ ， $d<0$ 时转过 $-\pi/4$ 。因为首项可能为零或多个贡献可能抵消，完整定理最好采用这种加法余项形式，而非无条件使用比值符号 $\sim$ 。

直觉

从二次相位找出常数

先看模型积分。对 $\eta>0$ ，Gaussian 阻尼使

\int_{\mathbb R}e^{-\eta u^2}e^{is u^2/2}\,du =\frac{\sqrt\pi}{\sqrt{\eta-is/2}},\qquad s\in\{1,-1\},

绝对收敛，其中平方根取实部为正的分支。令 $\eta\downarrow0$ 得

\int_{\mathbb R}^{\mathrm{osc}}e^{is u^2/2}du =\sqrt{2\pi}\,e^{is\pi/4}.

上标说明这是振荡积分，不能用绝对可积性解释：被积函数模长恒为 1。它也等于对称截断极限。对 $|u|\ge R$ 分部积分，利用相位导数 $su$ ，可把尾积分界为常数乘 $1/R$ ，从而证明截断收敛并与阻尼极限一致。

由Taylor 展开，原相位在驻点附近满足 $\phi(t)=\phi(t_0)+d(t-t_0)^2/2+\cdots$ 。换元 $u=\sqrt{\lambda|d|}(t-t_0)$ ，便产生宽度因子 $1/\sqrt{\lambda|d|}$ 和上面的相位常数。

为什么局部二次模型足够

先用一个在 $t_0$ 附近等于 1 的光滑截断函数，把积分分成近端与远端。远端的 $\phi'$ 在支撑上远离零，非驻相分部积分使其对每个固定 $N$ 都是 $O(\lambda^{-N})$ 。

近端可把相位精确变成二次式。令 $s=\operatorname{sgn}(d)$ ，定义

u=(t-t_0) \sqrt{\frac{2s[\phi(t)-\phi(t_0)]}{(t-t_0)^2}}.

根号内在 $t=t_0$ 光滑延拓为 $|d|>0$ 。因此这是局部可逆的光滑坐标，且 $\phi(t)=\phi(t_0)+s u^2/2$ 。变换后的光滑紧支撑振幅 $B(u)$ 满足 $B(0)=A(t_0)/\sqrt{|d|}$ 。

还需证明用 $B(0)$ 替代 $B(u)$ 的误差确为 $O(\lambda^{-3/2})$ 。一种直接核算方式使用既有 Fourier 变换的 $2\pi$ 约定。Gaussian 阻尼后换序，再让阻尼趋零，可得

\int B(u)e^{is\lambda u^2/2}du =\sqrt{\frac{2\pi}{\lambda}}e^{is\pi/4} \int\widehat B(\xi)e^{-is2\pi^2\xi^2/\lambda}\,d\xi.

光滑紧支撑保证 $\int\xi^2|\widehat B(\xi)|d\xi<\infty$ 。用 $|e^{iv}-1|\le|v|$ ，最后一个积分与 $\int\widehat B=B(0)$ 的差不超过 $2\pi^2\lambda^{-1}\int\xi^2|\widehat B|$ 。再乘外面的平方根因子，就得到声称的余项阶。这一步避免把一个在整个实轴上不绝对可积的常振幅模型直接放进控制收敛定理。

例子与边界

一个有精确答案的数值检查

取快速衰减的 Gaussian 振幅，计算

K(\lambda)=\int_{\mathbb R}e^{-t^2}e^{i\lambda t^2/2}dt =\frac{\sqrt\pi}{\sqrt{1-i\lambda/2}}.

这里振幅不是紧支撑，但其各阶导数都快速衰减，远端分部积分仍成立。也可以完全绕开推广定理，直接对右端精确式作展开：

K(\lambda)=e^{i\pi/4}\sqrt{\frac{2\pi}{\lambda}} \left(1-\frac{i}{\lambda}+O(\lambda^{-2})\right).

主项既有大小 $\sqrt{2\pi/\lambda}$ ，也有相位 $\pi/4$ 。在 $\lambda=8,32,128$ 时，仅用主项的相对误差约为 $0.124286,0.031239,0.007812$ ；参数每增大四倍，误差约缩小四倍，符合括号内 $O(\lambda^{-1})$ 的预测。

若把 $i\lambda$ 改成 $-i\lambda$ ，原积分是复共轭，主相位也应从 $\pi/4$ 变成 $-\pi/4$ 。这是检查符号最方便的对称性测试。

边界、退化与干涉

在有限区间积分时，如果振幅在端点不消失，端点还会留下 $O(\lambda^{-1})$ 项。它比内部驻点的 $\lambda^{-1/2}$ 主项小，却比这里紧支撑定理的 $\lambda^{-3/2}$ 余项大。因此不能把同一个余项阶直接搬到任意有限区间。

如果相位是 $t^4$ ，二阶导数在零点也消失，换元给出 $\lambda^{-1/4}$ 的宽度；二次驻相公式失效。若相位依赖另一参数并出现两个驻点合并， $|d|^{-1/2}$ 会发散，提醒我们需要统一的退化相位模型，而非把发散系数当成真实积分发散。

若有有限多个分离的非退化驻点 $t_j$ ，通过局部截断可将上述贡献相加，每项都要保留 $e^{i\lambda\phi(t_j)}$ 和曲率相位。某些参数处主项之和可能为零，此时整体绝对误差仍受控，但相对于主项之和的相对误差没有意义。

推论与应用

为新积分选择局部模型

鞍点法在复平面中通过轮廓变形寻找下降方向；本页在实轴上依赖振荡抵消。它们共享局部二次结构，却有不同的合法性条件。学习终点是对新积分分别列出驻点、端点、退化点及支撑，再决定每一部分使用哪一种估计，而不是看到指数就套同一个常数。

大参数积分复算脚本核对本例精确复数式、主项误差和前几页的正实指数积分，并输出每个断言的输入与残差。

参考资料

NIST Digital Library of Mathematical Functions, §2.3(iv) Method of Stationary Phase，驻点与端点贡献及相位常数。
John K. Hunter, Asymptotic Analysis and Singular Perturbation Theory, §3.3, Theorem 3.8，非退化驻点的展开；§3.4 讨论退化点。
Elias M. Stein, Harmonic Analysis: Real-Variable Methods, Orthogonality, and Oscillatory Integrals, Princeton University Press, 1993, Chapter VIII。

形式陈述 ​

一个没有端点干扰的定理 ​

从二次相位找出常数 ​

为什么局部二次模型足够 ​

一个有精确答案的数值检查 ​

边界、退化与干涉 ​

为新积分选择局部模型 ​

反向链接