|
|
问一个浅显问题,请教大家:
% N0 G' p. C1 p/ ] U# n
6 m2 [0 [2 y- |- a0 t对于服从对称三角概率密度分布的随机数,比如以(a,0),(b,0),((a+b)/2,2/(b-a))三点决定的三角分布,假设b>a>0,要找出这样的随机数系列。我用反变换法求得结果是:" s* Y& Y2 S8 ~8 Z
1)当a<=x<(a+b)/2时,x=a+(b-a)sqrt(u/2)* J5 M. J; V, o3 { D1 v4 @3 p! O! c
2)当(a+b)/2<=x<b时,x=b-(b-a)sqrt(u/2)# [# [# T; j; S8 n8 q
其中,u是服从(0,1)均匀分布的随机数
7 y6 X, L6 [9 y' ]% J0 X& ?' U0 u) g, E- C1 w/ l5 @! ?
能否把两者合起来?
8 r5 i( _! @" t% ~& A3 O我看书上讲可用如下公式求服从这种分布的随机数:x=a+(b-a)*(u1+u2)/2
9 f; F% u2 P, ~ ?! N- |* J$ p }$ Q- R2 H不知道以上结果是怎么来的? |
|