Probabilitas bahwa titik-titik acak seragam dalam sebuah persegi panjang memiliki jarak Euclidean kurang dari ambang yang diberikan

8

Asumsikan kita memiliki titik dalam persegi panjang dengan batas , dan titik-titik ini terdistribusi secara merata di bidang ini. (Saya tidak begitu familiar dengan statistik, jadi saya tidak tahu perbedaan antara memilih secara seragam sebuah simpul di area , atau secara seragam memilih -aksi dari dan -aksi dari secara independen).n[0,a]×[0,b][0,a]×[0,b]x[0,a]y[0,b]

Diberi ambang jarak , saya mungkin ingin tahu probabilitas bahwa jarak dua titik 'Euclidean kurang dari , atau lebih tepatnya, berapa banyak pasangan node' jarak akan kurang dari ?ddd


Mungkin uraian berikut tidak akan mendua.

Biarkan saya menentukan masalah ini. Diberikan node dan ambang batas . Ini poin didistribusikan secara merata dalam persegi panjang . Nyatakan variabel acak sebagai jumlah pasangan titik dalam jarak . Temukan .ndn[0,a]×[0,b]ξdE[ξ]

zhouzhuojie
sumber
Anda harus menelusuri pertanyaan-pertanyaan di math.SE juga, karena saya ingat beberapa yang terkait di sana. Mereka kemungkinan ditandai probability.
kardinal
1
Berikut adalah beberapa pertanyaan yang saya ingat pernah saya lihat di math.SE, tetapi tidak satu pun dari mereka yang Anda tanyakan: ( 1 ) math.stackexchange.com/questions/64028 ( 2 ) math.stackexchange.com/questions/66777 ( 3 ) math.stackexchange.com/questions/101692 ( 4 ) math.stackexchange.com/questions/50775
kardinal

Jawaban:

15

Kita dapat memecahkan masalah ini secara analitis menggunakan beberapa intuisi dan argumen geometris . Sayangnya, jawabannya cukup panjang dan agak berantakan.

Pengaturan dasar

Pertama, mari kita buat beberapa notasi. Asumsikan kita menggambar titik secara seragam secara acak dari kotak . Kami berasumsi tanpa kehilangan keumuman bahwa . Biarkan menjadi koordinat dari titik pertama dan menjadi koordinat dari titik kedua. Kemudian, , , , dan saling independen dengan didistribusikan secara seragam pada dan didistribusikan secara seragam pada .[0,a]×[0,b]0<b<a(X1,Y1)(X2,Y2)X1X2Y1Y2Xi[0,a]Yi[0,b]

Pertimbangkan jarak Euclidean antara dua titik. Ini adalah manadan.

D=(X1X2)2+(Y1Y2)2=:Z12+Z22,
Z1=|X1X2|Z2=|Y1Y2|

Distribusi segitiga

Karena dan adalah seragam independen, maka memiliki distribusi segitiga, di manamemiliki distribusi dengan fungsi kerapatan Fungsi distribusi yang sesuai adalah untuk . Demikian pula,memiliki kepadatan dan fungsi distribusi .X1X2X1X2Z1=|X1X2|

fa(z1)=2a2(az1),0<z1<a.
Fa(z1)=1(1z1/a)20z1aZ2=|Y1Y2|fb(z2)Fb(z2)

Perhatikan bahwa karena adalah fungsi hanya dari dua dan adalah fungsi hanya dari , maka dan independen. Jadi jarak antara titik adalah norma euclidean dari dua variabel acak independen (dengan distribusi yang berbeda).Z1XiZ2YiZ1Z2

Panel kiri gambar menunjukkan distribusi dan panel kanan menunjukkandi mana dalam contoh ini.X1X2Z1=|X1X2|a=5

Kepadatan segitiga

Beberapa probabilitas geometris

Jadi dan independen dan didukung pada danZ1Z2[0,a][0,b]masing-masing. Untuk diperbaikid, fungsi distribusi jarak euclidean adalah

P(Dd)={z12+z22d2}fa(z1)fb(z2)dz1dz2.

Kita dapat menganggap ini secara geometris memiliki distribusi pada persegi panjang dan mempertimbangkan seperempat lingkaran jari-jari . Kami ingin mengetahui probabilitas yang ada di dalam persimpangan dua wilayah ini. Ada tiga kemungkinan berbeda untuk dipertimbangkan:[0,a]×[0,b]d

Wilayah 1 (oranye): . Di sini lingkaran seperempat terletak sepenuhnya di dalam persegi panjang.0d<b

Wilayah 2 (merah): . Di sini lingkaran seperempat memotong persegi panjang di sepanjang tepi atas dan bawah.bda

Wilayah 3 (biru): . Lingkaran seperempat memotong persegi panjang di sepanjang tepi atas dan kanan.a<da2+b2

Berikut adalah gambar, di mana kita menggambar contoh radius dari masing-masing dari ketiga jenis. Persegi panjang didefinisikan oleh , . Peta panas skala abu-abu dalam persegi panjang menunjukkan kepadatan mana area gelap memiliki kepadatan lebih tinggi dan area yang lebih terang memiliki kepadatan lebih kecil. Mengklik pada gambar akan membuka versi yang lebih besar.a=5b=4fa(z1)fb(z2)dz1dz2

