randn.h

来自「实现矩阵加、减、乘和求逆的C++程序」· C头文件 代码 · 共 31 行

/*产生一个服从均值为mu,方差为sigma的正态分布随机数*/
#include <stdio.h>
#include <stdlib.h>   
    
double random(void)
{  
	int a;   
    double r;   
    
    a=rand()%32767;   
    r=(a+0.00)/32767.00;   
    
    return r;   
 }   
    
double randn(double mu,double sigma)
{   
    int i;   
    double r,sum=0.0;   
    
    if(sigma<=0.0)
	{
		printf("Sigma<=0.0!");exit(1);
	}   
    for(i=1;i<=12;i++)   
    sum=sum+random();   
    r=(sum-6.00)*sigma+mu;   
    
    return r;   
}