平穩(wěn)過(guò)程講學(xué)PPT課件

1、 此外當(dāng)我們知道一個(gè)隨機(jī)過(guò)程是平穩(wěn)過(guò)程時(shí)，它應(yīng)不隨時(shí)間的推移而變幻無(wú)常。例如當(dāng)我們要測(cè)定一個(gè)電阻的熱噪聲的統(tǒng)計(jì)特性，由于它是平穩(wěn)過(guò)程，因而我們?cè)谌魏螘r(shí)間進(jìn)行測(cè)試都能得到相同的結(jié)果。第1頁(yè)/共76頁(yè)4.1 4.1 定義和例子定義和例子 定義嚴(yán)平穩(wěn)隨機(jī)過(guò)程：對(duì)于任意的t，隨機(jī)過(guò)程X(t)的任意n維概密度都有 則稱X(t)為嚴(yán)平穩(wěn)隨機(jī)過(guò)程。 研究平穩(wěn)過(guò)程的意義在于：該過(guò)程在任何時(shí)刻計(jì)算它的統(tǒng)計(jì)結(jié)果都是相同的。由定義知平穩(wěn)隨機(jī)過(guò)程的n維概度密度函數(shù)不隨時(shí)間而變化，這一特性具體反映在隨機(jī)過(guò)程的一、二維概率密度及數(shù)字特征方面具有如下性質(zhì)：12121212( ,; , ,)( ,)XnnXnnPx xx t

2、 ttPx xx ttt第2頁(yè)/共76頁(yè) 性質(zhì)4.1 若X(t)為平衡過(guò)程，則它的一維概率密度與時(shí)間無(wú)關(guān) 證 設(shè)X(t)的一維概率密度函數(shù)為 ，由于X(t)為平穩(wěn)過(guò)程 令 則 由此我們可求平穩(wěn)過(guò)程X(t)的均值、均方值、方差。11( ; )XPx t1111( ; )( ;)XXPx tPx t1t 1111( ; )( ;0)()XXXPx tPxPx222111( )()XXE XtX Px dx 第3頁(yè)/共76頁(yè)2221112222( )( )()()( )( )XXXXXxD X tE X tMxMPxdxE XtEX tM 顯然，X(t)的均方值、方差都與時(shí)間t無(wú)關(guān) 。由此知，當(dāng)隨機(jī)

3、過(guò)程為平穩(wěn)過(guò)程時(shí)，該過(guò)程的所有樣本函數(shù)總是它們均值水平直線上下波動(dòng)，樣本曲線偏離水平直線的幅度正好是( )XD X x。第4頁(yè)/共76頁(yè)如圖4.1所示，圖中細(xì)實(shí)線表示隨機(jī)過(guò)程的樣本函數(shù)，粗實(shí)線表示隨機(jī)過(guò)程的數(shù)學(xué)期望，虛線表示隨機(jī)過(guò)程對(duì)數(shù)學(xué)期望的偏差。第5頁(yè)/共76頁(yè) 性質(zhì)4.2 平穩(wěn)過(guò)程X(t)的二維概率密度只與 的時(shí)間間隔有關(guān)，而與時(shí)間起點(diǎn)無(wú)關(guān)。證：設(shè)X(t)的二維概率密度函數(shù)為由于X(t)為平穩(wěn)過(guò)程，所以對(duì)任意 有若令 ，則而 正是隨機(jī)過(guò)程二維概率密度函數(shù)的時(shí)間間隔，令 ，則：12,t t1212( ,; , )XPx xt t12121212( ,; , )( ,;,)XXPx x t

4、tPx x tt1t 12121212( ,; , )( ,;0,)( ,;)XXXPx x t tPx xPx x12121221( ,; , )( ,;0,)XXPx x t tPx xtt21tt21tt第6頁(yè)/共76頁(yè) 此式表明，平穩(wěn)隨機(jī)過(guò)程的二維概率密度函數(shù)僅依賴于 ，而時(shí)間的個(gè)別值 無(wú)關(guān)。由此，我們可以進(jìn)一步來(lái)討論平穩(wěn)過(guò)程X(t)的協(xié)方差函數(shù)應(yīng)具有什么樣的表達(dá)形式。12,t t121212121212121212( , )( )( )(; , )(;)( )XXXXRt tE X t X tx x Px x t t dx dxx x Px xdx dxR 1212122( , )(

5、 , )( )( )( )( )( )XXXXXXXXXXCt tRt tMt MtRMMRMC又 2( )( )XXXCRM 第7頁(yè)/共76頁(yè) 順便指出，由一個(gè)隨機(jī)過(guò)程的平穩(wěn)性研究可推廣到關(guān)于兩個(gè)隨機(jī)過(guò)程的平穩(wěn)性研究，可以這樣說(shuō)，若兩個(gè)隨機(jī)過(guò)程的聯(lián)合概率密度函數(shù)不隨時(shí)間的平移而變化，與時(shí)間的起點(diǎn)無(wú)關(guān)，則可稱這兩個(gè)隨機(jī)過(guò)程是聯(lián)合平衡的，或稱平穩(wěn)相依。第8頁(yè)/共76頁(yè) 從上面介紹的嚴(yán)平穩(wěn)隨機(jī)過(guò)程的定義知，要判斷一個(gè)隨機(jī)過(guò)程是否是嚴(yán)平穩(wěn)，需要確定該隨機(jī)過(guò)程的任意n維概率密度函數(shù)族，它的變化是否與時(shí)間的平穩(wěn)無(wú)關(guān)，這本身就是一個(gè)十分困難的工作，然而在工程上根據(jù)實(shí)際需要，我們往往只在所謂的相關(guān)理論范圍內(nèi)

