特殊类型二项式和椭圆积分

吴炳洋

已于 2022-03-27 13:02:02 修改

阅读量495

点赞数

分类专栏：数学基础文章标签：数字通信统计学

于 2022-03-27 06:37:08 首次发布

本文链接：https://blog.csdn.net/wubyatseu/article/details/123767639

版权

数学基础专栏收录该内容

20 篇文章 1 订阅

订阅专栏

我们所熟知的二项式展开 $(a+b)^N=\sum_{j=0}^N\begin{pmatrix}N\\j\end{pmatrix}a^jb^{(N-j)}=\sum_{j=0}^N\begin{pmatrix}N\\j\end{pmatrix}a^{(N-j)}b^j \tag 1$
二项式展开的系数 $\begin{pmatrix}N\\j\end{pmatrix}=\frac{N!}{(N-j)!\ j!}=\frac{\Gamma(N+1)}{\Gamma(N-j+1)\Gamma(j+1)} \tag 2$
如果 $N$ 不是正整数，例如 $N = 1 / 2$ 或 $N = - 1 / 2$ ，这种特殊类型的二项式展开类似地可表示为
$\sqrt{1-x}=(1-x)^{1/2}=\sum_{j=0}^\infty\begin{pmatrix}1/2\\j\end{pmatrix}(-x)^j \tag 3$
${1\over\sqrt{1-x}}=(1-x)^{-1/2}=\sum_{j=0}^\infty\begin{pmatrix}-1/2\\j\end{pmatrix}(-x)^j \tag 4$
这里
$\begin{pmatrix}1/2\\j\end{pmatrix}=(-1)^{j-1}\frac{(2j-3)!!}{(2j)!!} \tag 5$
$\begin{pmatrix}-1/2\\j\end{pmatrix}=(-1)^{j}\frac{(2j-1)!!}{(2j)!!} \tag 6$
一般地
$\begin{pmatrix}x\\j\end{pmatrix}=\frac{\Gamma(x+1)}{\Gamma(x-j+1)\Gamma(j+1)} \tag 7$ 利用
$\Gamma(x)={\Gamma(x+1)\over x} \tag 8$ 可计算（7）。

利用这种特殊类型的二项式展开，可以计算如下的椭圆积分

The complete elliptic integral of the first kind is defined as follows:
$K(m)=F(\left. {\pi\over 2}\right|m)=\int_0^{\pi/2}{1\over \sqrt{1-m^2\sin^2(\theta)}}d\theta \tag 9$
The complete elliptic integral of the second kind is defined as follows:
$E(m)=E(\left. {\pi\over 2}\right|m)=\int_0^{\pi/2}{ \sqrt{1-m^2\sin^2(\theta)}}d\theta \tag {10}$
利用二项式展开可得到(Ref 1)
$\begin{aligned}K(m)&=\int_0^{\pi/2}\sum_{j=0}^\infty \begin{pmatrix}-1/2\\j\end{pmatrix}(-1)^j(m^2\sin^2 \theta)^j \\&={\pi\over 2}{}_2F_1({1\over 2}, {1\over 2};1;k^2), \quad (\textrm{also Ref 4, } 8.113)\end{aligned} \tag {11}$
$\begin{aligned}E(m)&=\int_0^{\pi/2}\sum_{j=0}^\infty \begin{pmatrix}1/2\\j\end{pmatrix}(-1)^j(m^2\sin^2 \theta)^j \\&={\pi\over 2}{}_2F_1(-{1\over 2}, {1\over 2};1;k^2),\quad (\textrm{also Ref 4, } 8.114)\end{aligned} \tag {12}$
这里 $_2F_1(\ ,\ ;\ ;\ )$ 是超几何函数，涉及到了The Pochhammer symbol （Ref 2），利用下面的公式可计算负数的对应结果
$(x)_n={\Gamma(x+n)\over\Gamma(x)} \tag {13}$ $\Gamma(x)={\Gamma(x+1)\over x} \tag {14}$ $(-x)_n=(-1)^n(x-n+1)_n \tag {15}$
$(-1/2)_n={\Gamma(-1/2+n)\over\Gamma(-1/2)} ={(-1/2)\Gamma(-1/2+n)\over\Gamma(1/2)}={-(2n-3)!!\over 2^n}$
在Ref 3 p170 习题12，利用椭圆积分计算窄带高斯信号包络的相关。

Ref

https://archive.lib.msu.edu/crcmath/math/math/e/e078.htm
https://mathworld.wolfram.com/PochhammerSymbol.html
W. B. Davenport, Jr., and W. L. Root, An Introduction to the Theory of Randum Signals and Noise, McGraw-Hill, New York, 1958.
Table of Integrals, Series, and Products