产生正态分布(高斯分布)随机数 C语言实现

正态分布的概率密度函数如下

f(x) = \frac{1}{\sqrt{2\pi} \sigma }*e^{-(y-\mu)^2/2{\sigma}^2}

通常使用N(\mu,{\sigma}^2)表示均值为\mu,方差为\sigma^2

产生正太分布的方法如下:

r_{1},r_{2},...r_{n}为(0,1)上n个相互独立的均匀分布的随机数,由于E(r_{i})=\frac{1}{2},D(r_i)=\frac{1}{12}.

根据中心极限定理可知,当n趋向于无穷时

x=\sqrt{\frac{12}{n}}(\sum_{r=1}^{n}r_{i}-\frac{n}{2})

的分布近似于正态分布N(0,1),通常取n=12,

最后通过变换y=\mu+\sigma x,便可以得到均值为\mu,方差为\sigma^2的正太分布随机数。

C语言程序代码如下:

1、guass.h头文件代码

#pragma once
#ifndef GUASS_H_
#define GUASS_H_

/*
函数功能:	产生正太分布随机数
输入参数说明:
mean		给定均值
sigma		给定方差
seed	长整型指针变量, *seed 为伪随机数的种子
*/
double guass_data(double mean, double sigma, long int * seed);


#endif

2、guass.c文件

#include "guass.h"

#include <stdint.h>
#include "string.h"
#include "stdio.h"
#include "uniform.h"

/*
函数功能:	产生正太分布随机数
输入参数说明:
mean		给定均值
sigma		给定方差
seed	长整型指针变量, *seed 为伪随机数的种子
*/
double guass_data(double mean, double sigma, long int * seed)
{
	int i;
	double x, y;

	for (x=0, i = 0; i < 12; i++)
	{
		x += uniform_data(0.0,1.0,seed);
	}
	x = x - 6;
	y = mean + x * sigma;
	return y;
}

3、主程序Main.代码

#include <stdio.h>
#include <string.h>
#include "uniform.h"
#include "guass.h"

int main()
{
	// 产生50个0-1之间的随机数
	
	int i, j;
	long int s;

	double x, mean, sigma;
	mean = 0.0; sigma = 1.0; s = 13579;

	for (i = 0; i < 10; i++)
	{
		for (j = 0; j < 5; j++)
		{
			x = guass_data(mean, sigma, &s);
			printf("%13.7f",x);
		}
		printf("\n");
	}
	
	getchar();		// 此行代码是为了保持输出窗口,按任意按键关闭弹窗
	return 0;
}

4、程序运行结果

注:程序中包含的头文件 uniform.h和函数文件uniform.c 见文章产生(a,b)区间上均匀分布的随机数 C语言实现_jk99528878769的博客-CSDN博客使用C语言实现 在设定区间内产生一个均匀分布的随机数https://blog.csdn.net/jk99528878769/article/details/124435414?spm=1001.2014.3001.5501

matlab 验证结果

待添加。。。

  • 2
    点赞
  • 39
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 3
    评论
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

jk99528878769

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值