|
|
问一个浅显问题,请教大家:
$ c9 g2 I5 T* L, j( k4 _$ y& ~5 ~% S
对于服从对称三角概率密度分布的随机数,比如以(a,0),(b,0),((a+b)/2,2/(b-a))三点决定的三角分布,假设b>a>0,要找出这样的随机数系列。我用反变换法求得结果是:
. }2 e% d! ^7 D& p1)当a<=x<(a+b)/2时,x=a+(b-a)sqrt(u/2), N) {/ B8 c' q& J& r# n6 N3 I
2)当(a+b)/2<=x<b时,x=b-(b-a)sqrt(u/2) ]* Z* ~+ V4 L% w3 U# i6 ~6 b
其中,u是服从(0,1)均匀分布的随机数+ j7 g# [8 I4 ?: O
0 S/ R+ w0 q" @( a% y! j% ?% x
能否把两者合起来?0 j% j' y+ U1 h/ L4 c
我看书上讲可用如下公式求服从这种分布的随机数:x=a+(b-a)*(u1+u2)/2
0 w2 g7 S7 I' V5 C不知道以上结果是怎么来的? |
|