6、考慮隨機(jī)過(guò)程的平穩(wěn)性問(wèn)題，這里所指的相關(guān)理論，就是指隨機(jī)過(guò)程的數(shù)字特征，即數(shù)學(xué)期望、相關(guān)函數(shù)和今后要介紹的功率普密度等。當(dāng)在相關(guān)理論又可指研究隨機(jī)過(guò)程的一、二階矩理論。 第9頁(yè)/共76頁(yè) 前面已經(jīng)介紹過(guò)，對(duì)于一個(gè)隨機(jī)過(guò)程X(t)，我們當(dāng)然希望能建立起它的多維分布函數(shù)，因?yàn)殡S機(jī)過(guò)程的多維分函數(shù)能較完整地描述隨機(jī)過(guò)程的統(tǒng)計(jì)特性，但是要建立多維分布函數(shù)往往很困難，因此我們一般在相關(guān)理論范圍內(nèi)也就是用數(shù)字特征來(lái)描述過(guò)程的重要特性，這種用數(shù)字特征來(lái)描述過(guò)程X(t)統(tǒng)計(jì)特性變化規(guī)律，對(duì)很多實(shí)際問(wèn)題往往已能獲得很好的效果，可以提取到所需的參數(shù)。 第10頁(yè)/共76頁(yè) 定義寬平穩(wěn)過(guò)程：給定隨機(jī)過(guò)程X(t)，如果

7、 常數(shù) 且則稱X(t)為寬平穩(wěn)過(guò)程（廣義平穩(wěn)過(guò)程）。顯然由寬平穩(wěn)定義可知，要求就要考慮X(t)的一維概率密度函數(shù) 和二維概率密度函數(shù) 。( )XE X tM21212( ),( , )( )( )( )XXE XtRt tE X t X tR 21tt12 ( ),( , )XE E tRt t11( , )XPx t1212( ,; , )XPx xt t第11頁(yè)/共76頁(yè) 下面我們來(lái)分析一下嚴(yán)平穩(wěn)和寬平穩(wěn)之間的關(guān)系。對(duì)于一個(gè)隨機(jī)過(guò)程X(t)，如果它是嚴(yán)平穩(wěn)的，且它的二階矩存在及均方有界 ，則由嚴(yán)平穩(wěn) 雙因嚴(yán)平穩(wěn)的一維概率密度與時(shí)間無(wú)關(guān)，即 常數(shù) 又因嚴(yán)平穩(wěn)的二維概率密度只與時(shí)間間隔有關(guān)，即

8、2( )E Xt 111212( ,; ,)( ,;,)XnnXPxxttPxxtt111( ; )( )XXPx tPx( )XE X tM121212( ,; , )( ,; )XXPx x t tPx x第12頁(yè)/共76頁(yè)21tt12( , )( )XXRt tR 222111( )()XXE Xtx Pxdx 綜上所述，嚴(yán)平穩(wěn)一定是寬平穩(wěn) 反之不一定成立，除非是高斯過(guò)程（正態(tài)過(guò)程）。類似地，我們還可以給出兩個(gè)隨機(jī)過(guò)程聯(lián)合寬平穩(wěn)定義。定義聯(lián)合寬平穩(wěn)：對(duì)于平穩(wěn)過(guò)程 若1( ), ( )X tY t1221( , )( ),XYXYRt tRtt 則稱( ), ( )X t Y t聯(lián)合寬平穩(wěn)

9、。 第13頁(yè)/共76頁(yè) 順便指出，今后凡提到“平穩(wěn)過(guò)程”，通常是指寬平穩(wěn)過(guò)程。 例4.1 設(shè)Y是隨機(jī)變量，試分別考慮隨機(jī)過(guò)程 的平穩(wěn)性。 解 Y是隨機(jī)變量， 這一過(guò)程是一個(gè)與時(shí)間無(wú)關(guān)的特殊的過(guò)程，它的任何n維概率密度函數(shù) 與時(shí)間無(wú)關(guān)，所以是一個(gè)嚴(yán)平穩(wěn)。 是嚴(yán)平穩(wěn) , 只要 則X1(t)是寬平穩(wěn)。對(duì)于 12( ),( )X tY XttY1( )X tY1(,)YnP yy1( )X tY221( )E XtE Y 2( ),XttY2( ) E E tE tYtE Y第14頁(yè)/共76頁(yè) 都與時(shí)間 有關(guān)，所以 為非平穩(wěn)。 例4.2 設(shè) 是一周期為T(mén)的函數(shù)， 是（0,T）上具有均勻分布的隨機(jī)變量，

10、稱為 隨機(jī)相位周期過(guò)程，試討論它的平穩(wěn)性。解 由題設(shè)知 的概率密度函數(shù)為22122122121 2( , )( )( )XRt tE Xt XtE tYt Yt t E Y12,t t2( )XttY( )S t( )()X tS t10( )0Tf tT其它第15頁(yè)/共76頁(yè) 要討論X(t)的平穩(wěn)性，由寬平穩(wěn)定義知，需要求 。 當(dāng)取定 為一隨機(jī)變量 的函數(shù) ，由求隨機(jī)變量函數(shù)的數(shù)學(xué)期望公式知 令 ，則12( ),( , )XE X tRt t,( )()tX tX t時(shí)()Yg X ( ) ( )E Yg x f x dx001( )()( )()TTE X ttfdS tdTt 011(

11、)( )( )t TTtE X tSdSdTT常數(shù) 第16頁(yè)/共76頁(yè) 又 令120( , )( ,)( )() ()()() () ( )XXTRt tRt tE X t X tE S tS tS tS tfd,t 01( )( ) ()1( ) ()( )ttTXE X tSSdTSSdRT第17頁(yè)/共76頁(yè)4.2 4.2 遍歷性定理遍歷性定理 1. 各態(tài)歷經(jīng)問(wèn)題的提出 對(duì)于一個(gè)隨機(jī)過(guò)程X(t)，我們當(dāng)然希望知道它們的分布函數(shù)，但很困難，于是我們退而求其次，考慮求它的數(shù)字特征即數(shù)學(xué)期望、相關(guān)函數(shù)等。但要求X(t)的數(shù)字特征，首先需要知道它的一、二維概率密度函數(shù)，即 這實(shí)際上又很難辦，進(jìn)而為

