![水灾害防治中的多变量概率问题](https://wfqqreader-1252317822.image.myqcloud.com/cover/642/37204642/b_37204642.jpg)
第3章 多变量联合概率分布
3.1 两变量联合概率分布
3.1.1 两变量联合概率分布的性质
水科学领域经常遇到的多变量水文事件的描述以及不同水文事件遭遇和联合分布概率常常属于两变量联合概率分布问题,其联合概率分布函数表述如下:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_1.jpg?sign=1739258923-w91Qv9zRZYdldNSU7YTmkSXhWOwOLaYV-0-40f54c0a0e03d936c5e871760b55ccc1)
式(3.1)所示的联合概率分布函数F(x,y)表示事件X≤x和Y≤y同时发生的概率。
当随机变量相互独立时,联合概率分布函数直接为变量边际分布函数的乘积,即:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_2.jpg?sign=1739258923-dYccPlq1zo6BLNLVmwC7MGonP1F05LLL-0-1cf14b55e554c03864647af7bfac2bc2)
设点(x0,y0)的概率值F(x0,y0)表示二维随机变量(X,Y)落在点(x0,y0)左下方矩形区域,如图3.1所示。根据两变量联合概率分布函数F(x,y)的几何涵义,可以得到二维随机变量(X,Y)落在矩形区域(x1<X≤x2,y1<Y≤y2)内的概率如式(3.3)所示,图3.2为事件(x1<X≤x2,y1<Y≤y2)的几何表示。
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_3.jpg?sign=1739258923-pl3LN6UPhHuEBIrF9mFH5dZPiqJGEj6q-0-0581d3abb636f65a3f1629632982d755)
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_4.jpg?sign=1739258923-YrDFTcFS7D5yi0PGsfRKLNPKPlZNyvU0-0-7dc584220a3cdd2401fa93388f7a8d00)
图3.1 F(x0,y0)联合概率示意图
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_5.jpg?sign=1739258923-LeO1KRPhTjnegK9zy68Wbp9ujoyEZCfQ-0-66b4bd33686d033963630d3773b4f553)
图3.2 F(x1<X≤x2,y1<Y≤y2)联合概率示意图
3.1.1.1 两变量联合概率的基本性质
两变量联合概率分布函数F(x,y)具有以下基本性质[49-50]:
(1)对一切x,y,有0≤F(x,y)≤1:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_6.jpg?sign=1739258923-SRmRweSnDgez9m1tdZPcunu6MigDn49e-0-562e7a9fedb808de8b5160d41c26d466)
对于任意固定的x或y,有:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_7.jpg?sign=1739258923-WFtPS4S443FK3TDwq7st9mrTBwWm0IPF-0-2cd9991063040f9a7879bc8a7a09faba)
(2)F(x,y)是变量x或y的单调不减函数,即对任意固定的y,当x1<x2时,F(x1,y)≤F(x2,y);对任意固定的x,当y1<y2时,F(x,y1)≤F(x,y2)。
(3)F(x,y)分别是x或y的右连续函数,即对任意固定的y,F(x+0,y)=F(x,y);即对任意固定的x,F(x,y+0)=F(x,y)。
(4)对于任意的x1<x2,y1<y2,都有:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_8.jpg?sign=1739258923-YHn4cu3GJTGtRZmeU8GLKJ4ThIlvj41Y-0-f0b5ca64799e1edf43f933677a8e4455)
上述的基本性质(1)~(4)是F(x,y)为概率分布函数的充要条件,即:如果一个二维单值实函数F(x,y)满足上述基本性质,则该函数为某个二维随机变量(X,Y)的分布函数[49]。
3.1.1.2 二维连续型随机变量的性质
对于二维连续型随机变量(X,Y),联合密度函数概率为f(x,y),则对任意实数x或y,(X,Y)的概率分布函数为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_9.jpg?sign=1739258923-sofooHDxuTUWQe1Q84RL2o0DXTUa8uQM-0-675acc994b1b1793c0a94f5532599222)
根据二维概率分布的基本性质,二维连续型随机变量概率分布函数具有以下性质:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_10.jpg?sign=1739258923-b6LWmGU2Tof3pRj9ToGJp3wKujSLZymS-0-2c2106623c27b9c51151e973d4b16c00)
式(3.10)和式(3.11)为f(x,y)为二维连续型随机变量概率密度函数的充要条件,可以证明,满足式(3.10)和式(3.11)的任意一个二维函数f(x,y),必可作为某个二维随机变量(X,Y)的分布函数的概率密度函数。
若f(x,y)在点(x,y)处连续,则有:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_11.jpg?sign=1739258923-RzuQjFT3HODTYnSCHJIDK3ebeVh2TsQK-0-98694cd405d0ecfcef752b00fca84a3f)
设D为平面上的任意区域,则点(X,Y)落在D内的概率为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_12.jpg?sign=1739258923-B1olZhsvTnMO4IhTG9kZUylQe9957Rq3-0-073846b0568c349c256e774b92d6cf09)
即(X,Y)落在区域D内的概率等于以D为底,以曲面Z=f(x,y)为顶面的柱体体积。
3.1.1.3 边际分布及条件分布的性质
设n维随机变量(X1,X2,…,Xn)的联合分布函数为F(x1,x2,…,xn),边际分布函数记为FX1(x1),FX2(x2),…,FXn(xn)。则联合分布函数和边缘分布函数的关系为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_13.jpg?sign=1739258923-A25rPzVToMRXdRHnoFQuARgsbpI3knqY-0-73c2386e86c44225f20a387007af9cde)
当n=2时,对二维随机变量则有:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_14.jpg?sign=1739258923-25vN5hsfdnEYK7lhEVtdf7XhYpRx4Fl7-0-68706ce03a48ae2a0d1082feddfb8a91)
此外,当(X1,X2,…,Xn)的联合分布函数已知时,可以确定n维随机变量(X1,X2,…,Xn)的k维边缘分布函数,比如(X1,X2,…,Xn)关于(X1,X2),(X1,X2,X3)的边缘分布函数为:[49]
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_15.jpg?sign=1739258923-ff6Cc1bfQXsJDgY8tq8N2fLTOXuwMwai-0-97f6bb2f0f0f7c17a3805c5278dd1b2b)
对二维连续型随机变量(X,Y),其分布函数与边际概率具有以下关系:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_16.jpg?sign=1739258923-jNcRCUWGmMnOLDn29v6Ny9tJHes9jm9k-0-ad7a34dd1d07ab1207704e8009eb34a9)
3.1.2 各种形式的两变量联合概率及重现期
3.1.2.1 两变量联合概率及重现期
单变量的频率分析可以通过变量观测值的频次分析直观表示,其累积不超过概率(CDF)可以用式(3.18)表示:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_17.jpg?sign=1739258923-Gc1nj7opi2vmBQ7C5lofOTxqAiSJjeaZ-0-42fb9a7f8b4d1d75cfaffdfca7ad49a5)
水科学领域常常关心的是变量超过某一值的概率,则超过概率为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_18.jpg?sign=1739258923-xo1HjChn7VTLye130kc7tNJkWeZYR0N4-0-7f902d34df641a559104e75ae04541d2)
相应的重现期表示为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_19.jpg?sign=1739258923-QOa2KvtQoxXdYbCTfFfsn5Aq7emMaip7-0-ada832b35f1a22fe4260f8a997555d08)
两变量联合不超过概率:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_20.jpg?sign=1739258923-CM5LTwLINUMGMbEYGtWLvsdYspky029p-0-d236a07119d7ab7b23b0427636c5c876)
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_21.jpg?sign=1739258923-IgAM5henQxcSWsi9x9e0ISnMQdxxjvXY-0-5813ff66f08e2c992d74ded9895a8fbd)
图3.3 两变量概率分布函数——X,Y的区域
对于两变量频率分析而言,其变量频次分析没有单变量直观,但两变量联合分布函数的内涵与单变量分布函数本质上相同。由于变量之间的相关性和实际应用的需要,常常需要分析各种形式的联合概率和重现期。Sheng Yue[51]系统总结了各种变量的联合概率表示。图3.3为两变量联合分布函数变量的范围,可以直观描述两变量概率分布的内涵。将x—y平面分为四个象限,落入不同象限的(X,Y)表示不同的联合概率分布函数。
(1)象限Ⅰ:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_22.jpg?sign=1739258923-6qD0sKns55HkWyVQAoe61wmEAO1fZbYC-0-188056233e666288553f300494760d97)
(2)象限Ⅱ:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_23.jpg?sign=1739258923-FE8CYKmyRQLCEUbbb4EC59eVcSTbB3G7-0-b7afb93b17bd11ad0087a28e9347986f)
(3)象限Ⅲ:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_24.jpg?sign=1739258923-n6ysqczjLnR30UUnheGFZmY07bDEvGWe-0-ace84506bb7a4ac7400ec9265affc0eb)
(4)象限Ⅳ:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_25.jpg?sign=1739258923-LKZVxlTIEV1jNNiEy9rYAz38x7ypIj1c-0-864db9b2bddd907591b2328d552b9921)
根据概率分布内涵可知:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_26.jpg?sign=1739258923-Gy6OFnXfGgwVOQA5sPhkK1T0Yz4CqiOb-0-5c8cba6beda100603963ac1df5ac0f84)
相应事件(X>x,或Y>y,或X>x且Y>y)的重现期为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_27.jpg?sign=1739258923-Xn7Jqhb68XzwBppwwSQi7ly94RXEilhb-0-da6cf1d7240436be0e8e62d54f4157af)
相应(X>x,Y>y)的重现期为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_28.jpg?sign=1739258923-81YEljEBwPAJrpSN8YTMSSeZ6D37bHv4-0-023400ad3e1e1dfc928dda4565307e13)
当随机变量X,Y相互独立时,其联合分布函数及重现期表示如下:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_29.jpg?sign=1739258923-rWnSzgPyhKhHBAFkUmdThJdQDUCJRTAn-0-43a354753e57a741ca9f02d4a7b9fc92)
由于TX≥1及TY≥1,则不等式(3.31)成立:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_30.jpg?sign=1739258923-zxDClCyiabBBDUTWpJE9dP0y0sJCcgsA-0-f752521783ff74a9f1a03c4aed667c75)
3.1.2.2 两变量条件分布及重现期
水科学领域中各种随机变量之间常常具有不同程度的相关关系。条件分布是分析研究随机变量相依关系的有力工具。根据联合分布函数F(x,y)、边际分布函数FX(x)和FY(y),可以推导出各种条件分布函数。
设f(x,y)为连续随机变量X,Y的联合概率密度函数(PDF),fX(x)和fY(y)分别为变量X,Y的边际概率密度函数。
FX|Y(x|y)为在条件Y=y下,随机变量X的条件分布函数,fX|Y(x|y)≥0为条件概率密度函数,则:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_31.jpg?sign=1739258923-tEi0zgTrVJouDh4gPolmOVTSbxKBrVTS-0-757c1d41d9d1ae288f3a8a3b710ae587)
即
相应地,可得到给定Y=y条件下,X>x重现期为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_32.jpg?sign=1739258923-hXH5lnnOwhIDMmpBa7a0sqnfX3CljyKE-0-262bf5be69b15e0dfd1149ce963830c0)
同理可以得到在给定X=x条件下,Y的条件累积概率(CCDF)、条件概率密度函数及条件重现期如式(3.35)~式(3.37)所示:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_33.jpg?sign=1739258923-plnzS7WFo98bBqnUzizJnKQ6C0h4JVND-0-e256eeec0d0630d7e8ae0abff50db905)
当变量相互独立时,f(x,y)=f(x)f(y),则FX|Y(x|y)=FX(x),FY|X(y|x)=FY(y)。此时其相应的条件重现期为T(x|y)=TX,T(y|x)=TY。表明当变量相互独立时,条件概率和条件重现期等于相应边际分布的概率和重现期。
实际进行工程规划、设计、管理时,经常关心的是当某一变量大于或小于某一阈值条件下,另一变量的概率分布和重现期大小。
在X≤x条件下,变量Y≤y的概率为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_34.jpg?sign=1739258923-sv9Vco7q7MMkUUFHY1jnD0EGQVBnFxni-0-9c26f6f420e7ad785b10cc01c628b83b)
在X≤x条件下,变量Y>y的概率为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_35.jpg?sign=1739258923-pB04uuJ4c38pyhwlSpJ68aTMhR1bXyGW-0-2a72d1dea683796c1e973b719abb3d30)
给定X≥x时,变量Y的分布函数为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_36.jpg?sign=1739258923-BN5eN2aZFQxKjkieF57PNv3lwNZTIZa9-0-9304b9c8183f06875e06d4ff97432fff)
给定Y≥y时,变量X的分布函数为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_37.jpg?sign=1739258923-ymaYaUFkKQvbe8odmvezFgId5jhIxFcd-0-477f51008311e4a14e86fbff3ebf454a)
在Y≥y条件下,变量X≥x的概率为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_38.jpg?sign=1739258923-UTIFngsdrNTAx6CyueS2ZSguiEP9ctqt-0-2f54b04d138146de4b4ecfedd0bdf34f)
在X≥x条件下,变量Y≥y的概率为:
![img](https://epubservercos.yuewen.com/DF778D/19720710008529106/epubprivate/OEBPS/Images/txt003_39.jpg?sign=1739258923-yokSjVTo5AuC0FGcRTa4tPoY7al3S15v-0-fdf6fb4b779ed503f22f6bd820d1ae8d)