主題
Search

橢圓積分

橢圓積分是積分形式

int(A(x)+B(x)sqrt(S(x)))/(C(x)+D(x)sqrt(S(x)))dx,
(1)

或

int(A(x)dx)/(B(x)sqrt(S(x))),
(2)

其中 A(x), B(x), C(x), 和 D(x) 是關於 x 的多項式，並且 S(x) 是 3 次或 4 次多項式。更簡單地說，橢圓積分是形式的積分

intR(w,x)dx,
(3)

其中 R(w,x) 是 x 和 w 的有理函式，w^2 是關於 x 的函式，它是關於 x 的三次或四次函式，R(w,x) 至少包含一個 w 的奇次冪，並且 w^2 沒有重複因子（Abramowitz 和 Stegun 1972, p. 589）。

橢圓積分可以被看作是反三角函式的推廣，併為更廣泛的問題類別提供解決方案。例如，雖然圓的弧長可以表示為引數的簡單函式，但計算橢圓的弧長需要使用橢圓積分。類似地，對於小角度振盪，擺的位置由作為時間函式的三角函式給出，但是對於任意大位移的完整解需要使用橢圓積分。電磁學和引力中的許多其他問題都可以透過橢圓積分來解決。

透過反演橢圓積分以獲得三角函式的推廣，可以得到一類非常有用的函式，稱為橢圓函式。橢圓函式（其中雅可比橢圓函式和魏爾斯特拉斯橢圓函式是兩種最常見的形式）為分析數論以及數學其他領域的許多深刻問題提供了強大的工具。

所有橢圓積分都可以用三種“標準”型別來表示。為了看到這一點，寫出

R(w,x)=(P(w,x))/(Q(w,x))
(4)
=(wP(w,x)Q(-w,x))/(wQ(w,x)Q(-w,x)).
(5)

但是由於 w^2=f(x),

Q(w,x)Q(-w,x)=Q_1(w,x)
(6)
=Q_1(-w,x),
(7)

那麼

wP(w,x)Q(-w,x)=A+Bx+Cw+Dx^2+Ewx+Fw^2+Gw^2x+Hw^3x
(8)
=(A+Bx+Dx^2+Fw^2+Gw^2x)+w(c+Ex+Hw^2x+...)
(9)
=P_1(x)+wP_2(x),
(10)

所以

R(w,x)=(P_1(x)+wP_2(x))/(wQ_1(w))
(11)
=(R_1(x))/w+R_2(x).
(12)

但是任何函式 intR_2(x)dx 都可以用初等函式來評估，所以唯一需要考慮的部分是

int(R_1(x))/wdx.
(13)

現在，任何四次式都可以表示為 S_1S_2，其中

S_1=a_1x^2+2b_1x+c_1
(14)
S_2=a_2x^2+2b_2x+c_2.
(15)

這裡的係數是實數，因為成對的複數根是共軛複數

[x-(R+Ii)][x-(R-Ii)]=x^2+x(-R+Ii-R-Ii)+(R^2-I^2i)
(16)
=x^2-2Rx+(R^2+I^2).
(17)

如果所有四個根都是實數，則必須將它們排列成不交錯的形式（Whittaker 和 Watson 1990, p. 514）。現在定義一個量 lambda，使得 S_1-lambdaS_2

(a_1-lambdaa_2)x^2-(2b_1-2b_2lambda)x+(c_1-lambdac_2)
(18)

是一個平方數，並且

2sqrt((a_1-lambdaa_2)(c_1-lambdac_2))=2(b_1-b_2lambda)
(19)
(a_1-lambdaa_2)(c_1-lambdac_2)-(b_1-lambdab_2)^2=0.
(20)

將此方程的根稱為 lambda_1 和 lambda_2，那麼

