求函数$f(L,\Delta)=2^{-L}\binom{L}{\frac{L-\Delta}{2}}$等势线的连续近似方法
嘿,我之前在处理组合数的渐近行为时也遇到过几乎一模一样的问题!你用斯特林公式的思路完全没问题,卡在那些指数项上太正常了——其实只要取个对数就能把复杂的指数项拆成可处理的形式,我来一步步给你捋清楚:
首先,我们先明确函数定义:$f(L,\Delta)=2^{-L}\binom{L}{\frac{L-\Delta}{2}}$,这里默认$\frac{L-\Delta}{2}$是整数(后续连续近似可以忽略奇偶性限制)。
用斯特林公式$n! \approx \sqrt{2\pi n}\left(\frac{n}{e}\right)n$展开组合数,然后对整个函数取自然对数——这是处理指数项的核心技巧,能把$(L-\Delta){L-\Delta}$这类项转化为线性的对数项:
$$
\ln f(L,\Delta) \approx -\frac{L}{2}\left[(1-x)\ln(1-x)+(1+x)\ln(1+x)\right] + \frac{1}{2}\ln\left(\frac{2}{\pi L(1-x^2)}\right)
$$
这里我提前做了变量替换:$x=\frac{\Delta}{L}$(也就是$\Delta=xL$),这个替换能把二元函数转化为单变量关系,大幅简化后续分析。
等势线的要求是$f(L,\Delta)=C$($C$为常数),对应$\ln f=\ln C$,我们可以根据$x$的大小(即$\Delta$与$L$的相对关系)分两种情况处理:
场景1:$\Delta$与$L$同阶($x$不趋近于0)
当$L$很大时,含$L$的主导项远大于后面的对数小项,我们可以先忽略对数项,得到近似关系:
$$
(1-x)\ln(1-x)+(1+x)\ln(1+x) \approx -\frac{2\ln C}{L}
$$
这个式子没有解析解,但可以用牛顿迭代法做数值求解,比如给定$L$求$\Delta$,或者反过来给定$\Delta$求$L$,收敛速度很快。
场景2:$\Delta \ll L$($x$很小)
这时候可以对左边的对数项做泰勒展开(奇函数项会抵消,只保留偶次幂):
$$
(1-x)\ln(1-x)+(1+x)\ln(1+x) = x^2 + \frac{x^4}{3} + \frac{x^6}{5} + \dots
$$
代入等势线条件,先取一阶近似(忽略高次项):
$$
-Lx^2 \approx 2\ln C
$$
注意$f(L,\Delta)$是概率类的函数,最大值约为$\sqrt{\frac{2}{\pi L}}<1$,所以$\ln C$是负数,两边都是正数。把$x=\frac{\Delta}{L}$代回去,得到:
$$
\Delta^2 \approx -2L\ln C \quad \Rightarrow \quad \Delta \approx \sqrt{-2L\ln C}
$$
这是一条抛物线型的等势线,$\Delta$和$\sqrt{L}$成正比。如果需要更精确的结果,可以加入$x^4$项修正:
$$
\Delta^2 \approx -2L\ln C - \frac{4(\ln C)^2}{3}
$$
如果$L$不是特别大,不能忽略对数小项,可以把完整的等式整理为:
$$
-L\left[(1-x)\ln(1-x)+(1+x)\ln(1+x)\right] = 2\ln C + \ln\left(\frac{\pi L(1-x^2)}{2}\right)
$$
再用数值方法求解,这个修正能让近似结果在中等大小的$L$下更准确。
内容的提问来源于stack exchange,提问作者Alice Schwarze

