|
|
问一个浅显问题,请教大家: N5 r7 g. U: K$ x% k ~/ P& {
X) f) e J( P, Z' R对于服从对称三角概率密度分布的随机数,比如以(a,0),(b,0),((a+b)/2,2/(b-a))三点决定的三角分布,假设b>a>0,要找出这样的随机数系列。我用反变换法求得结果是:, g* Z! d6 X+ g$ Z' R7 G# s% G
1)当a<=x<(a+b)/2时,x=a+(b-a)sqrt(u/2)6 g; {6 }( ^) ?0 _7 x+ y0 M
2)当(a+b)/2<=x<b时,x=b-(b-a)sqrt(u/2)
$ Q6 q$ z9 B6 l4 q+ x9 x% k其中,u是服从(0,1)均匀分布的随机数
* K( J7 y# Z* I5 ~9 d/ Y+ x# A9 i$ t1 S) J! P
能否把两者合起来?
" z! s. W' W9 R' l4 X我看书上讲可用如下公式求服从这种分布的随机数:x=a+(b-a)*(u1+u2)/2
3 I4 Z0 Q7 u: ^6 K t7 V4 |" y+ ?不知道以上结果是怎么来的? |
|