通常,它有助于使用累积分布函数。
第一,
F(x)=Pr((a−d)2≤x)=Pr(|a−d|≤x−−√)=1−(1−x−−√)2=2x−−√−x.
下一个,
G(y)=Pr(4bc≤y)=Pr(bc≤y4)=∫y/40dt+∫1y/4ydt4t=y4(1−log(y4)).
令介于(a - d )2 + 4 b c的最小(0)和最大(5)可能值之间。写入X = (一- d )2与CDF ˚F和ÿ = 4 b Ç与PDF 克= g ^ ',我们需要计算δ05(a−d)2+4bcx=(a−d)2Fy=4bcg=G′
H(δ)=Pr((a−d)2+4bc≤δ)=Pr(x≤δ−y)=∫40F(δ−y)g(y)dy.
我们可以预料这很讨厌- 均匀分布的PDF不连续,因此应该在的定义上产生破绽-因此Mathematica获得封闭形式(在此我不再赘述)有点令人惊讶。相对于δ求差得到所需的密度。在三个间隔内分段定义。在0 < δ < 1时Hδ0<δ<1
H′(δ)=h(δ)=18(8δ√+δ(−(2+log(16)))+2(δ−2δ√)log(δ)).
在中1<δ<4
h(δ)=14(−(δ+1)log(δ−1)+δlog(δ)−4δ√coth−1(δ√)+3+log(4)).
并且在,4<δ<5
h(δ)=14(δ−4δ−4−−−−√+(δ+1)log(4δ−1)+4δ√tanh−1((δ−4)δ−−−−−−√−δ√δ−δ−4−−−−√)−1).
该图将的曲线图覆盖在(a - d )2 + 4 b c的10 6个 iid实现的直方图中。两者几乎无法区分,这表明h公式的正确性。h106(a−d)2+4bch
以下是几乎没有头脑的,蛮力的Mathematica解决方案。它几乎可以自动执行有关计算的所有操作。例如,它甚至可以计算结果变量的范围:
ClearAll[ a, b, c, d, ff, gg, hh, g, h, x, y, z, zMin, zMax, assumptions];
assumptions = 0 <= a <= 1 && 0 <= b <= 1 && 0 <= c <= 1 && 0 <= d <= 1;
zMax = First@Maximize[{(a - d)^2 + 4 b c, assumptions}, {a, b, c, d}];
zMin = First@Minimize[{(a - d)^2 + 4 b c, assumptions}, {a, b, c, d}];
这就是所有的整合和差异化。(请耐心;计算需要花费几分钟。)H
ff[x_] := Evaluate@FullSimplify@Integrate[Boole[(a - d)^2 <= x], {a, 0, 1}, {d, 0, 1}];
gg[y_] := Evaluate@FullSimplify@Integrate[Boole[4 b c <= y], {b, 0, 1}, {c, 0, 1}];
g[y_] := Evaluate@FullSimplify@D[gg[y], y];
hh[z_] := Evaluate@FullSimplify@Integrate[ff[-y + z] g[y], {y, 0, 4},
Assumptions -> zMin <= z <= zMax];
h[z_] := Evaluate@FullSimplify@D[hh[z], z];
最后,模拟并与的图进行比较:h
x = RandomReal[{0, 1}, {4, 10^6}];
x = (x[[1, All]] - x[[4, All]])^2 + 4 x[[2, All]] x[[3, All]];
Show[Histogram[x, {.1}, "PDF"],
Plot[h[z], {z, zMin, zMax}, Exclusions -> {1, 4}],
AxesLabel -> {"\[Delta]", "Density"}, BaseStyle -> Medium,
Ticks -> {{{0, "0"}, {1, "1"}, {4, "4"}, {5, "5"}}, Automatic}]