第六章 平稳随机过程

更新时间:2023-12-17 14:02:01 阅读量: 教育文库 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

第六章 平稳随机过程

在自然科学与工程技术研究中遇到的随机过程有很多并不具有Markov性,这就是说从随机过程本身随时间的变化和互相关联来看,不仅它当前的状况,而且它过去的状况都对未来的状况有着不可忽略的影响,并且其统计特征不随时间推移而变化,这类随机过程称为平稳过程. 例如,恒温条件下热噪声电压X(t)是由于电路中电子的热扰动引起的,这种热扰动不随时间推移而改变;又如,连续测量飞机飞行速度产生的测量误差X(t),它有很多因素(如仪器振动,电磁波干扰与气候等)造成,但主要因素不随时间推移而改变.

平稳过程是一种特殊的二阶矩过程,其表现在过程的统计特性不随时间的推移而改变.用概率论语言来描述:相隔时间h的两个时刻t与t?h处随机过程所处的状态X(t)与

X(t?h)具有相同的概率分布.一般地,两个n维随机向量?X(t1),X(t2),?,X(tn)?与

?X(t1?h),X(t2?h),?,X(tn?h)?具有相同的概率分布. 这一思想抓住了没有固定时间

(空间)起点的物理系统中最自然现象的本质,因而平稳过程在通讯理论、天文学、生物学、生态学、和经济学个领域中有着十分广泛的应用.

6.1 随机微积分

在高等数学的微积分中,连续、导数和积分等概念都是建立在极限概念的基础上.对于随机过程的研究,也需要建立在随机过程的连续性、可导性和可积性等概念的基础上,这些内容形式上与高等数学极为相似,但实质不同,高等数学研究的对象是函数,随机微积分研究的对象是随机函数(即随机过程),有关这部分的内容统称为随机分析(stochastic

analysis).

在随机分析中,随机序列极限的定义有多种,下面我们简单介绍常用的定义.由于我们主要研究广义平稳过程(具体的定义将在第二节介绍),因此,以下的随机过程都假定为二阶矩过程.为了讨论的方便,我们约定:今后如不加说明,二阶矩过程{X(t),t?T}的均值函数mX(t)?EX(t)?0,自协方差函数CX(s,t)?E?X(s)X(t)?.

??6.1.1 均方收敛

定义6.1 称二阶矩随机序列{Xn(?)}以概率为1收敛于二阶矩随机变量X(?),若使

limXn(?)?X(?)成立集合的概率为1,即

n??P?:limXn(?)?X(?)?1

n????或称{Xn(?)}几乎处处收敛(almost everywhere converge)于X(?),记作Xn

a.e.??? X.

第六章《平稳随机过程》

定义6.2 称二阶矩随机序列{Xn(?)}以概率收敛(convergence in probability)于二阶矩随机变量X(?),若对于任意给定的??0,有

limP?|Xn(?)?X(?)|????0

n??p? X. 记作Xn ??定义6.3 若二阶矩随机序列{Xn(?)}和二阶矩随机变量X(?)满足

limE[|Xn?X|2]?0 (6.1)

n??m.s? X. 成立,则称Xn均方收敛(convergencein mean square)于X,记作Xn ??(6.1)式的极限常常写成l?i?mXn?X或l?i?mXn?X(.l?i?m是英文limit in mean的

n??缩写).

定义6.4称二阶矩随机序列{Xn(?)}依分布收敛(convergencein distribution)于二阶矩随机变量X(?). 若{Xn(?)}相应的分布函数列{Fn(x)},在X的分布函数每一个连续点处,有

Fnx(?)Fx( ) limn??d?X. 记作Xn ??对于以上四种收敛定义进行比较,有下列关系:

m.sp? X,则Xn ??? X; (1) 若Xn ??a.ep? X,则Xn ??? X; (2) 若Xn ??pd? X,则Xn ??? X. (3) 若Xn ??值得注意的是,在四种收敛定义中,均方收敛是最简单的收敛形式,它只涉及单独一个

序列.下面我们讨论随机序列的收敛性,都是指均方收敛.

定理6.1 二阶矩随机序列{Xn}收敛于二阶矩随机变量X的充要条件是

X|n?Xm limE[n,m??2|?] 0定理6.2 设{Xn},{Yn},{Zn}都是二阶矩随机序列,U为二阶矩随机变量,{cn}为常数序列,a,b,c为常数.令l?i?mXn?X,l?i?mYn?Y,l?i?mZn?Z,limcn?c,则

n??n??n??n??(1)l?i?mcn?limcn?c;

n??n?? 2

第六章《平稳随机过程》

(2)l?i?mU?U;

n??(3)l?i?m(cnU)?cU;

n??(4)l?i?m(aXn?bYn)?aX?bY;

n??(5)l?i?mEXn?EX?E?l?i?mXn?;