S_1-lambda_2S_2=[sqrt((a_1-lambda_2a_2)x^2)+sqrt(c_1-lambda_2c_2)]^2
(21)
=(a_1-lambda_2a_2)(x+sqrt((c_1-lambda_2c_2)/(a_1-lambda_2a_2)))
(22)
=(a_1-lambda_2a_2)(x-beta)^2
(23)
S_1-lambda_1S_2=[sqrt((a_1-lambda_1a_2)x^2)+sqrt(c_1-lambda_1c_2)]^2
(24)
=(a_1-lambda_1a_2)(x+sqrt((c_1-lambda_1c_2)/(a_1-lambda_1a_2)))
(25)
=(a_1-lambda_1a_2)(x-alpha)^2.
(26)

取 (25)-(26) 和 lambda_2(1)-lambda_1(2) 得到

S_2(lambda_2-lambda_1)=(a_1-lambda_1a_2)(x-alpha)^2-(a_1-lambda_2a_2)(x-beta)^2
(27)
S_1(lambda_2-lambda_1)=lambda_2(a_1-lambda_1a_2)(x-alpha)^2-lambda_1(a_1-lambda_2a_2)(x-beta)^2.
(28)

求解得到

S_1=(a_1-lambda_1a_2)/(lambda_2-lambda_1)(x-alpha)^2-(a_1-lambda_2a_2)/(lambda_2-lambda_1)(x-beta)^2
(29)
=A_1(x-alpha)^2+B_1(x-beta)^2
(30)
S_2=(lambda_2(a_1-lambda_1a_2))/(lambda_2-lambda_1)(x-alpha)^2-(lambda_1(a_1-lambda_2a_2))/(lambda_2-lambda_1)(x-beta)^2
(31)
=A_2(x-alpha)^2+B_2(x-beta)^2,
(32)

所以我們有

w^2=S_1S_2=[A_1(x-alpha)^2+B_1(x-beta)^2][A^2(x-alpha)^2+B^2(x-beta)^2].
(33)

現在令

t=(x-alpha)/(x-beta)
(34)
dt=[(x-beta)^(-1)-(x-alpha)(x-beta)^(-2)]dx
(35)
=((x-beta)-(x-alpha))/((x-beta)^2)dx
(36)
=(alpha-beta)/((x-beta)^2)dx,
(37)

所以

w^2=(x-beta)^4[A_1((x-alpha)/(x-beta))^2+B_1][A_2((x-alpha)/(x-beta))+B_2]
(38)
=(x-beta)^4(A_1t^2+B_1)(A_2t^2+B_2),
(39)

和

w=(x-beta)^2sqrt((A_1t^2+B_1)(A_2t^2+B_2))
(40)
(dx)/w=[((x-beta)^2)/(alpha-beta)dt]1/((x-beta)^2sqrt((A_1t^2+B_1)(A_2t^2+B_2)))
(41)
=(dt)/((alpha-beta)sqrt((A_1t^2+B_1)(A_2t^2+B_2))).
(42)

現在令

R_3(t)=(R_1(x))/(alpha-beta),
(43)

所以

int(R_1(x)dx)/w=int(R_3(t)dt)/(sqrt((A_1t^2+B_1)(A_2t^2+B_2))).
(44)

重寫偶部和奇部

R_3(t)+R_3(-t)=2R_4(t^2)
(45)
R_3(t)-R_3(-t)=2tR_5(t^2),
(46)

得到

R_3(t)=1/2(R_(even)-R_(odd))
(47)
=R_4(t^2)+tR_5(t^2),
(48)

所以我們有

int(R_1(x)dx)/w=int(R_4(t^2)dt)/(sqrt((A_1t^2+B_1)(A_2t^2+B_2)))+int(R_5(t^2)tdt)/(sqrt((A_1t^2+B_1)(A_2t^2+B_2))).
(49)

令

u=t^2
(50)
du=2tdt
(51)

將第二個積分簡化為

1/2int(R_5(u)du)/(sqrt((A_1u+B_1)(A_2u+B_2))),
(52)

