主題
Search

Watson 三重積分

DOWNLOAD Mathematica Notebook下載 Wolfram Notebook

Watson (1939) 考慮了以下三個三重積分，

I_1=1/(pi^3)int_0^piint_0^piint_0^pi(dudvdw)/(1-cosucosvcosw)
(1)
=(4[K(1/2sqrt(2))]^2)/(pi^2)
(2)
=(Gamma^4(1/4))/(4pi^3)
(3)
=1.39320393...
(4)
I_2=1/(pi^3)int_0^piint_0^piint_0^pi(dudvdw)/(3-cosvcosw-coswcosu-cosucosv)
(5)
=(sqrt(3)[K(1/4(sqrt(6)-sqrt(2)))]^2)/(pi^2)
(6)
=(3Gamma^6(1/3))/(2^(14/3)pi^4)
(7)
=0.448220394...
(8)
I_3=1/(pi^3)int_0^piint_0^piint_0^pi(dudvdw)/(3-cosu-cosv-cosw)
(9)
=(4(18+12sqrt(2)-10sqrt(3)-7sqrt(6))[K((2-sqrt(3))(sqrt(3)-sqrt(2)))]^2)/(pi^2)
(10)
=(18+12sqrt(2)-10sqrt(3)-7sqrt(6))[1+2sum_(k=1)^(infty)exp(-k^2pisqrt(6))]^4
(11)
=(18+12sqrt(2)-10sqrt(3)-7sqrt(6))theta_3^4(0,e^(-pisqrt(6)))
(12)
=(sqrt(6))/(96pi^3)Gamma(1/(24))Gamma(5/(24))Gamma(7/(24))Gamma((11)/(24))
(13)
=0.505462019...
(14)

(OEIS A091670, A091671, 和 A091672), 其中 K(k) 是第一類完全橢圓積分, theta_3(0,q) 是 雅可比 theta 函式, 並且 Gamma(z) 是 伽瑪函式。這些積分的解析計算相當具有挑戰性，特別是 I_2 和 I_3。

Watson (1939) 透過進行變換處理了所有三個積分

x=tan(1/2u)
(15)
y=tan(1/2v)
(16)
z=tan(1/2w),
(17)

關於 x, y, 和 z 作為 笛卡爾座標，並更改為 極座標，

x=rsinthetacosphi
(18)
y=rsinthetasinphi
(19)
z=rcostheta
(20)

在寫下 2phi=psi 之後。

對 I_1 執行此變換得到

I_1=8/(pi^3)int_0^inftyint_0^inftyint_0^infty(dxdydz)/((1+x^2)(1+y^2)(1+z^2)-(1-x^2)(1-y^2)(1-z^2))
(21)
=4/(pi^3)int_0^inftyint_0^inftyint_0^infty(dxdydz)/(x^2+y^2+z^2+x^2y^2z^2)
(22)
=4/(pi^3)int_0^(pi/2)int_0^(pi/2)int_0^infty(sinthetadrdthetadphi)/(1+r^4sin^4thetacos^2thetasin^2phicos^2phi)
(23)
=4/(pi^3)int_0^(pi/2)int_0^(pi/2)int_0^infty(sinthetadrdthetadpsi)/(1+1/4r^4sin^4thetacos^2thetasin^2psi).
(24)

I_1 然後可以直接使用計算機代數進行積分，儘管 Watson (1939) 使用了額外的變換

t=rsinthetasqrt(1/2costhetasinpsi)
(25)

將積分分離為

I_1=(4sqrt(2))/(pi^3)int_0^infty(dt)/(1+t^4)int_0^(pi/2)(dtheta)/(sqrt(costheta))int_0^(pi/2)(dpsi)/(sqrt(sinpsi))
(26)
=(4sqrt(2))/(pi^3)·pi/(2sqrt(2))·(Gamma^2(1/4))/(2sqrt(2pi))·(Gamma^2(1/4))/(2sqrt(2pi))
(27)
=(Gamma^4(1/4))/(4pi^3).
(28)

積分 I_1 也可以透過執行其中一個積分來完成

