C语言生成正态随机变量,用C语言程序生成符合正态分布的随机数列

2007-10-05 11:33

如何用C语言程序生成符合正态分布的随机数列？

一般有两种算法：       算法一产生12个(0，1)平均分布的随机函数，用大数定理可以模拟出正态分布。       算法二用到了数学中的雅可比变换，直接生成正态分布，但此算法在计算很大规模的数时       会出现溢出错误。       测试程序：       #include           #include           #include           #include           #include                double    _random(void)       {         int    a;         double    r;              a=rand()%32767;         r=(a+0.00)/32767.00;              return    r;            }            double    _sta(double    mu,double    sigma)       {         int    i;         double    r,sum=0.0;              if(sigma<=0.0)    {    printf("Sigma<=0.0    in    _sta!");    exit(1);    }         for(i=1;i<=12;i++)         sum    =    sum    +    _random();         r=(sum-6.00)*sigma+mu;              return    r;            }            double    _sta2(double    mu,double    sigma)       {         double    r1,r2;              r1=_random();         r2=_random();              return    sqrt(-2*log(r1))*cos(2*M_PI*r2)*sigma+mu    ;            }                 int    main()       {       int    i;       double    mu,sigma;            srand(    (unsigned)time(    NULL    )    );       mu=0.0;       sigma=1.0;       printf("Algorithm    1:\n");       for(i=0;i<10;i++)       printf("%lf\t",_sta(mu,sigma));       printf("Algorithm    2:\n");       for(i=0;i<10;i++)       printf("%lf\t",_sta2(mu,sigma));       return    0;            }