可以使用初等函式來評估它。然後，第一個積分可以透過分部積分法簡化為三種勒讓德橢圓積分之一（也稱為勒讓德-雅可比橢圓積分），稱為不完全第一類橢圓積分、第二類橢圓積分和第三類橢圓積分，分別表示為 F(phi,k)、E(phi,k) 和 Pi(n;phi,k)（von Kármán 和 Biot 1940, Whittaker 和 Watson 1990, p. 515）。如果 phi=pi/2，則這些積分稱為完全橢圓積分，表示為 K(k)、E(k)、Pi(n;k)。

不完全橢圓積分使用橢圓模 k、引數 m=k^2 或模角 alpha=sin^(-1)k 來表示。當使用引數時，橢圓積分寫為 I(phi|m)，當使用橢圓模時，寫為 I(phi,k)，當使用模角時，寫為 I(phi\alpha)。當 phi=pi/2 時，定義完全橢圓積分，並且可以使用展開式表示

(1-k^2sin^2theta)^(-1/2)=sum_(n=0)^infty((2n-1)!!)/((2n)!!)k^(2n)sin^(2n)theta.
(53)

標準形式的橢圓積分

int_a^x(dx)/(sqrt(f(x))),
(54)

其中

f(x)=a_4x^4+a_3x^3+a_2x^2+a_1x+a_0,
(55)

可以用魏爾斯特拉斯橢圓函式和不變數進行解析計算 (Whittaker 和 Watson 1990, p. 453)

g_2=a_0a_4-4a_1a_3+3a_2^2
(56)
g_3=a_0a_2a_4-2a_1a_2a_3-a_4a_1^2-a_3^2a_0.
(57)

如果 a=x_0 是 f(x)=0 的根，那麼解是

x=x_0+1/4f^'(x_0)[P(z;g_2,g_3)-1/(24)f^('')(x_0)]^(-1).
(58)

對於任意下限，