12、我們求數(shù)字特征又帶來(lái)困難。怎么解決這個(gè)問(wèn)題呢？實(shí)際上，在工程中，要求X(t)的數(shù)字特征，我們自先是通過(guò)試驗(yàn)來(lái)產(chǎn)生一族時(shí)間樣本函數(shù) 111212( ; ),( ,; , )XXPx tPx x t t1( ),( ),nx tx t第18頁(yè)/共76頁(yè)X(t)或者是做試驗(yàn)產(chǎn)生一個(gè)樣本函數(shù)x(t)，然后再對(duì)樣本函數(shù)x(t)取不同時(shí)刻，如 ，得所對(duì)應(yīng)的結(jié)果 ，即此時(shí)隨機(jī)過(guò)程可表示為 。 對(duì)任意指定時(shí)刻 的數(shù)學(xué)期望可近似表示為 協(xié)方差函數(shù)可近似表示為 來(lái)計(jì)算，顯然這種用近似計(jì)算的方法來(lái)估計(jì)隨機(jī)過(guò)程的數(shù)學(xué)期望及協(xié)方差函數(shù)要求n很大，即樣本函數(shù)xk(t)很多。但這在實(shí)際工程又常常又很難做到，于是人們自然想到

13、能不能夠通過(guò)測(cè)試一個(gè)樣本函數(shù)如 01,nt tt01( ), ( ),( )nnx tx tx t0( ) ( ), ( ),nX tx tx t1,( )t X t1111( )( )nkhE X tx tn121211( , )( )( )nXkkkRt tx t x tn( ),1,2,ix t i 第19頁(yè)/共76頁(yè) 用一個(gè)樣本函數(shù)xi(t)的均值和相關(guān)函數(shù)來(lái)近似隨機(jī)過(guò)程的均值和相關(guān)函數(shù)，如果能，這為我們求隨機(jī)過(guò)程的數(shù)學(xué)特征就帶來(lái)了很大方便。 這里提出一個(gè)問(wèn)題：怎樣表示一個(gè)樣本函數(shù)如x1(t)的均值呢？我們以下式來(lái)表示 顯然x1(t)不同其積分結(jié)果一般不同。 于是對(duì)一個(gè)隨機(jī)過(guò)程， ，其

14、樣本函數(shù)的積數(shù)結(jié)果可能不同。此時(shí)顯然用一個(gè)樣本函數(shù)的數(shù)字特征如 ，近似 是不正確的。但是如果當(dāng)時(shí)間區(qū)間T充分大時(shí)，如果X(t)的絕大多數(shù)樣本函數(shù)的均值 111( )2TxTMx t drT1( ) ( ),( ),nX tx tx t1xM( )E X t111( )lim( )2TTTx tx t dtT第20頁(yè)/共76頁(yè)都有則我們可用其中一個(gè)樣本函數(shù)的均值 作為 X(t)的近似，即 定義隨機(jī)過(guò)程的時(shí)間均值和時(shí)間相關(guān)函數(shù)：稱為隨機(jī)過(guò)程的時(shí)間相關(guān)函數(shù)。221( )lim( )2TTTx tx t dtT1( )lim( )2TnnTTx tx t dtT12( )( )( )nx tx tx

15、t( )nx t( )( ),1,2nx tE X tn1( )()lim( )()2TTTX t X tX t X tdtT第21頁(yè)/共76頁(yè) 注意：定義中 一般都是隨機(jī)變量（常數(shù)可看作特殊的隨機(jī)變量）。 由上述分析可知，是不是任何一個(gè)隨機(jī)過(guò)程 ，它的數(shù)學(xué)期望、相關(guān)函數(shù)都可用其中的一個(gè)樣本函數(shù)的均值和協(xié)方差函數(shù)來(lái)近似呢，顯然不一定，一個(gè)自然的問(wèn)題是X(t)在什么條件下可用一個(gè)樣本函數(shù)的均值和協(xié)方差函數(shù)作為整個(gè)過(guò)程X(t)的均值，協(xié)方差函數(shù)的近似呢？( ),( )()X tX t X t1( ) ( ),( )nX tx tx t第22頁(yè)/共76頁(yè) 2. 平均隨機(jī)過(guò)程的各態(tài)歷經(jīng)性 要回答上述的

16、問(wèn)題，我們?cè)O(shè)當(dāng)X(t)為平穩(wěn)過(guò)程且滿足一定條件時(shí)，可用一個(gè)樣本函數(shù)的均值和協(xié)方差函數(shù)作為過(guò)程X(t)的數(shù)字特征近似，為此我們給出如下定義：定義：設(shè)X(t)是一個(gè)平穩(wěn)過(guò)程 （1）若 以概率1成立，則稱隨機(jī)過(guò)程X(t)均值具有各態(tài)歷經(jīng)性這里依概率1成立是指對(duì)X(t)的所有樣本函數(shù)即( )( )XX tE X tM1( )( ),( )( )Xnxx tE X tMx tE X tM第23頁(yè)/共76頁(yè) 由此知，此時(shí)，我們可用一個(gè)樣本函數(shù)的均值如 的值作為 的近似值。反之，若已知X(t)的均值各態(tài)歷程，則可用一個(gè)樣本函數(shù)的均值作為過(guò)程X(t)的均值。（2）若 以概率1成立，則稱X(t)的協(xié)方差函數(shù)具有

