正态分布及matlab实现

来源：互联网发布：dnf一直网络连接中断编辑：程序博客网时间：2024/06/06 01:13

正态分布（Normal distribution）又名高斯分布（Gaussian distribution），是一个在数学、物理及工程等领域都非常重要的概率分布，在统计学的许多方面有着重大的影响力。

若随机变量 $X$ 服从一个数学期望为 $μ$ 、标准方差为 $σ 2$ 的高斯分布，记为：

X \sim N (μ,σ 2),

则其概率密度函数为

$f(x) = {1 \over \sigma\sqrt{2\pi} }\,e^{- {{(x-\mu )^2 \over 2\sigma^2}}}$

正态分布的期望值 $μ$ 决定了其位置，其标准差 $σ$ 决定了分布的幅度。因其曲线呈钟形，因此人们又经常称之为钟形曲线。我们通常所说的标准正态分布是 $μ = 0,σ = 1$ 的正态分布（见右图中绿色曲线）。

1 概要
- 1.1 历史
2 正态分布的定义
- 2.1 概率密度函数
- 2.2 累积分布函数
- 2.3 生成函数
  - 2.3.1 动差生成函数
  - 2.3.2 特征函数
3 性质
- 3.1 标准化正态随机变量
- 3.2 矩(英文:moment)
- 3.3 生成正态随机变量
- 3.4 中心极限定理
- 3.5 无限可分性
- 3.6 稳定性
- 3.7 标准偏差
4 正态测试
5 相关分布
6 参量估计
- 6.1 参数的极大似然估计
  - 6.1.1 概念一般化
- 6.2 参数的矩估计
7 常见实例
- 7.1 光子计数
- 7.2 计量误差
- 7.3 生物标本的物理特性
- 7.4 金融变量
- 7.5 寿命
- 7.6 测试和智力分布
8 计算统计应用
- 8.1 生成正态分布随机变量
9 参见
10 引用条目
11 外部连接

[编辑]概要

正态分布是自然科学与行为科学中的定量现象的一个方便模型。各种各样的心理学测试分数和物理现象比如光子计数都被发现近似地服从正态分布。尽管这些现象的根本原因经常是未知的，理论上可以证明如果把许多小作用加起来看做一个变量，那么这个变量服从正态分布(在R.N.Bracewell的Fourier transform and its application中可以找到一种简单的证明)。正态分布出现在许多区域统计:例如, 采样分布均值是近似地正态的，既使被采样的样本总体并不服从正态分布。另外，常态分布信息熵在所有的已知均值及方差的分布中最大，这使得它作为一种均值以及方差已知的分布的自然选择。正态分布是在统计以及许多统计测试中最广泛应用的一类分布。在概率论，正态分布是几种连续以及离散分布的极限分布。

[编辑]历史

常态分布最早是亚伯拉罕·棣莫弗在1734年发表的一篇关于二项分布文章中提出的。拉普拉斯在1812年发表的《分析概率论》（Theorie Analytique des Probabilites）中对棣莫佛的结论作了扩展。现在这一结论通常被称为棣莫佛－拉普拉斯定理。

拉普拉斯在误差分析试验中使用了正态分布。勒让德于1805年引入最小二乘法这一重要方法；而高斯则宣称他早在1794年就使用了该方法，并通过假设误差服从正态分布给出了严格的证明。

“钟形曲线”这个名字可以追溯到Jouffret他在1872年首次提出这个术语"钟形曲面"，用来指代二元正态分布（bivariate normal）。正态分布这个名字还被Charles S. Peirce、Francis Galton、Wilhelm Lexis在1875分布独立的使用。这个术语是不幸的，因为它反应和鼓励了一种谬误，即很多概率分布都是正态的。（请参考下面的“实例”）

这个分布被称为“正态”或者“高斯”正好是Stigler名字由来法则的一个例子，这个法则说“没有科学发现是以它最初的发现者命名的”。

[编辑]正态分布的定义

有几种不同的方法用来说明一个随机变量。最直观的方法是概率密度函数，这种方法能够表示随机变量每个取值有多大的可能性。累积分布函数是一种概率上更加清楚的方法，但是非专业人士看起来不直观（请看下边的例子）。还有一些其他的等价方法，例如cumulant、特征函数、动差生成函数以及cumulant-生成函数。这些方法中有一些对于理论工作非常有用，但是不够直观。请参考关于概率分布的讨论。

[编辑]概率密度函数

四个不同参数集的概率密度函数（绿色线代表标准正态分布）

正态分布的概率密度函数均值为 $μ$ 方差为 $σ 2$ (或标准差 $σ$ )是高斯函数的一个实例：

