椭圆积分及变换关系

吴炳洋

已于 2022-03-27 20:02:40 修改

阅读量2.5k

点赞数

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

于 2022-03-27 19:59:44 首次发布

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

版权

数学基础专栏收录该内容

20 篇文章 1 订阅

订阅专栏

一种一般化的椭圆积分可以表示为以下函数：
$\begin{aligned} T(a,b)&\equiv {2\over\pi} \int^{\pi/2}_0 {d\theta\over\sqrt{a^2\cos^2\theta+b^2\sin^2\theta}}\\ &= {2\over\pi} \int^{\pi/2}_0 {d\theta\over\cos\theta\sqrt{a^2+b^2\tan^2\theta}} \end{aligned} \tag 1$
使 $t\equiv b\tan\theta,\quad dt=b\sec^2\theta\,d\theta={d\theta \over \cos \theta }\sqrt{b^2+t^2}.$ 得到 $\begin{aligned} T(a,b)&=&{2\over \pi} \int_0^\infty {dt\over \sqrt{(a^2+t^2)(b^2+t^2)}}\\ &=&{1\over \pi} \int_{-\infty}^\infty {dt\over \sqrt{(a^2+t^2)(b^2+t^2)}}. \end{aligned}\tag 2$
进一步进行变量代换 $u\equiv {\textstyle{1\over 2}}(t-ab/t), \quad du={\textstyle{1\over 2}}(1+ab/t^2)\,dt$ 得到
$\begin{aligned} T(a,b)&=& {2\over \pi} \int_{-\infty}^\infty {du\over \sqrt{[4u^2+(a+b)^2](u^2+ab)}}\\ &=& {1\over \pi} \int_{-\infty}^\infty {du\over \sqrt{\left[{u^2+\left({a+b\over 2}\right)^2}\right](u^2+ab)}}. \end{aligned} \tag 3$
因此 $T({\textstyle{1\over 2}}(a+b),\sqrt{ab}). \tag 4$ 以算术平均和几何平均代替原变量得到相同的结果。

完全椭圆积分可以利用此算术平均和几何平均变换方法进行计算
$\begin{aligned} T(a,b)&\equiv& {2\over \pi} \int^{\pi/2}_0 {d\theta\over\sqrt{a^2\cos^2\theta+b^2\sin^2\theta}}\\ &=& {2\over\pi}\int^{\pi/2}_0 {d\theta \over a\sqrt{\cos ^2\theta +\left({b\over a}\right)^2 \sin^2 \theta}}\\ &=&{2\over a\pi}\int^{\pi/2}_0 {d\theta \over \sqrt{1-\left({1-{b^2\over a^2}}\right)^2\sin^2 \theta}}. \end{aligned} \tag 5$
因此
${2\over a\pi}K\left({1-{b^2\over a^2}}\right) \tag 6$

$K (k)$ 是第一类完全椭圆积分，令 $a\equiv a_0\equiv 1,\quad b\equiv b_0\equiv k', \quad k=\sqrt{1-k'^2},$ 等式 (6) 为
${2\over \pi}K(\sqrt{1-k'^2}\,) = {2\over \pi}K(k),$
定义
$K'(k)\equiv K(\sqrt{1-k^2}\,)=K(k'). \tag 7$
应用变量算术平均和几何平均进行变量代换
$\begin{aligned} {1\over a}K\left({\sqrt{1-{b^2\over a^2}}\,}\right)&=& {2\over a+b}K\left({\sqrt{1-{4ab\over (a+b)^2}}\,}\right)\\ &=& {2\over a+b} K\left({\sqrt{a^2+b^2-2ab\over(a+b)^2}\,}\right)\\ &=& {2\over a+b}K\left({a-b\over a+b}\right) \end{aligned} \tag 8$
得到
$K\left({\sqrt{1-{b^2\over a^2}}\,}\right)= {2\over 1+{b\over a}}K\left({1-{b\over a}\over 1+{b\over a}}\right). \tag 9$ 令
$k'\equiv {b\over a}, \quad k\equiv \sqrt{1-k'^2},$ 得到
$K(k)={2\over 1+k'} K\left({1-k'\over 1+k'}\right). \tag {10}$
令
$l\equiv (1-k')/(1+k'),$ 得 $k'={1-l\over 1+l}, \quad k=\sqrt{1-k'^2}= {2\sqrt{l}\over 1+l},$ 因此 $K(k)={1\over k+1} K\left({2\sqrt{k}\over 1+k}\right). \tag {11}$ 类似地可以得到
$E(k)={1+k\over 2}E\left({2\sqrt{k}\over 1+k}\right)+{k'^2\over 2}K(k) \tag {12}$ $E(k)=(1+k')E\left({1-k'\over 1+k'}\right)-k'K(k). \tag {13}$ 这里 $E (k)$ 为第二类完全椭圆积分。
$\begin{aligned} K'(k)&=&K(k') = {2\over 1+k} K\left({1-k\over 1+k}\right)\\ &=& {2\over 1+k}K'\left({\sqrt{1-\left({1-k\over 1+k}\right)^2}\,}\right)\\ &=& {2\over 1+k}K'\left({2\sqrt{k}\over 1+k}\right) \end{aligned} \tag {14}$
$K'(k)={1\over 1+k'} K\left({2\sqrt{k'}\over 1+k'}\right)= {1\over 1+k'} K'\left({1-k'\over 1+k'}\right), \tag {15}$
$E'(k)=(1+k)E'\left({2\sqrt{k}\over 1+k}\right)-kK'(k) \tag {16}$
$E'(k)=\left({1+k'\over 2}\right)E'\left({1-k'\over 1+k'}\right)+{k^2\over 2}K'(k).\tag {17}$
${K'(k)\over K(k)}=2 {K'\left({2\sqrt{k}\,\over 1+k}\right)\over K\left({2\sqrt{k}\,\over 1+k}\right)} = {1\over 2} {K'\left({1-k'\over 1+k'}\right)\over K\left({1-k'\over 1+k'}\right)} \tag {18}$