Biarkan koordinat Cartesian dari titik acak dipilih st .
Dengan demikian, jari-jari, , tidak merata seperti yang tersirat oleh 's pdf .
Meskipun demikian saya berharap hampir seragam, tidak termasuk artefak karena sisa 4 di tepi:
Berikut ini adalah fungsi kepadatan probabilitas yang dihitung secara grafis dari dan :
Sekarang jika saya membiarkan didistribusikan st maka tampaknya terdistribusi secara merata:x , y ∼ N ( 0 , 20 2 ) × N ( 0 , 20 2 ) θ
Mengapa tidak seragam ketika dan seragam ketika ?( x , y ) ∼ U ( - 10 , 10 ) × U ( - 10 , 10 ) x , y ∼ N ( 0 , 20 2 ) × N ( 0 , 20 2 )
Kode Matlab yang saya gunakan:
number_of_points = 100000;
rng('shuffle')
a = -10;
b = 10;
r = (b-a).*randn(2,number_of_points);
r = reshape(r, [2,number_of_points]);
I = eye(2);
e1 = I(:,1); e2 = I(:,2);
theta = inf*ones(1,number_of_points);
rho = inf*ones(1,number_of_points);
for i=1:length(r(1,:))
x = r(:,i);
[theta(i),rho(i)] = cart2pol(x(1),x(2));
end
figure
M=3;N=1; bins = 360;
subplot(M,N,1);
histogram(rad2deg(theta), bins)
title('Polar angle coordinate p.d.f');
subplot(M,N,2);
histogram(rho, bins);
title('Polar radius coordinate p.d.f');
subplot(M,N,3);
histogram(r(:));
title('The x-y cooridnates distrbution (p.d.f)');
Mengganti baris ke-3: r = (b-a).*randn(2,number_of_points);
dengan r = (b-a).*randn(2,number_of_points) +a ;
akan mengubah distribusi dari normal ke seragam.
Jawaban:
Anda mengacu pada transformasi dari sepasang varian independen ke representasi polar ( R , θ ) (jari-jari dan sudut), dan kemudian melihat distribusi marginal θ .(X,Y) (R,θ) θ
Saya akan menawarkan penjelasan yang agak intuitif (meskipun derivasi matematis dari kepadatan pada dasarnya melakukan apa yang saya uraikan secara informal).
Perhatikan bahwa jika Anda skala dua variabel, X dan Y oleh beberapa skala umum (misalnya pergi dari U (-1,1) ke U (-10,10) atau dari N (0,1) ke N (0,20) pada kedua variabel secara bersamaan) yang tidak membuat perbedaan pada distribusi sudut (hanya mempengaruhi skala distribusi jari-jari). Jadi mari kita pertimbangkan kasus unit saja.
Pertama-tama pertimbangkan apa yang terjadi dengan kasing seragam. Perhatikan bahwa distribusinya seragam di atas satuan kuadrat, sehingga kepadatan probabilitas di wilayah yang terkandung dalam sebanding dengan luas wilayah. Secara khusus, lihat kepadatan yang terkait dengan elemen sudut, d θ dekat horizontal (dekat sudut ) dan pada diagonal (dekat sudut ):[−1,1]2 dθ θ = π / 4θ=0 θ=π/4
Jelas elemen probabilitas (yaitu area) yang sesuai dengan elemen sudut ( ) lebih besar ketika sudutnya dekat dengan salah satu diagonal. Memang pertimbangkan untuk menuliskan lingkaran di dalam kotak; area yang direntang oleh sudut kecil yang diberikan di dalam lingkaran adalah konstan, dan kemudian bagian di luar lingkaran tumbuh ketika kita mendekati diagonal, di mana ia pada titik maksimum. d θdfθ dθ
Ini sepenuhnya menjelaskan pola yang Anda lihat dalam simulasi.
Memang, kita dapat melihat bahwa kerapatan harus proporsional dengan panjang segmen dari pusat persegi ke ujungnya; trigonometri sederhana sudah cukup untuk memperoleh kerapatan dari sana dan kemudian mudah untuk menemukan konstanta yang diperlukan untuk membuat kerapatan berintegrasi ke 1.
[Sunting: menambahkan bit ini selanjutnya untuk membahas radius, karena pertanyaan telah berubah sejak jawaban awal saya.]
Perhatikan bahwa jika kita memiliki distribusi yang seragam di atas lingkaran satuan (yaitu yang kita tulis di alun-alun sebelumnya) maka kepadatan jari-jari untuk itu akan sebanding dengan jari-jari (pertimbangkan luas elemen annular kecil lebar di radius - yaitu antara dan - memiliki luas proporsional dengan ). Kemudian ketika kita melewati di luar lingkaran, daerah annular baru dengan jari-jari lebih besar hanya mendapatkan kontribusi kepadatan dari bagian di alun-alun, sehingga kepadatan berkurang (awalnya cukup cepat, kemudian lebih lambat) antara dan . (Sekali lagi, pengertian geometris yang cukup sederhana cukup untuk mendapatkan bentuk fungsional kepadatan jika diperlukan.)dr r r r+dr r 1 2–√
Sebaliknya, jika distribusi sambungan simetris secara rotasi tentang asal maka elemen probabilitas pada beberapa sudut tidak tergantung pada sudut (ini pada dasarnya adalah sebuah tautologi!). Distribusi bivariat dari dua standar Gaussians secara simetris berputar tentang asal:
(kode untuk gambar ini berdasarkan kode Elan Cohen di sini tetapi ada alternatif yang bagus di sini , dan sesuatu di antara keduanya di sini )
Akibatnya volume yang terkandung dalam beberapa sudut adalah sama untuk setiap , sehingga kepadatan yang terkait dengan sudut adalah seragam pada .dθ θ [0,2π)
[Trik polar biasanya digunakan untuk mengintegrasikan kerapatan normal di atas garis nyata dapat digunakan untuk mengetahui bahwa kerapatan jari-jari kuadrat adalah eksponensial negatif, dan dari sana kerapatan jari-jari mudah diidentifikasi dengan argumen transformasi sederhana dari fungsi distribusi]
sumber
Saya akan menjawab pertanyaan tentang kasus normal yang mengarah ke distribusi seragam. Diketahui bahwa jika dan Y adalah independen dan terdistribusi normal, kontur kerapatan probabilitas konstan adalah sebuah lingkaran pada bidang x - y . Jari-jari R = √X Y x−y memilikidistribusi Rayleigh. Untuk diskusi yang baik tentang ini, artikel wikipedia berjudul Distribusi Rayleigh.R=X2+Y2−−−−−−−√
Sekarang mari kita lihat variabel acak dan Y menggunakan koordinat kutub.X Y
sumber
Wilayah yang menarik untuk pertanyaan kami adalah sektor merah pada gambar ini:
Verifikasi:
sumber