|
|
问一个浅显问题,请教大家:
7 g2 d, B( N# p" X/ _1 D8 f7 j% J) g8 T6 s- P z
对于服从对称三角概率密度分布的随机数,比如以(a,0),(b,0),((a+b)/2,2/(b-a))三点决定的三角分布,假设b>a>0,要找出这样的随机数系列。我用反变换法求得结果是:; D! j0 f3 @1 d- a0 E1 ?8 [$ v* G$ A
1)当a<=x<(a+b)/2时,x=a+(b-a)sqrt(u/2)1 R0 ]9 M& [3 Z( j$ Z/ A9 W
2)当(a+b)/2<=x<b时,x=b-(b-a)sqrt(u/2). @+ y/ m7 |6 A; q* {
其中,u是服从(0,1)均匀分布的随机数/ a% G; { A1 P
3 r) c6 U; Y! H/ Z1 x+ l% \9 `/ Z能否把两者合起来?1 X2 n& s0 q5 Q+ c
我看书上讲可用如下公式求服从这种分布的随机数:x=a+(b-a)*(u1+u2)/2! @/ X7 O+ p5 m' O" K: {/ ]$ Q
不知道以上结果是怎么来的? |
|