随机过程第四章和第五章

工程随机过程第四章平衡随机过程和各态历经过程在随机过程的大家族中，有一类随机过程，它的统计特性或者说统计变化规律与所选取的时间起点无关。或者说，整个随机过程的统计特性不随时间的推移而变化。例如，飞机在某一水平高度h上飞行时，由于受到气流的影响，实际飞行高度H(t)总是在理论设计高度h水平上下随机波动，此时飞机的实际飞行高度H(t)是一个随机过程，显然此过程可看作不随推移面变化的过程，这个随机过程，我们把它看作是平衡的随机过程。此外当我们知道一个随机过程是平稳过程时，它应不随时间的推移而变化。例如当我们要测定一个电阻的热噪声的统计特性，由于它是平稳过程，因而我们在任何时间进行测试都能得到相同的结果。§4.1严平稳随机过程及其数字特征

定义严平稳随机过程：对于任意的t，随机过程X(t)的任意n维概密度都有

则称X(t)为严平稳随机过程。研究平稳过程的意义在于：该过程在任何时刻计算它的统计结果都是相同的。由定义知平稳随机过程的n维概度密度函数不随时间而变化，这一特性具体反映在随机过程的一、二维概率密度及数字特征方面具有如下性质：性质4.1

若X(t)为平衡过程，则它的一维概率密度与时间无关。

性质4.1

若X(t)为平衡过程，则它的一维概率密度与时间无关证设X(t)的一维概率密度函数为，由于X(t)为平稳过程∴令则由此我们可求平稳过程X(t)的均值、均方值、方差。显然是一个数。这说明平稳随机过程的数学期望与时间无关。显然，X(t)的均方值、方差都与时间t无关。由此知，当随机过程为平稳过程时，该过程的所有样本函数总是它们均值——水平直线上下波动，样本曲线偏离水平直线的幅度正好是。如图4.1所示，图中细实线表示随机过程的样本函数，粗实线表示随机过程的数学期望，虚线表示随机过程对数学期望的偏差。性质4.2

平稳过程X(t)的二维概率密度只与的时间间隔有关，而与时间起点无关。性质4.2

平稳过程X(t)的二维概率密度只与的时间间隔有关，而与时间起点无关。证：设X(t)的二维概率密度函数为由于X(t)为平稳过程，所以对任意有若令，则而正是随机过程二维概率密度函数的时间间隔，令，则：此式表明，平稳随机过程的二维概率密度函数仅依赖于，而时间的个别值无关。由此，我们可以进一步来讨论平稳过程X(t)的自相关函数应具有什么样的表达形式。又∵

∴

顺便指出，由一个随机过程的平稳性研究可推广到关于两个随机过程的平稳性研究，可以这样说，若两个随机过程的联合概率密度函数不随时间的平移而变化，与时间的起点无关，则可称这两个随机过程是联合平衡的，或称平稳相依。§4.2宽平稳随机过程

从上面介绍的严平稳随机过程的定义知，要判断一个随机过程是否是严平稳，需要确定该随机过程的任意n维概率密度函数族，它的变化是否与时间的平稳无关，这本身就是一个十分困难的工作，然而在工程上根据实际需要，我们往往只在所谓的相关理论范围内考虑随机过程的平稳性问题，这里所指的相关理论，就是指随机过程的数字特征，即数学期望、相关函数和今后要介绍的功率普密度等。当在相关理论又可指研究随机过程的一、二阶矩理论。

前面已经介绍过，对于一个随机过程X(t)，我们当然希望能建立起它的多维分布函数，因为随机过程的多维分函数能较完整地描述随机过程的统计特性，但是要建立多维分布函数往往很困难，因此我们一般在相关理论范围内也就是用数字特征来描述过程的重要特性，这种用数字特征来描述过程X(t)统计特性变化规律，对很多实际问题往往已能获得很好的效果，可以提取到所需的参数。定义宽平稳过程：给定随机过程X(t)，如果常数且则称X(t)为宽平稳过程（广义平稳过程）。显然由宽平稳定义可知，要求就要考虑X(t)的一维概率密度函数和二维概率密度函数。下面我们来分析一下严平稳和宽平稳之间的关系。对于一个随机过程X(t)，如果它是严平稳的，且它的二阶矩存在及均方有界，则由严平稳双因严平稳的一维概率密度与时间无关，即∴常数又因严平稳的二维概率密度只与时间间隔有关，即∴综上所述，严平稳一定是宽平稳

