Tengo cuatro variables independientes uniformemente distribuidas , cada una en . Quiero calcular la distribución de . Calculé la distribución de para ser
Hice cuatro conjuntos independientes constaban de 10 6 números cada uno y dibujé un histograma de ( a - d ) 2 + 4 b c :
y dibujó una gráfica de :
Generalmente, la gráfica es similar al histograma, pero en el intervalo mayor parte es negativa (la raíz está en 2.27034). Y la integral de la parte positiva es ≈ 0.77 .
¿Dónde está el error? ¿O dónde me estoy perdiendo algo?
EDITAR: Escalé el histograma para mostrar el PDF.
EDIT 2: Creo que sé dónde está el problema en mi razonamiento, en los límites de integración. Debido y x - y ∈ ( 0 , 1 ] , no puedo simplemente ∫ x 0 El gráfico muestra la región tengo que integrar en:.
Esto significa que tengo para y ∈ ( 0 , 1 ] (es por eso que parte de mi f era correcta), ∫ x x - 1 en y ∈ ( 1 , 4 ] y ∫ 4 x - 1 en y ∈ ( 4 , 5. ] Desafortunadamente, Mathematica no puede calcular las dos últimas integrales (bueno, calcula la segunda, porque hay una unidad imaginaria en la salida que estropea todo ...).
EDITAR 3: Parece que Mathematica PUEDE calcular las últimas tres integrales con el siguiente código:
(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]
que da una respuesta correcta :)