|
|
问一个浅显问题,请教大家:
8 n' C. W, I- P: ^3 q }+ T$ s7 t3 N- M9 ^/ G
对于服从对称三角概率密度分布的随机数,比如以(a,0),(b,0),((a+b)/2,2/(b-a))三点决定的三角分布,假设b>a>0,要找出这样的随机数系列。我用反变换法求得结果是:
' I2 ]& d: Z& S& C/ {" R; o1)当a<=x<(a+b)/2时,x=a+(b-a)sqrt(u/2)
3 d0 c6 Q8 I8 D) s. \% [5 ?2)当(a+b)/2<=x<b时,x=b-(b-a)sqrt(u/2)- ]% a0 J* P! _; x" G) K2 J
其中,u是服从(0,1)均匀分布的随机数% n* b, `* B% v) p+ ^+ E- ^- j8 {
( A8 k' v1 D5 ^/ \7 I2 Y9 F能否把两者合起来?) }$ q# x/ ]6 f" Z" g' P6 h# ~
我看书上讲可用如下公式求服从这种分布的随机数:x=a+(b-a)*(u1+u2)/2
4 R$ s8 [* C& U: j: r0 V$ Z不知道以上结果是怎么来的? |
|