[求助]C语言中如何产生gauss分布的随机数

adminftp · 发表于 2005-9-5 17:25

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

C语言中如何产生gauss分布的随机数？

VibInfo · 发表于 2005-9-6 11:17

一理论上，计算任意f分布的随机数： 假如你要产生分布函数为f(x)的随机数序列,首先获得一个均匀分布序列A, 同时求出f(x)的反函数g(x), 然后将序列A的元素代入g(x)中, 然后计算出新的一个序列B. B序列就满足f(x)分布. 二 Box-Muller 变换方法 y1 = sqrt( - 2 ln(x1) ) cos( 2 pi x2 ) y2 = sqrt( - 2 ln(x1) ) sin( 2 pi x2 ) 如果只想得到一个则y2不需计算 说明：x1,x2是0到1的均匀分布随机数。 由于需要计算三角函数，所以速度很慢，而且当x1接近0时不稳定。最好采用Box-Muller的极坐标方法 三 Box-Muller的极坐标方法 float x1, x2, w, y1, y2; do { x1 = 2.0 * ranf() - 1.0; x2 = 2.0 * ranf() - 1.0; w = x1 * x1 + x2 * x2; } while ( w >= 1.0 ); w = sqrt( (-2.0 * ln( w ) ) / w ); y1 = x1 * w; y2 = x2 * w; 代码： 代码: #ifndef _M_GAUSSRANDOM_2004_6_28_ZHY #define _M_GAUSSRANDOM_2004_6_28_ZHY #pragma warning(disable:4786) #include <vector> #include <cmath> using namespace std; #include <time.h> namespace Statical { class CUniRandom { public: CUniRandom():s(65536.0),u(2053.0),v(13849.0){} void SetSeed(const double& seed) { r = seed; } void GetRandom(double& result) { r = u * r + v; m = (int)(r/s); r = r - m * s; result = r/s; } private: double s,u,v,r; int m; }; void Uniform_Random(double& r,vector<double>& vResult,int nNum) { vResult.resize(nNum); int i,m; double s = 65536.0,u = 2053.0,v = 13849.0; for(i = 0; i <= nNum-1; i++) { r = u * r + v; m = (int)(r/s); r = r - m * s; vResult = r/s; } } void Gauss_Random(double& y) { float x1, x2, w; CUniRandom randContainer; randContainer.SetSeed(time(NULL)); do { double r; randContainer.GetRandom(r); x1 = 2.0 * r - 1.0; randContainer.GetRandom(r); x2 = 2.0 * r - 1.0; w = x1 * x1 + x2 * x2; } while ( w >= 1.0 ); w = sqrt( (-2.0 * log( w ) ) / w ); y = x1 * w; } void Gauss_Random(double& y1,double& y2) { float x1, x2, w; CUniRandom randContainer; randContainer.SetSeed(time(NULL)); do { double r; randContainer.GetRandom(r); x1 = 2.0 * r - 1.0; randContainer.GetRandom(r); x2 = 2.0 * r - 1.0; w = x1 * x1 + x2 * x2; } while ( w >= 1.0 ); w = sqrt( (-2.0 * log( w ) ) / w ); y1 = x1 * w; y2 = x2 * w; } void Gauss_Random(vector<double>& result,int nNum) { result.resize(nNum); float x1, x2, w; CUniRandom randContainer; randContainer.SetSeed(time(NULL)); for(int i = 0; i < nNum; i++) { do { double r; randContainer.GetRandom(r); x1 = 2.0 * r - 1.0; randContainer.GetRandom(r); x2 = 2.0 * r - 1.0; w = x1 * x1 + x2 * x2; } while ( w >= 1.0 ); w = sqrt( (-2.0 * log( w ) ) / w ); result = x1 * w; } } }; #endif

[此贴子已经被作者于2005-9-6 11:18:05编辑过]

VibInfo · 发表于 2005-9-6 11:18

上面的F分布的，改动一下就可以了

账号		自动登录	找回密码
密码			我要加入

[C/C++] [求助]C语言中如何产生gauss分布的随机数

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（beyond）[求助]C语言中如何产生gauss分布的随...

回复：（beyond）[求助]C语言中如何产生gauss分布的随...