Distribusi terinduksi: titik-temu

Beberapa kalkulus jelek

Untuk menghitung probabilitas, kita perlu melakukan beberapa kalkulus. Mari kita pertimbangkan masing-masing daerah secara bergantian dan kita akan melihat bahwa integral yang sama akan muncul. Integral ini memiliki bentuk tertutup, meskipun tidak terlalu cantik.

Wilayah 1 : .0d<b

P(Dd)=0d0d2y2fb(y)fa(x)dxdy=0dfb(y)0d2y2fa(x)dxdy.

Sekarang, integral integral menghasilkan . Jadi, kita dibiarkan menghitung integral dari bentuk di mana dalam kasus ini bunga . Penangkal integrand adalah 1a2d2y2(2ad2y2)

G(c)G(0)=0c(by)d2y2(2ad2y2)dy,
c=d
G(y)=(by)d2y2(2ad2y2)dy=a3d2y2(y(3b2y)+2d2)+abd2tan1(yd2y2)bd2y+by33+(dy)22y44.

Dari sini kita dapatkan .P(Dd)=2a2b2(G(d)G(0))

Wilayah 2 : .bda

P(Dd)=2a2b2(G(b)G(0)),
dengan alasan yang sama seperti untuk Wilayah 1, kecuali sekarang kita harus mengintegrasikan sepanjang sumbu sampai ke bukan hanya .ybd

Wilayah 3 : . a<da2+b2

P(Dd)=0d2a2fb(y)dy+d2a2bfb(y)0d2y2fa(x)dxdy=Fb(d2a2)+2a2b2(G(b)G(d2a2))

Di bawah ini adalah simulasi 20.000 titik di mana kami memplot distribusi empiris sebagai titik abu-abu dan distribusi teoretis sebagai garis, diwarnai sesuai dengan wilayah tertentu yang berlaku.

Cdf empiris dan teoretis

Dari simulasi yang sama, di bawah ini kami plot 100 pasang poin pertama dan menarik garis di antara mereka. Masing-masing diwarnai sesuai dengan jarak antara sepasang titik dan wilayah mana dari jarak ini.

Sampel acak poin

Jumlah pasangan poin yang diharapkan dalam jarak hanyalah dengan linearitas harapan.d

E[ξ]=(n2)P(Dd),
kardinal
sumber
3
+1. Kerja bagus! Akan sangat bagus untuk melihat jawaban yang diungkapkan dalam hal sifat geometris intrinsik dari persegi panjang: ia harus bergantung pada hal-hal seperti area, perimeter, dan konfigurasi keempat sudutnya. (Literatur - yang saya lihat direferensikan tetapi belum memiliki akses - tampaknya berfokus pada domain dengan batas-batas yang halus.)
whuber
Terima kasih. Itu saran yang bagus. Saya akan mencoba membuat penyederhanaan dan reformulasi tersebut.
kardinal
@ kardinal Kerja yang sangat bagus! Saya terkejut bahwa Anda benar-benar menjawab masalah bahkan dengan cdf terperinci. Terima kasih.
zhouzhuojie
0

Jika titik benar-benar terdistribusi secara seragam, yaitu dalam pola yang diketahui tetap, maka untuk jarak apa pun d, Anda dapat dengan mudah mengulang semua pasangan dan menghitung yang berada dalam jarak tersebut. Probabilitas Anda adalah (angka itu / n).

Jika Anda memiliki kebebasan tambahan untuk memilih bagaimana n poin didistribusikan / dipetik, maka ini adalah versi persegi panjang dari paradoks Bertrand . Halaman itu menunjukkan sejumlah cara untuk menjawab pertanyaan ini berdasarkan pada bagaimana Anda mendistribusikan poin Anda.

cape1232
sumber
Pertanyaannya adalah tentang distribusi untuk titik-titik yang didistribusikan secara seragam: ini adalah variabel acak, bukan "pola yang diketahui pasti", dan seseorang tidak bisa begitu saja mengulangi pasangannya!
whuber
Saya pikir Anda mungkin salah memahami pertanyaan OP. Juga, distribusi yang diinginkan didefinisikan secara jelas dalam pertanyaan. Komentar saya kepada OP mengisyaratkan bahwa sudah ada solusi pada jaringan SE untuk pertanyaan ini, maka yang ini kemungkinan besar bisa ditutup. :)
kardinal
Apakah Anda yakin ada solusi pada matematika. SE, kardinal? Ini adalah masalah yang sulit karena efek tepi. Mungkin ada solusi pada torus yang rata.
whuber
@whuber: Solusi? Tidak. Tapi, saya hampir yakin pertanyaan ini muncul. :) Saya akan melihat apakah saya dapat menemukannya. Bagaimanapun, saya tidak yakin masalah ini sangat sulit, bahkan dalam kasus ini. Saya yakin Anda dapat menggunakan terjemahan invarian untuk menyederhanakannya. Tapi, saya belum mengerjakan detailnya.
kardinal
1
@ kardinal Terima kasih. Sebenarnya saya telah memeriksa semua pertanyaan di Math.SE, tetapi saya masih belum dapat menemukan beberapa masalah yang dekat dengannya.
zhouzhuojie