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

2007-10-05 11:33

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

一般有两种算法：   
   算法一产生12个（0，1）平均分布的随机函数，用大数定理可以模拟出正态分布。   
   算法二用到了数学中的雅可比变换，直接生成正态分布，但此算法在计算很大规模的数时   
   会出现溢出错误。   
   测试程序：   
   #include    <math.h>   
   #include    <stdio.h>   
   #include    <conio.h>   
   #include    <stdlib.h>   
   #include    <time.h>   
    
   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;   
    
   }