|
|
问一个浅显问题,请教大家:2 k+ `/ }) z5 C, R# K2 ?3 \
6 j+ \( w/ e. z5 m n" f: w R对于服从对称三角概率密度分布的随机数,比如以(a,0),(b,0),((a+b)/2,2/(b-a))三点决定的三角分布,假设b>a>0,要找出这样的随机数系列。我用反变换法求得结果是:
4 U0 d* L7 F* \$ S# t4 g1)当a<=x<(a+b)/2时,x=a+(b-a)sqrt(u/2)
1 s( }* `9 ^: Y0 ]/ }$ L1 Y; t2)当(a+b)/2<=x<b时,x=b-(b-a)sqrt(u/2)
P; z1 A, g4 H* {9 ]其中,u是服从(0,1)均匀分布的随机数9 Q- j) @0 p. ^2 p& Z4 T! @/ }; `
7 {1 `1 F7 g/ m) a5 b5 g, ]; ]2 h
能否把两者合起来?
0 P7 B/ n% n6 N6 H C5 a# \我看书上讲可用如下公式求服从这种分布的随机数:x=a+(b-a)*(u1+u2)/2
3 X& ?3 T5 F/ q6 b不知道以上结果是怎么来的? |
|