C++生成随机数:高斯/正态分布(gaussian/normal distribution) |
您所在的位置:网站首页 › 生成随机数的方法有哪些 › C++生成随机数:高斯/正态分布(gaussian/normal distribution) |
高斯分布也称为正态分布(normal distribution) 常用的成熟的生成高斯分布随机数序列的方法由Marsaglia和Bray在1964年提出,C++版本如下: #include #include double gaussrand() { static double V1, V2, S; static int phase = 0; double X; if ( phase == 0 ) { do { double U1 = (double)rand() / RAND_MAX; double U2 = (double)rand() / RAND_MAX; V1 = 2 * U1 - 1; V2 = 2 * U2 - 1; S = V1 * V1 + V2 * V2; } while(S >= 1 || S == 0); X = V1 * sqrt(-2 * log(S) / S); } else X = V2 * sqrt(-2 * log(S) / S); phase = 1 - phase; return X; } 这样生成的高斯分布随机数序列的期望为0.0,方差为1.0。若指定期望为E,方差为V,则只需增加: X = X * V + E; 期望E= wiki: http://zh.wikipedia.org/wiki/%E9%AB%98%E6%96%AF%E5%88%86%E5%B8%83 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |