主題
Search

超立方體線段選取

下載 Mathematica 筆記本下載 Wolfram 筆記本

設兩個點 x 和 y 從一個單位 n 維 超立方體 中隨機選取。點 Delta(n) 之間的期望距離，即 平均線段長度，則為

Delta(n)=int_0^1...int_0^1_()_(2n)sqrt((x_1-y_1)^2+(x_2-y_2)^2+...+(x_n-y_n)^2)dx_1...dx_ndy_1...dy_n.
(1)

這個多重積分僅針對 n 的小值進行了分析計算。情況 Delta(1) 對應於區間 [0,1] 中兩個隨機點之間的 線段選取。

HypercubeLinePlot

以下表格給出了 Delta(n) 的前幾個值。

nOEISDelta(n)
1--0.3333333333...
2A0915050.5214054331...
3A0730120.6617071822...
4A1039830.7776656535...
5A1039840.8785309152...
6A1039850.9689420830...
7A1039861.0515838734...
8A1039871.1281653402...

函式 Delta(n) 滿足

1/3n^(1/2)<=Delta(n)<=(1/6n)^(1/2)sqrt(1/3[1+2(1-3/(5n))^(1/2)])
(2)

(Anderssen等人，1976)，與實際值一起繪製在上方。

HypercubeLinePickingIntegrands

M. Trott（私人通訊，2005年2月23日）設計了一種巧妙的演算法，用於將 2n 維積分簡化為一維被積函式 I_n(x) 上的積分，使得

Delta(n)=int_0^inftyI_n(x)dx.
(3)

前幾個值為

I_1=(2e^(-x^2))/(3sqrt(pi))
(4)
I_2=(e^(-x^2)erf(x))/(3x)+(4e^(-2x^2))/(15sqrt(pi))+(4e^(-x^2))/(15sqrt(pi))
(5)
I_3=-2/5e^(-x^2)sqrt(pi)erf^2(x)+(4e^(-2x^2)erf(x))/(5x)+(e^(-x^2)erf(x))/(5x)-(12e^(-3x^2))/(35sqrt(pi))+(68e^(-2x^2))/(105sqrt(pi))+(8e^(-x^2))/(105sqrt(pi))
(6)
I_4=-(2e^(-x^2)pierf^3(x))/(15x)-(136)/(105)e^(-2x^2)sqrt(pi)erf^2(x)-(32)/(105)e^(-x^2)sqrt(pi)erf^2(x)+(197e^(-3x^2)erf(x))/(210x)+(104e^(-2x^2)erf(x))/(105x)+(e^(-x^2)erf(x))/(14x)-(676e^(-4x^2))/(945sqrt(pi))+(16e^(-3x^2))/(35sqrt(pi))+(146e^(-2x^2))/(315sqrt(pi))+(16e^(-x^2))/(945sqrt(pi)).
(7)

在 x->0 的極限情況下，這些值對於 n=1, 2, ... 由 1/sqrt(pi) 乘以 2/3, 6/5, 50/21, 38/9, 74/11, ... 給出（OEIS A103990 和 A103991）。

這等價於計算 箱積分

Delta_n(s)=s/(Gamma(1-1/2s))int_0^infty(1-[d(u)]^n)/(u^(s+1))du
(8)

其中

d(u)=int_0^1int_0^1e^(-u^2(x-y)^2)dxdy
(9)
=int_0^1int_0^1(e^(-u^2)-1+sqrt(pi)uerf(u))/(u^2)du
(10)

(Bailey等人，2006)。

這些給出了 n=1, 2, 3 和 4 的閉合形式結果

