码迷,mamicode.com
首页 > 其他好文 > 详细

使用均匀分布随机数产生器构造任意分布的随机数产生器

时间:2015-08-16 13:51:41      阅读:255      评论:0      收藏:0      [点我收藏+]

标签:均匀分布

一.算法原理

     在概率论中,有下面重要的结论:设随机变量Y的概率分布函数F(y)为连续函数,而X是在(0,1)上均匀分布的随机变量,设技术分享技术分享,则Z与Y有相同的分布。
    因此,理论上,我们可以用(0,1)上的均匀分布来得到任意分布的随机数
    设随机变量Y的概率密度函数为f(y),则,我们可以按照下面的步骤来获得满足其分布的随机数:
1.根据给定的概率密度函数(PDF)f(y),通过积分算出概率分布函数F(y);
2.计算概率分布函数F(y)的反函数技术分享
3.生成一个(0,1)上均匀分布的随机数变量X,
4.那么技术分享就是满足要求的分布的随机数。



二.使用均匀分布产生指数分布

1.指数分布是概率密度函数(PDF)为
技术分享
  通过计算可以得到,指数分布的概率分布函数为
技术分享

2.指数分布的概率分布函数的反函数为
技术分享

3.获得一个在(0,1)之间均匀分布的随机数x,则通过下面的式子计算得到的y满足指数分布
技术分享



2.MATLAB代码
 
%% clear
clc;
clear;
 
%% uniform distribution
N = 100000000;
U1 = rand(1, N);
 
%% inverse function
lamda = 1;
Z1 = -1 / lamda * log(1 - U1);
 
 
 
%% visualisztion of PDF
subplot(1, 2, 1), hist(U1, 100), title( ‘The PDF of U1‘);
subplot(1, 2, 2), hist(Z1, 100), title( ‘The PDF of Z1‘);
  


3.实验结果
技术分享








参考文献
徐全智,数学建模


版权所有,欢迎转载,转载请注明出处,谢谢技术分享


版权声明:本文为博主原创文章,未经博主允许不得转载。

使用均匀分布随机数产生器构造任意分布的随机数产生器

标签:均匀分布

原文地址:http://blog.csdn.net/sheng_ai/article/details/47701523

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!