反之不一定成立，除非是高斯过程（正态过程）。类似地，我们还可以给出两个随机过程联合宽平稳定义。定义联合宽平稳：对于平稳过程若则称联合宽平稳。

顺便指出，今后凡提到“平稳过程”，通常是指宽平稳过程。例4.1设Y是随机变量，试分别考虑随机过程的平稳性。

顺便指出，今后凡提到“平稳过程”，通常是指宽平稳过程。例4.1设Y是随机变量，试分别考虑随机过程的平稳性。解∵Y是随机变量，∵这一过程是一个与时间无关的特殊的过程，它的任何n维概率密度函数与时间无关，所以是一个严平稳。∵是严平稳,∴只要则X1(t)是宽平稳。对于都与时间有关，所以为非平稳。例4.2设是一周期为T的函数，是（0,T）上具有均匀分布的随机变量，称为随机相位周期过程，试讨论它的平稳性。解由题设知的概率密度函数为要讨论X(t)的平稳性，由宽平稳定义知，需要求。当取定为一随机变量的函数，由求随机变量函数的数学期望公式知∵令，则常数又∵令§4.3各态历经过程

1.各态历经问题的提出

对于一个随机过程X(t)，我们当然希望知道它们的分布函数，但很困难，于是我们退而求其次，考虑求它的数字特征即数学期望、相关函数等。但要求X(t)的数字特征，首先需要知道它的一、二维概率密度函数，即这实际上又很难办，进而为我们求数字特征又带来困难。怎么解决这个问题呢？实际上，在工程中，要求X(t)的数字特征，我们首先是通过试验来产生一族时间样本函数或者是做试验产生一个样本函数x(t)，然后再对样本函数x(t)取不同时刻，如，得所对应的结果，即此时随机过程可表示为。对任意指定时刻的数学期望可近似表示为自相关函数可近似表示为来计算，显然这种用近似计算的方法来估计随机过程的数学期望及自相关函数要求n很大，即样本函数xk(t)很多。但这在实际工程又常常又很难做到，于是人们自然想到能不能够通过测试一个样本函数如用一个样本函数xi(t)的均值和相关函数来近似随机过程的均值和相关函数，如果能，这为我们求随机过程的数学特征就带来了很大方便。这里提出一个问题：怎样表示一个样本函数如x1(t)的均值呢？我们以下式来表示显然x1(t)不同其积分结果一般不同。

于是对一个随机过程，，其样本函数的积数结果可能不同。此时显然用一个样本函数的数字特征如，近似是不正确的。但是如果当时间区间T充分大时，如果X(t)的绝大多数样本函数的均值都有则我们可用其中一个样本函数的均值作为[X(t)]的近似，即定义随机过程的时间均值和时间相关函数：称为随机过程的时间相关函数。注意：定义中一般都是随机变量（常数可看作特殊的随机变量）。由上述分析可知，是不是任何一个随机过程，它的数学期望、相关函数都可用其中的一个样本函数的均值和自相关函数来近似呢，显然不一定，一个自然的问题是X(t)在什么条件下可用一个样本函数的均值和自相关函数作为整个过程X(t)的均值，自相关函数的近似呢？2.平稳随机过程的各态历经性

要回答上述的问题，我们设当X(t)为平稳过程且满足一定条件时，可用一个样本函数的均值和自相关函数作为过程X(t)的数字特征近似，为此我们给出如下定义：定义：设X(t)是一个平稳过程（1）若

以概率1成立，则称随机过程X(t)均值具有各态历经性这里依概率1成立是指对X(t)的所有样本函数即由此知，此时，我们可用一个样本函数的均值如的值作为的近似值。反之，若已知X(t)的均值各态历程，则可用一个样本函数的均值作为过程X(t)的均值。（2）若

以概率1成立，则称X(t)的自相关函数具有各态历经性。

这里若X(t)的自相关函数各态历经，就是指我们可用过程X(t)的一个样本函数、xn(t)的时间相关函数

