2.6.2 連續(xù)隨機(jī)變量函數(shù)的分布(Distribution of Functions of Random Variables)

參考:概率論與數(shù)理統(tǒng)計(jì)教程(第二版) 茆詩(shī)松

個(gè)人覺得“連續(xù)隨機(jī)變量函數(shù)的分布”這個(gè)表述有點(diǎn)繞裸弦,遠(yuǎn)不如英語(yǔ)的“Distribution of Functions of Random Variables”,所以加了個(gè)英文的標(biāo)題

幾個(gè)定理的證明的練習(xí)和筆記
先總結(jié)下思路脈絡(luò):

當(dāng)g(x)為嚴(yán)格單調(diào)時(shí)
  • 定理2.6.1是重點(diǎn)狰住,后面的定理2.6.2~定理2.6.4都是基于定理2.6.1推導(dǎo)

  • 定理2.6.1給出了y=g(x)嚴(yán)格單調(diào)時(shí),隨機(jī)變量Y=g(X)的通用PDF公式
  • 定理2.6.2給出了初始隨機(jī)變量為正態(tài)分布辅斟,F(xiàn)unctions of (正態(tài)分布隨機(jī)變量)的Function為一次線性函數(shù)時(shí)的分布
  • 定理2.6.3給出了初始隨機(jī)變量為正態(tài)分布转晰,F(xiàn)unctions of (正態(tài)分布隨機(jī)變量)的Function為對(duì)數(shù)函數(shù)時(shí)的分布,并引出了對(duì)數(shù)正態(tài)分布
  • 定理2.6.4給出了初始隨機(jī)變量為伽馬分布士飒,F(xiàn)unctions of (正態(tài)分布隨機(jī)變量)的Function為kX時(shí)的分布查邢,并引出了其應(yīng)用——任一伽馬分布轉(zhuǎn)化為卡方分布

  • 定理2.6.5給出了特定條件下F_{X}(X)服從(0,1)上的均勻分布U(0,1)這么一個(gè)結(jié)論,并推演出了各種分布隨機(jī)數(shù)的獲取方法——隨機(jī)模擬法(蒙特卡羅法)

當(dāng)g(x)為其他形式時(shí)
  • 給出了g(X)的通用求法思路酵幕,結(jié)合例題思考


下面進(jìn)入正題



當(dāng)g(x)為嚴(yán)格單調(diào)時(shí)


定理2.6.1