17、各態(tài)歷經(jīng)性。11( )lim( )2TTtx tx t dtT( )E X t1( )lim( ),1,2,2TnTTE X tx t dt nT( )()( )()( )XX t X tE X t X tR第24頁(yè)/共76頁(yè) 這里若X(t)的協(xié)方差函數(shù)各態(tài)歷經(jīng)，就是指我們可用過(guò)程X(t)的一個(gè)樣本函數(shù)、xn(t)的時(shí)間相關(guān)函數(shù) 即 作為過(guò)程的相關(guān)函數(shù)。 （3）若X(t)的均值和協(xié)方差函數(shù)都具有各態(tài)歷經(jīng)性，則稱X(t)是寬各態(tài)歷經(jīng)過(guò)程，簡(jiǎn)稱X(t)為各態(tài)歷經(jīng)過(guò)程。 綜上所述，如果X(t)是各態(tài)歷經(jīng)過(guò)程，則必為平穩(wěn)過(guò)程，此時(shí)可用過(guò)程的一個(gè)樣本函數(shù)的數(shù)字特征作為過(guò)程的數(shù)字特征近似。( )()( )

18、()( )nnXx t x tE X t X tR1( )lim( )()2TXnnTTRx t x tdtT第25頁(yè)/共76頁(yè) 例4.3 設(shè)隨機(jī)過(guò)程 式中 為參數(shù)，是（0.2, ）上均勻分布隨機(jī)變量。 求證X(t)是寬平穩(wěn)過(guò)程； 該過(guò)程是否是各態(tài)歷經(jīng)過(guò)程。解 0( )cos()X tAt0,A2001( )cos()02E X tAtd第26頁(yè)/共76頁(yè)000200022000020( )()cos()cos()coscos(22 )21coscos(22 )22cos( )2XE X t X tE AtAtAEtAdAR X(t)為一寬平穩(wěn)過(guò)程。 01( )limcos()2TTTX tA

19、tdtT00cossinlim0TATT( ,)XRt t第27頁(yè)/共76頁(yè)顯然由、結(jié)果再由隨機(jī)過(guò)程各態(tài)歷經(jīng)定義知 X(t)為寬各態(tài)歷經(jīng)過(guò)程。00020cos() cos()( )()limcos2TTTAtdrX t X tA ( )( )0X tE X t0( )()( )()( )2XX t X tE X t X tRA 第28頁(yè)/共76頁(yè) 如果兩個(gè)隨機(jī)過(guò)程X(t)，Y(t)，當(dāng)它們各自都是各態(tài)歷經(jīng)時(shí)，并且時(shí)間互相關(guān)函數(shù)與統(tǒng)計(jì)相關(guān)函數(shù)以概率1相等時(shí)，我們有如下定義： 定義兩個(gè)隨機(jī)過(guò)程聯(lián)合各態(tài)歷經(jīng)： 設(shè)X(t)，Y(t)各自都各態(tài)歷經(jīng) 則稱X(t)，Y(t)為聯(lián)合各態(tài)歷經(jīng)過(guò)程。 同理當(dāng)X(

20、t)，Y(t)聯(lián)合各態(tài)歷經(jīng)時(shí)，可用它們的一對(duì)樣本函數(shù)的數(shù)字特征作為X(t)，Y(t)的數(shù)字特征近似。( ) ()( ) ()( )XYX t Y tE X t Y tR第29頁(yè)/共76頁(yè) 3. 隨機(jī)過(guò)程成為各態(tài)歷經(jīng)過(guò)程的判定 從前面的分析知，如果一個(gè)隨機(jī)過(guò)程能成為一個(gè)平衡過(guò)程，這對(duì)我們研究各態(tài)歷經(jīng)，則該過(guò)程一定是平衡過(guò)程，反之則一定成立，于是很自然提出這樣一個(gè)問(wèn)題，能不能給出一些判定定理，使其可以很方便地判定一個(gè)平穩(wěn)過(guò)程成為各態(tài)歷經(jīng)過(guò)程。通過(guò)對(duì)平穩(wěn)過(guò)程的分析研究，我們給出如下的幾個(gè)判定定理。 性質(zhì)4.3 平穩(wěn)過(guò)程X(t)的均值具有各態(tài)歷經(jīng)性的充要條件是 式中： 為平穩(wěn)過(guò)程的協(xié)方差函數(shù)； 為平穩(wěn)

21、過(guò)程的數(shù)學(xué)期望。2201lim1( )02TXXTRMdTT( )XRXM第30頁(yè)/共76頁(yè) 例4.4 已知隨機(jī)電報(bào)信號(hào)X(t)，它的 ， ，問(wèn)X(t)是否均值各態(tài)歷經(jīng)。解 ( )0E X t( )XRe 2201lim1( )2TXXTRMdrTT2| |02| |0222201lim1021lim12111lim02TTTTTTTedrTTedTTeTTT X(t)是均值各態(tài)歷經(jīng)的。第31頁(yè)/共76頁(yè) 性質(zhì)4.4 平衡過(guò)程X(t)的協(xié)方差函數(shù)具備各態(tài)歷經(jīng)性的充要條件是 式中 平穩(wěn)過(guò)程X(t)和Y(t)的互相關(guān)函數(shù)具有聯(lián)合各態(tài)歷經(jīng)性的充要條件（4.7）式相似，只是將（4.7）式中相應(yīng)的協(xié)方差

