主題
Search

圓盤線段拾取

DOWNLOAD Mathematica Notebook下載 Wolfram 筆記本
Disk line picking

使用圓盤點拾取，

x=sqrt(r)costheta
(1)
y=sqrt(r)sintheta
(2)

對於 r in [0,1], theta in [0,2pi), 在單位圓盤中隨機選擇兩個點，並找到兩點之間距離 s 的分佈。不失一般性，取第一個點為 (r,theta)=(r_1,0)，第二個點為 (r_2,theta)。然後

s^_=(int_0^1int_0^1int_0^pisqrt(r_1+r_2-2sqrt(r_1r_2)costheta)dr_1dr_2dtheta)/(int_0^1int_0^1int_0^pidr_1dr_2dtheta)
(3)
=1/piint_0^1int_0^1int_0^pisqrt(r_1+r_2-2sqrt(r_1r_2)costheta)dr_1dr_2dtheta
(4)
=(128)/(45pi)
(5)
=0.905414787...
(6)

(OEIS A093070; Uspensky 1937, p. 258; Solomon 1978, p. 36)。

DiskLinePickingP

這是球體線段拾取在 n=2 時的特殊情況，因此半徑為 R 的圓盤的完整機率函式為

P_2(s,R)=(4s)/(piR^2)cos^(-1)(s/(2R))-(2s^2)/(piR^3)sqrt(1-(s^2)/(4R^2))
(7)

(Solomon 1978, p. 129; Mathai 1999, p. 204)。

線段長度分佈的原點矩由下式給出

mu_n^'=int_0^2s^nP_2(s,1)ds
(8)
=(2Gamma(n+3))/((n+2)Gamma(2+1/2n)Gamma(3+1/2n)),
(9)

其中 Gamma(x) 是伽瑪函式，n>-2。 1/r 的期望值由 n=-1 給出，得到

mu_(-1)=(16)/(3pi)
(10)

(Solomon 1978, p. 36; Pure et al. )。前幾個矩為

mu_0^'=1
(11)
mu_1^'=(128)/(45pi)
(12)
mu_2^'=1
(13)
mu_3^'=(2048)/(525pi)
(14)
mu_4^'=5/3
(15)
mu_5^'=(16384)/(2205pi)
(16)

(OEIS A093526 和 A093527 和 OEIS A093528 和 A093529)。矩 mu_(2n-2)^' 為整數的情況發生在 n=1, 2, 6, 15, 20, 28, 42, 45, 66, ... (OEIS A014847)，這非常令人驚訝地正是 n 的值，使得 n|C_n，其中 C_n 是卡塔蘭數 (E. Weisstein, 3月 30, 2004)。

另請參閱

球體線段拾取, 圓線段拾取, 扇形線段拾取, 圓盤三角形拾取, 直線線段拾取

使用 探索

參考文獻

Sloane, N. J. A. Sequences A014847, A093070, A093526, A093527, A093528, 和 A093529 in "整數數列線上大全。"Mathai, A. M. 幾何機率導論：具有應用的分佈方面。 Amsterdam, Netherlands: Gordon and Breach, 1999。Pure, R.; Durran, S.; Tong, F.; Pan, J. "任意多邊形中兩個隨機點之間的距離分佈。" 即將發表於 Math. Meth. Appl. Sci.Solomon, H. 幾何機率。 Philadelphia, PA: SIAM, 1978。Uspensky, J. V. Ch. 12, Problem 5 in 數學機率導論。 New York: McGraw-Hill, pp. 257-258, 1937。

在 上引用

圓盤線段拾取

請引用本文為

Weisstein, Eric W. "圓盤線段拾取。" 來自 Web 資源。 https://mathworld.tw/DiskLinePicking.html

主題分類