【C++基于中心极限定理生成高斯随机数（高斯白噪声）】-程序员宅基地

C++基于中心极限定理生成高斯随机数（高斯白噪声）

前言

在c++的标准头文件中，‘math.h’并没有直接提供高斯随机数的生成函数。‘math.h’头文件主要包含了数学函数和常量的声明，例如三角函数、对数函数和指数函数等。要生成高斯随机数，你可以使用c++标准库中一个名为‘rand()’的函数，可以用于生成伪随机数，这个函数在‘<stdlib.h>’头文件中声明。生成的伪随机数服从均匀分布，通过采样足够多的样本，依据中心极限定理便可得到高斯随机数。

中心极限定理

中心极限定理是概率论中的一个重要的定理，它描述了当从任意分布中抽取大量样本时，样本均值的分布将近似服从于正态分布。

生成算法

高斯随机数又可以称为高斯白噪声，它是一种最为随机的噪声，存在于许多自然或工程领域中。首先定义rand(A)函数表示为能够产生在区间 $[0, A]$ 上的一个随机值的随机数发生器，它在该区间上满足均匀分布。一般区间右边界 $A$ 的值可取 $2^{15}-1$ （单精度浮点数）或 $2^{31}-1$ （双精度浮点数）。假设 $X$ 是一个随机数，将随机变量 $X$ 进行归一化可以构造出随机变量 $Z$ ：
$Z=\frac{X-E(X)}{\sqrt{Var(X)}}\qquad(1)$
式(1)中， $E (X)$ 为 $X$ 的数学期望， $Va r (X)$ 为 $X$ 的方差。假设对任意 $i = 1, ..., n$ ， $X_i$ 均在区间 $[0, A]$ 上服从均匀分布，因此有：
$E(X_i)=\frac{1}{2}A\qquad(2)$
$Var(X_i)=\frac{1}{12}A^2\qquad(3)$
$E(\sum_{i=1}^nX_i)=\frac{n}{2}A\qquad(4)$
$Var(\sum_{i=1}^nX_i)=\frac{n}{2}A^2\qquad(5)$

通过概率论中心极限定理生成服从高斯噪声的随机数。生成过程如下：对于公式(1)，假设 $X$ 是由 $n$ 个具有相同分布的随机变量相加得到的结果，记为 $X=\sum_{i=1}^nX_i$ ，当随机变量的样本足够大时，即 $n\rightarrow\infty$ ，公式(1)中随机变量 $Z$ 服从均值为0，方差为1的高斯分布。把服从高斯分布的随机变量记为 $Z_n$ ,
$Z_n=\frac{\sum_{i=1}^nX_i-\frac{1}{2}A}{\sqrt{\frac{n}{12}}A}=\frac{1}{A}\sqrt{\frac{12}{n}}\sum_{i=1}^{n}X_i-\sqrt{3n}\qquad(6)$

c++代码实现

代码需要包含如下头文件

#include <iostream>
#include <cstdlib>
#include<vector>
using namespace std;

生成高斯噪声代码

double generateGaussianNoise() {
    
	double Sum = 0.0;
	int numSample = 1000;  // 采样个数

	for (int i = 0; i < numSample; i++) {
    
		Sum += static_cast<double>(rand());
	}
	// 公式(6)
	double Res = (Sum - numSample / 2.0 * RAND_MAX) / (sqrt(numSample / 12.0) * RAND_MAX);
	return Res;
}

主函数

int main() {
    
	for (int i = 0; i < 1000; ++i) {
    
		double sample = generateGaussianNoise();
		cout << "Sample " << i << ": " << sample << endl;
	}
	return 0;
}

验证结果

利用generateGaussianNoise()生成多个随机数，求取其均值和方差，看其是否为0和1。
计算均值

double cal_mean(vector<double> val) {
    
	double length = static_cast<double>(val.size());
	double Sum = 0.0;
	for (double num : val)
		Sum += num;
	return Sum / length;
}

计算方差

double cal_Var(vector<double> val) {
    
	double length = static_cast<double>(val.size());
	double mean = cal_mean(val);
	double Sum = 0.0;
	for (double num : val)
		Sum += pow(num - mean, 2);
	return Sum / length;
}

主函数

int main() {
    
	vector<double> Res;
	for (int i = 0; i < 1000; ++i) {
    
		double sample = generateGaussianNoise();
		Res.push_back(sample);
		cout << "Sample " << i << ": " << sample << endl;
	}
	double mean = cal_mean(Res);
	double var = cal_Var(Res);
	cout << "mean: " << mean << endl;
	cout << "var: " << var << endl;
	return 0;
}

打印结果，其均值接近0，方差接近1，验证正确。
在这里插入图片描述

本文链接：https://blog.csdn.net/weixin_40603846/article/details/137114690

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

51单片机的中断系统_51单片机中断篇-程序员宅基地

