PS還沒試過 有問題再幫我補
1 先準備一個 pdf的副程式 anypdf(x_0.005,x_0.995,x)
x0.005是0.005的quantile當作下界 另一個是上界 x是等會要丟的變數
2 對pdf做積分 找到 F(x<X_0.1)....F(x<X_0.9) 的 X_%
3 生成 1個 u=unif(0,1)
4 找到一個max X_%1 < u
x_%
5 對pdf積分從X_%1 積到x_ans s.t F(x_ans) ~=u
6 goto 3 until n r.v generated
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.114.48
※ 編輯: shyfang 來自: 140.113.114.46 (10/17 16:39)