即作为过程的相关函数。（3）若X(t)的均值和自相关函数都具有各态历经性，则称X(t)是宽各态历经过程，简称X(t)为各态历经过程。综上所述，如果X(t)是各态历经过程，则必为平稳过程，此时可用过程的一个样本函数的数字特征作为过程的数字特征近似。例4.3设随机过程式中为参数，是（0.2,）上均匀分布随机变量。①求证X(t)是宽平稳过程；②该过程是否是各态历经过程。例4.3设随机过程式中为参数，是（0.2,）上均匀分布随机变量。①求证X(t)是宽平稳过程；②该过程是否是各态历经过程。解①∴X(t)为一宽平稳过程。②显然由①、②结果再由随机过程各态历经定义知∴X(t)为宽各态历经过程。如果两个随机过程X(t)，Y(t)，当它们各自都是各态历经时，并且时间互相关函数与统计相关函数以概率1相等时，我们有如下定义：定义两个随机过程联合各态历经：设X(t)，Y(t)各自都各态历经则称X(t)，Y(t)为联合各态历经过程。同理当X(t)，Y(t)联合各态历经时，可用它们的一对样本函数的数字特征作为X(t)，Y(t)的数字特征近似。3.

随机过程成为各态历经过程的判定

从前面的分析知，如果一个随机过程能成为一个平衡过程，这对我们研究各态历经非常方便。各态历经则该过程一定是平衡过程，反之则不一定成立，于是很自然提出这样一个问题，能不能给出一些判定定理，使其可以很方便地判定一个平稳过程成为各态历经过程。通过对平稳过程的分析研究，我们给出如下的几个判定定理。

性质4.3平稳过程X(t)的均值具有各态历经性的充要条件是

式中：为平稳过程的自相关函数；为平稳过程的数学期望。例4.4已知随机电报信号X(t)，它的，，问X(t)是否均值各态历经。例4.4已知随机电报信号X(t)，它的，，问X(t)是否均值各态历经。解∵∴X(t)是均值各态历经的。性质4.4平衡过程X(t)的自相关函数具备各态历经性的充要条件是式中平稳过程X(t)和Y(t)的互相关函数具有联合各态历经性的充要条件（4.7）式相似，只是将（4.7）式中相应的自相关函数改为互相关函数即可。性质4.5对于高斯平稳过程，如果它的均值为零，自相关函数连续，则该过程各态历经的一个充分条件是

综上所述，对一个平稳随机过程X(t)通过性质1、2判定以后，如果X(t)各态历经了，则对于该过程的数字特征，即求，我们可用也就是当X(t)各态历经时，我们可用一个样本函数的时间均值和时间自相关函数作为过程X(t)的数学期望、自相关函数的近似。最后顺便说明，对于许多实际问题，如果要从理论上判定一个过程是否为各态历经过程，往往是比较困难。因此工程上经常都是凭经验把各态历经性作为一种假设，然后根据实验来检验这个假设是否合理。在实际应用中一般不可能给出随机过程X(t)的样本函数x(t)的表达式，因此，确定各态历经过程的数学期望、自相关函数，有两种方法：第一种方法用模拟自相关分析仪，自动画出自相关曲线。这种仪器的功能是当输入样本函数时，X-Y记录仪自动描绘出自相关函数的曲线。它的方框图如图4.2所示。图4.2图4.3

第二种方法用数字处理方法（即近似计算方法）。如图4.3把[0，T]等分为N个长为的小区间，再在时刻，取样，得N个函数值。于是再把积分表达式表示为基本区间上的和，就有数字估计式（4.8）。类似可以写出在时的自相关函数估计式（4.9）式，由这个估计式可算出自相关函数的一系列近似值，从而可作出自相关函数的近似图形，见图4.4。图4.4最后指出，工程上遇到的很多平稳过程，我们一般都把它看成各态历经过程，然后用各态历经的方法来确定过程X(t)的统计特性，看处理出来的结果是否与实际相符合，如果不相符合，再对过程的假设作修改。

4.平稳过程自相关函数的性质对于一个随机过程，它的基本数字特征是数学期望和相关函数，但是当随机过程为平稳过程时，它的数学期望是一个常数，经中心化后可以变为零，所以当过程X(t)平稳后其基本数字特征实际上就是相关函数。此外，相关函数不仅可向我们提供随机过程各状态间的关联特性的信息，而且也是求取随机过程的功率谱密度以及从噪声中提取有用信息工具。为此下面我们专门研究一下平稳过程相关函数的性质。性质4.6证∵当，即平稳过程的均方值可自由相关函数得到。性质4.7，即平稳过程的自相关函数为偶函数。同理证∵性质4.8平稳过程X(t)自相关函数的最大点在处，即

性质4.8平稳过程X(t)自相关函数的最大点在处，即证∵任何非负函数的数学期望恒为非负值，的平方均值，即∴∴∴又∵X(t)平稳，∴∴∵平稳随机过程的一维率密度函数不随时间的平缓而变化，即∴∴∴同理对于