n???n???(6)l?i?mE?XnYm??EXY?E?l?i?mXn?l?i?mYm?;

n,m???????n????m?????222特别地,有l?i?mE[|Xn|]?E|X|?E?|l?i?mXn|?

n???n??? 证明 (1),(2),(3),(4)由均方收敛的定义可以得证. (5)由Schwartz不等式 E|XY|?将X取为Xn?X,Y取为1,则有

E|X|2?E|Y|2 0?|EXn?EX|2?|E[Xn?X]|2?E|Xn?X|2?0 (n??)

因此 l?i?mEXn?EX?E?l?i?mXn?

n???n??? (6)由Schwartz不等式

|E[XnYm]?E[XY]|?|E[XnYm?XY]|

?E[(Xn?X)(Ym?Y)?XnY?XYm?2XY]

?E[(Xn?X)(Ym?Y)]?E[(Xn?X)Y]?E[(Ym?Y)X]

??????E??(Xn?X)(Ym?Y)??E?[(Xn?X)Y]??E?(Ym?Y)X?

?E|Xn?X|2E|Ym?Y|2?E|Xn?X|2E|Y|2 ?E|Ym?Y|2E|X|2 ?0

因此 l?i?mE[XnYm]?E[XY].

n,m??(5)式和(6)式表明:极限运算和求数学期望运算可以交换顺序.

定理6.3 二阶矩随机序列{Xn}均方收敛的充要条件是limE?XnXm?=c(常数)

n,m????证明 必要性由定理6.2之(6)易知,下证充分性. 设limE?XnXm??E|X|?c,由

n,m????2E|Xn?Xm|2?E[|Xn|2?XnXm?XnXm?|Xm|2]

?E|Xn|2?E[XnXm]?E[XnXm]?E|Xm|2

3

第六章《平稳随机过程》

因此 limE|Xn?Xm|?c?2c?c?0.

n,m??2定理6.3给出了判定二阶矩随机序列{Xn}均方收敛的方法,该条件称为洛弗(Loeve)准则.

6.1.2 均方连续