22、函數(shù)改為互相關(guān)函數(shù)即可。2211101lim1 ( )( )02XTTBRdT111( )()()()( )BE X tX tX tX t第32頁(yè)/共76頁(yè) 性質(zhì)4.5 對(duì)于高斯平穩(wěn)過(guò)程，如果它的均值為零，協(xié)方差函數(shù)連續(xù)，則該過(guò)程各態(tài)歷經(jīng)的一個(gè)充分條件是 綜上所述，對(duì)一個(gè)平穩(wěn)隨機(jī)過(guò)程X(t)通過(guò)性質(zhì)1、2判定以后，如果X(t) 各態(tài)歷經(jīng)了，則對(duì)于該過(guò)程的數(shù)字特征，即求 ，我們可用0|( )|XRdr ( ),( )XE X tR11( )( )lim( )2TTTX tE X tx t dtT111( )()( )lim()( )2TXTTRX tX tx tx t dtT第33頁(yè)/共76頁(yè)

23、也就是當(dāng)X(t)各態(tài)歷經(jīng)時(shí)，我們可用一個(gè)樣本函數(shù)的時(shí)間均值和時(shí)間協(xié)方差函數(shù)作為過(guò)程X(t)的數(shù)學(xué)期望、協(xié)方差函數(shù)的近似。 最后順便說(shuō)明，對(duì)于許多實(shí)際問(wèn)題，如果要從理論上判定一個(gè)過(guò)程是否為各態(tài)歷經(jīng)過(guò)程，往往是比較困難。因此工程上經(jīng)常都是憑經(jīng)驗(yàn)把各態(tài)歷經(jīng)性作為一種假設(shè)，在后根據(jù)實(shí)驗(yàn)來(lái)檢驗(yàn)這個(gè)假設(shè)是否合理。 在實(shí)際應(yīng)用一般不可能給出隨機(jī)過(guò)程X(t)的樣本函數(shù)x(t)的表達(dá)式，因此，確定各態(tài)歷經(jīng)過(guò)程的數(shù)學(xué)期望、協(xié)方差函數(shù)，有兩種方法：第34頁(yè)/共76頁(yè) 第一種方法用模擬協(xié)方差分析儀，自動(dòng)畫(huà)出協(xié)方差曲線。 這種儀器的功能是當(dāng)輸入樣本函數(shù)時(shí)，X-Y記錄儀自動(dòng)描繪出協(xié)方差函數(shù)的曲線。它的工作原理如圖4.2所

24、示。圖4.2第35頁(yè)/共76頁(yè) 第二種方法用數(shù)字處理方法（即近似計(jì)算方法）。 如圖4.2把0，T等分為N個(gè)長(zhǎng)為 的小區(qū)間，再在時(shí)刻 ， 取樣，得N個(gè)函數(shù)值 。于是再把積分表過(guò)式表示為基本區(qū)間上的和，就有數(shù)字估計(jì)式（4.8）。類似可以寫(xiě)出在 時(shí)的協(xié)方差函數(shù)估計(jì)式（4.9）式，由這個(gè)估計(jì)式可算出協(xié)方差函數(shù)的一系列近似值，從而可作出協(xié)方差函數(shù)的近似圖形，見(jiàn)圖4.3。 TtN 1()2ktkt1,2,( )kNx t對(duì)( ),1,2,kkxx tkNrr t 第36頁(yè)/共76頁(yè) 最后指出，工程上遇到的很多平穩(wěn)過(guò)程，我們一般都把它看成各態(tài)歷經(jīng)過(guò)程，然后用各態(tài)歷經(jīng)的方法來(lái)確定過(guò)程X(t)的統(tǒng)計(jì)特性，看處理

25、出來(lái)的結(jié)果是否與實(shí)際相符合，如果不相符合，再對(duì)過(guò)程的假設(shè)作修改。 111( )( )lim( )2NTkTTkE X tx tX t dtxTN1( )()( )lim( )()2TXTTRX tX tX t X tdrT11(),0,1,2,()N rXrkk rkRx xNrrm mn第37頁(yè)/共76頁(yè) 1. 平穩(wěn)過(guò)程協(xié)方差函數(shù)的性質(zhì) 對(duì)于一個(gè)隨機(jī)過(guò)程，它的基本數(shù)字特征是數(shù)學(xué)期望和相關(guān)函數(shù)，但是當(dāng)隨機(jī)過(guò)程為平穩(wěn)過(guò)程時(shí)，它的數(shù)學(xué)期望是一個(gè)常數(shù)，經(jīng)中心化后可以變?yōu)榱?，所以?dāng)過(guò)程X(t)平穩(wěn)后其基本數(shù)字特征實(shí)際上就是相關(guān)函數(shù)。此外，相關(guān)函數(shù)不僅可向我們提供隨機(jī)過(guò)程各狀態(tài)間的關(guān)聯(lián)特性的信息，而且也

26、是求取隨機(jī)過(guò)程的功率譜密度以及從噪聲中提取有用信息工具。為此下面我們專門(mén)研究一下平穩(wěn)過(guò)程相關(guān)函數(shù)的性質(zhì)。第38頁(yè)/共76頁(yè) 性質(zhì) 4.6 證 當(dāng) ， 即平穩(wěn)過(guò)程的均方值可自由相關(guān)函數(shù) 得到。 性質(zhì) 4.7 ，即平穩(wěn)過(guò)程的協(xié)方差函數(shù)為偶函數(shù)。 同理 證 2(0)( )XRE Xt( )( )()XRE X t X t2(0)( )XRE Xt0( )XR令 =0( )()XXRR( )()XXCC( )( )(),()( )(),XXRE X t X tRE X t X tut 第39頁(yè)/共76頁(yè) 性質(zhì)4.8 平穩(wěn)過(guò)程X(t)協(xié)方差函數(shù)的最大點(diǎn)在 處 ，即 證 任何非負(fù)函數(shù)的數(shù)學(xué)期望恒為非負(fù)值，