$f(x;\mu,\sigma)=\frac{1}{\sigma\sqrt{2\pi}} \, \exp \left( -\frac{(x- \mu)^2}{2\sigma^2} \right)$ 。

(请看指数函数以及 $π$ .)

如果一个随机变量 $X$ 服从这个分布，我们写作 $X$ ~ $N (μ,σ 2)$ . 如果 $μ = 0$ 并且 $σ = 1$ ，这个分布被称为标准正态分布，这个分布能够简化为

$f(x) = \frac{1}{\sqrt{2\pi}} \, \exp\left(-\frac{x^2}{2} \right)$ 。

右边是给出了不同参数的正态分布的函数图。

正态分布中一些值得注意的量：

密度函数关于平均值对称
平均值是它的众数（statistical mode）以及中位数（median）
函数曲线下68.268949%的面积在平均值左右的一个标准差范围内
95.449974%的面积在平均值左右两个标准差 $2σ$ 的范围内
99.730020%的面积在平均值左右三个标准差 $3σ$ 的范围内
99.993666%的面积在平均值左右四个标准差 $4σ$ 的范围内
反曲点（inflection point）在离平均值的距离为标准差之处

[编辑]累积分布函数

上图所示的概率密度函数的累积分布函数

累积分布函数是指随机变量 $X$ 小于或等于 $x$ 的概率，用密度函数表示为

$F(x;\mu,\sigma)=\frac{1}{\sigma\sqrt{2\pi}}\int_{-\infty}^x \exp \left( -\frac{(x - \mu)^2}{2\sigma^2}\ \right)\, dx.$

正态分布的累积分布函数能够由一个叫做误差函数的特殊函数表示：

$\Phi(z)=\frac12 \left[1 + \mathrm{erf}\,(\frac{z-\mu}{\sigma\sqrt2})\right] .$

标准正态分布的累积分布函数习惯上记为 $Φ$ ，它仅仅是指 $μ = 0$ ， $σ = 1$ 时的值，

$\Phi(x)=F(x;0,1)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^x\exp\left(-\frac{x^2}{2}\right)\, dx.$

将一般正态分布用误差函数表示的公式简化，可得：

$\Phi(z)=\frac{1}{2} \left[ 1 + \operatorname{erf} \left( \frac{z}{\sqrt{2}} \right) \right].$

它的反函数被称为反误差函数，为：

$\Phi^{-1}(p)=\sqrt2\;\operatorname{erf}^{-1} \left(2p - 1 \right).$

该分位数函数有时也被称为probit函数。probit函数已被证明没有初等原函数。

正态分布的分布函数Φ(x)没有解析表达式，它的值可以通过数值积分、泰勒级数或者渐进序列近似得到。

[编辑]生成函数

[编辑]动差生成函数

动差生成函数被定义为 $exp(tX)$ 的期望值。

正态分布的矩生成函数如下：

$M_X(t)\,$ $=\mathrm{E}\left( e^{tX}\right)$ $=\int_{-\infty}^{\infty} \frac {1} {\sigma \sqrt{2\pi} } e^{\left( -\frac{(x - \mu)^2}{2 \sigma^2} \right)} e^{tx}\, dx$ $=e^{\left( \mu t + \frac{\sigma^2 t^2}{2}\right)}$

可以通过在指数函数内配平方得到。

[编辑]特征函数

特征函数被定义为 $exp(itX)$ 的期望值，其中 $i$ 是虚数单位. 对于一个正态分布来讲，特征函数是：

$\phi_X(t;\mu,\sigma)\!$ $=\mathrm{E}\left[ \exp(i t X)\right]$ $=\int_{-\infty}^{\infty} \frac{1}{\sigma \sqrt{2\pi}} \exp \left(- \frac{(x - \mu)^2}{2\sigma^2} \right) \exp(i t x)\, dx$ $=\exp\left( i \mu t - \frac{\sigma^2 t^2}{2}\right).$

把矩生成函数中的 $t$ 换成 $it$ 就能得到特征函数。

[编辑]性质

正态分布的一些性质:

