作换元 `t = pi/2 - x`, 显然有 `int_0^(pi/2) f(sin t) dt = int_0^(pi/2) f(cos x) dx`, 这一结论在下文直接应用, 不再说明.
第一类椭圆积分 完全椭圆积分: `K(k) := int_0^(pi/2) ("d"theta)/sqrt(1-k^2 sin^2 theta)`, `quad "Re"(k) lt 1`. (Legendre 形式) 其中变元 `k` 称为模. 令 `t = sin theta` 得到: `K(k) = int_0^1 dt/(sqrt(1-t^2) sqrt(1-k^2 t^2))`. (Jacobi 形式) 利用下文的 Landen 变换得到: `K(z) = int_0^(pi/2) ("d"theta)/sqrt(1+z^2 +- 2z cos 2 theta)`. (第二形式) 不完全椭圆积分 (Legendre 积分): `F(k, varphi) := int_0^varphi ("d"theta)/sqrt(1-k^2 sin^2 theta)`.
Euler 变换 `K("i"x) = 1/sqrt(1+x^2) K(sqrt(x^2/(1+x^2)))`.
右边 `= 1/sqrt(1+x^2) int_0^(pi/2) ("d"theta)/sqrt(1-x^2/(1+x^2) cos^2 theta)` `= int_0^(pi/2) ("d"theta)/sqrt(1+x^2-x^2 cos^2 theta)` `= int_0^(pi/2) ("d"theta)/sqrt(1+x^2 sin^2 theta)` `=` 左边.
Landen 变换 取共轭模 `k' = sqrt(1-k^2)`. 令 `z = (1-k')/(1+k')`, 则 `K(z) = (1+k')/2 K(k)` `= 1/(1+z) K((2sqrt z)/(1+z))`.
令
`t = (1+k') (s sqrt(1-s^2))/sqrt(1-k^2 s^2)`,
求微分,
`dt = (1+k')/sqrt(1-s^2)
(1-2s^2 + k^2 s^4)/(1-k^2 s^2)^(3//2) "d"s`,
从 Jacobi 形式出发, 令 `z = (1-k')/(1+k')`,
由
`1-k^2 s^2 - (1+-k')^2 s^2 (1-s^2)`
`= 1 - k^2 s^2 + (s^4 - s^2) (1 + {:k':}^2 +- 2 k')`
`= 1 - 2(1+-k') s^2 + (1+-k')^2 s^4`
`= [1 - (1+-k') s^2]^2`,
有
`sqrt(1-t^2) = (1 - (1+k') s^2)/sqrt(1-k^2 s^2)`,
`quad sqrt(1-z^2 t^2) = (1 - (1-k') s^2)/sqrt(1-k^2 s^2)`.
当 `t = 0` 时, 相应的 `s = 0` 或 `1`;
当 `t = 1` 时, 解二次方程知, 相应的 `s = 1/sqrt(1+k')`.
`K(z)`
`= int_0^1 dt/(sqrt(1-t^2) sqrt(1-z^2 t^2))`
`= (1+k') int_0^(1//sqrt(1+k')) ("d"s)/(sqrt(1-s^2) sqrt(1-k^2 s^2))`
`= (1+k') int_(1//sqrt(1+k'))^1 ("d"s)/(sqrt(1-s^2) sqrt(1-k^2 s^2))`
`= (1+k')/2 int_0^1 ("d"s)/(sqrt(1-s^2) sqrt(1-k^2 s^2))`.
Gauss 公式 定义 `M(a, b) := int_0^(pi/2)("d"theta)/sqrt(a^2 cos^2 theta+b^2 sin^2 theta)`, 则 `M((a+b)/2, sqrt(a b)) = 2/(a+b) K((a-b)/(a+b)) = M(a, b)`. 取 `a, b = 1 +- k`: `M(1, k') = M(1, sqrt(1-k^2)) = K(k) = M(1+k, 1-k)`.
首先 `a^2 cos^2 theta + b^2 sin^2 theta` `= a^2 - (a^2-b^2) sin^2 theta`, 故 `((a+b)/2)^2 cos^2 theta + a b sin^2 theta` `= ((a+b)/2)^2 - ((a-b)/2)^2 sin^2 theta`. 于是 `M((a+b)/2, sqrt(a b))` `= int_0^(pi/2) ("d"theta) / sqrt(((a+b)/2)^2 - ((a-b)/2)^2 sin^2 theta)` `= 2/(a+b) K((a-b)/(a+b))`. 这就证明了第一个等号. 接下来利用 `K(z)` 的第二形式, `2/(a+b) K((a-b)/(a+b))` `= 1/(a+b) int_0^pi [1 + ((a-b)/(a+b))^2 - 2(a-b)/(a+b) cos theta]^(-1/2) "d"theta` `= int_0^pi [(a+b)^2 + (a-b)^2 - 2(a-b)(a+b) cos theta]^(-1/2) "d"theta` `= int_0^pi [a^2(1-cos theta) + b^2(1+cos theta)]^(-1/2) ("d"theta)/sqrt 2` `= int_0^pi (a^2 sin^2{:theta/2:} + b^2 cos^2{:theta/2:})^(-1/2) ("d"theta)/2` `= M(a, b)`. 这就证明了第二个等号.
[来自 菲赫金哥尔茨 第二卷第 315 目] 令
`sin theta = (2 a sin varphi)/((a+b) + (a-b)sin^2 varphi)`,
注意到 `varphi` 从 `0` 到 `pi//2` 时, `theta` 也在同一区间变化.
`cos^2 theta = 1 - sin^2 theta`
`= 1/()^2 [(a+b)^2 + (a-b)^2 sin^4 varphi + 2(a^2-b^2)sin^2 varphi
- 4 a^2 sin^2 varphi]`
`= 1/()^2 [(a+b)^2 (cos^2 varphi + sin^2 varphi):}
`quad + {:(a-b)^2 sin^2 varphi (1 - cos^2 varphi) - 2(a^2+b^2) sin^2 varphi]`
`= (cos^2 varphi)/()^2 [(a+b)^2 - (a-b)^2 sin^2 varphi]`,
我们得到
`cos theta = sqrt((a+b)^2 - (a-b)^2 sin^2 varphi)/((a+b) + (a-b) sin^2
varphi) cos varphi`.
另一方面,
`a^2 cos^2 theta + b^2 sin^2 theta`
`= a^2/()^2 (((a+b)^2 - (a-b)^2 sin^2 varphi) cos^2 varphi + 4 b^2
sin^2 varphi)`
`= a^2/()^2 ((a+b)^2 (1-sin^2 varphi) - (a-b)^2 sin^2 varphi (1-sin^2
varphi) 1 + 4 b^2 sin^2 varphi)`
`= a^2/()^2 ((a+b)^2 - 2 (a^2-b^2) sin^2 varphi + (a-b)^2 sin^4
varphi)`
被积函数的分母可以有理化:
`sqrt(a^2 cos^2 theta + b^2 sin^2 theta)`
`= a ((a+b) - (a-b) sin^2 varphi)/((a+b) + (a-b) sin^2 varphi`.
最后由 微分得
`cos theta "d"theta`
`= 2a cos varphi ((a+b) - (a-b) sin^2 varphi)/((a+b) + (a-b) sin^2
varphi)^2 "d"varphi`,
于是
`("d"theta)/sqrt(a^2 cos^2 theta + b^2 sin^2 theta)`
`= (2 "d"varphi)/sqrt((a+b)^2 - (a-b)^2 sin^2 varphi)`
`= ("d"varphi)/sqrt(((a+b)/2)^2 - ((a-b)/2)^2 sin^2 varphi)`
`= ("d"varphi)/sqrt(((a+b)/2)^2 cos^2 varphi + a b sin^2 varphi)`.
即 `M(a, b) = M((a+b)/2, sqrt(a b))`.
(Gauss, 1799) 设数列 `{a_n}`, `{b_n}` 满足: `0 lt b_0 lt a_0`, 且 `a_(n+1) = (a_n + b_n)/2`, `quad b_(n+1) = sqrt(a_n b_n)`, `quad n ge 0`, 则两数列收敛到同一极限. 此极限称为 `a_0` 与 `b_0` 的算术几何平均值 `"agm"(a_0, b_0)`, 我们有 `"agm"(a_0, b_0) M(a_0, b_0) = pi/2`. 从而得到高效求 `K(k)` 的公式 `K(k) = pi/2 1/("agm"(1, sqrt(1-k^2)))`.
`K(k) = (1 + k_1) K(k_1)`, 其中 `k_1 = (1-k')/(1+k')`, `k' = sqrt(1-k^2)`. 反复应用此公式得 `K(k) = pi/2 prod_(n ge 1) (1+k_n)`.
Gauss 常数 `G = "agm"(1, sqrt 2)^-1 = (Gamma^2(1/4))/(2pi)^(3/2)`.
`pi/2 G` `= pi/2 "agm"(1, sqrt 2)^-1` `= K("i")` `= int_0^1 dt/(sqrt(1-t^2) sqrt(1+t^2))` `= 1/4 B(1/4, 1/2)` `= (sqrt pi Gamma(1/4)^2 sin{:pi/4:})/(4 pi)`.
一道经典例题: `K(z)` 与单摆周期.
第二类椭圆积分 `E(z) = int_0^(pi/2) sqrt(1-z^2 sin^2 theta) "d"theta` `= int_0^1 sqrt((1-z^2 t^2)/(1-t^2)) dt`.
Euler 变换 `E("i"x) = sqrt(1+x^2) E(sqrt(x^2/(1+x^2)))`.
椭圆的周长 记椭圆的离心率为 `k = sqrt(a^2-b^2) // a`, 则椭圆周长 `L = 4 int_0^(pi/2) sqrt(a^2 cos^2 theta + b^2 sin^2 theta) "d"theta` `= 4 a int_0^(pi/2) sqrt (1 - k^2 sin^2 theta) "d"theta` `= 4 a E(k)`.
求 `int_0^1 sqrt((1+x^2)/(1-x^2)) dx`.
原式等于
`int_0^1 (1+x^2)/sqrt(1-x^4) dx`
`= int_0^1 1/sqrt(1-x^4) dx + int_0^1 x^2/sqrt(1-x^4) dx`
`= I_1 + I_2`.
令 `t = x^4`, `dt = 4x^3 dx`, `dx = dt/(4 t^(3//4))`,
`I_1 = 1/4 int_0^1 t^(-3//4)(1-t)^(-1//2) dt = 1/4 B(1/4, 1/2)`.
类似 `I_2 = 1/4 B(3/4, 1/2)`.
另一方面, 令 `x = sin theta`, 原式等于
`int_0^(pi/2) sqrt(1+sin^2 theta) "d"theta`
`= E("i") = sqrt 2 E(1/sqrt 2)`.
`E` 和 `K` 的导数 [来自 菲赫金哥尔茨 第二册511目] `("d"E)/("d"k) = (E-K)/k`, `quad ("d"K)/("d"k) = (E//{:k':}^2 - K)/k`. 其中共轭模 `k'` 满足 `k^2 + {:k':}^2 = 1`. 进一步可得 `E` 满足的微分方程 `("d"^2 E)/{:"d"k:}^2 + 1/k ("d"E)/("d"k) + E/{:k':}^2 = 0`.
直接计算: `("d"E)/("d"k) = "d"/("d"k) int_0^(pi/2) sqrt(1-k^2 sin^2 t) dt` `= int_0^(pi/2) (-k sin^2 t)/sqrt(1-k^2 sin^2 t) dt` `= (E-K)/k`. `K` 的导数则不太简单: 先验证等式 `"d"/dt (sin t cos t)/sqrt(1-k^2 sin^2 t)` `= 1/k^2 [(1-k^2 sin^2 t)^(1/2) - (1-k^2)/k^2 (1-k^2 sin^2 t)^(-3/2)]`, 于是 `int_0^(pi/2) (1-k^2 sin^2 t)^(-3/2) dt` `= E//{:k':}^2`. 从而 `("d"K)/("d"k) = "d"/("d"k) int_0^(pi/2) dt/sqrt(1-k^2 sin^2 t)` `= int_0^(pi/2) k sin^2 t (1-k^2 sin^2 t)^(-3/2) dt` `= 1/k int_0^(pi/2) [(1-k^2 sin^2 t)^(-3/2) - (1-k^2 sin^2 t)^(-1/2)] dt` `= (E//{:k':}^2 - K)/k`.
Legendre 关系式 [菲赫金哥尔茨 第二册 511 目 12)] `E K' + E' K - K K' = pi/2`. 其中 `K = K(k)`, `K' = K(k')`, `E = E(k)`, `E' = E(k')`.
由 `("d"k')/("d"k) = "d"/("d"k) sqrt(1-k^2) = -k/sqrt(1-k^2) = -k/k'`, 将 `("d"E)/("d"k) = (E-K)/k`, `("d"K)/("d"k) = (E//{:k':}^2-K)/k`, `("d"E')/("d"k) = (-k/k')(E'-K')/k'`, `("d"K')/("d"k) = (-k/k')(E'//k^2 - K')/k'` 代入整理得 `"d"/("d"k) (E K' + E' K - K K') = 0`. 即 `E K' + E' K - K K'` 等于一个常数 `c`. 现在令 `k to 0`, `k' to 1` 求极限, 来求出 `c`. 事实上 `lim_(k to 0) E' K = 1 * pi/2 = pi/2`, 又 `|K'(E - K)|` `= |int_0^(pi/2) dt/sqrt(1-{:k':}^2 sin^2 t)| |int_0^(pi/2) (1/sqrt(1-k^2 sin^2 t) - sqrt(1-k^2 sin^2 t)) dt|` `le pi/(2 k) * |int_0^(pi/2) (k^2 sin^2 t dt)/sqrt(1-k^2 sin^2 t)|` `le pi/(2 k) * pi/2 * k^2/(k') to 0`, 所以极限为 `pi/2`.
只证第一式. 利用幂级数 `1/sqrt(1-z) = sum_(n ge 0) ((2n-1)!!)/((2n)!!) z^n`, `quad |z| lt 1` 和积分 `int_0^(pi/2) sin^(2n) theta "d"theta = pi/2 ((2n-1)!!)/((2n)!!)` 对 `K(z)` 逐项积分即可.
第三类椭圆积分 `int dz/((1+h z^2) sqrt((1-z^2)(1-k^2 z^2)))` `= int ("d"varphi)/((1 + h sin^2 varphi)sqrt(1 - k^2 sin^2 varphi))`.
考虑 `u(varphi) = F(k, varphi)` `= int_0^varphi ("d"theta)/sqrt(1-k^2 sin^2 theta)` 的反演. 其中 `0 lt k lt 1` 为常数. 由于被积函数恒为正, `u(varphi)` 是从 `RR` 到 `RR` 的单射, 将它的反函数记为 `varphi = "am" u`, 称为振幅 (amplitude) 函数. `sin "am" u := "sn" u` 称为辐角正弦或椭圆正弦. 这些名称来自单摆 (数学摆) 的振幅公式 `sin {:theta/2:} = sin{:theta_0/2:} "sn" sqrt(g/l) t`, 其中 `theta` 为 `t` 时刻与竖直方向的夹角, `theta_0` 为初始角, `l` 为摆长, `g` 为重力加速度.
超几何函数 (超几何级数) 定义为 `{::}_2 F_1(a,b";"c";"z)` `= sum_(n ge 0) (a^(bar n) b^(bar n))/c^(bar n) z^n/n!` `= 1 + (a b)/c z + (a(a+1)b(b+1))/(c(c+1)) z^2/(2!) + cdots`. 它可以用来表示许多初等/非初等的函数. 观察得到: 超几何函数中, `a, b` 的地位对称; 当 `a` 或 `b` 为 `0` 时, 函数值为 `1`; 当 `a` 或 `b` 为负整数时, 级数只有有限项.