性质4.9周期平稳过程X(t)的自相关函数是周期函数，且与周期平稳地程的周期相同，即证设∵性质4.10非周期平稳过程X(t)的自相关函数满足例4.6非周期平稳过程X(t)的自相关函数。求。解∵∴∴为了方便表征随机过程在两个不同时刻状态之间的线性关联程度，我们给出自相关系数定义：定义自相关系数：

特别取，一般有显然，自相关系数越接近1，状态之间的关联程度越高。也可以说，当状态与状态之间的时间间隔越小，状态之间的关系越高。因此相关系数可直观地说明随机过程不同两个状态的自相关程度的强弱或随机过程起伏的快慢。相关时间是另一个表示随机过程相关程度的量，它是利用相关系数来定义的。一般相关时间的定义有两种，一种是把满足时的作为相关时间。其物理意义为：若随机过程X(t)的相关时间为，则认为随机过程的时间间隔大于的两个时刻的取值不相关。另一种定义相关时间间隔大于的两相时刻的取值不相关。另一种定义相关时间方法是将曲线在之间的面积等效成的矩形，如图4.3所示。因此有图4.3自相关系数5.设随机过程互相关函数性质设为两个平稳过程。性质4.11一般情况下，互相关函数是非奇非偶函数，同理，互协方差函数也是非奇非偶函数。性质4.12互相关函数的幅度平方满足同理，互协方差函数满足性质4.13互相关函数和互协方差函数的幅度满足同理性质4.14互相关系数为了研究两个平稳过程的相互关联程度，我们引入互相关系定义定义互相关系数：可以证明，且当时互不相关。习题四1.考虑一个具有随机相位的余弦波，它由如下定义的随机过程描述：，其中是常数，服从上的均匀分布，证明X(t)是宽平稳过程。2.考虑一个具有随机振幅的正弦波，它由如下定义的随机过程描述其中，A、B为两个随机变量，且满足，，度X(t)为宽平稳过程。3.设随机过程是方差不为零的随机变量，试讨论其各态历经性。4.设X(t)是雷达的发射信号，遇到目标后返回接收机的微弱信号是是信号返回时间，由于接收到的信号总是伴有噪声，记噪声为，于是接收机收到的全信号。①若X(t)和Y(t)是联合平稳，求互相关函数。②在①的条件下，假如N(t)的均值为零，且X(t)是相互独立，求（这是利用互相关函数从全信号中检测小信号的接收法）。5.设有随机过程，其中A是具有瑞利分布的随机变量，其概率密度为是在（0,2）上具有均匀分布且与A相互独立的随机变量，是一个常数，问X(t)是否是宽平稳过程。第五章随机过程的功率谱密度

当我们在时间域内研究某一函数的特性时，如果确定起来不方便，在数学上我们可以考虑将此函数通过某种变换将它变换到另一区域，比如说频率域内进行研究，最终目的是使问题简化。傅里叶变换提供了一种方法，就是如何将时间域的问题转换到频率域，进而使问题简化。在频率域内，频率意味着信息变化的速度。即，如果一个信号有“高”频成分，我们在频率域内就可以看到“快”的变化。这方面的应用在数字信号分析和电路理论等方面应用极广。是不是任何一个时间函数都可以将其通过傅氏变换变到频率域去研究呢？我们说当时间函数满足绝对可积条件时可以。然而，随机过程的样本函数，即一般不满足绝对条件，因此随机过程不能直接进行付氏变换。此外，很多随机过程的样本函数不规则，无法用方程描述。这样，若想直接对随机过程进行谱分解，显然也不行。但是，对随机过程进行某种处理后，同样可对随机过程施行傅里叶变换。§5.1功率谱密度

为了研究随机信号的傅氏变换，我们首先简单复习一下确定信号的频谱、能谱密度及能量概念，然后再引入随机过程的功率谱密度概念。定理5.1设S(t)是一个确定信号，且在上，则S(t)的傅氏变换存在，或者说具有频谱记为一般频谱是一个复数，且有，*表示共轭。我们知道，对于复数有∴对于定理的物理解释是，或代表电流或电压，则定理条件要求，即是要求的总能量必须有限。

