“数值计算时,可以先提出 $e^{ \lambda p(a)}$,只计算量级适中的缩放积分。如果这个指数已经下溢为零,应报告对数值或缩放值,而不是宣称原积分精确为零。内部 Laplace 方法…”
向一个窄谷里倒入指数权重
形式陈述
定理:曲率进入分母
设
则
若
注意
直觉
为什么宽度是平方根尺度
在谷底附近,Taylor 展开给出
要让指数变化保持在数量级 1,需要
固定
邻域外由唯一极小性和紧性得到正相位间隙,因此贡献至多是
得到主项。Gaussian 常数可以把积分平方,化成平面上的极坐标积分求得,并不是任意选取的归一化。
例子与边界
具体计算与可证明误差
取
此时
本例可以给两项近似
而原区间外有
两界相加就是原积分的证书。
适用边界:四次谷底和多个谷底
若
峰变平后,重要区间更宽,积分衰减更慢;非驻定端点的宽度则为
若有几个分离的非退化全局极小点,应将各个 Gaussian 贡献相加;只保留其中一个会漏掉同阶量。如果点的位置或曲率随另一个参数变化,当两个点靠近或曲率趋零时,各固定参数的公式可能不再一致。对这样的过渡问题,先重新找缩放变量,比继续增加普通展开项更重要。
推论与应用
再算一项,哪些导数会出现
若进一步有
其中
公式可在缩放坐标中逐项核对。振幅二次项给
把这些项分别除以 Gaussian 总质量,对应的偶次矩是
参考资料
- N. G. de Bruijn, Asymptotic Methods in Analysis, 3rd ed., North-Holland, 1970, Chapter 4,内部极值处的 Laplace 方法。
- NIST Digital Library of Mathematical Functions, §2.3(iii) Laplace’s Method 及 §2.3(v) Coalescing Peak and Endpoint,局部化与非一致过渡的边界。
- John K. Hunter, Asymptotic Analysis and Singular Perturbation Theory, §3.2,Gaussian 权重的展开与积分余项控制。