对数正态分布的概率密度函数为
对数正态分布均值为,方差为
首先产生正态分布的随机变量 y ,然后通过变换产生对数正态分布的随机变量 x 。具体方法如下:
(1)产生正态分布的随机数 y ,即
(2)计算
1、lognorm.h
#pragma once
#ifndef LOGNORM_H_
#define LOGNORM_H_
/*
函数功能: 产生对数的随机数
输入参数说明:
u 对数正太分布的参数
sigma 对数正太分布的参数
seed 长整型指针变量, *seed 为伪随机数的种子
*/
double lognorm_data(double u, double sigma, long int * seed);
#endif // !LOGNORM_H_
2、lognrm.c 文件
#include "lognorm.h"
#include <stdint.h>
#include "string.h"
#include "stdio.h"
#include <math.h>
#include "guass.h"
/*
函数功能: 产生对数的随机数
输入参数说明:
u 对数正太分布的参数
sigma 对数正太分布的参数
seed 长整型指针变量, *seed 为伪随机数的种子
*/
double lognorm_data(double u, double sigma, long int * seed)
{
double x, y;
y = guass_data(u, sigma, seed);
x = exp(y);
return x;
}
高斯分布随机数的产生有关代码见以下链接
3、main.c
#include <stdio.h>
#include <string.h>
#include "lognorm.h"
int main()
{
// 产生50个对数正太分布的随机数
int i, j;
long int s;
double u, x, sigma;
u = 0.0; sigma = 0.5; s = 13579;
for (i = 0; i < 10; i++)
{
for (j = 0; j < 5; j++)
{
x = lognorm_data(u, sigma, &s);
printf("%13.7f",x);
}
printf("\n");
}
getchar(); // 此行代码是为了保持输出窗口,按任意按键关闭弹窗
return 0;
}
4、运行输出结果