?t0) ,即定义6.5 设{X(t),t?T}是二阶矩过程,若对t0?T,有l?i?mX(t)?X(t?t0 limE?|X(t)?Xt(0)2|????0t?t0则称{X(t),t?T}在t0点均方连续(continuity in mean square). 如果{X(t),t?T}在t?T每点都均方连续,则称{X(t)}在T上均方连续.

定理6.4 (均方连续准则)二阶矩过程{X(t),t?T}在t点均方连续的充要条件是相关函数RX(t1,t2)在点(t,t)处连续.

证明 必要性:若l?i?mX(t?h)?X(t),由定理6.2中(6),得到

h?0limRX(t1,t2)?limE[X(t1)X(t2)]?E[X(t)X(t)]?RX(t,t)

t1?tt2?tt1?tt2?t充分性:若RX(t1,t2)在点(t,t)处连续,考虑到

E[|X(t?h)?X(t)|2]?RX(t?h,t?h)?RX(t,t?h)?RX(t?h,t)?RX(t,t)

令h?0取极限可得.

推论6.4.1 若相关函数RX(t1,t2)在{(t,t),t?T}上连续,则它在T?T上连续. 证明 若RX(t1,t2)在{(t,t),t?T}上连续,由定理6.4知X(t)在上均方连续,因此

l?i?mX(s)?X(t1),l?i?mX(s)?X(t2)

s?t1s?t2再由定理6.2中(6),得

limRX(t1,t2)?limE[X(s)X(t)]?E[X(t1)X(t2)]?RX(t1,t2)

s?t1t?t2s?t1t?t2知RX(t1,t2)在T?T上连续.

推论6.4.2 如果{X(t),t?T}是平稳过程,则X(t)在T上均方连续的充分必要条件是

X(t)的相关函数RX(?)在??0处连续,并且此时RX(?)是连续函数.

4

第六章《平稳随机过程》

证明:由于平稳过程的相关函数RX(?)本质上是RX(t,t??),所证结论很显然. 定理6.4表明:对于一般二阶矩过程在T上均方连续性与它的相关函数(作为二元函数)在T?T上连续性等价,而相关函数在T?T上的连续性又等价于它在第一、三象限平分线

{(t,t),t?T}上的连续性;对于平稳随机过程,均方连续等价于相关函数(作为一元函数)

在原点的连续性.

6.1.3 均方导数

定义6.6 设{X(t),t?T}是二阶矩过程,若存在另一随机过程X'(t),满足

X(t?h)?X(t)limE?X'(t)?0 h??h则称X(t)在t点均方可微(differentiability in mean square),记作

2X'(t)?dX(t)X(t?h)?X(t)?l?i?m h?0dth称X'(t)为X(t)在t点的均方导数.若X(t)在每点t都均方可微,则称它在T上均方可微.

类似地,若随机过程{X'(t),t?T}在t点均方可微,则称X(t)在t点二次均方可微,

d2X记为X''(t)或,称它为二阶矩过程X(t)的二阶均方导数.同理可定义高阶均方导数.

dt2定理6.5 (均方可导准则)二阶矩过程{X(t),t?T}在t点均方可微的充要条件是相关函数RX(t1,t2)在点(t,t)处广义二阶导数存在.

证明 由定理6.3知,X(t)在t点均方可微的充要条件为

?X(t?h1)?X(t)??X(t?h2)?X(t)?limE?? ??h1?0hh?1??2?h2?0存在,将其展开得

?R(t?h1,t?h2)?RX(t?h1,t)?RX(t,t?h2)?RX(t,t)?lim?X? h1?0hh12?h2?0?上式极限存在的充要条件是相关函数RX(t1,t2)在点(t,t)处广义二阶导数存在. 6.1.4 均方积分

设{X(t),t?T}是二阶矩过程,f(t)为普通函数,其中T?[a,b],用一组分点将T划分如下:a?t0?t1?????tn?b,记max{ti?ti?1}??n,作和式

1?i?n 5

第六章《平稳随机过程》

Sn??f(ti?)X(ti?)(ti?ti?1)

i?1n其中ti?1?ti??ti(i?1,2,?,n).

定义6.7 如果当?n?0时,Sn均方收敛于S,即

?n?0limE|Sn?S|2?0

则称f(t)X(t)在区间[a,b]上均方可积(integral in mean square ),并记

S??f(t)X(t)dt? l?i?m?f(ti?)X(ti?)(ti?ti?1) (6.2)

a?n?0i?1bn称(6.2)式为f(t)X(t)在区间[a,b]上的(Riemann)均方积分. 需要说明的是:均方积分

b?baf(t)X(t)dt是一个随机变量,而不是一个随机过程.,当

nf(t)?1时,?X(t)dt? l?i?m?X(ti?)(ti?ti?1)

a?n?0i?1定理6.6(均方可积准则)f(t)X(t)在区间[a,b]上均方可积的充要条件是

??abbaf(t1)f(t2)RX(t1,t2)dt1dt2

存在,特别地,二阶矩过程X(t)区间[a,b]上均方可积的充要条件是RX(t1,t2)在

[a,b]?[a,b]上可积.

定理6.7 (数学期望与积分交换次序)f(t)X(t)在区间[a,b]上均方可积,则有

bbbb????(1)E?f(t)X(t)dt??f(t)E[X(t)]dt,特别地E?X(t)dt??E[X(t)]dt; ?????a?a?a?a(2)E?bbb?b??f(t1)f(t2)RX(t1,t2)dt1dt2 f(t)X(t)dtf(t)X(t)dt111?a222???aa?a??2特别地,E?baX(t)dt??ba?baRX(t1,t2)dt1dt2.

证明 由定理6.2之(5),有

E?ban??f(t)X(t)dt?E?l?i?m?f(ti?)X(ti?)(ti?ti?1)?

??n?0i?1? 6

第六章《平稳随机过程》

?n? ?limE??f(ti?)X(ti?)(ti?ti?1)?

?n?0?i?1? ?lim类似地,可证明(2)式.

均方积分有类似于普通函数积分的许多性质,如X(t)均方连续,则它均方可积;均方积分唯一性;对于a?c?b,有

?n?0?f(ti?)E?X(ti?)(ti?ti?1)???f(t)E[X(t)]dt

i?1nba?baf(t)X(t)d?t?caf(t)X(t)?d?tbc(f)tX(;)t若dtX(t),Y(t在区间)[a,b]上均方连续,则

?[?X(t)??Y(t)]dt???abbaX(t)dt???Y(t)dt

ab其中?,?为常数,等等.

定理6.8 二阶矩过程{X(t),t?T}在区间[a,b]上均方连续,则

Y(t)??X(?)d?, (a?t?b)

at在均方意义下存在,且随机过程{Y(t),t?T}在[a,b]上均方可微,且有Y?(t)?X(t). 推论 设X(t)均方可微,且X?(t)均方连续,则

X(t)?X(a)??X?(t)dt (6 .3)

at特别地,X(b)?X(a)??baX?(t)dt

上式相当于普通积分中的Newton?Leibniz公式.

最后,对本节的内容作一些说明.

(1)均方积分可以把区间[a,b]推广到无穷区间上,得到广义均方积分;

(2) 均方连续、均方导数、均方可积对复随机过程依然适应,但要把前面的绝对值理解为复数的模;

(3)均方连续、均方可导、均方可积都取决于相关函数的性质;

(4)在计算均方导数与均方积分时,可以把随机过程当成普通函数来处理; (5)均方导数是随机过程,均方极限与均方积分都是随机变量.

6.2 平稳过程及其相关函数

平稳过程作为特殊的二阶矩过程在工程技术中有着广泛的应用.

定义6.8 设{X(t),t?T}是随机过程,如果对任意常数?和正整数n,t1,t2,?,tn?T,

7

第六章《平稳随机过程》

t1??,t2??,?,tn???T,

?X(t1),X(t2),?,X(tn)?与?X(t1??),X(t2??),?,X(tn??)?

有相同的联合分布,则称{X(t),t?T}为严平稳过程,也称狭义平稳过程.

定义6.9 设{X(t),t?T}是随机过程,如果 (1){X(t),t?T}是二阶矩过程; (2)对任意t?T,mX(t)?EX(t)?常数;

(3)对任意s,t?T,RX(s,t)?E[X(s)X(t)]?RX(s?t),则称{X(t),t?T}为广义平稳过程,也称平稳过程(stationary process).

若T为离散集,则称平稳过程{X(t),t?T}为平稳序列(stationary sequence). 比较两种定义:广义平稳过程对时间推移的不变性是表现在统计平均的一阶矩、二阶矩上,而严平稳过程对时间推移的不变性是表现在概率分布上. 两者的要求是不一样的,一般来说,严平稳过程要求的条件比广义平稳过程要求的条件要严格得多. 显然,广义平稳过程不一定是严平稳过程;反之,严平稳过程只有当二阶矩存在时为广义平稳过程. 值得注意的是对于正态过程来说,二者是一样的.

例6.1 设随机过程X(t)?Ycos(?t)?Zsin(?t),t?0.

其中,Y,Z是相互独立的随机变量,且EY?EZ?0,DY?DZ??,则

2EX(t)?EYcos(?t)?EZsin(?t)?0

RX(s,t)?E[X(s)X(t)]?E[Ycos(?s)?Zsin(?s)][Ycos(?t)?Zsin(?t)]

=cos(?s)cos(?t)EY?sin(?s)sin(?t)EZ ??cos[(t?s)?] 因此,{X(t),t?0}为广义平稳过程.

例6.2 (随机电报信号过程)设随机过程{N(t),t?0}是具有参数为?的Poisson过程,随机过程{X(t),t?0}定义为:若随机点在[0,t]内出现偶数次,则X(t)?1;若出现奇数次,则X(t)??1.(1)讨论随机过程X(t)的平稳性;(2)设随机过程V具有概率分布

222P{V?1}?P{V??1}?12

8

第六章《平稳随机过程》

且V与X(t)独立,令Y(t)?VX(t),试讨论随机过程Y(t)的平稳性.

解 (1)由于随机点N(t)是具有参数为?的Poisson过程,因此,在[0,t]内随机点出

现k次的概率 Pk(t)?e??t(?t)k,k?0,1,2,? k!因此 P{X(t)?1}?P0(t)?P2(t)?P4(t)????

(?t)2(?t)4?e[1??????]?e??tch(?t)

2!4!??tP{X(t)??1}?P1(t)?P3(t)?P5(t)????

(?t)3(?t)5?????]?e??tsh(?t) ?e[?t?3!5!??t于是 mX(t)?EX(t)?1?e??tch(?t)?1?e??tsh(?t)

?e??t[ch(?t)?sh(?t)]?e??t?e??t?e?2?t.

为了求X(t)的相关函数,先求X(t1),X(t2)的联合分布

P{X(t1)?x1,X(t2)?x2}?P{X(t2)?x2|X(t1)?x1}P{X(t1)?x1}

其中xi??1或1(i?1,2).

设t2?t1,令??t2?t1,因为事件{X(t1)?1,X(t2)?1}等价于事件{X(t1)?1,且在

(t1,t2]内随机点出现偶数次}.由假设知,在X(t1)?1的条件下,在区间(t1,t2]内随机点出现

偶数次的概率与在区间(0,?]内随机出现偶数次的概率相等,故

P{X(t2)?1|X(t1)?1}?e???ch(??)

由于 P{X(t1)?1}?e??t1ch(?t1)

??t1所以 P{X(t1)?1,X(t2)?1}?ech(?t1)e???ch(??). sh(?t1)e???ch(??);

类似可得 P{X(t1)??1,X(t2)??1}?e??t1P{X(t1)??1,X(t2)?1}?e??t1sh(?t1)e???sh(??); P{X(t1)?1,X(t2)??1}?e??t1ch(?t1)e???sh(??)

9

第六章《平稳随机过程》

因此 RX(t )]?E[X((t1,t2)1t)X2?1?1?e??t1ch(?t1)e???ch(??)+(?1)?(?1)?e??t1sh(?t1)e???ch(??) ?(?1)?1?e??t1sh(?t1)e???sh(??)?1?(?1)?e??t1ch(?t1)e???sh(??) ?e??(t1??)[ch?(??t1)?sh?(??t1)]

?e??(t1??)e??(??t1)?e?2???e?2?(t2?t1).

当t2?t1,同理可得

RX(t1,t2)?e2?(t2?t1)?e2??

因此,对于任意t1,t2,有 RX(t1,t2)?e?2?|t2?t1||| ?e?2?? 由于mX(t)?e?2?t与时间t有关,故X(t)不是平稳随机过程,值得注意的是非平稳过程相关函数也可以与时间起点无关.

(2)由于EV?0,EV2?1,由V与X(t)独立知

EY(t)?EVEX(t)?0

||RY(t,t??)?EV2E[X(t)X(t??)]?e?2???RY(?)

所以,Y(t)是平稳过程.

例6.3 设X(t)?Xf(t)为复随机过程,其中X是均值为0的实随机变量,f(t)是确定函数. 证明X(t)是平稳过程的充要条件是f(t)?cei(?t??).其中i??1,c,?,?为常数.

证明 充分性:若f(t)?cei(?t??),记DX??,因此

2mX(t)?EX(t)?E[Xf(t)]?0

RX(t,t??)?E[X(t)X(t??)]?EX2c2ei(?t??)e?i[?(t??)??]?c2?2ei??

所以,X(t)是平稳过程.

必要性:若X(t)是平稳过程,则

RX(t,t??)?E[X(t)X(t??)]?EX2f(t)f(t??)

10

第六章《平稳随机过程》

上式必须与t无关,取??0,有|f(t)|2?c2(常数) 因此,f(t)?cei?(t),其中?(t)为实函数,于是 f(t)f(t??)?c2exp{i[?(t)??(t??)]} 上式应与t无关,因此

d[?(t)??(t??)]?0 dt即

d?(t)d?(t??)?对一切?成立,于是 dtdt?(t)??t??.

故 f(t)?cei(?t??).

例6.3显示了相关函数在平稳过程中的重要性,平稳过程的统计特性往往通过相关函数来表现.

例6.4 (随机相位过程)给定随机相位过程X(t)??(???),其中?(t)时周期为l的函数,?是服从(0,l)上均匀分布的随机变量,讨论其平稳性.

l11t?l1l??(t??)?d???(s)ds??(s)ds 解 mX(t)?EX(t)?E?(t??)???0t0lll与t无关;

l1?(t??)?(t????)?d ??)?E(?t??)(??t????0l1l?t1l ???(s)?(s??)ds???(s)?(s??)ds

lll0与t无关. 因此,随机相位周期过程是平稳过程.

RX(t,t??)?EX(t)X(??t 下面我们来讨论联合平稳过程及相关函数的性质.

定义6.10 设{X(t),t?T}和{Y(t),t?T}是两个平稳随机过程,若它们的相关函数

???E??X(t)Y(t??)?及E?Y(t)X(t??)?仅与?有关,而与t无关,则称X(t)和Y(t)是联合

平稳随机过程. 由定义有

?RXY(t,t??)?E??X(t)Y(t??)??RXY(?) ?RYX(t,t??)?E??Y(t)X(t??)??RYX(?)

当两个平稳过程X(t)和Y(t)是联合平稳随机过程时,则它们的和W(t)是平稳过程,此时有

11

第六章《平稳随机过程》

?E??W(t)W(t??)??RX(?)?RY(?)?RXY(?)?RYX(?)?RW(?)

定理6.9 (相关函数的性质) 设{X(t),t?T}是平稳过程,则其相关函数RX(?)具有下列性质:

(1)RX(0)?0; (2)RX(?)?RX(??); (3)|RX(?)|?RX(0);

(4)(非负定性)对于任意实数t1,t2,???,tn及复数?1,?2,???,?n,有

i,j?1?RnX(ti,tj)?i?j?0

(5)若X(t)是周期为T的周期函数,即X(t)?X(t?T),则 RX(?)?RX(??T)

(6)若X(t)是不含周期分量的非周期过程,当|?|??时,X(t)与X(t??)相互独立,则limRX(?)?mXmX

|?|??证明 由平稳过程相关函数的定义,得

(1)RX(0)?E?X(t)X(t)??E|X(t)|?0;

??2(2)RX(?)?E?X(t)X(t??)??E?X(t??)X(t)??RX(??); ??????对于实平稳过程,由于RX(?)为实数,因此,RX(??)?RX(?),即实平稳过程的相关函数为偶函数.

(3)由Schwartz不等式有

???EX(t)X(t??)??E|X(t)|2EX(t??) E?X(t)X(t??)????即 |RX(?)|?[RX(0)],因此|RX(?)|?RX(0);

(4) 显然;

22222???(5) RX(??T)?E??X(t)X(t???T)??E?X(t)X(t??)??RX(?);

(6) limRX(?)?limE?X(t)X(t??)??limEX(t)EX(t??)]?mXmX.

|?|??|?|????|?|??类似地,联合平稳过程X(t)和Y(t)的互相关函数具有下列性质:

12

第六章《平稳随机过程》

(1)|RXY(?)|2?RX(0)RY(0), |RYX(?)|2?RX(0)RY(0); (2)RXY(??)?RYX(?). 证明 (1)由Schwartz不等式,

|RXY(?)|2?|E[X(t)Y(t??)]|2?[E|X(t)Y(t??)|]2

?E|X(t)|2E|Y(t??)|2?RX(0)RY(0);

(2)RXY(??)?E[X(t??)Y(t)]?E[Y(t)X(t??)]?RYX(?).

当X(t)和Y(t)是实联合平稳过程时,(2)式变成RXY(??)?RYX(?).,这表明RXY(?)与RYX(?)在一般情况下是不相等的,且它们不是?的偶函数.

例6.5 设X(t)?Asin(?t??), Y(t)?Bsin(?t????)是两个平稳过程,其中

A,B,?为常数,?在(0,2?)上服从均匀分布,求RXY(?)和RYX(?).

解 RXY(?)?E[Xt(Y)t?(??)E[Asin?(t??B)s?int(???????

1d? 2????2?0ABsin(?t??)sin(?t???????)AB2?sin(?t??)[sin(?t??)cos(????) ?cos(?t??)sin(????)]d? ?02?1?ABcos(????). 2同理可得

RYX(?)?

1ABcos(????). 26.3 平稳过程的各态历经性

平稳随机过程的统计特征完全由前二阶矩函数确定,为了研究平稳过程的相关理论,必

须先明确均值函数与相关函数.但在实际应用中,随机过程的均值函数与相关函数一般是未知的,需要先通过大量的观察试验获得样本函数,然后用数理统计的点估计理论作出估计,其要求是很高的.为了提高估计的精度,需要做出多次试验,以获得许多样本函数.限于人力和财力,更限于试验周期等原因,这是不现实的.然而,对于平稳过程,它的均值函数是常数,相关函数只与时间间隔有关,它们都与起始时刻无关,也就是说,平稳过程的统计特性不随时间推移而改变,这就提供了一个是否在较宽的条件下,用样本函数估计平稳过程均值与相关函数的方法,它需要平稳过程具有各态历经性,即遍历性.

各态历经性的理论依据是大数定律.大数定律表明:随时间n的无限增大,随机过程的样本函数按时间平均以越来越大的概率近似于过程的统计平均.也就是说,时间平均与状态平均殊途同归,它的直观含义是:只要观测的时间足够长,随机过程的每一个样本函数都能

13

第六章《平稳随机过程》

够“遍历”各个可能状态.

定义6.11 设{X(t),???t??}为均方连续的平稳过程,称

X(t)?l?i?mT??12T?T?TX(t)dt (6.4)

T为该过程的时间均值;称

1X(t)X(t??)?l?i?mT??2T为时间相关函数.

??TX(t)X(t??)dt (6.5)

定义6.12 设{X(t),???t??}为均方连续的平稳过程,若?X(t)??EX(t),a.s.即

l?i?mT??12T?T?TX(t)dt?mX (6.6)

以概率为1成立,则称该平稳过程的均值具有各态历经性.

若X(t)X(t??)?E?X(t)X(t??)?,即

??l?i?mT??12T?T?TX(t)X(t??)dt?RX(?) (6.7)

则称该平稳过程的相关函数具有各态历经性.

如果均方连续平稳过程的均值和相关函数都具有各态历经性,则称该平稳过程具有各态历经性或遍历性(ergodicity),或称X(t)是各态历经过程(ergodic process).

由上面的讨论知,如果X(t)是各态历经过程,则X(t)和X(t)X(t??)不再依赖

??,而是以概率为1分别等于EX(t)和E??X(t)X(t??)?,这一方面表明各态历经过程各

样本函数的时间平均实际上可以认为是相同的,于是,对随机过程的时间平均也可以用样本函数的时间平均来表示,且可以用任一个样本函数的时间平均代替随机过程的统计平均;另一方面也表明EX(t)和E?X(t)X(t??)?必定与时间t无关,即各态历经过程必定是平稳过

??程.但是平稳过程只有在一定的条件下才是各态历经过程.

例6.6 随机相位正弦波X(t)?acos(?t??),???t??具有各态历经性,其中?是

(0,2?)上均匀分布的随机变量.

a2cos??,于是X(t)的时间平均为 容易求得mX?0,RX(?)?21Tacos(?t??)dt

T??2T??TaT?m?stc?o?s?stin??sdi nt ?l?i?coT??2T??TTaacos?si?nT?mco?s?c?otdst?l?i?m? 0 ?l?i?TT??2TT???TX(t)?l?i?m

14

第六章《平稳随机过程》

X(t)的时间相关函数为

X(t)X(t??)?l?i?mT??12T?T?Ta2cos(?t??)cos(?(t??)??)dt

a2?l?i?mT??4Ta2??T??cos(2?t????2?)?cos???dt?2cos??

T上述结果表明:随机相位正弦波X(t)的均值与相关函数都具有各态历经性,从而X(t)具有各态历经性.

下面我们讨论平稳过程具有遍历性的条件.

定理6.10 设{X(t),???t??}为均方连续的平稳过程,则它的均值具有各态历经性的充要条件是

1T??2Tlim?|?|?2??1?R(?)?|m|d??0 (6.8) XX????2T???2T?2T 证明 因X(t)是随机变量,先求它的期望与方差

1?EX(t)?E?l?i?m?T??2T1??lim X(t)dt??T?T??2T?T?T?TE[X(t)]dt?mX

因此,随机变量X(t)的均值函数为常数EX(t)?mX,由方差的性质,若能证明

DX(t)?0,则X(t)依概率为1等于EX(t).因此,要证明X(t)的均值具有各态历经

性等价于证明D?X(t)??0,由于

D?X(t)??E|?X(t)?|2?|mX|2 (6.9)

1而 E|X(t)|2?El?i?mT??2T?1?limE?2T???4T1T??4T21?lim2T??4T?limT?TT?T2?TX(t)dt

T? X(t)dtX(t)dt??T22??T11??TT?????T?TT?E??X(t2)X(t1)?dt1dt2 RX(t2?t1)dt1dt2

?(t1,t2)1?

?(?1,?2)2?T作变换?1?t1?t2,?2?t2?t1,变换的雅可比式为于是

15

第六章《平稳随机过程》

1 E|?X(t)?|?limT??4T221??2T??2T?|?2|2RX(?2)d?1d?2

2T2T??2| ?lim1T??2T?|?2|?R(?)1??d?2 (6.10) ??2TX2?2T??2T1又因为

2T?|?2|?1??d?2?1 ??2T?2T??2T故 |mX|2?12T|?2|?2? |m|1??d?2 (6.11)??2TX?2T??2T将(6.10) 和(6.11)代人(6.9),得

1 D?X(t)??limT??2T?|?|?21???RX(?)?|mX|?d? (6.12) ??2T??2T?2T(6.12)式等于0就是?X(t)?以概率1等于EX(t)?mX的充要条件,证毕.

当X(t)是实均方连续平稳过程时,RX(?)为偶函数,过程X(t)的均值各态历经性的充要条件可以写成

12T???2??(6.13) 1?R(?)?mX?d??0 ???XT??T?0?2T?lim由于CX(?)?RX(?)?|mX|2,因此,(6.8)式等价于

1limT??2T相应地,(6.13)式等价于

?|?|? 1??CX(?)d??0 (6.14)??2T?2T??2T12T??? 1???CX(?)d??0 (6.15)T??T?0?2T?lim定理6.11 设{X(t),???t??}为均方连续的平稳过程,则它的相关函数具有各态历经性的充要条件是

1T??2Tlim?|?1|?2??1?C(?)?|R(?)|1X1????2T?2T???d?1?0 (6.16)

2T其中CX(?1)?E?X(t)X(t??)X(t??1)X(t????1)? (6.17)

????证明 对于固定的?,记Y(t)?X(t)X(t??),则Y(t)为均方连续的平稳过程,且

?mY?EY(t)?E??X(t)X(t??)??RX(?)

16

第六章《平稳随机过程》

因此,RX(?)的各态历经性相当于EY(t)的各态历经性,由于

?X(t)X(t??)X(t??)X(t????)??C(?) ?RY(?1)?E?111??Y(t)Y(t??1)??E???由定理6.10得定理6.11成立.

定理6.12 对于均方连续平稳过程{X(t),0?t???},等式

1l?i?mT??T以概率1成立的充要条件为

?T0X(t)dt?mX (6.18 )

1 limT??2T|??|? (6.19 ) ?1)? 0?CX?(d???T?T??T若X(t)为实随机过程,则上式变为 lim1T????1??CX(?)d??0. T??T?0?T?定理6.13 对于均方连续平稳过程{X(t),0?t???},等式

l?i?mT??1TX(t)X(t??)d??RX(?) (6.20) T?0以概率1成立的充要条件为

1T?|?1|?2?? 1?C(?)?|R(?)|d?1?0 (6.21)1X????T??T??TT??lim1T??1?2若X(t)为实随机过程,则上式变为 lim??1???C(?)?R(?)?1X??d?1?0. T??T0T?? 例6.7 (续例6.2)考虑例6.2中随机电报信号过程Y(t)均值的各态历经性.

||因为它是实平稳过程,且EY(t)?0,RY(?)?e?2??,因此

12T????2?|?|lim??1??0????e?d??0 T??T02T??由(6.13知,Y(t)是均值具有各态历经性的平稳过程.

例6.8 讨论随机过程X(t)?Y的各态历经性,其中Y是方差不为0的随机变量. 解 容易知道X(t)?Y是平稳过程,事实上,EX(t)?EY?mX(常数),

2(与t无关),但此过程不具有各态历经性,因为 RX(t,t??)?EY2?DY?mXX(t)?l?i?mT??12T?T?TYdt?Y

17

第六章《平稳随机过程》

Y不是常数,不等于EX(t),因此,X(t)?Y的均值不具有各态历经性.类似地可证相关函

数也不具有各态历经性.

实际应用中,要严格验证平稳过程是否满足各态历经性条件是比较困难的,但各态历经性定理的条件较宽,工程中所遇到的平稳过程大多数都能满足. 因此,通常的处理方法是:先假设平稳过程是各态历经过程,然后由此假定出发,对各种数据进行分析处理,在实践中考察是否会产生较大的偏差,如果偏差较大,便认为该平稳过程不具有各态历经性.

各态历经性定理的重要意义在于它从理论上给出了如下的结论:一个实平稳过程,如果它是各态历经的,则可用任意一个样本函数的时间平均代替平稳过程的统计平均,即

mX?l?i?mT??1T1Tx(t)dtR(?)?l?i?mx(t)x(t??)dt ;X??00T??TT若样本函数x(t)只在有限区间[0,T]上给出,则对于实平稳过程有下面的估计式

?X? mX?m?(?)?RX(?)?RX1T?T0x(t)dt (6.22) x(t)x(t??)dt. (6.23)

1T???T??0(6.23) 式取积分区间[0,T??]是因为x(t??)只对t???T为已知,即0?t?T??.

习 题 六

6.1 设X1,X2,?是独立同分布随机变量,证明:随机序列{Xn,n?1}是严平稳时间序列.

6.2 设随机过程X(t)?Ucost?Vsint,???t??,其中U与V相互独立,且都服从

N(0,1).

(1) X(t)是平稳过程吗?为什么? (2) X(t)是严平稳过程吗?为什么?

6.3 设随机过程X(t)?Acos(?t??),???t??,其中,?为正常数,随机变量A与

aa???2exp{?2?},a?0?相互独立,且A的密度函数为f(a)??,?服从区间[0,2?]上的

其它??0,22均匀分布,求X(t)的均值函数与相关函数,并由此证明X(t)是平稳过程.

6.4设随机过程X(t)?sinUt,t?T,其中U服从区间[0,2?]上的均匀分布.

(1)如果T?{0,1,2,?},试求X(t)的均值函数与相关函数,并由此证明X(t)是平

18

第六章《平稳随机过程》

稳时间序列.

(2)如果T?[0,??],试求X(t)的均值函数,并由此证明X(t)不是平稳过程. 6.5 在习题6.2中,试求?X(t)?与?X(t)X(t??)?,并由此证明平稳过程X(t)的均值具有各态历经性,但相关函数不具有各态历经性.

6.6在习题6.3中,试求?X(t)?与?X(t)X(t??)?,并由此证明平稳过程X(t)的均值具有各态历经性,但相关函数不具有各态历经性.

6.7 证明相位周期过程X(t)??(t??)是各态历经过程,其中,?是有界函数.[提示:利用高等数学中周期函数的积分性质计算?X(t)?与?X(t)X(t??)?.

6.8 设平稳过程{X(t),???t??}的均值具有各态历经性,记随机过程

Y(t)?X(t)?U,其中,U是与X(t)不相关的随机变量,且EU?c,DU?1.

(1) 试求Y(t)函数与协方差函数,并由此证明Y(t)是平稳过程; (2) Y(t)函数是否具有各态历经性?为什么?

6.9 设有随机过程X(t)和Y(t)都不是平稳过程,且X(t)?A(t)cost,Y(t)?B(t)sint,其中A(t)和B(t)是均值为0的相互独立的平稳过程,它们有相同的相关函数,求证:

Z(t)?X(t)?Y(t)是平稳过程.

6.10 设X1(t),X2(t),Y1(t),Y2(t)都是均值为0的实随机过程,定义复随机过程

Z1(t)?X1(t)?iY1(t),Z2(t)?X2(t)?iY2(t)

求在下列情况下Z1(t)和Z2(t)的互相关函数.

(1) 所有实随机过程是相关的; (2) 所有实随机过程互不相关.

6.11 设X(t)是具有相关函数为RX(?)的平稳过程,令Y?实数,证明:E|Y|?2?a?TaX(t)dt,其中T?0,a为

?T?T(T?|?|)RX(?)d?.

2

6.12 设有随机过程X(t)?Asin(?t)?Bcos(?t),其中A,B是均值为0,方差为?的相互独立的正态随机变量.问:

(1)X(t)的均值是否具有各态历经性?

19

第六章《平稳随机过程》

(2)X(t)的均方值是否具有各态历经性?

(3)若A??2?sin?,B?2?cos?, ?是(0,2?)上均匀分布的随机变量,此时

E[X(t)]2是否具有各态历经性?

20

本文来源:https://www.bwwdw.com/article/rn25.html

Top