Ini adalah kelanjutan langsung dari pertanyaan terakhir saya . Hal yang sebenarnya ingin saya dapatkan adalah distribusi , di mana seragam dalam . Sekarang, distribusi berhasil dihitung di utas yang disebutkan , dan sebut saja . Distribusi hanyalah . Langkah terakhir adalah menghitung distribusi jumlah dan dengan cara yang mirip dengan yang sebelumnya , tetapi dan tidak mandiri, dan sekarang saya terjebak dan bahkan tidak tahu harus mulai dari mana.
Mungkin berguna untuk mencatat bahwa dan pada komponen terakhir komponen di bawah root (yaitu, dan ) mudah untuk dihitung. Kemudian, saya tertarik pada distribusi , mengetahui distribusi dan .
Saya tidak melihat perubahan variabel yang berguna. Saya berpikir tentang menggunakan probabilitas bersyarat, tetapi bagaimana saya bisa menemukan ? Saya mungkin terlalu maju dan mungkin harus mundur beberapa langkah.
Apakah mungkin menghitung sesuatu seperti ini?
Distribusi yang dihasilkan akan terlihat seperti ini:
EDIT: Jawaban yang diterima memberikan solusi yang saya cari, namun, saya masih penasaran bagaimana mendapatkannya secara analitis. Maksud saya, dalam pertanyaan saya sebelumnya CDF diberikan sebagai satu kesatuan:
dengan dan diberikan oleh fungsi sederhana. Secara teoritis, itu bisa diintegrasikan menggunakan pena dan kertas. Tentu saja menggunakan perangkat lunak itu wajar. Namun, saya masih penasaran bagaimana memberikan jawaban dalam bentuk tertutup di sini. serigala menjawab dering bel, tapi ... Konvolusi dari tiga pdf dari fungsi yang (relatif) rumit seperti itu?
sumber
Tepat setelah membaca jawaban serigala, saya mengerti bahwa saya dapat menghitung distribusi akhir dari awal tanpa semua langkah tengah:
M[x_] := M[x] = Evaluate@FullSimplify@ Integrate[ Boole[a + d + Sqrt[(a - d)^2 + 4 b c] <= x], {a, 0, 1}, {b, 0, 1}, {c, 0, 1}, {d, 0, 1}]
memberikan CDF danm[x_] := m[x] = Evaluate@FullSimplify@D[M[x], x]
memberikan PDF yang berfungsi sempurna dengan simulasi saya:Ini menggunakan pendekatan jawaban untuk pertanyaan saya sebelumnya.
sumber
Integrate[ Boole[(a-d)^2 + 4 b c < x], {a,0,1}, {b,0,1}, {c,0,1}, {d,0,1}]
mengembalikan integral yang tidak dievaluasi.