产生满足正态分布的随机数

一般有两种算法：

算法一产生12个（0，1）平均分布的随机函数，用大数定理可以模拟出正态分布。

算法二用到了数学中的雅可比变换，直接生成正态分布，但此算法在计算很大规模的数时 会出现溢出错误。

附加代码：

<!--

Code highlighting produced by Actipro CodeHighlighter (freeware)
http://www.CodeHighlighter.com/
-->

#include <math.h>


#include <stdio.h>


#include <conio.h>


#include <stdlib.h>


#include <time.h>




#define M_PI 3.14159265358979323846




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;




}