如果 $X \sim N(\mu, \sigma^2) \,$ 且 $a$ 与 $b$ 是实数，那么 $aX + b \sim N (a μ + b,(a σ) 2)$ (参见期望值和方差).
如果与是统计独立的正态随机变量，那么:
- 它们的和也满足正态分布 $U = X + Y \sim N(\mu_X + \mu_Y, \sigma^2_X + \sigma^2_Y)$ (proof).
- 它们的差也满足正态分布 $V = X - Y \sim N(\mu_X - \mu_Y, \sigma^2_X + \sigma^2_Y)$ .
- $U$ 与 $V$ 两者是相互独立的。
如果和是独立正态随机变量，那么:
- 它们的积 $XY$ 服从概率密度函数为 $p$ 的分布
  $p(z) = \frac{1}{\pi\,\sigma_X\,\sigma_Y} \; K_0\left(\frac{|z|}{\sigma_X\,\sigma_Y}\right),$ 其中 $K 0$ 是贝塞尔函数（modified Bessel function）
- 它们的比符合柯西分布，满足 $X / Y \simCauchy(0,σ X / σ Y)$ .
如果 $X_1, \cdots, X_n$ 为独立标准正态随机变量，那么 $X_1^2 + \cdots + X_n^2$ 服从自由度为n的卡方分布。

[编辑]标准化正态随机变量

[编辑]矩(英文:moment)

一些正态分布的一阶动差如下：

阶数原点矩中心矩累积量010 1

μ

μ

μ 2 + σ 2

σ 2

σ 2

μ 3 + 3μσ 2

004

μ 4 + 6μ 2 σ 2 + 3σ 4

3σ 4

正态分布的所有二阶以上的累积量为零。

[编辑]生成正态随机变量

[编辑]中心极限定理

主条目：中心极限定理

正态分布的概率密度函数，参数为μ = 12，σ = 3，趋近于n = 48、p = 1/4的二项分布的概率质量函数。

正态分布有一个非常重要的性质：在特定条件下，大量统计独立的随机变量的和的分布趋于正态分布，这就是中心极限定理。中心极限定理的重要意义在于，根据这一定理的结论，其他概率分布可以用正态分布作为近似。

参数为 $n$ 和 $p$ 的二项分布，在 $n$ 相当大而且 $p$ 不接近1或者0时近似于正态分布（有的参考书建议仅在 $np$ 与 $n (1 - p)$ 至少为5时才能使用这一近似）。

近似正态分布平均数为 $μ = np$ 且方差为 $σ 2 = np (1 - p)$ .

一泊松分布带有参数 $λ$ 当取样样本数很大时将近似正态分布 $λ$ .

近似正态分布平均数为 $μ = λ$ 且方差为 $σ 2 = λ$ .

这些近似值是否完全充分正确取决于使用者的使用需求

[编辑]无限可分性

正态分布是无限可分的概率分布。

[编辑]稳定性

正态分布是严格稳定的概率分布。

[编辑]标准偏差

深蓝色区域是距平均值小于一个标准差之内的数值范围。在正态分布中，此范围所占比率为全部数值之68%。根据正态分布，两个标准差之内（蓝，棕）的比率合起来为95%。根据正态分布，三个标准差之内（深蓝，橙，黄）的比率合起来为99%。

在实际应用上，常考虑一组数据具有近似于正态分布的概率分布。若其假设正确，则约68%数值分布在距离平均值有1个标准差之内的范围，约95%数值分布在距离平均值有2个标准差之内的范围，以及约99.7%数值分布在距离平均值有3个标准差之内的范围。称为"68-95-99.7法则"或"经验法则".

[编辑]正态测试

[编辑]相关分布

$R \simRayleigh(σ)$ 是瑞利分布，如果 $R = \sqrt{X^2 + Y^2}$ ，这里 $X \sim N (0,σ 2)$ 和 $Y \sim N (0,σ 2)$ 是两个独立正态分布。
$Y \sim \chi_{\nu}^2$ 是卡方分布具有 $ν$ 自由度，如果 $Y = \sum_{k=1}^{\nu} X_k^2$ 这里 $X k \sim N (0,1)$ 其中 $k=1,\dots,\nu$ 是独立的。
$Y \simCauchy(μ = 0,θ = 1)$ 是柯西分布，如果 $Y = X 1 / X 2$ ，其中 $X 1 \sim N (0,1)$ 并且 $X 2 \sim N (0,1)$ 是两个独立的正态分布。
$Y \simLog-N(μ,σ 2)$ 是对数正态分布如果 $Y = e X$ 并且 $X \sim N (μ,σ 2)$ .
与Lévy skew alpha-stable分布相关：如果 $X\sim \textrm{Levy-S}\alpha\textrm{S}(2,\beta,\sigma/\sqrt{2},\mu)$ 因而 $X \sim N(\mu,\sigma^2)$ .

