正态分布函数实现

//产生(min,max)之间均匀分布的随机数

double AverageRandom(double min,double max)                        

{

    int MINnteger = (int)(min*10000);

    int MAXnteger = (int)(max*10000);

    int randInteger = rand()*rand();

    int diffInteger = MAXnteger - MINnteger;

    int resultInteger = randInteger % diffInteger + MINnteger;

    return resultInteger/10000.0;

}

//正态分布概率密度函数

double Normal(double x,double miu,double sigma)                        

{

    return 1.0/(sqrt(2*3.14)*sigma) * exp(-1*(x-miu)*(x-miu)/(2*sigma*sigma));

}

//产生正态分布随机数

double Random_Normal(double miu,double sigma,double min,double max)    

{

    double x;

    double dScope;

    double y;

    do

    {

        x = AverageRandom(min,max);

        y = Normal(x, miu, sigma);

        dScope = AverageRandom(0, Normal(miu,miu,sigma));

    }while( dScope > y);

    return x;

}