24小时热门版块排行榜

返回列表

dubo

金虫 (著名写手)

[交流] 产生高斯分布随机数的原理

for (i=0;iatoms_num;i=i+2)
{
for (a=0;a<3;a++)
{
do{
u1=(real)rand()/RAND_MAX;
u2=(real)rand()/RAND_MAX;
v1=2.0*u1-1.0;
v2=2.0*u2-1.0;
s=v1*v1+v2*v2;
}while(s>=1.0);
r=-2.0*log(s)/s;
work->wx[a]=v1*sqrt(r);
work->wx[i+1][a]=v2*sqrt(r);
}
}

回复此楼

» 猜你喜欢

听说今天filecode变了已经有25人回复
基金中了已经有4人回复
filecode 已经有10人回复
filecode 已经有12人回复
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑已经有5人回复
一个有机合成实验室都需要哪些设备？已经有12人回复
求各位大神看下已经有33人回复
咨询面上基金已经有7人回复
影响面上的因素已经有11人回复
8月时间戳变的，举个手。玩一下，释放压力已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

高斯批处理已经有8人回复
高斯计算的问题？？已经有5人回复
请教向量Q满足高斯分布，用matlab如何表示？已经有11人回复
MATLAB中产生高斯白噪声，涉及到awgn和wgn函数已经有3人回复
dmol 中的GGA-RPBE相当于高斯中的什么计算方法？已经有5人回复
【求助】pdos高斯展开作图已经有4人回复
【求助】Ms如何导入通过高斯优化好的结构已经有9人回复
【求助】高斯 ESP 已经有5人回复
求助高斯分布、对数正态分布及威布尔（weibull）分布模型表达式及采用什么软件已经有5人回复
洛仑兹线性分布和高斯线性分布有啥不同？已经有8人回复
【求助】如何用Matlab产生高斯白噪声已经有6人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-05-31 14:59:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ryhzd

木虫 (正式写手)

程序强帖: 1
应助: 2 (幼儿园)
金币: 3244.8
帖子: 546
在线: 173.4小时
虫号: 235381

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+2): 鼓励交流！这个不要太认真哈！版主会视情况编辑的！ 2011-05-31 20:20:52
dubo(金币+5): 谢谢，哈，金币给你吧 2011-06-06 16:24:41

这个是The Box-Muller transformation的改进方法，称为Polar Method，while循环里面取代了Box-Muller方法中的sin和cos函数，从而提高了速度。
ps：很不明白，很多人明明是求助帖，却写的是交流贴，如果是交流贴，只是应该有你给大家分享的东西。

赞一下(3人)

回复此楼

2楼2011-05-31 18:42:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

gwyy

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2201.7
帖子: 109
在线: 63.8小时
虫号: 115170

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
jjdg(金币+1): 感谢参与 2011-06-01 01:39:42
dubo(金币+5): 谢谢，哈，金币给你吧 2011-06-06 16:24:48

在Java的Random类中就是采用这种方法生成正态（高斯）分布随机数的，在JFC帮助文档中有说明。该方法的详细原理可以参看Knuth所著的“计算机程序设计艺术”一书的数值算法卷。生成正态随机数有很多不同方法，这种方法最大好处就是速度快，对于较大规模的随机抽样、仿真和模拟有很大用处。学物理的经常用Monte Carlo，应该很熟。

赞一下(2人)

回复此楼

3楼2011-05-31 20:40:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