截断正态分布.如果 $X \sim N(\mu, \sigma^2),\!$ ，在 $A$ 以下和 $B$ 以上截取X 将产生一个平均值 $E(X)=\mu + \frac{\sigma(\varphi_1-\varphi_2)}{T},\!$ 这里 $T=\Phi\left(\frac{B-\mu}{\sigma}\right)-\Phi\left(\frac{A-\mu}{\sigma}\right), \; \varphi_1 = \varphi\left(\frac{A-\mu}{\sigma}\right), \; \varphi_2 = \varphi\left(\frac{B-\mu}{\sigma}\right)$ ， $φ$ 是一个标准正态随机变量的密度函数

如果 $X$ 是一个正态分布的随机变量, $Y = | X |$ ，那么 $Y$ 具有折叠正态分布.

[编辑]参量估计

[编辑]参数的极大似然估计

[编辑]概念一般化

多元正态分布的协方差矩阵的估计的推导是比较难于理解的。它需要了解谱原理（spectral theorem）以及为什么把一个标量看做一个1×1 matrix的trace而不仅仅是一个标量更合理的原因。请参考协方差矩阵的估计（estimation of covariance matrices）.

[编辑]参数的矩估计

[编辑]常见实例

[编辑]光子计数

[编辑]计量误差

《饮料装填量不足与超量的概率》

某饮料公司装瓶流程严谨，每罐饮料装填量符合平均600毫升，标准差3毫升的常态分配法则。随机选取一罐，容量超过605毫升的概率？容量小于590毫升的概率

容量超过605毫升的概率 = p ( X > 605)= p ( ((X-μ) /σ) > ( (605 – 600) / 3) )= p ( Z > 5/3) = p( Z > 1.67) = 0.9525

容量小于590毫升的概率 = p (X < 590) = p ( ((X-μ) /σ) < ( (590 – 600) / 3) )= p ( Z < -10/3) = p( Z < -3.33) = 0.0004

《6-标准差(6-sigma或6-σ)的品质管制标准》

6-标准差(6-sigma或6-σ)，是制造业流行的品质管制标准。在这个标准之下，一个标准常态分配的变量值出现在正负三个标准差之外，只有2* 0.0013= 0.0026 (p (Z < -3) = 0.0013以及p(Z > 3) = 0.0013)。也就是说，这种品质管制标准的产品不良率只有万分之二十六。假设例3-16的饮料公司装瓶流程采用这个标准，而每罐饮料装填量符合平均600毫升，标准差3毫升的常态分配法则。预期装填容量的范围应该多少？ 6-标准差的范围 = p ( -3 < Z < 3)= p ( - 3 < (X-μ) /σ < 3) = p ( -3 < (X- 600) / 3 < 3)= p ( -9 < X – 600 < 9) = p (591 < X < 609) 因此，预期装填容量应该介于591至609毫升之间。

[编辑]生物标本的物理特性

[编辑]金融变量

[编辑]寿命

[编辑]测试和智力分布

《计算学生智商高低的概率》

假设某校入学新生的智力测验平均分数与方差分别为100与12。那么随机抽取50个学生，他们智力测验平均分数大于105的概率？小于90的概率？

本例没有常态分配的假设，还好中心极限定理提供一个可行解，那就是当随机样本长度超过30，样本平均数xbar近似于一个常态变量，因此标准常态变量Z = (xbar –μ) /σ/ √n。

平均分数大于105的概率 = p(Z> (105 – 100) / (12 /√50))= p(Z> 5/1.7) = p( Z > 2.94) = 0.0016

平均分数小于90的概率 = p(Z< (90 – 100) / (12 /√50))= p(Z < 5.88) = 0.0000

[编辑]计算统计应用

[编辑]生成正态分布随机变量

在计算机模拟中，经常需要生成正态分布的数值。最基本的一个方法是使用标准的正态累积分布函数的反函数。除此之外还有其他更加高效的方法，Box-Muller变换就是其中之一。另一个更加快捷的方法是ziggurat算法。下面将介绍这两种方法。一个简单可行的并且容易编程的方法是：求12个在（0,1）上均匀分布的和，然后减6(12的一半)。这种方法可以用在很多应用中。这12个数的和是Irwin-Hall分布；选择一个方差12。这个随即推导的结果限制在（-6,6）之间，并且密度为12，是用11次多项式估计正态分布。

Box-Muller方法是以两组独立的随机数U和V，这两组数在(0,1]上均匀分布，用U和V生成两组独立的标准正态分布随即变量X和Y:

$X = \sqrt{- 2 \ln U} \, \cos(2 \pi V) ,$

$Y = \sqrt{- 2 \ln U} \, \sin(2 \pi V)$ 。

