高斯型积分的微扰展开(二)
By 苏剑林 | 2015-03-07 | 22840位读者 |为什么第二篇姗姗来迟? #
其实要写这系列之前,我已经构思好了接下来几篇的内容,本来想要自信地介绍自己想到的一些积分展开的技巧;而且摄动法我本身就比较熟悉,所以正常来说不会这么迟才有第二篇。然而,在我写完第一篇,准备写第二篇的期间,我看到了知乎上的这篇回复:
http://www.zhihu.com/question/24735673
这篇文章大大地拓展了我对级数的认识。里边谈及到了积分的展开是一个渐近级数。这让我犹豫了,怀疑这系列有没有价值,因为渐近级数意味着不管怎样的展开技巧,得到的级数收敛半径都是0。
后来再想想,就算是渐近级数,也有改进的空间,有加速收敛的方法,所以我想我这几篇文章,应该还有一点点意义吧,还可以顺便介绍一下渐近级数和奇点的相关理论。嗯,就这么办吧。
在指数中逐阶展开 #
前一篇文章已经提到的技巧就是对常系数进行变易,通过调整变量使得一阶项为0,该方法主要问题就是不能逐阶地求出高阶近似。然而,它的思想是很不错的,即引入可调整的变量,让展开的各阶项为0,而不要纯粹地做幂级数展开。下面,我们先来试着把展开式放到指数中去,得到一个可以逐级求近似的方法,作为我们以后要进一步发展的方法的原始模型。
总的来说,我们要求积分
$$\int_{-\infty}^{+\infty} e^{-a x^2-\varepsilon x^4} dx$$
上一篇文章中,我们希望用积分$\int_{-\infty}^{+\infty} e^{-A x^2} dx$来逼近它,逼近的思路是通过调整$A$来使得一阶项为0,但这样子我们只能做到这里的,因为只有一个可调整变量,一般来讲只能让一项为0。为了实现逐阶展开,即引入多个可调整的量,我们设
$$A=a+a_1 \varepsilon + a_2\varepsilon^2 + a_3 \varepsilon^3 + \dots$$
然后考虑
$$\int_{-\infty}^{+\infty} \left[e^{-a x^2-\varepsilon x^4}-e^{-(a+a_1 \varepsilon + a_2\varepsilon^2 + a_3 \varepsilon^3 + \dots)x^2}\right] dx$$
将括号内的被积函数按$\varepsilon$的阶展开,就有
$$\begin{aligned}x^2 e^{-a x^2} \left(a_1-x^2\right)\varepsilon &+ \frac{1}{2} x^2 e^{-a x^2} \left(-a_1^2 x^2+2 a_2+x^6\right)\varepsilon^2 \\&- \frac{1}{6} x^2 e^{-a x^2} \left(-a_1^3 x^4+6 a_1 a_2 x^2-6 a_3+x^{10}\right)\varepsilon+\dots\end{aligned}$$
逐项积分,就有
$$\begin{aligned}\sqrt{\frac{\pi}{a} }&\left[\frac{1}{4}(2 a a_1 -3)\left(\frac{\varepsilon}{a^2} \right)+\frac{1}{32}\left(4 a^2 \left(4 a a_2-3 a_1^2\right)+105\right)\left(\frac{\varepsilon}{a^2} \right)^2\right.\\
&\left.+\frac{1}{128}\left(3465-8 a^3 \left(8 a^2 a_3-12 a a_1 a_2+5 a_1^3\right)\right)\left(\frac{\varepsilon}{a^2} \right)^3+\dots\right]\end{aligned}$$
让每项为0,得到方程组
$$\left\{\begin{aligned}&2 a a_1 -3=0\\
&4 a^2 \left(4 a a_2-3 a_1^2\right)+105=0\\
&3465-8 a^3 \left(8 a^2 a_3-12 a a_1 a_2+5 a_1^3\right)=0\\
&\dots\end{aligned}\right.$$
解得
$$a_1=\frac{3}{2 a},\,a_2=-\frac{39}{8 a^3},\,a_3=\frac{657}{16 a^5},\,\dots$$
所以我们有级数
$$\begin{aligned}&\int_{-\infty}^{+\infty} e^{-a x^2-\varepsilon x^4} dx \\
=&\int_{-\infty}^{+\infty} e^{-\left(a+\frac{3}{2}\frac{\varepsilon}{a}-\frac{39}{8}\frac{\varepsilon^2}{a^3}+\frac{657}{16}\frac{\varepsilon^3}{a^5}\dots\right)x^2} dx\\=&\sqrt{\frac{\pi}{a+\frac{3}{2}\frac{\varepsilon}{a}-\frac{39}{8}\frac{\varepsilon^2}{a^3}+\frac{657}{16}\frac{\varepsilon^3}{a^5}\dots}}\end{aligned}\tag{5}$$
这是关于原积分的有一个渐近表达式。$(5)$式效果比$(3)$式好一些;相比于$(4)$式,当$\varepsilon/a^2$很小的时候,截断$(5)$的前几项进行计算,能够得到比$(4)$式更多位有效数字。当然,它在$a\to 0$或者$\varepsilon\to \infty$时则失效了——不过就算是失效,它也只是得到一个有限的0,而不是其它不确定的结果。这也是这种展开方式的好处之一。
渐近级数的意义 #
已经说过,取无穷多项时,渐近级数的收敛半径实际是0,那么渐近级数还有什么意义呢?事实上,首先,对于近似计算来说,可以截断渐近级数的前几项,从而对于某个范围的自变量得到一个比较不错的近似值。其次,对于某些渐近级数,我们可以通过再求和的技巧,来还原它的原始级数,从而得到精确解——量子场论中有大量的例子!
所以,我们需要渐近级数,同时,还可以考虑改进渐近级数的方式(即降低发散速度)。有各种各样的技巧可以做这样的事情,只要你有足够多的创意。在下一篇文章中,我们会再展示一次。
转载到请包括本文地址:https://kexue.fm/archives/3241
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Mar. 07, 2015). 《高斯型积分的微扰展开(二) 》[Blog post]. Retrieved from https://kexue.fm/archives/3241
@online{kexuefm-3241,
title={高斯型积分的微扰展开(二)},
author={苏剑林},
year={2015},
month={Mar},
url={\url{https://kexue.fm/archives/3241}},
}
最近评论