1-3、常用概率分布与随机数生成

杨利霞

1-3、常用概率分布与随机数生成本文首先介绍了概率论里在数学建模中的常用分布，包括连续型分布中的均匀分布、指数分布、正态分布以及离散型分布中的二项分布、泊松分布，并针对各种分布列举了一些对应的遵循该分布的例子。然后列举了十多种Matlab里按照各种要求生成随机数、随机序列的方法，并简要介绍了excel里生成随机数的方法。最后列举了一个需要用到生成随机数知识的例子。
4 y! F2 v* H  s3 n0 ?& [2 c
/ \; K' e4 f" n- Z4 _一、常用概率分布与服从该分布的事件举例

在实际生活中，有一些事情的发生遵循某些概率分布。因此可以用某些概率分布模型来刻画某些事件。本节简单介绍了一些概率论基础里的概率分布，与可用其刻画的一些事件。数学理论部分不详述，建议参考概率论教科书。

9 k/ ^  s4 H9 l0 t5 Z' t1、连续型分布

（1）、均匀分布(Uniform)

均匀分布是概率统计中的重要分布之一。顾名思义，均匀，表示可能性相等的含义。在实际问题中，当我们无法区分在区间[a,b]内取值的随机变量X取不同值的可能性有何不同时，我们就可以假定X服从[a,b]上的均匀分布。
6 n* U" E+ {$ E- {
. a& _' A3 Z" W! x/ C) c# w/ ?1 o概率密度函数：
! M2 R0 ^3 x/ ]4 M# v6 @f(x)={1b−a,0,a<x<botherwise
f(x)={1b−a,a<x<b0,otherwise
1 y# h$ }0 p5 ^4 Q7 J' ~* x' X
; v1 Z0 x* ]1 W: t- a+ Z0 a4 k
分布函数：
" B0 C# u# L6 EF(x)=⎧⎩⎨⎪⎪0,x−ab−a,1,x≤aa<x≤bx>b
F(x)={0,x≤ax−ab−a,a<x≤b1,x>b
  g: ~) P( U7 A, w& y8 h" c
" D& X3 s7 o+ `  A, v
$ j1 r5 Y7 T* y适用模型：

* m. ]: H: |2 q/ l* l3 c* Q* H在某区间内某事件发生的概率相等的问题。
  e4 F* @" w% x( m. [1
