java正态分布随机数_正态分布的随机数

一、功能

产生正态分布\(N(\mu, \ \sigma^2)\)。

二、方法简介

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

\[f(x)=\frac{1}{\sqrt{2\pi}\sigma}e^{-(x-\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}}\left ( \sum_{n}^{i=1}r_{i}-\frac{n}{2} \right )

\]

的分布近似正态分布\(N(0, \ 1)\)。通常取\(n=12\)，此时有

\[x=\sum_{12}^{i=1}r_{i}-6

\]

最后，再通过变换\(y=\mu+\sigma x\)，便可得到均值\(\mu\)、方差为\(\sigma^2\)的正态分布随机数\(y\)。

三、使用说明

使用C语言编程生成正态分布函数\(N(0, \ 1)\)

/************************************

a---给定区间下限

b---给定区间上限

seed---随机数种子

************************************/

#include "uniform.c"

double gauss(double mean, double sigma, long int *s)

{

int i;

double x;

double y;

for(x = 0, i = 0; i < 12; i++){

x += uniform(0.0, 1.0, s);

}

x = x - 6.0;

y = mean + x * sigma;

return(y);

}

uniform.c文件参见均匀分布的随机数