Skip to content

方法Method

内部 Laplace 方法

Interior Laplace method · Laplace method at a nondegenerate minimum

将实指数积分的非退化内部极小点放大成高斯轮廓,以局部曲率确定主贡献的宽度和常数。

向一个窄谷里倒入指数权重 e−λp(t),p 越低的位置权重越大。若谷底 t0 位于区间内部,而且底部近似抛物线,就会出现 Gaussian 积分。曲率决定谷有多窄,振幅决定谷底有多高,两者一起决定积分主项。

形式陈述 ​

定理:曲率进入分母 ​

设 a<t0<b,p∈C2([a,b]),q 在 [a,b] 连续。假设 t0 是 p 唯一的全局极小点,且其一、二阶导数满足

p′(t0)=0,d=p″(t0)>0.

则

∫abe−λp(t)q(t)dt=e−λp(t0)λ−1/2[q(t0)2πd+o(1)].

若 q(t0)≠0,可把它写成以首项为分母的比值渐近;渐近展开的加法余项形式则不要求首系数非零。若 q(t0)=0,上面的加法陈述仍成立,但首个非零阶需要另算。

注意 d 在平方根内部。把相位从 t2/2 换为 2t2,曲率从 1 变成 4,峰宽减半,积分主常数也应减半。

直觉

为什么宽度是平方根尺度 ​

在谷底附近,Taylor 展开给出

p(t)−p(t0)=d2(t−t0)2+o((t−t0)2).

要让指数变化保持在数量级 1,需要 λ(t−t0)2≈1,故令 t=t0+u/λ。提出最小相位后,局部积分变成

e−λp(t0)λ∫e−λ[p(t0+u/λ)−p(t0)]q(t0+u/λ)du.

固定 u 时,被积函数趋向 e−du2/2q(t0)。仅有逐点收敛还不够,因为积分区间正在伸长。由 d>0,可选固定小邻域,使 p(t)−p(t0)≥c(t−t0)2,其中 c>0。又因 q 有界,缩放后的局部被积函数被 Me−cu2 控制;这个函数在实轴可积,控制收敛定理便允许取极限。

邻域外由唯一极小性和紧性得到正相位间隙,因此贡献至多是 e−λp(t0)O(e−ηλ)。最后利用

∫−∞∞e−du2/2du=2π/d

得到主项。Gaussian 常数可以把积分平方,化成平面上的极坐标积分求得,并不是任意选取的归一化。

相位阶数决定贡献区宽度
例子与边界

具体计算与可证明误差 ​

取

J(λ)=∫−11e−λ(t2/2+t4/4)dt.

此时 d=1,p3=0,p4=6,q=1,所以

J(λ)=2πλ[1−34λ+O(λ−2)].

本例可以给两项近似 A2 一个明确界。先把积分延到实轴,换元 u=λt。对 v≥0,|e−v−1+v|≤v2/2;令 v=u4/(4λ),积分后得到

|JR−A2|≤2πλ10532λ2.

而原区间外有 t2/2+t4/4≥t2/2,所以

0≤JR−J≤2∫1∞e−λt2/2dt≤2e−λ/2λ.

两界相加就是原积分的证书。λ=32 时,A2≈0.4327279909,积分约为 0.4338768415,绝对误差约 0.0011488506,小于上述总界。这里上界的常数也有推导,可以和近似一起复算。

适用边界:四次谷底和多个谷底 ​

若 p(t)=t4,则 p″(0)=0,平方根公式不能用。直接换元 u=λ1/4t,得到

∫−11e−λt4dt∼Γ(1/4)2λ−1/4.

峰变平后,重要区间更宽,积分衰减更慢;非驻定端点的宽度则为 λ−1,也不能混用内部二次公式。这里 Γ(s)=∫0∞e−uus−1du,常数由 v=u4 算出。

若有几个分离的非退化全局极小点,应将各个 Gaussian 贡献相加;只保留其中一个会漏掉同阶量。如果点的位置或曲率随另一个参数变化,当两个点靠近或曲率趋零时,各固定参数的公式可能不再一致。对这样的过渡问题,先重新找缩放变量,比继续增加普通展开项更重要。

推论与应用

再算一项,哪些导数会出现 ​

若进一步有 p∈C4、q∈C2,写 pj=p(j)(t0)、qj=q(j)(t0),则

I(λ)=e−λp02πλd[q0+C1λ+o(λ−1)],

其中

C1=q22d−q1p32d2−q0p48d2+5q0p3224d3.

公式可在缩放坐标中逐项核对。振幅二次项给 q2u2/(2λ);相位四次项给 −q0p4u4/(24λ);两个三次相位因子相乘给 q0p32u6/(72λ);振幅一次项和相位三次项交叉给 −q1p3u4/(6λ)。奇函数项在对称局部区间积分为零。

把这些项分别除以 Gaussian 总质量,对应的偶次矩是 1/d,3/d2,15/d3,就得到 C1。这些矩由分部积分递推获得。为控制余项,可以先限制到随 λ 缓慢扩大的 u 区间,使有限 Taylor 余项一致小,再用 Gaussian 尾部控制其余部分;不能在整个实轴上无条件宣称 Taylor 余项一致。

参考资料
关系图谱11 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系