Saya memiliki empat variabel bebas terdistribusi seragam , masing-masing dalam . Saya ingin menghitung distribusi . Saya menghitung distribusi menjadi (maka ), dan dari menjadi Sekarang, distribusi jumlah u_1 + u_2 adalah ( u_1, \, u_2 juga independen) f_ {u_1 + u_2} (x) = \ int _ {- \ infty} ^ {+ \ infty} f_1 (xy) f_2 (y) dy = - \ frac {1} {4} \ int_0 ^ 4 \ frac {1- \ sqrt {xy}} {\ sqrt {xy}} \ cdot \ ln \ frac {y} {4} dy, karena y \ in (0,4]
Saya membuat empat set independen yang terdiri dari nomor masing-masing dan menarik histogram :
dan menggambar plot :
Secara umum, plotnya mirip dengan histogram, tetapi pada interval sebagian besar negatif (root berada di 2,27034). Dan integral dari bagian positif adalah .
Dimana kesalahannya? Atau di mana saya melewatkan sesuatu?
EDIT: Saya memperbesar histogram untuk menampilkan PDF.
EDIT 2: Saya pikir saya tahu di mana masalah dalam alasan saya - dalam batas integrasi. Karena dan , saya tidak bisa hanya . Plot menunjukkan wilayah yang harus saya integrasikan di:x - y ∈ ( 0 , 1 ] ∫ x 0
Ini berarti saya memiliki untuk (itu sebabnya bagian dari saya benar), dalam , dan in . Sayangnya, Mathematica gagal menghitung dua integral terakhir (well, ia menghitung yang kedua, dengan ada unit imajiner dalam output yang merusak segalanya ... ). y ∈ ( 0 , 1 ] f ∫ x x - 1 y ∈ ( 1 , 4 ] ∫ 4 x - 1 y ∈ ( 4 , 5 ]
EDIT 3: Tampaknya Mathematica CAN dapat menghitung tiga integral terakhir dengan kode berikut:
(1/4)*Integrate[((1-Sqrt[u1-u2])*Log[4/u2])/Sqrt[u1-u2],{u2,0,u1},
Assumptions ->0 <= u2 <= u1 && u1 > 0]
(1/4)*Integrate[((1-Sqrt[u1-u2])*Log[4/u2])/Sqrt[u1-u2],{u2,u1-1,u1},
Assumptions -> 1 <= u2 <= 3 && u1 > 0]
(1/4)*Integrate[((1-Sqrt[u1-u2])*Log[4/u2])/Sqrt[u1-u2],{u2,u1-1,4},
Assumptions -> 4 <= u2 <= 4 && u1 > 0]
yang memberikan jawaban yang benar :)