主題
Search

單位正方形積分

DOWNLOAD Mathematica Notebook下載 Wolfram Notebook

幾何機率中出現的單位正方形上的積分是

int_0^1int_0^1sqrt(x^2+y^2)dxdy=1/3[sqrt(2)+sinh^(-1)(1)] int_0^1int_0^1sqrt((x-1/2)^2+(y-1/2)^2)dxdy =1/6[sqrt(2)+sinh^(-1)(1)],
(1)

這些積分分別給出了在正方形點選取中，從單位正方形內隨機選取一點到角和中心的平均距離。

涉及絕對值的單位正方形積分由下式給出

int_0^1int_0^1|x-y|^ndxdy=2/((n+1)(n+2))
(2)
int_0^1int_0^1|x+y|^ndxdy=(2(2^(n+1)-1))/((n+1)(n+2)),
(3)

對於 R[n]>-1 和 R[n]>-2，分別。

另一個簡單的積分由下式給出

int_0^1int_0^1(dxdy)/(sqrt(1+x^2+y^2))=4ln(2+sqrt(3))-2/3pi
(4)

(Bailey et al. 2007, p. 67)。 分母平方得到

int_0^1int_0^1(dxdy)/(x^2+y^2+1)=int_0^1(tan^(-1)(1/(sqrt(1+y^2))))/(sqrt(1+y^2))dy
(5)
=int_0^11/(sqrt(4x^2-1))[tan^(-1)(3/(sqrt(4x^2-1)))-tan^(-1)((1+2x^2)/(sqrt(4x^2-1)))]dx
(6)
=int_0^(pi/4)(tan^(-1)(costheta))/(costheta)dtheta
(7)
=1/2pisinh^(-1)1-K+1/6_3F_2(1/2,1,1;3/2,3/2;1/9)
(8)
=0.63951...
(9)

(OEIS A093754; M. Trott, pers. comm.)，其中 K 是 卡塔蘭常數，_3F_2(a,b,c;d,e;z) 是廣義超幾何函式。 一個相關的積分由下式給出

int_0^1int_0^1(dxdy)/(x^2+y^2),
(10)

這在黎曼意義上是發散的，可以透過轉換為極座標快速看出。 然而，取Hadamard 積分來丟棄單位圓內的發散部分得到

Hint_0^1int_0^1(dxdy)/(x^2+y^2)=intint_(x^2+y^2>1; 0<x<1; 0<y<1)(dxdy)/(x^2+y^2)
(11)
=int_0^11/x[tan^(-1)(1/x)-tan^(-1)((sqrt(1-x^2))/x)]dx
(12)
=1/2piln2-K
(13)
=0.172827...
(14)

(OEIS A093753)，其中 K 是 卡塔蘭常數。

Guillera 和 Sondow (2005) 給出了單位正方形上的一系列優美的積分，它們來自一般積分

int_0^1int_0^1(x^(u-1)y^(v-1))/(1-xyz)[-ln(xy)]^sdxdy=Gamma(s+1)(Phi(z,s+1,v)-Phi(z,s+1,u))/(u-v)
(15)
int_0^1int_0^1((xy)^(u-1))/(1-xyz)[-ln(xy)]^sdxdy=Gamma(s+2)Phi(z,s+2,u),
(16)

對於 u,v>0, R[s]>-2 如果 z in C-[1,infty), 以及 R[s]>-1 如果 z=1，其中 Gamma(s) 是 伽馬函式，Phi(z,s,a) 是 萊奇超越函式。 在 (15) 中，為了處理 u=v 的情況，取極限 v->u，這得到 (16)。

另一個結果是

int_0^1int_0^1(1-x)/((1-xy)(-lnxy))(xy)^(u-1)dxdy=lnu-psi_0(u)
(17)

(Guillera 和 Sondow 2005)，對於 u>0 且其中 psi_0(z) 是 雙伽馬函式。

Guillera 和 Sondow (2005) 也給出了

int_0^1int_0^1([-ln(xy)]^s)/(1-xy)dxdy=Gamma(s+2)zeta(s+2)
(18)
int_0^1int_0^1([-ln(xy)]^s)/(1+xy)dxdy=Gamma(s+2)eta(s+2)
(19)
int_0^1int_0^1([-ln(xy)]^s)/(1+x^2y^2)dxdy=Gamma(s+2)beta(s+2),
(20)

