连续型分布随机数产生
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
连续型分布随机数的产生（反函数法）
一、定理：
二、例
（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<-c(rep(0,N))
for(k in 1:N){
u<-runif(1)
x<--log(1-u)/lamda
m
M[k]<-x
}
M
hist(M) ### 画直方图 ### mean(M) ### 均值 ###
var(M) ### 方差 ###