文章浏览阅读3.3k次，点赞7次，收藏39次。CPU 执行现行程序的过程中，出现某些急需处理的异常情况或特殊请求，CPU暂时中止现行程序，而转去对异常情况或特殊请求进行处理，处理完毕后再返回现行程序断点处，继续执行原程序。void 函数名(void) interrupt n using m {中断函数内容 //尽量精简 }编译器会把该函数转化为中断函数，表示中断源编号为n，中断源对应一个中断入口地址，而中断入口地址的内容为跳转指令，转入本函数。using m用于指定本函数内部使用的工作寄存器组，m取值为0~3。该修饰符可省略，由编译器自动分配。_51单片机中断篇

oracle项目经验求职,网络工程师简历中的项目经验怎么写-程序员宅基地

文章浏览阅读396次。项目经验(案例一)项目时间：2009-10 - 2009-12项目名称：中驰别克信息化管理整改完善项目描述：项目介绍一，建立中驰别克硬件档案(PC,服务器，网络设备，办公设备等)二，建立中驰别克软件档案(每台PC安装的软件，财务，HR，OA，专用系统等)三，能过建立的档案对中驰别克信息化办公环境优化(合理使用ADSL宽带资源，对域进行调整，对文件服务器进行优化，对共享打印机进行调整)四，优化完成后..._网络工程师项目经历

LVS四层负载均衡集群-程序员宅基地

文章浏览阅读1k次，点赞31次，收藏30次。LVS：Linux Virtual Server，负载调度器，内核集成，阿里的四层SLB(Server Load Balance)是基于LVS+keepalived实现。NATTUNDR优点端口转换WAN性能最好缺点性能瓶颈服务器支持隧道模式不支持跨网段真实服务器要求anyTunneling支持网络private（私网）LAN/WAN（私网/公网）LAN（私网）真实服务器数量High (100)High (100)真实服务器网关lvs内网地址。

「技术综述」一文道尽传统图像降噪方法_噪声很大的图片可以降噪吗-程序员宅基地

文章浏览阅读899次。https://www.toutiao.com/a6713171323893318151/作者 | 黄小邪/言有三编辑 | 黄小邪/言有三图像预处理算法的好坏直接关系到后续图像处理的效果，如图像分割、目标识别、边缘提取等，为了获取高质量的数字图像，很多时候都需要对图像进行降噪处理，尽可能的保持原始信息完整性（即主要特征）的同时，又能够去除信号中无用的信息。并且，降噪还引出了一..._噪声很大的图片可以降噪吗

Effective Java 【对于所有对象都通用的方法】第13条谨慎地覆盖clone_为继承设计类有两种选择,但无论选择其中的-程序员宅基地

文章浏览阅读152次。目录谨慎地覆盖cloneCloneable接口并没有包含任何方法，那么它到底有什么作用呢？Object类中的clone()方法如何重写好一个clone()方法1.对于数组类型我可以采用clone()方法的递归2.如果对象是非数组，建议提供拷贝构造器（copy constructor）或者拷贝工厂（copy factory）3.如果为线程安全的类重写clone()方法4.如果为需要被继承的类重写clone()方法总结谨慎地覆盖cloneCloneable接口地目的是作为对象的一个mixin接口（详见第20_为继承设计类有两种选择,但无论选择其中的

毕业设计基于协同过滤的电影推荐系统-程序员宅基地

文章浏览阅读958次，点赞21次，收藏24次。今天学长向大家分享一个毕业设计项目基于协同过滤的电影推荐系统项目运行效果：项目获取：https://gitee.com/assistant-a/project-sharing21世纪是信息化时代，随着信息技术和网络技术的发展，信息化已经渗透到人们日常生活的各个方面，人们可以随时随地浏览到海量信息，但是这些大量信息千差万别，需要费事费力的筛选、甄别自己喜欢或者感兴趣的数据。对网络电影服务来说，需要用到优秀的协同过滤推荐功能去辅助整个系统。系统基于Python技术，使用UML建模，采用Django框架组合进行设

随便推点

你想要的10G SFP+光模块大全都在这里-程序员宅基地

文章浏览阅读614次。10G SFP+光模块被广泛应用于10G以太网中，在下一代移动网络、固定接入网、城域网、以及数据中心等领域非常常见。下面易天光通信（ETU-LINK）就为大家一一盘点下10G SFP+光模块都有哪些吧。一、10G SFP+双纤光模块10G SFP+双纤光模块是一种常规的光模块，有两个LC光纤接口，传输距离最远可达100公里，常用的10G SFP+双纤光模块有10G SFP+ SR、10G SFP+ LR，其中10G SFP+ SR的传输距离为300米，10G SFP+ LR的传输距离为10公里。_10g sfp+

计算机毕业设计Node.js+Vue基于Web美食网站设计(程序+源码+LW+部署)_基于vue美食网站源码-程序员宅基地

文章浏览阅读239次。该项目含有源码、文档、程序、数据库、配套开发软件、软件安装教程。欢迎交流项目运行环境配置：项目技术：Express框架 + Node.js+ Vue 等等组成，B/S模式 +Vscode管理+前后端分离等等。环境需要1.运行环境：最好是Nodejs最新版，我们在这个版本上开发的。其他版本理论上也可以。2.开发环境：Vscode或HbuilderX都可以。推荐HbuilderX;3.mysql环境：建议是用5.7版本均可4.硬件环境：windows 7/8/10 1G内存以上；_基于vue美食网站源码