其中第一個式子適用於 R[s]>-1，第二個和第三個式子適用於 R[s]>-2，zeta(s) 是 黎曼 zeta 函式，eta(s) 是 狄利克雷 eta 函式，beta(s) 是 狄利克雷 beta 函式。 (19) 是 Hadjicostas (2002) 為整數 s>=0 發現的。 公式 (18) 和 (19) 是透過設定 u=1 ，然後分別取 z=1 和 z=-1 得到的 (16) 的特殊情況。

優美的公式

int_0^1int_0^1(dxdy)/(1-xy)=zeta(2)
(21)
int_0^1int_0^1(-ln(xy))/(1-xy)dxdy=2zeta(3)
(22)

是 Beukers (1979) 給出的。 這些積分是透過取 s=0 和 1 分別得到的 (19) 的特殊情況。 涉及卡塔蘭常數 K 的類似公式由下式給出

int_0^1int_0^1(dxdy)/((1-xy)sqrt(x(1-y)))=8K
(23)

(Zudilin 2003)。

其他與 Hadjicostas 公式相關的優美積分由下式給出

int_0^1int_0^1(x-1)/((1-xy)ln(xy))dxdy=gamma
(24)
int_0^1int_0^1(x-1)/((1+xy)ln(xy))dxdy=ln(4/pi)
(25)

(Sondow 2003, 2005; Borwein et al. 2004, p. 49)，其中 gamma 是 尤拉-馬歇羅尼常數。

其他特殊情況的集合 (Guillera 和 Sondow 2005) 包括

int_0^1int_0^1(-xln(xy))/(1+x^2y^2)dxdy=K-1/(48)pi^2
(26)
int_0^1int_0^1(-xln(xy))/(1-x^2y^2)dxdy=1/(12)pi^2
(27)
int_0^1int_0^1(-ln(xy))/(1+x^2y^2)dxdy=1/(16)pi^3
(28)
int_0^1int_0^1(dxdy)/(1+x^2y^2)=K
(29)
int_0^1int_0^1(-dxdy)/((2-xy)ln(xy))=ln2
(30)
int_0^1int_0^1(dxdy)/(2-xy)=1/(12)pi^2-1/2ln^22
(31)
int_0^1int_0^1(-ln(xy))/(2-xy)dxdy=7/4zeta(3)-1/6pi^2ln2+1/3ln^32
(32)
int_0^1int_0^1(-x)/((2-xy)ln(xy))dxdy=lnsigma
(33)
int_0^1int_0^1(-dxdy)/((phi-xy)ln(xy))=2lnphi
(34)
int_0^1int_0^1(-dxdy)/((phi^2-xy)ln(xy))=lnphi
(35)
int_0^1int_0^1(dxdy)/(phi-xy)=1/(10)pi^2-ln^2phi
(36)
int_0^1int_0^1(dxdy)/(phi^2-xy)=1/(15)pi^2-ln^2phi
(37)
int_0^1int_0^1(-ln(xy))/(phi^2-xy)dxdy=8/5zeta(3)-4/(15)pi^2lnphi+4/3ln^3phi
(38)
int_0^1int_0^1(-x)/((1+xy)ln(xy))dxdy=ln(1/2pi)
(39)
int_0^1int_0^1(-x)/((1+x^2y^2)ln(xy))dxdy=ln[(sqrt(2pi))/(Gamma^2(3/4))]
(40)
int_0^1int_0^1(-1)/((1+x^2y^2)ln(xy))dxdy=1/4pi
(41)
int_0^1int_0^1x/(1-x^3y^3)dxdy=pi/(3sqrt(3))
(42)
int_0^1int_0^1(dxdy)/(1-x^2y^2)=1/8pi^2
(43)
int_0^1int_0^1(ln(2-x))/(1-xy)dxdy=1/4pi^2ln2-zeta(3)
(44)
int_0^1int_0^1(ln(2-xy))/(1-xy)dxdy=5/8zeta(3)
(45)
int_0^1int_0^1(ln(1+x))/(1-xy)dxdy=5/8zeta(3)
(46)
int_0^1int_0^1(ln(1+xy))/(1-xy)dxdy=1/4pi^2ln2-zeta(3)
(47)
int_0^1int_0^1(-ln(1-x))/(1-xy)dxdy=2zeta(3)
(48)
int_0^1int_0^1(-ln(1-xy))/(1-xy)dxdy=zeta(3)
(49)
int_0^1int_0^1(xdxdy)/([-ln(xy)]^(3/2))=2(sqrt(2)-1)sqrt(pi)
(50)
int_0^1int_0^1(dxdy)/([-ln(xy)]^(3/2))=sqrt(pi)
(51)
int_0^1int_0^1(dxdy)/([-ln(xy)]^(5/4))=Gamma(3/4)
(52)
int_0^1int_0^1(xln^2(xy))/((1+x^2y^2)^2)dxdy=K-1/(48)pi^2+1/(32)pi^3
(53)
int_0^1int_0^1(ln^4(xy))/((1+xy)^2)dxdy=(225)/2zeta(5)
(54)
int_0^1int_0^1(-dxdy)/((1+x^2y^2)^2ln(xy))=1/8(pi+2)
(55)
int_0^1int_0^1(-xdxdy)/((1+xy)^2ln(xy))=ln((A^6)/(2^(1/6)sqrt(pie)))
(56)
int_0^1int_0^1(1-x)/((1+xy)ln^2(xy))dxdy=ln((pi^(1/2)A^6)/(2^(7/6)e))
(57)
int_0^1int_0^1(1-x^2)/((1+x^2y^2)ln^2(xy))dxdy=ln[(Gamma(1/4))/(2Gamma(3/4))]+(2K)/pi-1/2
(58)
int_0^1int_0^1(1-x)/([-ln(xy)]^(5/2))dxdy=1/3sqrt(pi)(8sqrt(2)-10),
(59)

