作者shyfang (超值水果)
看板NCTU-STAT94G
标题[情报] 统算第二题的演算法
时间Mon Oct 17 09:33:09 2005
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)