• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    迪恩网络公众号

逆分布函数法生成随机数(指数分布) R语言实现

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

 

先说明一下符号:U(0,1)-均匀分布,”~“表示服从xxx分布,F(x),为需要生成的随机数的分布函数,invF(x)表示逆分布函数,那么算法步骤如下:

step 1: 产生 u~U(0,1)

step 2:计算X=invF(u)

那么X就是服从分布函数为F(x)的随机数,一个简单的证明:

        证明: P(X<y)=P(invF(u)<y)=P(u<F(y))=F(y); 证毕。

你会感觉好简单啊,对是很简,下面就给出常用的指数分布
       首先,指数分布的分布函数:

F(x)=1-exp(-lamda*x);x>0( lamda 为参数)

        因此逆分布函数:

invF(x)=-(1/lamda)*ln(1-x);0<x<1

       所以随机数数产生过程:

step 1:产生u~U(0,1)

step 2: 计算X=-ln(1-u)

得到随机数

 

 

 

x <- runif(100)
lamda <- 0.1
U <- runif(100)
lamda <- 0.1
X <- -1/lamda * log(U)
hist(X,prob=T,col=gray(0.9),
main="exp from uniform")
curve(dexp(x,lamda),add=T,col="red")

  


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
R语言笔记-set.seed()函数发布时间:2022-07-18
下一篇:
R语言runif生成均匀随机数有问题??发布时间:2022-07-18
热门推荐
热门话题
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap