![水文水资源系统风险分析(普通高等教育“十二五”规划教材 全国水利行业规划教材)](https://wfqqreader-1252317822.image.myqcloud.com/cover/756/23370756/b_23370756.jpg)
2.5 风险分析中常用的概率分布
2.5.1 几种离散型概率分布
1.(0-1)分布(又称两点分布)
(0-1)分布的分布列如下:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0038_0003.jpg?sign=1739671198-XWr39vvDKECMIwTYgtjWoaEazKsnEBD1-0-056506a8433e6771b113369727f3bb54)
若X服从(0-1)分布,则
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0038_0004.jpg?sign=1739671198-idee1SI1Ag2smL04OsWkomj7Cy8LuO9K-0-6aea257e7062248a6c96441720c8b01b)
2.二项分布
设试验只有两种可能结果:A、Ā,且每次试验中A发生的概率为p, Ā发生的概率为q=1-p,将试验重复n次,则称这种试验为伯努利试验。
在n次伯努利试验中,事件A发生的次数X是一个随机变量。可以证明,在n次试验中,事件A恰好发生k次的概率Pn(X=k)可用下式求得
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0038_0005.jpg?sign=1739671198-C8p5Y0NFpxf0zDSqmISKJIFKkhw4uive-0-e46f206973456ab31412d19138a0ad2c)
由于式(2-52)的右侧是二项式(p+q)n展开式的第k+1项,所以称此随机变量X服从二项分布。二项分布记为B(n, p)。
若X服从二项分布,则
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0038_0006.jpg?sign=1739671198-eaAZFN8NWQKqGDENW2Cb347FK2N2Zq44-0-2554cf45b2ee5b3c20924c8727521387)
3.泊松分布
如果随机变量X的可能取值为0,1,2, …,而(X=k)的概率为
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0038_0007.jpg?sign=1739671198-SYLekN2YOiha00JlAF8VbCSwJg9AW7Sl-0-668a47bc40b3483ef215a433c5084384)
则称X服从泊松分布,其中参数λ>0。泊松分布记为Pλ(k)。
若X服从泊松分布,则
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0039_0001.jpg?sign=1739671198-zvkJcE1fPHwQcsM0AMQdVgfjdtJ4snuV-0-4e3a0673c738ab5f16ce81220fea78ba)
可以证明,随着n的增大,二项分布与以λ=np为参数的泊松分布趋于接近,因此,当n很大时,可用泊松分布近似求出二项分布的值。
2.5.2 几种连续型概率分布
1.均匀分布U(a, b)
(1)密度函数:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0039_0002.jpg?sign=1739671198-yyeT34wwqdW9Gh6dXvntVgj5c0pZlTv0-0-0400c26361ee6a7fd58c2f3046b5ef58)
均匀分布的密度曲线见图2-3。
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0039_0003.jpg?sign=1739671198-nLAnxt8R5j2eFC91Ig5dnNQgTyYysppR-0-4bc07c8997f441401f0b0d7902957046)
图2-3 均匀分布密度曲线
(2)均值:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0039_0004.jpg?sign=1739671198-huNXB1V4HSyPsxlqEt65CS3ItFgZ6gPY-0-98538e32ce82336ca826819e24e31628)
(3)方差:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0039_0005.jpg?sign=1739671198-u4aa7ybVfyn6faEOKYJfU9BajXb5a2cy-0-0e556bd70c5b979a5413f39c341c3fa2)
若X在区间 [a, b]上服从均匀分布,则 X取值于[a, b]中任一小区间内的概率只与该小区间的长度成正比,而与小区间的具体位置无关。因此,X取值于 [a, b]内任意等大小区间内的概率都是相等的,这也是一种等可能性的意思。
均匀分布是一种常用的分布,它在统计仿真等方法中占有重要地位。
2.正态分布N(a, σ2)
(1)密度函数:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0039_0006.jpg?sign=1739671198-AwEML1AUv82bbY6C9ByGMsl2ACE3xz7G-0-36d906b4f3ce3921048f7bc6305630f0)
正态分布密度曲线见图2-4。
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0039_0007.jpg?sign=1739671198-NgD6J6VADaKjhHhmuzSc51zfeIzGxmEm-0-2b3e2380d4747f62d757fe12d508abdc)
图2-4 正态分布密度曲线
(2)均值:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0040_0001.jpg?sign=1739671198-du18FWVfhFPAWKZD9fXvYBVRPTcu4tt6-0-03b654c2cd66b777a8be7d5b5cafba85)
(3)方差:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0040_0002.jpg?sign=1739671198-MUx40Es8elnDdL5jKiJZKr0iQHiSOB1b-0-f0b71ec35e834e4d3153a84853b418c3)
正态分布中的参数a和σ2分别为均值和方差。
可以证明,若X服从正态分布N(a, σ2),则X的取值落在 [a-3σ, a+3σ]区间的概率为99.7%,即
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0040_0003.jpg?sign=1739671198-HD8LOyohtBJYMGXTvt2tDEcaC5dMkxpl-0-b235ce5ceae8bc31970eeb6030346658)
这就是重要的“3σ法则”。
特别地,当a=0, σ2=1时,N(a, σ2)称为标准化正态分布,记为N(0,1)。标准化正态分布的密度函数为:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0040_0004.jpg?sign=1739671198-KJs6PW5Oe79v6KqAZrv4PZYrT2sMdZvW-0-8f25df109ac49a401abff261bec81f9d)
标准化正态分布的分布函数为
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0040_0005.jpg?sign=1739671198-VZn41h0Qq8qKOlMEr7xBQFBqXG0Fu6Uq-0-c1fe93831b2b8722668667e94f421034)
式中:Φ(x)称为拉普拉斯函数,可由正态分布表查得。
一般正态分布可通过变换转化为标准化正态分布。引进变量代换:,则
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0040_0007.jpg?sign=1739671198-RxpZUQMTGPAUr006QksaFqGGBj7ByYbI-0-5e6a7f183f839f663c9f707ca4c7e099)
上式是一重要的关系式。要计算F(x),只要将其转化为,便可查标准化正态分布表了。
无论在理论上,还是在实际中,正态分布都起着重要的作用。许多随机变量的概率分布都可以用正态分布来描述。若一个随机变量受大量作用微小且相互独立的因素综合影响,那么它将服从或近似服从正态分布。
3.对数正态分布
(1)两参数对数正态分布ln(a, σ2)。
1)密度函数。若Y=lnX服从正态分布N(ay,),则X服从对数正态分布,其密度函数为:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0040_0010.jpg?sign=1739671198-Bk4sOcXgX86P6cJUGrqyyJGdFzWRHmPC-0-19e3133c1c05219be5b76e6ab8adf3c7)
式中:ay、σy分别为Y的均值、均方差。
对数正态分布密度函数曲线见图2-5。
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0001.jpg?sign=1739671198-TDN4KcqfoHzWxbFgA7gkAnAI8DPwjdHv-0-84c42ad442cc73dadb2dae8f6179910d)
图2-5 对数正态分布密度函数
2)均值:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0002.jpg?sign=1739671198-NNLEMaAbOUcAuxnMzz1XEChIoQAtH5Gd-0-b69d2d92a42670055ae9890cd633a939)
3)方差:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0003.jpg?sign=1739671198-vfmIHzqHW9Kgma3ya4tkLmV0Mv8PdfbY-0-d3b26fb499c5056b301ea483fb9d82e4)
当对数分布应用于经济资料时,常称为Cobb-Douglas分布。
(2)三参数对数正态分布ln(ay,, b)。
1)密度函数:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0005.jpg?sign=1739671198-Mf0kxHCwMpSYGrSolDczAQT1QUkn8JS9-0-526e0b94ac371c251dfd9a063529296a)
2)均值:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0006.jpg?sign=1739671198-I06cGPEN9Lb1HT2n0AfwI92X0HUxCUd8-0-d6cffca86db995dbad3ce1336310c0a2)
3)方差:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0007.jpg?sign=1739671198-9MgVTNhQGLWNCw9i4v2DyV4VrgNZLV5B-0-3ee480ebc05b6307d6fba34ef7031462)
对两参数的对数正态分布Y,采用平移变换X=Y+b,就可以得到三参数的对数正态分布。即两参数正态分布ln(ay,)可视为b=0的三参数对数正态分布的特例。
4.指数分布 e(λ)
(1)密度函数:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0009.jpg?sign=1739671198-QpYRQWppKMY2pk0wFtjIxrTpJN4dYR5A-0-5880fd3178693879e7db9ad9df1e5598)
式中参数λ>0。
(2)均值:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0010.jpg?sign=1739671198-To5vvItDlXjStv8RegWEBfzh1hFyl85i-0-c7be0a19cbe2f49c5bf5b6387631f363)
(3)方差:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0041_0011.jpg?sign=1739671198-uUgeqry5Vp5qgZDxMQLdEX0w6vqOqMV8-0-14071025f13012542a5a0f8637be6ebf)
5.三角分布Triangular(a, b, c)
(1)密度函数:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0001.jpg?sign=1739671198-4MuaoktVEZqaXjOZq5mfFsYvCdO6rZk5-0-acffabafaf42b23bf0608d2acca48b4a)
三角分布密度函数曲线见图2-6。
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0002.jpg?sign=1739671198-lRjeUTRmmfIdJBA41XXBDmIycnerRSMo-0-a0c219c9b48349efd50e2707151b1b15)
图2-6 三角分布密度函数
a<c<b,其中,a为位置参数;b-a为比例参数;c为形状参数。当c=b时,为右三角分布;当c=a时,为左三角分布。
(2)均值:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0003.jpg?sign=1739671198-U96QxSTU2ZZvPcP0VHp018bTgAC42Nle-0-e7f93474a5cc8dca1c5534fb457cfc45)
(3)方差:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0004.jpg?sign=1739671198-0yrsRTtKCsC2K4dais0kpHnI1vdGMWoS-0-f833ba51afa746245bc621736d6c2271)
三角分布在风险管理中也经常使用。
6.极值分布(Ⅰ型)G(u, α)
(1)密度函数:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0005.jpg?sign=1739671198-Mxx5XXAzhZOtv1xJuxW2V833TS73R3fS-0-4cce689490ef9be4c752235c0af9056a)
式中:u为位置参数;α>0为比例参数。
分布函数
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0006.jpg?sign=1739671198-bHnFZr7SeL0kDN6NoL1fQ8SGq6AwHP3A-0-3bc3be46ebf3ee38721bef28df87425e)
(2)均值:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0007.jpg?sign=1739671198-9sQNW2czywHEmyEjC1Vcl4dKvEU7bJWi-0-5ce8b33d5fb7af3ffed5f09a98e98372)
(3)方差:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0008.jpg?sign=1739671198-wGBomydNC6dennBdhLFW96Pt0RpIMMhc-0-1241235b480337efe9183a61dd9d2561)
极值分布指的是n次观测中的极大值或极小值的概率分布。理论上,极值分布有3种可能的渐近极值分布,此处介绍的极值分布(Ⅰ型)为指数原型极值分布,其在气象、水文和地震的风险分析中用得较为广泛。
7.P-Ⅲ型分布
(1)密度函数:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0009.jpg?sign=1739671198-b2TDNZqSF6QmO2128uIzxac94AcbbTOt-0-bb7eff8c6c18d1495c58549efcb9d32f)
式中:参数α>0, x>a0。
P-Ⅲ型分布密度曲线如图2-7所示。
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0001.jpg?sign=1739671198-cQqJ1nO9S3BdNBNSqScbT2DQV1k2R3dk-0-dd75a5163d9d77237909735148452cda)
图2-7 P-Ⅲ型分布密度函数
(2)均值:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0042_0010.jpg?sign=1739671198-Gur7pdA6k5vXsdhz1XkC0wEqILF4qJoE-0-525d15f1b321bec46f296cd1cf023daa)
(3)方差:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0002.jpg?sign=1739671198-W4gMzp8GSIRnVHShuTPGKRxwVdbJN1TG-0-7612dc9b68d18bacf3daf1946f2d445f)
在实际工作中,常要计算指定概率p所相应的随机变量取值xp,即求出满足下述等式的xp:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0003.jpg?sign=1739671198-ApSumxknS24qe83btRY3YavEwJgrCFdC-0-08ae58042e9a48235d84c5eb231b20b0)
亦即
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0004.jpg?sign=1739671198-yyTbJ6VP5L9fgLIfLbE6StqZkj2W6WC9-0-5595590dc7470e1675debbec648a1eff)
由上式可知,当α、β、a0已知时,xp只取决于p。α、β与a0与数字特征E(X)、Cv、Cs有下列关系:
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0005.jpg?sign=1739671198-MUiBKHDnFI1lIWKJtPkD06Iaz0T2k0A6-0-d5099bc64a8ba1c9e52e20a22d39250d)
因此只要E(X)、Cv与Cs一经确定,xp仅与p有关。但是直接由积分式来计算是非常困难的,实用上,通过查算已制成的专用表可以使计算工作大大简化。
令
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0006.jpg?sign=1739671198-sH7vwqVTEgjGuBpkMPM59h6oHkdwavqp-0-533be4d7382d952353cf0fd9088aa110)
由上式可得
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0007.jpg?sign=1739671198-PCD6ziXgTP7RPqZs5KrzmAPdqQBraf59-0-4261b9c781c096ff54f45fb6e4ea20fd)
将之代入式(2-86),并同时将式(2-86)中的α, β, a0以相应的E(X), Cv, Cs来表示,通过简化可得
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0008.jpg?sign=1739671198-w3iI1D35SuSnhMBKXdjfZ2cUKYQTPLjW-0-2dccb188ac2c2c23ff9c46f1d5aae46c)
上式等号右边与参数α有关而由式(2-87)可知,α是Cs的函数,因此,可对若干给定的Cs值,编制φp和p的对应数值表。此表已先后由美国工程师福斯特和苏联工程师雷布京制订出来。常称为P-Ⅲ型分布离均系数φ值表。
例2设某地年雨量X的分布符合P-Ⅲ型分布,且E(X)=1000mm, Cv=0.5, Cs=1.0,试求该地百年一遇的年降雨量。
解:所谓百年一遇的年降雨量,即求满足下式的xp。
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0009.jpg?sign=1739671198-ZKcSRic48sSgnMtOcDV1FDt7loeMAmPl-0-80daa6bd413c76b8d790f0a64d2c8792)
由Cs=1.0, p=0.01查P-Ⅲ型分布离均系数表,得
φp=3.02
所以
![](https://epubservercos.yuewen.com/4652D5/11989104804194506/epubprivate/OEBPS/Images/figure_0043_0010.jpg?sign=1739671198-BkRDPne53ovZhNpgo79yaEd35PswftcF-0-fc272ebbb5b6c92e591c0655add9bc9a)
P-Ⅲ型分布适应性较强,计算比较简便。1927年福斯特首先将它用于水文现象,以后得到很多国家水文学者的广泛研究,也是我国水利水电工程水文计算规范中推荐采用的概率分布。