由积分变换的巴塞伐能量公式有下面我们来解释一下公式的物理含义等式左边表示在上的总能量，而右边积分中被积函数相应地称为能谱密度。巴塞伐公式理解为时间域上的总能量可用频率域上的频谱能量表示。然而，工程技术上有许多重要的时间函数总能量是无限的，不能满足傅氏变换绝对可积条件，如正弦就是。我们要研究的随机过程，由于持续时间是无限的，所以其总能量也是无限的，即所以随机过程的频谱不存在。那么该如何应用傅氏变换工具来对随机过程进行化简研究呢？我们是这样考虑的，一个随机过程，尽管它的样本函数总能量是无限的，但它的平均功率是有限的，即这是随机过程的样本函数在时间域上的平均功率表示。这样，对随机过程的样本函数而方，虽然研究它的频谱没有意义，但研究它的平均功率确有意义。图5.1及其截取函数

怎样具体表示随机过程一个样本函数的平均功率呢，我们是这样操作的：首先定义的一个样本函数，不妨设为，再次将样本函数任意截取一段，长度为2T，并记为。称为原样本函数的截取函数，如图5.1所示。用公式表示即为于是满足绝对可积条件。∴存在付氏变换，即这里称为的频谱函数。又由于随机过程在随机试验中取哪一个样本函数具有不确定性。因此，不同的试验结果，就意味着随机过程可能取不同的样本函数，亦即样本函数与试验结果有关，为此，可将样本函数进一步表示为，当然该样本函数的截取函数也可相应表示为，显然它的傅氏变换也可表示为。又∵由于引入随机过程样本函数的截取函数定义，所以又可给出上式随机过程的样本函数平均功率在频率域的表示形式。在上式中，令则称（5.1）式为随机过程X(t)的样本函数的功率谱密度函数。定义样本函数的功率谱密度式中，为截取函数的频谱。又∵随机过程是由一族样本函数组成，即显然对每一个样本函数，按照上面类似的方法都呆以求出它的一个样本函数的功率谱密度，于是对所有的样本函数取统计平均就可给出随机过程的功率谱密度定义。定义随机过程的功率谱密度：随机过程的一个样本函数的平均功率的表示形式，有两种类似的，可求出X(t)的所有样本函数的平均功率表示形式，然后取统计平均，则可以给出随机过程的平均功率定义，定义随机过程的平均功率：由随机过程平均功率定义可知，要求随机过程的平均功率可用两种方法，一种方法是求出，即过程的功率谱密度，然后再积分，另一种方法是先求出过程的均方值，再积分。特别地，当我们研究的随机过程是平稳过程时，此时的平稳过程平均功率可表示为：∵X(t)平稳∴∴该式说明：平稳过程的平均功率等于该过程的均方值，也可由随机过程的功率谱密度在全频域上积分得到。若随机过程再各态历经，则各态历经过程的功率谱密度可用一个样本函数的功率谱密度来表示：例5.1随机过程式中，是常数，上均匀分布随机变量，求的平均功率。例5.1随机过程式中，是常数，上均匀分布随机变量，求的平均功率。解∵∴显然该过程不平稳。∴§5.2功率谱密度与自相关函数的关系通过对随机过程的分析，我们知道随机过程的相关函数是从时间有度描述了过程的重要统计特性，而随机过程的功率谱密度是从频率角度描述了过程的统计特性，二者是异曲同工，研究的都是一个对象，于是人们自然提出一个问题，随机过程的相关函数和它的功率谱密度之间是否存在一定关系，我们说当随机过程平稳且满足一定条件时，它们之间存在一定关系。定理5.2如果平稳过程X(t)的相关函数绝对可积，即则过程X(t)的相关函数和功率谱密度之间存在付氏变换，即例5.2设X(t)是平稳过程，共相关函数，其中、是正数，求X(t)的谱密度。例5.2设X(t)是平稳过程，自相关函数，其中、是正数，求X(t)的谱密度。解最后需要指出，在实际问题中常常碰到一些平稳过程，它的自相关函数或功率谱密度在通常意义可能付氏变换不存在。如果我们允许自相关函数或功率谱密度函数含有函数定义下自相关函数与功率谱密度存在付氏变换。定义函数：如果函数满足则称函数为函数。性质5.1若为无穷次可微函数，则或例5.3求的付氏变换。例5.4求的付氏逆变换。解∵又∵

∴例5.5若随机过程X(t)的自相关函数为求谱密度。例5.5若随机过程X(t)的自相关函数为解§5.3功率谱密度性质

由随机过程的功率谱密度定义，即可得如下几个常用的性质。性质5.2

∵

式中，必为非负。∴∴性质5.3是实函数。证∵是实函数∴是实函数∴是实函数性质5.4是偶函数。证∵∴

而∴∴∴是偶函数。性质5.5设X(t)为一个随机过程，为随机过程的导数，则证