产生随机数的代码

1.产生介于MIN与MIX之间的符合正态分布的代码:

/****************************************************
    产生N=100个在(0,50)区间内满足对数正态分布的随机数
/*****************************************************/

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <time.h>
#define N 100 
#define MAX 50
#define MIN 0
#define MIU 2
#define SIGMA 1
#define PI 3.14159
 
double AverageRandom(double min,double max)//产生(min,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 LogNormal(double x,double miu,double sigma) //对数正态分布概率密度函数
{
        return 1.0/(x*sqrt(2*PI)*sigma) * exp(-1*(log(x)-miu)*(log(x)-miu)/(2*sigma*sigma));
}
double Random_LogNormal(double miu,double sigma,double min,double max)//产生对数正态分布随机数
{
        double x;
        double dScope;
        double y;
        do
        {
        x = AverageRandom(min,max); 
        y = LogNormal(x, miu, sigma);
        dScope = AverageRandom(0, LogNormal(miu,miu,sigma));
        }while( dScope > y);
        return x;
}
void main(void)
{
    srand(time(NULL));
    int i,j;
    for(i=0,j=0;i<N;i++)
    {
        printf("%14f", Random_LogNormal(MIU,SIGMA,MIN,MAX));
        j++;
        if(j==5)//if 代码的意义是:使每行显示5个数
        {
            printf("\n");    
            j=0;
        }
 
    }
	system("pause");
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值