其中 zeta(n) 是 黎曼 zeta 函式，zeta(3) 是 Apéry 常數，phi 是 黃金比例，sigma 是 Somos 二次遞推常數，A 是 Glaisher-Kinkelin 常數。 方程 (57) 出現在 Sondow (2005) 中，但它是 Guillera 和 Sondow (2005) 考慮的型別的一個特例。

透過變數替換 x=X/Y, y=Y 可以找到這些積分中大多數積分在 [0,1] 上的對應單積分。 雅可比行列式給出 dxdy=Y^(-1)dXdY，新的積分限是 {X,0,1}, {Y,X,1}。 然後對 Y 積分得到 [0,1] 上的一個一維積分。 詳情請參見 Guillera-Sondow 定理 3.1 證明的第一部分。

另請參閱

雙重積分, Hadjicostas 公式, 超立方體線段選取, 正方形點選取, 單位圓盤積分, 單位正方形

使用 探索

參考文獻

Bailey, D. H.; Borwein, J. M.; Calkin, N. J.; Girgensohn, R.; Luke, D. R.; and Moll, V. H. 實驗數學在行動。 Wellesley, MA: A K Peters, 2007.Beukers, F. "關於 zeta(2) 和 zeta(3) 的無理性證明的註記。" Bull. London Math. Soc. 11, 268-272, 1979.Borwein, J.; Bailey, D.; and Girgensohn, R. 數學實驗：通往發現的計算路徑。 Wellesley, MA: A K Peters, 2004.Guillera, J. and Sondow, J. "透過 Lerch 超越函式的解析延拓得到的若干經典常數的雙重積分和無窮乘積。" 16 June 2005. http://arxiv.org/abs/math.NT/0506319.Hadjicostas, P. "Beukers 積分的一些推廣。" Kyungpook Math. J. 42, 399-416, 2002.Sloane, N. J. A. 整數序列線上百科全書中的序列 A093753 和 A093754。"Sondow, J. "尤拉常數的無理性判據。" Proc. Amer. Math. Soc. 131, 3335-3344, 2003. http://arxiv.org/abs/math.NT/0209070.Sondow, J. "尤拉常數和 ln(4/pi) 的雙重積分以及 Hadjicostas 公式的一個類似物。" Amer. Math. Monthly 112, 61-65, 2005.Zudilin, W. "卡塔蘭常數的類 Apéry 差分方程。" Electronic J. Combinatorics 10, No. 1, R14, 1-10, 2003. http://www.combinatorics.org/Volume_10/Abstracts/v10i1r14.html.

在 中被引用

單位正方形積分

請引用為

Weisstein, Eric W. "單位正方形積分。" 來自 -- 資源. https://mathworld.tw/UnitSquareIntegral.html

主題分類