对数正态分布的随机数产生 C语言实现

对数正态分布的概率密度函数为

$f(x)=\left\{\begin{matrix} \frac{1}{x \sqrt{2\pi} \sigma}exp(-\frac{(lnx - \mu)^2}{2\sigma ^2}) & x>0\\ 0 & x \leq0 \end{matrix}\right.$

对数正态分布均值为 $e^{\mu + \sigma^2/2}$ ，方差为 $(e^{\sigma^2} -1)e^{2\mu + \sigma^2}$

首先产生正态分布的随机变量 y ，然后通过变换 $x=e^y$ 产生对数正态分布的随机变量 x 。具体方法如下：
(1）产生正态分布的随机数 y ，即 $y\sim N(\mu, \sigma)$
(2）计算 $x=e^y$

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;
}

高斯分布随机数的产生有关代码见以下链接

产生正态分布（高斯分布）随机数 C语言实现_jk99528878769的博客-CSDN博客C语言实现产生正态分布（高斯分布）随机数https://blog.csdn.net/jk99528878769/article/details/124457234?spm=1001.2014.3001.5501

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、运行输出结果

对数正态分布的随机数产生 C语言实现相关推荐

根据对数正态分布产生随机数
系统自带的对数正态分布函数lognrnd example m = 1; % 分布的均值 v = 2; % 分布的方差 mu = log((m^2)/sqrt(v+m^2)); % 公式中的均值参数 s ...
R语言对数正态分布分布函数Log Normal Distribution（dlnorm, plnorm, qlnorm rlnorm）实战
R语言对数正态分布函数Log Normal Distribution(dlnorm, plnorm, qlnorm & rlnorm)实战目录 R语言对数正态分布函数Log Normal D ...
R语言作业一：矩估计、极大似然估计、拟合、对数正态分布、泊松分布、负二项分布
一.矩估计.极大似然估计.拟合.对数正态分布 ##导入数据 setwd("C:/Users/chang/Documents/SRM-PA/R简介/上课练习数据集") healthe ...
如何产生正态分布的随机数？
如何产生正态分布的随机数? 添加评论分享按投票排序按时间排序 28 个回答 86赞同反对,不会显示你的姓名 Milo Yip,游戏程序员.<游戏引擎架构>译者 Porphyah.熊墩 ...
对应于对数正态分布的拉丁超立方抽样——Python版
拉丁超立方抽样-对数正态分布 0.拉丁超立方抽样的理论基础 0.1.概况 0.2.基本原理 0.3.基本步骤 1.导入库和基本准备 2.生成(具有对数正态分布的随机变量)参数的随机数 3.将生成的随机 ...
使用 C / C++ 产生符合正态分布的随机数
使用 C / C++ 产生符合正态分布的随机数利用有box 和 muller 提供的,在 knuth的网上讨论过的方法 (比较常用的方法) 目的:Box-Muller,一般是要得到服从正态分布的随机 ...
R可视化绘制对数正态分布（Log Normal Distribution）
R可视化绘制对数正态分布(Log Normal Distribution) 为了绘制R中对数正态分布的概率密度函数,我们可以使用以下函数: dlnorm(x, meanlog = 0, sdlog = ...
对数正态分布_正态分布与幂律分布
1.正态分布英国生物统计学家法兰西斯·高尔顿做了一个实验.他在一块木板上画了一块等腰三角形,并在三角形区域内钉上n+1层钉子.第1层钉2个钉子,第2层钉3个钉子,下面每一层都比上一层增加一个钉子,上 ...
[数学]二维对数正态分布的概率分布，期望，方差和相关系数
最近遇到了一个联合对数正态分布的相关系数的问题,搜遍全网无果,索性自己动手.本文借鉴了这个知乎回答首先我们有二维正态分布: X , Y ∼ B V N ( μ x , μ y , σ x 2 , σ ...
机器学习小组知识点17：对数正态分布（Logarithmic Normal Distribution）
这块儿我是真的没听说过,所以直接抄了维基百科,维基万岁! 概率密度函数在概率论与统计学中,对数正态分布是对数为正态分布的任意随机变量的概率分布.如果 X X 是正态分布的随机变量,则 exp(X)e ...

对数正态分布的随机数产生 C语言实现

对数正态分布的随机数产生 C语言实现相关推荐

最新文章

热门文章