这个方程的提出是因为二自由度的卡方分布（见性质4）很容易由指数随机变量（方程中的lnU）生成。因而通过随机变量V可以选择一个均匀环绕圆圈的角度，用指数分布选择半径然后变换成（正态分布的）x,y坐标。

正态分布的matlab实现：

1）使用MatLab画出正态分布的概率密度函数图像。
x=[-10:0.01:10];
y=normpdf(x,0,1);%正态分布函数。
figure;
axes1=axes('Pos',[0.1 0.1 0.85 0.85]);
plot(x,y);
set(axes1,'YLim',[-0.01 0.43],'XLim',[-3 3]);
图1：

2）验证概率密度函数在区间(-∞,∞)上的积分为1。
这里取参数mu=3，sigma=5（注：下文全用这两个参数）。
y='exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5)';
s=int(y,-inf,inf) %int积分函数（inf代表无穷大）。
输出：s=1

3）验证x=mu时取最大值。
思路：求解函数一阶导数为零的点。

* 求一阶导数
y='exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5)';
d=diff(y);%微分函数。
sd=simplify(d)
输出结果：sd = -1/250*(x-3)*exp(-1/50*(x-3)^2)*2^(1/2)/pi^(1/2)

* 通过图像判断解的位置
x=[0:0.001:40];
sd=-1/250.*(x-3).*exp(-1/50.*(x-3).^2).*2.^(1/2)./pi.^(1/2);
axes1=axes('Pos',[0.1 0.1 0.85 0.85]);
plot(x,sd);
set(axes1,'YLim',[-0.01 0.01],'XLim',[0 40]);
图2：

从图中可以看出在3附近有解。

* 定义函数并求解
function y=f(x)
y=-1/250.*(x-3).*exp(-1/50.*(x-3).^2).*2.^(1/2)./pi.^(1/2);

r=fzero('f',3)
输出：r=3

从图上看在x > 20以后，几乎是一条直线，若用20：
r=fzero('f',20)
输出：r=3

这说明是无限趋近于0。

* 进一步说明该点为最大值点
该概率密度函数一阶导数为0的解为3，此值正好为mu，再取x=1，x=4与x=3时的函数值比较。
>> normpdf(1,3,5)
ans = 0.0737

>> normpdf(4,3,5)
ans = 0.0782

>> normpdf(3,3,5)
ans = 0.0798

显然在x=3的两边函数值都比x=3小，说明该点为极大值点。根据正态分布函数的图像特点可知该点是最大值点。

4）验证x=mu +- sigma （即8或-2）处曲线有拐点。
思路：求二阶导数为零的点。

* 先求二阶微分
y='exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5)';
d=diff(y,2);%微分函数。
sd=simplify(d)
输出：sd = 1/6250*exp(-1/50*(x-3)^2)*2^(1/2)*(-16+x^2-6*x)/pi^(1/2)

* 通过图像判断解的位置
x=[-20:0.001:20];
sd=1./6250.*exp(-1/50.*(x-3).^2).*2.^(1/2).*(-16+x.^2-6.*x)/pi.^(1/2);
axes1=axes('Pos',[0.1 0.1 0.85 0.85]);
plot(x,sd);
set(axes1,'YLim',[-0.005 0.005],'XLim',[-20 20]);
图3：
从上图可以看出，曲线在(-5,0)和(5,10)之间分别都与y=0有交点，因此有两个解。

* 定义函数并求解
function y=f(x)
y=1./6250.*exp(-1/50.*(x-3).^2).*2.^(1/2).*(-16+x.^2-6.*x)/pi.^(1/2);

r=fzero('f',-5)
r = -2

>> r=fzero('f',5)
r = 8.0000

从而得到了两个拐点x=8和x=-2，也即mu +- sigma。

5）验证曲线以x轴为渐近线渐近线求解：
A 垂直渐近线 x=a是y=f(x)的渐近线<==>lim f(x)=∞或lim f(x)=∞
x->a+0 x->a-0
其中a在间断点中找——∞型第二类间断点B 水平渐近线
x→+∞（-∞）时，y=b是y=f(x)的渐近线<==>lim f(x)=b （或lim f(x)=b)
x->+∞ x->-∞

求其一阶倒数在x趋向于无穷大时的极限值b，若存在，即有水平渐近线y=b。
* 先定义函数：
function y=f(x)
syms x; %定义符号变量。
y=exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5);

* 求x趋向无穷大时一阶导数的极限
limit(f,inf)
ans = 0

6）验证 3 sigma 法则
思路：求解概率密度函数在[mu-3*sigma,mu+3*sigma]区间上的积分。
y='exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5)';
double(int(y,-12,18))
ans = 0.9973

0 0