27、 的平方均值，即 又 X(t)平穩(wěn)， 0(0)|( )|XXRR( )()X tX t2|( )()| 0EX tX t22( )2( )()()E XtX t X tXt022( )2( )()()0E XtE X t X tE Xt22( )2( )()()0E XtE X t X tE Xt( )()( )XE X t X tR22( )2( )()0 xE XtRE Xt第40頁(yè)/共76頁(yè) 平穩(wěn)隨機(jī)過(guò)程的一維率密度函數(shù)不隨時(shí)間的平緩而變化，即 同理 22()( ,)XE Xtx Px tdx( ; )( ;)XXPx tPx t22( )()(0)XE XtE XtR2(0)2( )0

28、XXRR(0)XXRR|( )|(0)XXCC|( )|22( )( , )XE Xtx Px t dxv對(duì)于 第41頁(yè)/共76頁(yè) 性質(zhì)4.9 周期平穩(wěn)過(guò)程X(t)的協(xié)方差函數(shù)是周期函數(shù)，且與周期平穩(wěn)地程的周期相同，即 證 設(shè) 性質(zhì)4.10 非周期平穩(wěn)過(guò)程X(t)的協(xié)方差函數(shù)滿足()( )XXRtTR( )()X tX tT()( )()( )()( )xXRTE X t X tTE X t X tR2lim( )( )XXXRRM 222( )(0)( )( )( )XXXD X tRRE XtEX t第42頁(yè)/共76頁(yè) 例4.6 非周期平穩(wěn)過(guò)程X(t)的協(xié)方差函數(shù)。 求 。 解 為了方便表

29、征隨機(jī)過(guò)程在兩個(gè)不同時(shí)刻狀態(tài)之間的線性關(guān)聯(lián)程度，我們給出協(xié)方差系定義： 定義協(xié)方差系數(shù)： 29( )161 3XR( ),( )E X tD X t229lim( )lim16161XXMR4XM 2(0)( )25 169XXXRR 222( )( )( )XXXXXXCRM第43頁(yè)/共76頁(yè) 特別取 ， 一般有 顯然，協(xié)方差系數(shù)越接近1，狀態(tài)之間的關(guān)聯(lián)程度越高。也可以說(shuō)，當(dāng)狀態(tài)與狀態(tài)之間的時(shí)間間隔越小，狀態(tài)之間的關(guān)系越高。因此相關(guān)系數(shù)可直觀地說(shuō)明隨機(jī)過(guò)程不同兩個(gè)狀態(tài)的協(xié)方差程度的強(qiáng)弱或隨機(jī)過(guò)程起伏的快慢。 0222(0)(0)(0)0XXXXXXCRM22222( )1XXXXE XtM|

30、( )|1X第44頁(yè)/共76頁(yè) 相關(guān)時(shí)間 是另一個(gè)表示隨機(jī)過(guò)程相關(guān)程度的量，它是利用相關(guān)系數(shù)來(lái)定義的。 一般相關(guān)時(shí)間的定義有兩種，一種是把滿足時(shí)的 作為相關(guān)時(shí)間 。其物理意義為：若隨機(jī)過(guò)程X(t)的相關(guān)時(shí)間為 ，則認(rèn)為隨機(jī)過(guò)程的時(shí)間間隔大于 的兩個(gè)時(shí)刻的取值不相關(guān)。另一種定義相關(guān)時(shí)間間隔大于的兩相時(shí)刻的取值不相關(guān)。 0()0.05X0000第45頁(yè)/共76頁(yè) 另一種定義相關(guān)時(shí)間方法是將 曲線在 之間的面積等效成 的矩形，如圖4.3所示。因此有( )X1,0(0)X00( )Xd圖4.3 協(xié)方差系數(shù)第46頁(yè)/共76頁(yè)2. 隨機(jī)過(guò)程協(xié)方差函數(shù)性質(zhì) 設(shè) 為兩個(gè)平穩(wěn)過(guò)程。 性質(zhì)4.11 一般情況下，互

31、相關(guān)函數(shù) 是非奇非偶函數(shù)，同理，互協(xié)方差函數(shù) 也是非奇非偶函數(shù)。 性質(zhì)4.12 互相關(guān)函數(shù)的幅度平方滿足 同理，互協(xié)方差函數(shù)滿足( ), ( )X t Y t( )XYR( )XYC2( )(0)(0)XYXYRRR222( )(0)(0)XYXYXYCCC第47頁(yè)/共76頁(yè) 性質(zhì)4.13 互相關(guān)函數(shù)和互協(xié)方差函數(shù)的幅度滿足 同理 性質(zhì)4.14 互相關(guān)系數(shù) 為了研究?jī)蓚€(gè)平穩(wěn)過(guò)程的相互關(guān)聯(lián)程度，我們引入互相關(guān)系定義 定義互相關(guān)系數(shù)： 可以證明 ，且當(dāng) 時(shí) 互不相關(guān)。1( )(0)(0)2XYXYRRR2211( )(0)(0)22XYXYXYCCC( )( )(0)(0)XYXYXYCCC( )