設(shè)X是連續(xù)隨機(jī)變量扰藕,其密度函數(shù)為p_{X}(x)Y=g(X)是另一個(gè)隨機(jī)變量芳撒。
y=g(x)嚴(yán)格單調(diào)邓深,其反函數(shù)h(y)有連續(xù)導(dǎo)函數(shù),則Y=g(X)密度函數(shù)
p_{Y}(y)=\left\{\begin{matrix} p_{X}[h(y)]|{h}'(y)|,a<y<b \\ 0,\quad \quad \quad \quad \quad \quad \quad others \end{matrix}\right.
其中a=min\left\{g(-\infty,g(\infty))\right \},b=max\left\{g(-\infty,g(\infty))\right \}

證明:

設(shè)g(x)是嚴(yán)格單調(diào)增函數(shù)笔刹,這時(shí)它的反函數(shù)h(y)也嚴(yán)格單調(diào)增函數(shù)芥备。
{h}'(y)>0。記a=g(-\infty),b=g(\infty)舌菜,這意味著y=g(x)僅在區(qū)間(a,b)取值萌壳,于是y的CDF F_{Y}(y)

  • 當(dāng)y<a時(shí)
    F_{Y}(y)=P(Y\leqslant y)=0
  • 當(dāng)y>b時(shí)
    F_{Y}(y)=P(Y\leqslant y)=1
  • 當(dāng)a\leqslant y\leqslant b時(shí)
    F_{Y}(y)=P(Y\leqslant y)=P(g(X)\leqslant y)
    \qquad \qquad \qquad \quad \quad =P(X\leqslant h(y))=\int_{-\infty}^{h(y)}p_{X}(x)dx
    由此得Y的PDF為
    p_{Y}(y)=\left\{\begin{matrix} p_{X}[h(y)]|{h}'(y)|,a<y<b \\ 0,\quad \quad \quad \quad \quad \quad \quad others \end{matrix}\right.
    同理可證當(dāng)g(x)是嚴(yán)格單調(diào)減函數(shù)時(shí),結(jié)論也成立。但要注意{h}'(y)<0袱瓮,故要加絕對(duì)值符號(hào)缤骨,這時(shí)a=g(\infty),b=g(-\infty),綜上所述尺借,定理得證绊起。

定理2.6.2

設(shè)隨機(jī)變量X服從正態(tài)分布N(\mu,\sigma ^2),則當(dāng)a\neq 0時(shí)燎斩,有Y=aX+b \sim N(a\mu+b,a^2\sigma ^2)

證明:

當(dāng)a>0時(shí)虱歪,Y=aX+b是嚴(yán)格增函數(shù),仍在(-\infty,\infty)上取值瘫里,其反函數(shù)h(y)X=(Y-b)/a实蔽,{h}'(y)=\frac{1}{a}定理2.6.1可得y的PDF
\quad p_{Y}(y)
=p_{X}[h(y)]|{h}'(y)|
=p_{X}(\frac{y-b}{a})\frac{1}{a}
=\frac{1}{\sqrt{2\pi}\sigma}exp\left \{ -\frac{1}{2\sigma^ 2}(\frac{y-b}{a}-\mu)^2\right \}\frac{1}{a}
=\frac{1}{\sqrt{2\pi}(a\sigma)}exp\left\{ -\frac{(y-a\mu-b)^2}{2a^2\sigma ^2}\right\}
這就是正態(tài)分布N(a\mu+b,a^2\sigma^2)的PDF。

當(dāng)a<0時(shí)谨读,Y=aX+b是嚴(yán)格減函數(shù),仍在(-\infty,\infty)上取值坛吁,其反函數(shù)h(y)X=(Y-b)/a劳殖,由定理2.6.1可得y的PDF
\quad p_{Y}(y)
=\frac{1}{\sqrt{2\pi}|a|\sigma}exp\left\{ -\frac{(y-a\mu-b)^2}{2 a^{2}\sigma^{2}} \right\}
這是正態(tài)分布N(a\mu+b,a^2\sigma^2)的PDF。
Q.E.D.


定理2.6.3(對(duì)數(shù)正態(tài)分布)

設(shè)隨機(jī)變量X\sim N(\mu,\sigma^2)拨脉,則Y=e^{X}的概率密度函數(shù)為
p_{Y}(y)=\left \{\begin{matrix} \frac{1}{\sqrt{2\pi}y\sigma}exp \left \{-\frac{(lny-\mu)^2}{2\sigma^2}\right \}, y>0 \\ 0 \qquad \qquad\qquad \qquad \quad, y\leqslant 0 \end{matrix} \right.

證明:

y=e^{x}是嚴(yán)格增函數(shù)哆姻,它僅在(0,\infty)上取值,其反函數(shù)h(y)x=\ln y玫膀,由定理2.6.1可得

  • 當(dāng)y \leqslant 0時(shí)矛缨,F_{Y}(y)=0,從而p_{Y}(y)=0.
  • 當(dāng)y>0時(shí)帖旨,Y的PDF為
    \quad p_{Y}(y)
    =p_{X}[h(y)]|{h}'(y)|
    =p_{X}(\ln y)|{\ln}'y|
    =\frac{1}{\sqrt{2\pi}\sigma}exp\left \{ -\frac{(\ln y-\mu)^2}{2\sigma^2}\right \}\frac{1}{y}
    =\frac{1}{\sqrt{2\pi}y\sigma }exp\left \{ -\frac{(\ln y-\mu)^2}{2\sigma^2}\right \}

對(duì)數(shù)正態(tài)分布

這個(gè)分布被稱為“對(duì)數(shù)正態(tài)分布”箕昭,記為LN(\mu,\sigma^2),其中\mu稱為對(duì)數(shù)均值解阅,\sigma^2稱為對(duì)數(shù)方差落竹。對(duì)數(shù)正態(tài)分布LN(\mu,\sigma^2)是一個(gè)偏態(tài)分布,也是一個(gè)常用分布货抄,實(shí)際中有不少隨機(jī)變量服從對(duì)數(shù)正態(tài)分布述召,譬如

  • 絕緣材料的壽命服從對(duì)數(shù)正態(tài)分布
  • 設(shè)備故障的維修時(shí)間服從對(duì)數(shù)正態(tài)分布
  • 家中僅有兩個(gè)小孩的年齡差服從對(duì)數(shù)正態(tài)分布
對(duì)數(shù)正態(tài)分布

定理2.6.4

設(shè)隨機(jī)變量X \sim Ga(\alpha,\lambda),則當(dāng)k>0時(shí)蟹地,有Y=kX\sim Ga(\alpha,\lambda/k)

證明:

x>0時(shí)积暖,Ga(\alpha,\lambda)PDF:
p(x)=\frac{(\lambda)^{\alpha}}{\Gamma (\alpha)}x^{\alpha-1}exp\left \{-\lambda x \right \}

因?yàn)?img class="math-inline" src="https://math.jianshu.com/math?formula=k%3E0" alt="k>0" mathimg="1">,所以y=kx是嚴(yán)格增函數(shù)怪与,它仍在(0,\infty)上取值夺刑,其反函數(shù)為x=y/k,由定理2.6.1可得

  • 當(dāng)y<0時(shí),p_{Y}(y)=0
  • 當(dāng)y\geqslant 0時(shí)性誉,
    \quad p_{Y}(y)=
    =p_{X}[h(y)]|{h}'(y)|
    =p_{X}(\frac{y}{k})\frac{1}{k}
    =\frac{\lambda^{\alpha}}{k\Gamma (\alpha)}(\frac{y}{k})^{\alpha-1}exp\left \{-\lambda \frac{y}{k} \right \}
    =\frac{(\lambda/k)^{\alpha}}{\Gamma (\alpha)}y^{\alpha-1}exp\left \{-\lambda \frac{y}{k} \right \}
    此即Ga(\alpha,\lambda/k)的PDF
    Q.E.D.

用途:

將任一伽馬分布轉(zhuǎn)化為\chi ^2分布窿吩,如
當(dāng)X \sim Ga(\alpha,\lambda),則2\lambda X\sim Ga(\alpha,\lambda/2\lambda)=Ga(\alpha,1/2)=\chi^2(2\alpha)


定理2.6.5

若隨機(jī)變量的分布函數(shù)F_{X}(x)為嚴(yán)格單調(diào)遞增的連續(xù)函數(shù),其反函數(shù)F_{X}^{-1}(y)存在狞洋,則Y=F_{X}(X)服從(0,1)上的均勻分布U(0,1)

證明:

下求Y=F_{X}(X)的分布函數(shù)锅论。由于分布函數(shù)F_{X}(X)僅在[0,1]區(qū)間上取值,故

  • 當(dāng)y<0時(shí)轧邪,因?yàn)?img class="math-inline" src="https://math.jianshu.com/math?formula=%5C%7BF_%7BX%7D(X)%5Cleqslant%20y%5C%7D" alt="\{F_{X}(X)\leqslant y\}" mathimg="1">是不可能事件,所以
    F_{Y}(y)=P(Y\leqslant y)=P(F_{X}(X)\leqslant y)=0
  • 當(dāng)0\leqslant y \leqslant 1時(shí)羞海,有
    F_{Y}(y)=P(Y\leqslant y)=P(F_{X}(X)\leqslant y)
    \qquad \qquad \qquad \qquad =P(X\leqslant F_{X}^{-1}(y))=F_{X}(F_{X}^{-1}(y))=y
  • 當(dāng)y\geqslant 1時(shí)忌愚,因?yàn)?img class="math-inline" src="https://math.jianshu.com/math?formula=%5C%7BF_%7BX%7D(X)%5Cleqslant%20y%5C%7D" alt="\{F_{X}(X)\leqslant y\}" mathimg="1">是必然事件,所以
    F_{Y}(y)=P(Y\leqslant y)=P(F_{X}(X)\leqslant y)=1