int_0^pi(du)/(1-ccosu)=pi/(sqrt(1-c^2))
(29)

其中 c=cosvcosw 得到

I_1=1/(pi^2)int_0^piint_0^pi(dvdw)/(sqrt(1-cos^2vcos^2w)).
(30)

使用二項式級數展開

(1-c)^(-1/2)=sum_(n=0)^(infty)((1/2)_n)/(n!)c^n
(31)
=sum_(n=0)^(infty)a_nc^n
(32)

其中 (z)_n 是 Pochhammer 符號，並且

a_n=((2n-1)!!)/((2n)!!).
(33)

積分得到

I_1=1/(pi^2)int_0^piint_0^pisum_(n=0)^(infty)a_ncos^(2n)vcos^(2n)wdvdw
(34)
=1/(pi^2)sum_(n=0)^(infty)a_n(int_0^picos^(2n)vdv)^2
(35)
=1/(pi^2)sum_(n=0)^(infty)a_n(pia_n)^2
(36)
=sum_(n=0)^(infty)a_n^3.
(37)

現在，由於第一類完全橢圓積分 K(k) 的驚人恆等式

[K(k)]^2=1/4pi^2sum_(n=0)^inftya_n^3(2kk^')^(2n),
(38)

其中 k^' 是互補模量，並且 0<k<=1/sqrt(2) (Watson 1908, Watson 1939), 立即得出結論，當 k=k^'=1/sqrt(2) (即，k=k_1，第一個奇異值) 時，

K(1/2sqrt(2))=1/4pi^2sum_(n=0)^inftya_n^3=1/4pi^2I_1,
(39)

因此

I_1=(4[K(1/2sqrt(2))]^2)/(pi^2)
(40)
=(Gamma^4(1/4))/(4pi^3).
(41)

I_2 可以使用相同的規則進行變換，得到

I_2=2/(pi^3)int_0^inftyint_0^inftyint_0^infty(dxdydz)/(3product(1+x^2)-sum(1-y^2)(1-z^2)(1+x^2))
(42)
=2/(pi^3)int_0^inftyint_0^inftyint_0^infty(dxdydz)/(x^2+y^2+z^2+y^2z^2+z^2x^2+x^2y^2)
(43)
=2/(pi^3)int_0^(pi/2)int_0^(pi/2)int_0^infty(sinthetadrdthetadphi)/(1+r^2sin^2theta(cos^2theta+sin^2psisin^2phicos^2phi))
(44)
=2/(pi^3)int_0^(pi/2)int_0^(pi/2)int_0^infty(sinthetadrdthetadpsi)/(1+r^2sin^2theta(cos^2theta+1/4sin^2thetasin^2psi))
(45)
=1/(pi^2)int_0^(pi/2)int_0^(pi/2)(dthetadpsi)/(sqrt(cos^2theta+1/4sin^2thetasin^2psi))
(46)
=1/(pi^2)int_0^(pi/2)int_0^infty(dtdpsi)/(sqrt((1+t^2)(1+1/4t^2sin^2psi))),
(47)

其中在最後一步中進行了替換 t=tantheta。計算機代數可以以 Meijer G-函式 的形式返回此積分

I_2=1/(2pi^(5/2))G_(3,3)^(3,2)(4|1/2,1/2,1/2; 0,0,0),
(48)

但是，需要更巧妙的處理才能以更漂亮的形式獲得它。例如，Watson (1939) 指出

K^'(k^')=int_0^infty(dt)/(sqrt((1+t^2)(1+k^('2)t^2)))
(49)

立即給出

I_2=1/(pi^2)int_0^piK^'(1/2sinpsi)dpsi.
(50)

然而，這個積分的求積需要非常巧妙地使用 K(k) 的複雜級數恆等式，以獲得可以重新組合為被識別為的逐項積分

I_2=(K(k_1)K(k_1^'))/(pi^2)
(51)

(Watson 1939)。

對於 I_3，只能解析地完成單個積分，即

int_(-pi)^pi1/(3-cosx-cosy-cosz)dz =-(2pi)/(sqrt((cosx+cosy-2)(cosx+cosy-4))).
(52)

透過定義 w=(cosx+cosy+cosz)/3 並使用 二項式級數 展開來寫，它可以簡化為單個無限和

1/w=1/3sum_(k=0)^inftyw^k=sum_(k=0)^infty1/(3^(k+1))(cosx+cosy+cosz)^k.
(53)

但是，這可以寫成 多項式級數 並插回以獲得

I_3=1/(pi^3)int_(-pi)^piint_(-pi)^piint_(-pi)^pisum_(k=0)^infty1/(3^(k+1))×sum_(n_1,n_2,n_3>=0; n_1+n_2+n_3=k)(k!)/(n_1!n_2!n_3!)cos^(n_1)xcos^(n_2)ycos^(n_3)zdxdydz.
(54)

交換積分和求和的順序允許完成積分，從而得到

I_3=(pi^3)/3sum_(k=0)^infty1/(3^k)sum_(n_1,n_2,n_3>=0; n_1+n_2+n_3=k)(k!)/(n_1!n_2!n_3!) ×((-1)^(n_1+n_2+n_3)2^(n_1+n_2+n_3))/(Gamma^2(1/2(1-n_1))Gamma^2(1/2(1-n_2))Gamma^2(1/2(1-n_3))).
(55)

令人驚訝的是，n_i 上的和可以用閉合形式完成，產生

I_3=1/3sum_(n=0)^infty([(2n)!]^2)/(36^n(n!)^3)_3F_2(-n,-n,-n;1,1/2-n;1/4),
(56)

其中 _3F_2(a,b,c;d,e;z) 是廣義超幾何函式。然而，這個和不能以閉合形式完成。

Watson (1939) 將積分轉換為

I_3=(2sqrt(2))/piint_0^(pi/2)int_0^infty(dtdpsi)/(sqrt((1+4t^2+3t^4sin^2psi)(1+t^2sin^2psi))).
(57)

然而，為了獲得完全閉合的形式，有必要進行一些分析魔法（詳見 Watson 1939）。因此，對於這個積分來說，閉合形式的存在本身就非常令人驚訝。

另請參閱

Pólya 隨機遊走常數, Watson 公式, Watson 恆等式

使用 探索

參考文獻

Bailey, D. H.; Borwein, J. M.; Kapoor, V.; 和 Weisstein, E. W. "Ten Problems in Experimental Mathematics." Amer. Math. Monthly 113, 481-509, 2006b.Borwein, J. 和 Bailey, D. Mathematics by Experiment: Plausible Reasoning in the 21st Century. Wellesley, MA: A K Peters, 2003.Domb, C. "On Multiple Returns in the Random-Walk Problem." Proc. Cambridge Philos. Soc. 50, 586-591, 1954.Glasser, M. L. 和 Zucker, I. J. "Extended Watson Integrals for the Cubic Lattices." Proc. Nat. Acad. Sci. U.S.A. 74, 1800-1801, 1977.Joyce, G. 和 Zucker, I. J. "On the Evaluation of Generalized Watson Integrals." Proc. Amer. Math. Soc. 133, 71-81, 2005.McCrea, W. H. 和 Whipple, F. J. W. "Random Paths in Two and Three Dimensions." Proc. Roy. Soc. Edinburgh 60, 281-298, 1940.Sloane, N. J. A. Sequences A091670, A091671, 和 A091672 in "The On-Line Encyclopedia of Integer Sequences."Watson G. N. "The Expansion of Products of Hypergeometric Functions." Quart. J. Pure Appl. Math. 39, 27-51, 1907.Watson G. N. "A Series for the Square of the Hypergeometric Function." Quart. J. Pure Appl. Math. 40, 46-57, 1908.Watson, G. N. "Three Triple Integrals." Quart. J. Math., Oxford Ser. 2 10, 266-276, 1939.

在 上被引用

Watson 三重積分

請引用為

Weisstein, Eric W. "Watson 三重積分。" 來自 Web 資源。 https://mathworld.tw/WatsonsTripleIntegrals.html

主題分類