连续型分布随机数的产生.docxVIP

连续型分布随机数的产生（反函数法） 一、定理：设随机变量的分布函数为，其为严格单调的连续函数，其反函数存在，则服从上的均匀分布。 也即，若，则 二、例 设，分布函数为 则，，或若， 产生分布随机数的算法 （1）产生 （2）计算 即为分布的随机数。 如：产生10000个分布随机数的Ｒ程序代码为 lamda-1/4 N-10000 M-c(rep(0,N)) for (k in 1:N) { u-runif(1) x--log(1-u)/lamda M[k]-x } M hist(M) mean(M) var(M) 实验六 分布随机数的产生 编制程序产生10000个分布的随机数，并完成 （1）绘制该10000个随机数的直方图； （2）求该10000个随机数的样本均值及样本方差。 ############卡方分布随机数############## lamda-1/2 N-10000 M-c(rep(0,N)) for(K in 1:N){ {n-5 m-c(rep(0,n)) for(k in 1:n){ u-runif(1) x--log(1-u)/lamda m[k]-x m1-sum(m[1:k]) } m1 } M[K]-m1 } M hist(M) mean(M) var(M) ############指数分布随机数########### lamda-1/2 N-10000 M