綜上所述却邓,Y=F_{X}(X)的分布函數(shù)為
F_{Y}(y)=\left\{ \begin{matrix} 0,y<0 \qquad\\ y,0\leqslant y <1\\ 1,y\geqslant 1 \qquad\\ \end{matrix} \right.
這正是(0,1)上均勻分布的CDF硕糊,所以Y\sim U(0,1)

意義:

任一個(gè)隨機(jī)變量X都可以通過(guò)其分布函數(shù)F(X)與均勻分布隨機(jī)變量U發(fā)生關(guān)系。譬如

  • X服從指數(shù)分布Exp(\lambda)腊徙,
    其分布函數(shù)為F(x)=1-e^{-\lambda x}简十,
    當(dāng)x換為X后,有
    U=1-e^{-\lambda X} \quadX=\frac{1}{\lambda}\ln \frac{1}{1-U}
    后一式表明:由均勻分布U(0,1)的隨機(jī)數(shù)(偽觀察值)u_{i}可得指數(shù)分布Exp(\lambda)的隨機(jī)數(shù)x_{i}=\frac{1}{\lambda}\ln \frac{1}{1-u_{i}},i=1,2,…,n,…撬腾。
    而均勻分布隨機(jī)數(shù)在任一個(gè)統(tǒng)計(jì)軟件都可產(chǎn)生螟蝙,從而指數(shù)分布(繼而其他分布)隨機(jī)數(shù)也可以獲得。而各種分布隨機(jī)數(shù)的獲得是進(jìn)行隨機(jī)模擬法(又稱蒙特卡羅法)的基礎(chǔ)


當(dāng)g(x)為其他形式時(shí)

可直接由Y的分布函數(shù)F_{Y}(y)=P(g(X)\leqslant y)出發(fā)民傻,按函數(shù)g(x)的特點(diǎn)做個(gè)案處理

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末胰默,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子漓踢,更是在濱河造成了極大的恐慌牵署,老刑警劉巖,帶你破解...
    沈念sama閱讀 221,888評(píng)論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件彭雾,死亡現(xiàn)場(chǎng)離奇詭異碟刺,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī)薯酝,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,677評(píng)論 3 399
  • 文/潘曉璐 我一進(jìn)店門半沽,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人吴菠,你說(shuō)我怎么就攤上這事者填。” “怎么了做葵?”我有些...
    開封第一講書人閱讀 168,386評(píng)論 0 360
  • 文/不壞的土叔 我叫張陵占哟,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我,道長(zhǎng)榨乎,這世上最難降的妖魔是什么怎燥? 我笑而不...
    開封第一講書人閱讀 59,726評(píng)論 1 297
  • 正文 為了忘掉前任,我火速辦了婚禮蜜暑,結(jié)果婚禮上铐姚,老公的妹妹穿的比我還像新娘。我一直安慰自己肛捍,他們只是感情好隐绵,可當(dāng)我...
    茶點(diǎn)故事閱讀 68,729評(píng)論 6 397
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著拙毫,像睡著了一般依许。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上缀蹄,一...
    開封第一講書人閱讀 52,337評(píng)論 1 310
  • 那天峭跳,我揣著相機(jī)與錄音,去河邊找鬼缺前。 笑死坦康,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的诡延。 我是一名探鬼主播,決...
    沈念sama閱讀 40,902評(píng)論 3 421
  • 文/蒼蘭香墨 我猛地睜開眼古胆,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼肆良!你這毒婦竟也來(lái)了?” 一聲冷哼從身側(cè)響起逸绎,我...
    開封第一講書人閱讀 39,807評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤惹恃,失蹤者是張志新(化名)和其女友劉穎,沒(méi)想到半個(gè)月后棺牧,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體巫糙,經(jīng)...
    沈念sama閱讀 46,349評(píng)論 1 318
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,439評(píng)論 3 340
  • 正文 我和宋清朗相戀三年颊乘,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了参淹。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 40,567評(píng)論 1 352
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡乏悄,死狀恐怖浙值,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情檩小,我是刑警寧澤开呐,帶...
    沈念sama閱讀 36,242評(píng)論 5 350
  • 正文 年R本政府宣布,位于F島的核電站,受9級(jí)特大地震影響筐付,放射性物質(zhì)發(fā)生泄漏卵惦。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,933評(píng)論 3 334
  • 文/蒙蒙 一瓦戚、第九天 我趴在偏房一處隱蔽的房頂上張望沮尿。 院中可真熱鬧,春花似錦伤极、人聲如沸蛹找。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,420評(píng)論 0 24
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)庸疾。三九已至,卻和暖如春当编,著一層夾襖步出監(jiān)牢的瞬間届慈,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,531評(píng)論 1 272
  • 我被黑心中介騙來(lái)泰國(guó)打工忿偷, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留金顿,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 48,995評(píng)論 3 377
  • 正文 我出身青樓鲤桥,卻偏偏與公主長(zhǎng)得像揍拆,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子茶凳,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,585評(píng)論 2 359