32、1XY( )0XY( )( )X tY t和第48頁(yè)/共76頁(yè)習(xí)題四 1. 考慮一個(gè)具有隨機(jī)相位的余弦波，它由如下定義的隨機(jī)過(guò)程描述： ，其中 是常數(shù)， 服從 上的均勻分布，證明X(t)是寬平穩(wěn)過(guò)程。 2. 考慮一個(gè)具有隨機(jī)振幅的正弦波，它由如下定義的隨機(jī)過(guò)程描述 其中，A、B為兩個(gè)隨機(jī)變量，且滿足 ， ，度X(t)為寬平穩(wěn)過(guò)程。 ( )cos()X tt(, ) ( )2sin2X tAxostt( )( )0,( )( )1E AE BD AD B()0E AB 第49頁(yè)/共76頁(yè) 3. 設(shè)隨機(jī)過(guò)程 是方差不為零的隨機(jī)變量，試討論其各態(tài)歷經(jīng)性。 4. 設(shè)X(t)是雷達(dá)的發(fā)射信號(hào)，遇到目標(biāo)后

33、返回接收機(jī)的微弱信號(hào)是 是信號(hào)返回時(shí)間，由于接收到的信號(hào)總是伴有噪聲，記噪聲為 ，于是接收機(jī)收到的全信號(hào) 。若X(t)和Y(t)是聯(lián)合平穩(wěn)，求互相關(guān)函數(shù) 。在的條件下，假如N(t)的均值為零，且X(t)是相互獨(dú)立，求 （這是利用互相關(guān)函數(shù)從全信號(hào)中檢測(cè)小信號(hào)的接收法）。( ),X tY Y111,(),aX ta( )N t1( )()( )Y taX tN t( )XYR( )XYR第50頁(yè)/共76頁(yè) 5. 設(shè)有隨機(jī)過(guò)程 ，其中A是具有瑞利分布的隨機(jī)變量，其概率密度為 是在（0, 2 ）上具有均勻分布且與A相互獨(dú)立的隨機(jī)變量， 是一個(gè)常數(shù)，問(wèn)X(t)是否是寬平穩(wěn)過(guò)程。( )cos ()X t

34、At2220( )200aaeat aa第51頁(yè)/共76頁(yè)3.3.功率譜密度功率譜密度 當(dāng)我們?cè)跁r(shí)間域內(nèi)研究某一函數(shù)的特性時(shí)，如果確定起來(lái)不方便，在數(shù)學(xué)上我們可以考慮將此函數(shù)通過(guò)某種變換將它變換到另一區(qū)域，比如說(shuō)頻率域內(nèi)進(jìn)行研究，最終目的是使問(wèn)題簡(jiǎn)化。傅里葉變換提供了一種方法，就是如何將時(shí)間域的問(wèn)題轉(zhuǎn)換到頻率域，進(jìn)而使問(wèn)題簡(jiǎn)化。在頻率域內(nèi)，頻率意味著信息變化的速度。即，如果一個(gè)信號(hào)有“高”頻成分，我們?cè)陬l率域內(nèi)就可以看到“快”的變化。這方面的應(yīng)用在數(shù)字信號(hào)分析和電路理論等方面應(yīng)用極廣。第52頁(yè)/共76頁(yè) 是不是任何一個(gè)時(shí)間函數(shù)都可以將其通過(guò)傅氏變換變到頻率域去研究呢？我們說(shuō)當(dāng)時(shí)間函數(shù) 滿足絕對(duì)

35、可積條件時(shí)可以。 然而，隨機(jī)過(guò)程的樣本函數(shù)，即 一般不滿足絕對(duì)條件，因此隨機(jī)過(guò)程不能直接進(jìn)行付氏變換。此外，很多隨要過(guò)程的樣本函數(shù)極不規(guī)則，無(wú)法用方程描述。這樣，若想直接對(duì)隨要過(guò)程進(jìn)行譜分解，顯然也不行。但是，對(duì)隨機(jī)過(guò)程進(jìn)行某種處理后，同樣可對(duì)隨機(jī)過(guò)程施行傅里葉變換。( )()X tt ( )x t dt 1( ) ( ),( ), nX tx tx t1( ),( )nx tx t第53頁(yè)/共76頁(yè)3.1 3.1 功率譜密度功率譜密度 為了研究隨機(jī)信號(hào)的傅氏變換，我們首先簡(jiǎn)單復(fù)習(xí)一下確定信號(hào) 的頻譜、能譜密度及能量概念，然后再引入隨機(jī)過(guò)程的功率譜密度概念。 定理3.1 設(shè)S(t)是一個(gè)確定信

36、號(hào)，且在 上，則S(t)的傅氏變換存在，或者說(shuō)具有頻譜 記為( )S t(,) ( )S( )( )j tSS t edt1( )( )2j tS tSed1( )( )FFS ts ( )( )F S tS第54頁(yè)/共76頁(yè) 一般頻譜 是一個(gè)復(fù)數(shù)，且有 ，*表示共軛。我們知道，對(duì)于復(fù)數(shù)有 對(duì)于定理的物理解釋是，或 代表電流或電壓，則定理?xiàng)l件要求 ，即是要求 的總能量必須有限。 由積分變換的巴塞伐能量公式有1 ( )( )FSS t( )S*( )()SS22*;,Zabi zabzabi2*2222(1)()( )z zbia biabiabz( )S t( )s t dt ( )S t22

37、1( )( )2St dtSd第55頁(yè)/共76頁(yè) 下面我們來(lái)解釋一下公式的物理含義 等式左邊表示 在 上的總能量，而右邊積分中被積函數(shù) 相應(yīng)地稱為能譜密度。巴塞伐公式理解為時(shí)間域上的總能量可用頻率域上的頻譜能量表示。 然而，工程技術(shù)上有許多重要的時(shí)間函數(shù)總能量是無(wú)限的，不能滿足傅氏變換絕對(duì)可積條件，如正弦 就是。我們要研究的隨機(jī)過(guò)程，由于持續(xù)時(shí)間是無(wú)限的，所以其總能量也是無(wú)限的，即 所以隨機(jī)過(guò)程的頻譜不存在。 ( )S t(,) 2( )Ssin( ) t dt ( ),1,2,nx tdtn第56頁(yè)/共76頁(yè) 那么該如何應(yīng)用傅氏變換工具來(lái)對(duì)隨機(jī)過(guò)程進(jìn)行化簡(jiǎn)研究呢？我們是這樣考慮的，一個(gè)隨機(jī)過(guò)