# j$ T" g; V9 z3 D' z4 F5 z（2）、指数分布
2 O5 j, S: e/ Q8 c7 N
; c) ]/ @" W( J: Q( i在概率理论和统计学中，指数分布（也称为负指数分布）是描述泊松过程中的事件之间的时间的概率分布，即事件以恒定平均速率连续且独立地发生的过程。 这是伽马分布的一个特殊情况。它是几何分布的连续模拟，它具有无记忆的关键性质。除了用于分析泊松过程外，还可以在其他各种环境中找到。
+ {2 t4 F6 R2 F) h/ i# s& J- F  z
概率密度函数：
f(x)={λe−λx,0,x>0x≤0
# W& N3 c: o7 @4 v8 q& E" df(x)={λe−λx,x>00,x≤0
- F( y# W9 E4 Y( ]4 b+ p

分布函数：

F(x)={1−e−λx,0,x≥0x<0
F(x)={1−e−λx,x≥00,x<0
9 c; c& ^) w. F

8 F$ `* F; D% u& H/ q6 D适用模型：
) Q: z$ [0 k* }% M- n" ~
在电子元器件的可靠性研究中，通常用于描述对发生的缺陷数或系统故障数的测量结果。这种分布表现为均值越小，分布偏斜的越厉害。可以近似地作为
& ^% `( H- |- K4 P; {% T* g- ~高可靠性的复杂部件、机器或系统的失效分布模型，特别是在部件或机器的整机试验中得到广泛的应用。
0 |: a, x. z- j/ b6 b2 [* [1
3 \- B" |. }3 z# T/ n( P# ~- E2
; A6 ^4 |( f0 N( j（3）、正态分布(Normal)（又称高斯分布）
) G: [- M, w+ z, @+ z- S
' m% _  B* n- J- Z4 Q& V6 Q/ ]正态分布是一个在数学、物理及工程等领域都非常重要的概率分布，在统计学的许多方面有着重大的影响力。正态曲线呈钟型，两头低，中间高，左右对称。因其曲线呈钟形，因此人们又经常称之为钟形曲线。
, J2 c9 q& R" I5 {7 Z- s, R& r若随机变量X服从一个数学期望为μ、方差为σ^2的正态分布，记为N(μ，σ^2)
( A( f: Q. E# K7 O" e8 ?/ |0 B& |。其概率密度函数为正态分布的期望值μ决定了其位置，其标准差σ决定了分布的幅度。当μ = 0,σ = 1时的正态分布是标准正态分布。

概率密度函数：
7 |8 D5 O9 O0 n" Sf(x)=12π−−√σe−(x−μ)22σ2
6 D4 K" I6 F  U  S! mf(x)=12πσe−(x−μ)22σ2
4 X+ U2 `/ v2 n
0 L$ Y# a' h8 N6 T
# T' Q. T+ m9 d/ \分布函数：
F(x)=12π−−√σ∫x−∞e−(t−μ)22σ2dt
* ~5 Q2 j3 G: U1 @) g% y' jF(x)=12πσ∫−∞xe−(t−μ)22σ2dt

% J" O& r& x, {. x8 \( V/ g$ g
适用模型：
0 y( P' N' Z: e! p) d+ _, _' s
7 \4 L" P& c4 N2 g+ E5 `8 n* \社会人群生活水平
; {) b; s8 l7 e; Z$ q) ]/ |人群身高分布
. y0 H5 R$ f( H通货膨胀率和能源价格
考试成绩及学生综合素质研究（教育统计学统计规律表明，学生的智力水平，包括学习能力，实际动手能力等呈正态分布）
8 {) [3 z; |5 @; F医学参考值（某些医学现象，如同质群体的身高、红细胞数、血红蛋白量，以及实验中的随机误差，呈现为正态或近似正态分布；有些指标（变量）虽服从偏态分布，但经数据转换后的新变量可服从正态或近似正态分布，可按正态分布规律处理）
8 l3 l9 B; e( P1
- j1 ?" S% c; ~8 a2
! b+ X. |. p2 ^3
4
3 F& M% ~) l  A3 M5
- T2 O$ t7 X) D+ }2、离散型分布
* E  B/ d  o6 A( M
6 ?# b3 X# D9 Y% M5 K6 k（1）、二项分布（伯努利概型）

二项分布就是重复n次独立的伯努利试验。在每次试验中只有两种可能的结果，而且两种结果发生与否互相对立，并且相互独立，与其它各次试验结果无关，事件发生与否的概率在每一次独立试验中都保持不变，则这一系列试验总称为n重伯努利实验，当试验次数为1时，二项分布服从0-1分布。
在n重伯努利试验中，设事件A发生的概率为p，事件A发生的次数是随机变量，设为X，则X可能取值为0,1,2，…，n

分布律：
P{X=k}=Cknpk(1−p)(n−k)(k=0,1,2,...,n)
& r: W- t# U% f5 vP{X=k}=Cnkpk(1−p)(n−k)(k=0,1,2,...,n)
n=20，p=0.7时分布图像如下：
. @6 C* O9 j, Q- z1 p9 p) u
# B, X$ l- o3 e( I! a9 L) \: B1 Z适用模型：
0 J( z) c. A+ s: x& A! K
打枪、投篮问题（实验n次发生k次）
8 v  n( U' _! |0 x$ H, b3 n设备使用设备故障等确定基数下发生或不发生问题
% T- X7 k  g( R7 w1
% Z6 p: y- t; n4 p7 L" k2
（2）、泊松分布（n趋于无穷时二项分布的近似）

3 H' }; ~2 s2 i* o1 T1 j: p+ e7 }+ d9 x泊松分布的参数λ是单位时间(或单位面积)内随机事件的平均发生率。 泊松分布适合于描述单位时间内随机事件发生的次数。可参考《随机过程》中的泊松过程进一步学习。

分布律：
5 k: j' E; U9 z8 s$ Y+ p4 O4 @P{X=k}=λkk!e−λ(k=0,1,2,...)
P{X=k}=λkk!e−λ(k=0,1,2,...)
. N; `+ d' a) ~; n" b% I
8 ~* S0 ?2 \( [8 q
+ R+ R; E/ N3 b  d
: Q1 G# H1 Y( M5 h% T' O) o2 ]. }