x=a+(sqrt(f(a))P^'(z)+1/2f^'(a)[P(z)-1/(24)f^('')(a)]+1/(24)f(a)f^(''')(a))/(2[P(z)-1/(24)f^('')(a)]^2-1/(48)f(a)f^((iv))(a)),
(59)

其中 P(z)=P(z;g_2,g_3) 是魏爾斯特拉斯橢圓函式 (Whittaker 和 Watson 1990, p. 454)。

廣義橢圓積分可以透過函式定義

T(a,b)=2/piint_0^(pi/2)(dtheta)/(sqrt(a^2cos^2theta+b^2sin^2theta))
(60)
=2/piint_0^(pi/2)(dtheta)/(costhetasqrt(a^2+b^2tan^2theta))
(61)

(Borwein 和 Borwein 1987)。現在令

t=btantheta
(62)
dt=bsec^2thetadtheta.
(63)

但是

sectheta=sqrt(1+tan^2theta),
(64)

所以

dt=b/(costheta)secthetadtheta
(65)
=b/(costheta)sqrt(1+tan^2theta)dtheta
(66)
=b/(costheta)sqrt(1+(t/b)^2)dtheta
(67)
=(dtheta)/(costheta)sqrt(b^2+t^2),
(68)

和

(dtheta)/(costheta)=(dt)/(sqrt(b^2+t^2)),
(69)

方程變為

T(a,b)=2/piint_0^infty(dt)/(sqrt((a^2+t^2)(b^2+t^2)))
(70)
=1/piint_(-infty)^infty(dt)/(sqrt((a^2+t^2)(b^2+t^2))).
(71)

現在我們進行進一步替換 u=1/2(t-ab/t)。微分變為

du=1/2(1+ab/t^2)dt,
(72)

但是 2u=t-ab/t，所以

2u/t=1-ab/t^2
(73)
ab/t^2=1-2u/t
(74)

和

1+ab/t^2=2-2u/t=2(1-u/t).
(75)

但是，左側始終為正，因此

1+ab/t^2=2-2u/t=2|1-u/t|
(76)

微分是

dt=(du)/(|1-u/t|).
(77)

我們需要注意積分的極限。將 (◇) 寫成

int_(-infty)^inftyf(t)dt=int_(-infty)^(0^-)f(t)dt+int_(0^+)^inftyf(t)dt.
(78)

現在將極限更改為適合 u 積分的極限

int_(-infty)^inftyg(u)du+int_(-infty)^inftyg(u)du=2int_(-infty)^inftyg(u)du,
(79)

所以我們提取了一個因子 2，必須包含在內。使用這個事實並將 (◇) 代入 (◇) 因此得到

T(a,b)=2/piint_(-infty)^infty(du)/(|1-u/t|sqrt(a^2b^2+(a^2+b^2)t^2+t^4)).
(80)

現在注意

u^2=(t^4-2abt^2+a^2b^2)/(4t^2)
(81)
4u^2t^2=t^4-2abt^2+a^2b^2
(82)
a^2b^2+t^4=4u^2t^2+2abt^2.
(83)

將 (◇) 代入 (◇) 得到

T(a,b)=2/piint_(-infty)^infty(du)/(|1-u/t|sqrt(4u^2t^2+2abt^2+(a^2+b^2)t^2))
(84)
=2/piint_(-infty)^infty(du)/(|t-u|sqrt(4u^2+(a+b)^2)).
(85)

但是

2ut=t^2-ab
(86)
t^2-2ut-ab=0
(87)
t=1/2(2u+/-sqrt(4u^2+4ab))
(88)
=u+/-sqrt(u^2+ab),
(89)

所以

t-u=+/-sqrt(u^2+ab),
(90)

並且 (◇) 變為

T(a,b)=2/piint_(-infty)^infty(du)/(sqrt([4u^2+(a+b)^2](u^2+ab)))
(91)
=1/piint_(-infty)^infty(du)/(sqrt([u^2+((a+b)/2)^2](u^2+ab))).
(92)

因此，我們證明了

T(a,b)=T(1/2(a+b),sqrt(ab)).
(93)

我們可以迭代

a_(i+1)=1/2(a_i+b_i)
(94)
b_(i+1)=sqrt(a_ib_i),
(95)

任意多次，而不會更改積分的值。但是此迭代與算術-幾何平均值相同並且因此收斂到它，因此迭代在 a_i=b_i=M(a_0,b_0) 處終止，並且我們有

T(a_0,b_0)=T(M(a_0,b_0),M(a_0,b_0))
(96)
=1/piint_(-infty)^infty(dt)/(M^2(a_0,b_0)+t^2)
(97)
=1/(piM(a_0,b_0))[tan^(-1)(t/(M(a_0,b_0)))]_(-infty)^infty
(98)
=1/(piM(a_0,b_0))[pi/2-(-pi/2)]
(99)
=1/(M(a_0,b_0)).
(100)

完全橢圓積分出現在求橢圓的弧長和擺的週期中。它們也自然地從 theta 函數理論中產生。可以使用涉及算術-幾何平均值的過程來計算完全橢圓積分。請注意

T(a,b)=2/piint_0^(pi/2)(dtheta)/(sqrt(a^2cos^2theta+b^2sin^2theta))
(101)
=2/piint_0^(pi/2)(dtheta)/(asqrt(cos^2theta+(b/a)^2sin^2theta))
(102)
=2/(api)int_0^(pi/2)(dtheta)/(sqrt(1-(1-(b^2)/(a^2))sin^2theta)).
(103)

所以我們有

T(a,b)=2/(api)K(sqrt(1-(b^2)/(a^2)))
(104)
=1/(M(a,b)),
(105)

其中 K(k) 是第一類完全橢圓積分。我們可以自由地令 a=a_0=1 和 b=b_0=k^'，所以

2/piK(sqrt(1-k^('2)))=2/piK(k)=1/(M(1,k^')),
(106)

由於 k=sqrt(1-k^('2))，所以

K(k)=pi/(2M(1,k^')).
(107)

但是算術-幾何平均值由下式定義

a_i=1/2(a_(i-1)+b_(i-1))
(108)
b_i=sqrt(a_(i-1)b_(i-1))
(109)
c_i={1/2(a_(i-1)-b_(i-1)) i>0; sqrt(a_0^2-b_0^2) i=0,
(110)

其中

c_(n-1)=1/2a_n-b_n=(c_n^2)/(4a_(n+1))<=(c_n^2)/(4M(a_0,b_0)),
(111)

所以我們有

K(k)=pi/(2a_N),
(112)

其中 a_N 是 a_n 收斂到的值。類似地，取 a_0^'=1 和 b_0^'=k 得到

K^'(k)=pi/(2a_N^').
(113)

Borwein 和 Borwein (1987) 也表明，定義

U(a,b)=pi/2int_0^(pi/2)sqrt(a^2cos^2theta+b^2sin^2theta)dtheta
(114)
=aE^'(b/a)
(115)

導致

2U(a_(n+1),b_(n+1))-U(a_n,b_n)=a_nb_nT(a_n,b_n),
(116)

所以

(K(k)-E(k))/(K(k))=1/2(c_0^2+2c_1^2+2^2c_2^2+...+2^nc_n^2)
(117)

對於 a_0=1 和 b_0=k^'，並且

(K^'(k)-E^'(k))/(K^'(k))=1/2(c_0^'^2+2c_1^'^2+2^2c_2^'^2+...+2^nc_n^'^2).
(118)

橢圓積分滿足大量恆等式。互補函式和模數定義為

K^'(k)=K(sqrt(1-k^2))=K(k^').
(119)

使用廣義橢圓積分的恆等式

T(a,b)=T(1/2(a+b),sqrt(ab))
(120)

寫出

1/aK(sqrt(1-(b^2)/(a^2)))=2/(a+b)K(sqrt(1-(4ab)/((a+b)^2)))
(121)
=2/(a+b)K(sqrt((a^2+b^2-2ab)/((a+b)^2)))
(122)
=2/(a+b)K((a-b)/(a+b))
(123)
K(sqrt(1-(b^2)/(a^2)))=2/(1+b/a)K((1-b/a)/(1+b/a)).
(124)

定義

k^'=b/a,
(125)

並使用

k=sqrt(1-k^('2)),
(126)

所以

K(k)=2/(1+k^')K((1-k^')/(1+k^')).
(127)

現在令 l=(1-k^')/(1+k^') 得到

l(1+k^')=1-k^'=>k^'(l+1)=1-l
(128)
k^'=(1-l)/(1+l)
(129)
k=sqrt(1-k^('2))
(130)
=sqrt(1-((1-l)/(1+l))^2)
(131)
=sqrt(((1+l)^2-(1-l)^2)/((1+l)^2))
(132)
=(2sqrt(l))/(1+l),
(133)

和

1/2(1+k^')=1/2(1+(1-l)/(1+l))
(134)
=1/2[((1+l)+(1-l))/(1+l)]
(135)
=1/(1+l).
(136)

寫 k 而不是 l,

K(k)=1/(k+1)K((2sqrt(k))/(1+k)).
(137)

類似地，根據 Borwein 和 Borwein (1987)，

E(k)=(1+k)/2E((2sqrt(k))/(1+k))+(k^('2))/2K(k)
(138)
E(k)=(1+k^')E((1-k^')/(1+k^'))-k^'K(k).
(139)

用互補函式表示的表示式可以從交換 (◇)、(◇)、(◇) 和 (◇) 中的模數及其補數匯出。

K^'(k)=K(k^')
(140)
=2/(1+k)K((1-k)/(1+k))
(141)
=2/(1+k)K^'(sqrt(1-((1-k)/(1+k))^2))
(142)
=2/(1+k)K^'((2sqrt(k))/(1+k))
(143)
=1/(1+k^')K((2sqrt(k^'))/(1+k^'))
(144)
=1/(1+k^')K^'((1-k^')/(1+k^')),
(145)

和

E^'(k)=(1+k)E^'((2sqrt(k))/(1+k))-kK^'(k)
(146)
E^'(k)=((1+k^')/2)E^'((1-k^')/(1+k^'))+(k^2)/2K^'(k).
(147)

取比率

(K^'(k))/(K(k))=2(K^'((2sqrt(k))/(1+k)))/(K((2sqrt(k))/(1+k)))=1/2(K^'((1-k^')/(1+k^')))/(K((1-k^')/(1+k^')))
(148)

得到 2 次模方程。同樣正確的是

K(x)=4/((1+sqrt(x^'))^2)K([(1-RadicalBox[{1, -, {x, ^, 4}}, 4])/(1+RadicalBox[{1, -, {x, ^, 4}}, 4])]^2).
(149)

另見

阿貝爾積分, 卡爾森橢圓積分, 第一類完全橢圓積分, 第二類完全橢圓積分, 第三類完全橢圓積分, Delta 幅度, 橢圓幅角, 橢圓特徵, 橢圓函式, 第一類橢圓積分, 第二類橢圓積分, 第三類橢圓積分, 橢圓積分奇異值, 橢圓模, 休曼Lambda函式, 雅可比幅度, 雅可比橢圓函式, 雅可比 Zeta 函式, 模角, Nome, 引數, 魏爾斯特拉斯橢圓函式

使用 探索

參考文獻

Abramowitz, M. 和 Stegun, I. A. (Eds.). "Elliptic Integrals." Ch. 17 in Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th printing. New York: Dover, pp. 587-607, 1972.Arfken, G. "Elliptic Integrals." §5.8 in Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 321-327, 1985.Borwein, J. M. 和 Borwein, P. B. Pi & the AGM: A Study in Analytic Number Theory and Computational Complexity. New York: Wiley, 1987.Hancock, H. Elliptic Integrals. New York: Wiley, 1917.Kármán, T. von 和 Biot, M. A. Mathematical Methods in Engineering: An Introduction to the Mathematical Treatment of Engineering Problems. New York: McGraw-Hill, p. 121, 1940.King, L. V. The Direct Numerical Calculation of Elliptic Functions and Integrals. London: Cambridge University Press, 1924.Prasolov, V. 和 Solovyev, Y. Elliptic Functions and Elliptic Integrals. Providence, RI: Amer. Math. Soc., 1997.Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; 和 Vetterling, W. T. "Elliptic Integrals and Jacobi Elliptic Functions." §6.11 in Numerical Recipes in FORTRAN: The Art of Scientific Computing, 2nd ed. Cambridge, England: Cambridge University Press, pp. 254-263, 1992.Prudnikov, A. P.; Brychkov, Yu. A.; 和 Marichev, O. I. Integrals and Series, Vol. 1: Elementary Functions. New York: Gordon & Breach, 1986.Timofeev, A. F. Integration of Functions. Moscow and Leningrad: GTTI, 1948.Weisstein, E. W. "Books about Elliptic Integrals." http://www.ericweisstein.com/encyclopedias/books/EllipticIntegrals.html.Whittaker, E. T. 和 Watson, G. N. A Course in Modern Analysis, 4th ed. Cambridge, England: Cambridge University Press, 1990.Woods, F. S. "Elliptic Integrals." Ch. 16 in Advanced Calculus: A Course Arranged with Special Reference to the Needs of Students of Applied Mathematics. Boston, MA: Ginn, pp. 365-386, 1926.

在 上被引用

橢圓積分

請引用本文為

Weisstein, Eric W. "橢圓積分。" 來自 -- 資源。 https://mathworld.tw/EllipticIntegral.html

學科分類