Delta(1)=1/3
(11)
Delta(2)=1/(15)[sqrt(2)+2+5ln(1+sqrt(2))]
(12)
Delta(3)=1/(105)[4+17sqrt(2)-6sqrt(3)+21ln(1+sqrt(2))+42ln(2+sqrt(3))-7pi]
(13)
Delta(4)=(136)/(105)sqrt(2)tan^(-1)(1/2sqrt(2))-(34)/(105)pisqrt(2)+8/(105)sqrt(3)+(73)/(630)sqrt(2)+4/5Cl_2(alpha)-4/5Cl_2(alpha+1/2pi)+(197)/(420)ln3+1/(14)ln(1+sqrt(2))-4/5alphaln(1+sqrt(2))-1/5piln(1+sqrt(2))+(52)/(105)ln(2+sqrt(3))-(23)/(135)-(16)/(315)pi+(26)/(15)K
(14)
Delta(5)=(65)/(42)K-(380)/(6237)sqrt(5)+(568)/(3465)sqrt(3)-4/(189)pi-(449)/(3465)-(73)/(63)sqrt(2)tan^(-1)(1/4sqrt(2))-(184)/(189)ln2+(64)/(189)ln(sqrt(5)+1)+1/(54)ln(1+sqrt(2))+(40)/(63)ln(sqrt(2)+sqrt(6))-5/(28)piln(1+sqrt(2))+(52)/(63)piln2+(295)/(252)ln3+4/(215)pi^2+(3239)/(62370)sqrt(2)-8/(21)sqrt(3)cot^(-1)(sqrt(15))-(52)/(63)piln(sqrt(2)+sqrt(6))-5/7alpha+5/7Cl_2(alpha)-5/7Cl_2(alpha+1/2pi)+(52)/(63)K_1,
(15)

其中 Cl_2(z) 是一個 克勞森函式，K 是 卡塔蘭常數，並且

alpha=sin^(-1)((sqrt(2))/6-2/3).
(16)

上面的 n=4 情況首次發表於這項工作；上面給出的簡化形式歸功於 Bailey等人（2007）。嘗試將 Delta(5) 簡化為求積法，得到了閉合形式的部分，但以下單個部分除外

K_1=pi[1/2ln(2+sqrt(3))-int_0^infty(e^(-2x^2)erf^3(x))/xdx]
(17)
=int_3^4(sec^(-1)x)/(sqrt((x-3)(x-1)))dx
(18)
=int_0^(I[cos^(-1)2])sec^(-1)(2+coshtheta)dtheta
(19)

這部分似乎難以以閉合形式積分（Bailey等人，2007，第272頁）。

為 立方體線段選取 獲得的 Delta(3) 值有時被稱為 羅賓斯常數。

另請參閱

立方體線段選取, 超立方體點選取, 平均線段長度, 羅賓斯常數, 正方形線段選取, 正方形點選取, 正方形三角形選取, 單位正方形積分

使用 探索

參考文獻

Anderssen, R. S.; Brent, R. P.; Daley, D. J.; 和 Moran, A. P. "關於 int_0^1...int_0^1sqrt(x_1^2+...+x_k^2)dx_1...dx_k 和泰勒級數方法。" SIAM J. Appl. Math. 30, 22-30, 1976.Bailey, D. H.; Borwein, J. M.; 和 Crandall, R. E. "箱積分。" J. Comput. Appl. Math. 206, 196-208, 2007.Bailey, D. H.; Borwein, J. M.; 和 Crandall, R. E. "箱積分理論進展。" Math. Comput. 79, 1839-1866, 2010.Bailey, D. H.; Borwein, J. M.; Calkin, N. J.; Girgensohn, R.; Luke, D. R.; 和 Moll, V. H. 實驗數學行動。 Wellesley, MA: A K Peters, p. 272, 2007.Finch, S. R. "幾何機率常數。" §8.1 在 數學常數。 Cambridge, England: Cambridge University Press, pp. 479-484, 2003.Le Lionnais, F. 卓越數。 Paris: Hermann, p. 30, 1983.Robbins, D. "盒子中兩點之間的平均距離。" Amer. Math. Monthly 85, 278, 1978.Sloane, N. J. A. 序列 A073012, A091505, A103983, A103984, A103985, A103986, A103987, A103988, A103989, A103990, 和 A103991 在 "整數序列線上百科全書" 中。Trott, M. "隨機三角形的面積。" Mathematica J. 7, 189-198, 1998.

在 上被引用

超立方體線段選取

引用為

Weisstein, Eric W. "超立方體線段選取。" 來源 Web 資源。 https://mathworld.tw/HypercubeLinePicking.html

主題分類