|
|
问一个浅显问题,请教大家:
c! [+ ~ Y; t; Z$ q4 X3 L: G$ E# E
对于服从对称三角概率密度分布的随机数,比如以(a,0),(b,0),((a+b)/2,2/(b-a))三点决定的三角分布,假设b>a>0,要找出这样的随机数系列。我用反变换法求得结果是:
( ^9 o1 i: k: E# R, C% V& ?1)当a<=x<(a+b)/2时,x=a+(b-a)sqrt(u/2)* p* m1 x# ^* w; Y! s# j9 G! P
2)当(a+b)/2<=x<b时,x=b-(b-a)sqrt(u/2)
1 p# ?+ e" d* R2 V其中,u是服从(0,1)均匀分布的随机数) ^5 J/ r1 \( |1 l2 M$ w( K
5 x" N/ ~6 t" ~5 A. [- G! B0 z
能否把两者合起来?& |" }1 l; V' n
我看书上讲可用如下公式求服从这种分布的随机数:x=a+(b-a)*(u1+u2)/25 y' C- c7 q2 {5 {6 y3 q2 [+ S8 Q
不知道以上结果是怎么来的? |
|