38、程 ，盡管它的樣本函數(shù)總能量是無(wú)限的，但它的平均功率是有限的，即 這是隨機(jī)過(guò)程的樣本函數(shù)在時(shí)間域上的平均功率表示。 這樣，對(duì)隨機(jī)過(guò)程的樣本函數(shù)而方，雖然研究它的頻譜沒(méi)有意義，但研究它的平均功率確有意義。1( ) ( ),( ),nX xx tx t21lim( )2TeTTWx tT 第57頁(yè)/共76頁(yè)圖3.1 .1 及其截取函數(shù) 怎樣具體表示隨機(jī)過(guò)程一個(gè)樣本函數(shù)的平均功率呢，我們是這樣操作的：首先定義 的一個(gè)樣本函數(shù)，不妨設(shè)為 ，再次地樣本函數(shù) 任意截取一段，長(zhǎng)度為2T，并記為 。稱 為原樣本函數(shù) 的截取函數(shù)，如圖5.1所示。( )x t( )Tx t( )X t( )x t( )x t(

39、)Tx t( )Tx t( )x t第58頁(yè)/共76頁(yè) 用公式表示即為 于是 滿足絕對(duì)可積條件。 存在付氏變換，即 這里 稱 為的頻譜函數(shù)。( )| |( )0| |Tx ttTx ttT( )Tx t dt ( )Tx t( )( )( )j tTj tTTTTxx t edtx t edt1( )()2jtTTxtxedt( )Tx( )Tx t第59頁(yè)/共76頁(yè) 又由于隨機(jī)過(guò)程 在隨機(jī)試驗(yàn)中取哪一個(gè)樣本函數(shù)具有不確定性。因此，不同的試驗(yàn)結(jié)果，就意味著隨機(jī)過(guò)程可能取不同的樣本函數(shù)，亦即樣本函數(shù)與試驗(yàn)結(jié)果有關(guān)，為此，可將樣本函數(shù)進(jìn)一步表示為 ，當(dāng)然該樣本函數(shù)的截取函數(shù)也可相應(yīng)表示為 ，顯然它

40、的傅氏變換也可表示為 。 又 1( ) ( ),( ),nX tx tx t( , )x t e( , )Tx t e( , )Tx t e21lim( )2TeTTWx tdtT21lim( , )2TTTx t edtT第60頁(yè)/共76頁(yè)11lim( , )( , )22Tj tTTTTx t exe eddtT()11lim( , )( , )22TjtTTTTxex t e edt dT11lim( , )(, )22TTtxe xe dT*11lim( , )( , )22TTtxe xe dT211lim( , )22TTxedT1( , )2XGe d第61頁(yè)/共76頁(yè) 由于引入

41、隨機(jī)過(guò)程樣本函數(shù)的截取函數(shù)定義，所以又可給出上式隨機(jī)過(guò)程的樣本函數(shù)平均功率在頻率域的表示形式。 在上式中，令 則稱（3.1）式為隨機(jī)過(guò)程X(t)的樣本函數(shù)的功率譜密度函數(shù)。 定義樣本函數(shù)的功率譜密度 式中， 為截取函數(shù) 的頻譜。21( , )lim( , )2XTTGexeT21( , )lim( , )2XTTGexeT( , )Txe( , )Tx t e第62頁(yè)/共76頁(yè) 又 隨機(jī)過(guò)程是由一族樣本函數(shù)組成，即 顯然對(duì)每一個(gè)樣本函數(shù)，按照上面類似的方法都呆以求出它的一個(gè)樣本函數(shù)的功率譜密度，于是對(duì)所有的樣本函數(shù)取統(tǒng)計(jì)平均就可給出隨機(jī)過(guò)程的功率譜密度定義。 定義隨機(jī)過(guò)程的功率譜密度：1( )

42、 ( ),( ),nX tx tx t221( )( , )lim( , )21lim( , )2XXTTTTGE GeExeTEXeT第63頁(yè)/共76頁(yè) 隨機(jī)過(guò)程的一個(gè)樣本函數(shù)的平均功率的表示形式，有兩種 類似的，可求出X(t)的所有樣本函數(shù)的平均功率表示形式，然后取統(tǒng)計(jì)平均，則可以給出隨機(jī)過(guò)程的平均功率定義，定義隨機(jī)過(guò)程的平均功率：11lim( , )22eTTWxe dT21lim( , )2TeTTTWx t edtT211lim( , )22eTTWE WExedT第64頁(yè)/共76頁(yè)211lim( , )22TTE xedT1( )2XGd21lim( , )2TeTTWE WEx t edtT221lim| ( , )| 21lim( )2TTTTTTEx t edtTE x tdtT第65頁(yè)/共76頁(yè) 由隨機(jī)過(guò)程平均功率定義可知，要求隨機(jī)過(guò)程的平均功率可用兩種方法，一種方法是求出 ，即過(guò)程的功率譜密度，然后再積分，另一種方法是先求出過(guò)程的增方值 ，再積分。 特別地，當(dāng)我們研究的隨機(jī)過(guò)程是平穩(wěn)過(guò)程時(shí)，此時(shí)的平穩(wěn)過(guò)程平均功率可表示為： X(t)平穩(wěn) ( )xG2( )E Xt22( )(0)XXE XtR2221lim( )21lim2TT