《第2章随机过程精.ppt》由会员分享,可在线阅读,更多相关《第2章随机过程精.ppt(85页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第2章随机过程第1页,本讲稿共85页第第 2 章章 随机过程随机过程 2.1 随机过程的基本概念和统计特性随机过程的基本概念和统计特性 2.1.1随机过程随机过程 自然界中事物的变化过程可以大致分成为两类。一类是其变化过程具有确定的形式,或者说具有必然的变化规律,用数学语言来说,其变化过程可以用一个或几个时间t的确定函数来描述,这类过程称为确定性过程。例如,电容器通过电阻放电时,电容两端的电位差随时间的变化就是一个确定性函数。而另一类过程没有确定的变化形式,也就是说,每次对它的测量结果没有一个确定的变化规律,用数学语言来说,这类事物变化的过程不可能用一个或几个时间t的确定函数来描述,这类过程称
2、为随机过程。下面我们给出一个例子:第2页,本讲稿共85页 设有n台性能完全相同的接收机。我们在相同的工作环境和测试条件下记录各台接收机的输出噪声波形(这也可以理解为对一台接收机在一段时间内持续地进行n次观测)。测试结果将表明,尽管设备和测试条件相同,记录的n条曲线中找不到两个完全相同的波形。这就是说,接收机输出的噪声电压随时间的变化是不可预知的,因而它是一个随机过程。由此我们给随机过程下一个更为严格的定义:设Sk(k=1,2,)是随机试验。每一次试验都有一条时间波形(称为样本函数或实现),记作xi(t),所有可能出现的结果的总体x1(t),x2(t),,xn(t),就构成一随机过程,记作(t)
3、。简言之,无穷多个样本函数的总体叫做随机过程,如图 2-1 所示。第3页,本讲稿共85页图 2-1样本函数的总体 第4页,本讲稿共85页 显然,上例中接收机的输出噪声波形也可用图 2-1 表示。我们把对接收机输出噪声波形的观测可看作是进行一次随机试验,每次试验之后,(t)取图 2-1 所示的样本空间中的某一样本函数,至于是空间中哪一个样本,在进行观测前是无法预知的,这正是随机过程随机性的具体表现。其基本特征体现在两个方面:其一,它是一个时间函数;其二,在固定的某一观察时刻t1,全体样本在t1时刻的取值(t1)是一个不含t变化的随机变量。因此,我们又可以把随机过程看成依赖时间参数的一族随机变量。
4、可见,随机过程具有随机变量和时间函数的特点。下面将会看到,在研究随机过程时正是利用了这两个特点。第5页,本讲稿共85页 2.1.2随机过程的统计特性随机过程的统计特性 随机过程的两重性使我们可以用与描述随机变量相似的方法,来描述它的统计特性。设(t)表示一个随机过程,在任意给定的时刻t1T,其取值(t1)是一个一维随机变量。而随机变量的统计特性可以用分布函数或概率密度函数来描述。我们把随机变量(t1)小于或等于某一数值x1的概率P(t1)x1,简记为F1(x1,t1),即F1(x1,t1)=P(t1)x1 (2.1-1)式(2.1-1)称为随机过程(t)的一维分布函数。如果F1(x1,t1)对
5、x1的偏导数存在,即有第6页,本讲稿共85页 则称f1(x1,t1)为(t)的一维概率密度函数。显然,随机过程的一维分布函数或一维概率密度函数仅仅描述了随机过程在各个孤立时刻的统计特性,而没有说明随机过程在不同时刻取值之间的内在联系,为此需要进一步引入二维分布函数。任给两个时刻t1,t2T,则随机变量(t1)和(t2)构成一个二元随机变量(t1),(t2),称F2(x1,x2;t1,t2)=P(t1)x1,(t2)x2 (2.1-3)为随机过程(t)的二维分布函数。如果存在第7页,本讲稿共85页 则称f2(x1,x2;t1,t2)为(t)的二维概率密度函数。同理,任给t1,t2,tnT,则(t
6、)的n维分布函数被定义为Fn(x1,x2,xn;t1,t2,tn)=P(t1)x1,(t2)x2,(tn)xn 第8页,本讲稿共85页 则称fn(x1,x2,xn;t1,t2,tn)为(t)的n维概率密度函数。显然,n越大,对随机过程统计特性的描述就越充分,但问题的复杂性也随之增加。在一般实际问题中,掌握二维分布函数就已经足够了。第9页,本讲稿共85页 2.1.3随机过程的数字特征随机过程的数字特征 分布函数或概率密度函数虽然能够较全面地描述随机过程的统计特性,但在实际工作中,有时不易或不需求出分布函数和概率密度函数,而用随机过程的数字特征来描述随机过程的统计特性,更简单直观。1.数学期望 设
7、随机过程(t)在任意给定时刻t1的取值(t1)是一个随机变量,其概率密度函数为f1(x1,t1),则(t1)的数学期望为第10页,本讲稿共85页 注意,这里t1是任取的,所以可以把t1直接写为t,x1改为x,这时上式就变为随机过程在任意时刻的数学期望,记作a(t),于是a(t)是时间t的函数,它表示随机过程的n个样本函数曲线的摆动中心。2.方差方差E第11页,本讲稿共85页 D(t)常记为2(t)。可见方差等于均方值与数学期望平方之差。它表示随机过程在时刻t对于均值a(t)的偏离程度。均值和方差都只与随机过程的一维概率密度函数有关,因而它们描述了随机过程在各个孤立时刻的特征。为了描述随机过程在
8、两个不同时刻状态之间的联系,还需利用二维概率密度引入新的数字特征。3.相关函数相关函数 衡量随机过程在任意两个时刻获得的随机变量之间的关联程度时,常用协方差函数B(t1,t2)和相关函数R(t1,t2)来表示。协方差函数定义为第12页,本讲稿共85页 B(t1,t2)=E(t1)-a(t1)(t2)-a(t2)=f2(x1,x2;t1,t2)dx1dx2 式中,t1与t2是任取的两个时刻;a(t1)与a(t2)为在t1及t2时刻得到的数学期望;f2(x1,x2;t1,t2)为二维概率密度函数。相关函数定义为 B(t1,t2)=R(t1,t2)-a(t1)a(t2)第13页,本讲稿共85页二者关
9、系为B(t1,t2)=R(t1,t2)-a(t1)a(t2)(2.1-10)若a(t1)=0或a(t2)=0,则B(t1,t2)=R(t1,t2)。若t2t1,并令t2=t1+,则R(t1,t2)可表示为R(t1,t1+)。这说明,相关函数依赖于起始时刻t1及t2与t1之间的时间间隔,即相关函数是t1和的函数。由于B(t1,t2)和R(t1,t2)是衡量同一过程的相关程度的,因此,它们又常分别称为自协方差函数和自相关函数。对于两个或更多个随机过程,可引入互协方差及互相关函数。设(t)和(t)分别表示两个随机过程,则互协方差函数定义为第14页,本讲稿共85页 B(t1,t2)=E(t1)-a(t
10、1)(t2)-a(t2)(2.1-11)而互相关函数定义为 R(t1,t2)=E(t1)(t2)(2.1-12)第15页,本讲稿共85页2.2平稳随机过程平稳随机过程 2.2.1定义定义 所谓平稳随机过程,是指它的统计特性不随时间的推移而变化。设随机过程(t),tT,若对于任意n和任意选定t1t2tn,tkT,k=1,2,n,以及h为任意值,且x1,x2,xnR,有fn(x1,x2,xn;t1,t2,tn)=fn(x1,x2,xn;t1+h,t2+h,tn+h)(2.2-1)则称(t)是平稳随机过程。该定义说明,当取样点在时间轴上作任意平移时,随机过程的所有有限维分布函数是不变的,具体到它的一
11、维分布,则与时间t无关,而二维分布只与时间间隔有关,即有第16页,本讲稿共85页 f1(x1,t1)=f1(x1)(2.2-2)和 f2(x1,x2;t1,t2)=f2(x1,x2;)(2.2-3)以上两式可由式(2.2-1)分别令n=1和n=2,并取h=-t1得证。于是,平稳随机过程(t)的均值 为一常数,这表示平稳随机过程的各样本函数围绕着一水平线起伏。同样,可以证明平稳随机过程的方差2(t)=2=常数,表示它的起伏偏离数学期望的程度也是常数。而平稳随机过程(t)的自相关函数第17页,本讲稿共85页R(t1,t2)=E(t1)(t1+)=仅是时间间隔=t2-t1的函数,而不再是t1和t2的
12、二维函数。以上表明,平稳随机过程(t)具有“平稳”的数字特征:它的均值与时间无关;它的自相关函数只与时间间隔有关,即 R(t1,t1+)=R()注意到式(2.2-1)定义的平稳随机过程对于一切n都成立,这在实际应用上很复杂。但仅仅由一个随机过程的均值是常数,自相关函数是的函数还不能充分说明它符合平稳条件,为此引入另一种平稳随机过程的定义:第18页,本讲稿共85页 设有一个二阶矩随机过程(t),它的均值为常数,自相关函数仅是的函数,则称它为宽平稳随机过程或广义平稳随机过程。相应地,称按式(2.2-1)定义的过程为严平稳随机过程或狭义平稳随机过程。因为广义平稳随机过程的定义只涉及与一维、二维概率密
13、度有关的数字特征,所以一个严平稳随机过程只要它的均方值E2(t)有界,则它必定是广义平稳随机过程,但反过来一般不成立。通信系统中所遇到的信号及噪声,大多数可视为平稳的随机过程。以后讨论的随机过程除特殊说明外,均假定是平稳的,且均指广义平稳随机过程,简称平稳过程。第19页,本讲稿共85页 2.2.2各态历经性各态历经性 平稳随机过程在满足一定条件下有一个有趣而又非常有用的特性,称为“各态历经性”。这种平稳随机过程,它的数字特征(均为统计平均)完全可由随机过程中的任一实现的数字特征(均为时间平均)来替代。也就是说,假设x(t)是平稳随机过程(t)的任意一个实现,它的时间均值和时间相关函数分别为第2
14、0页,本讲稿共85页如果平稳随机过程依概率1使下式成立:则称该平稳随机过程具有各态历经性。“各态历经”的含义:随机过程中的任一实现都经历了随机过程的所有可能状态。因此,我们无需(实际中也不可能)获得大量用来计算统计平均的样本函数,而只需从任意一个随机过程的样本函数中就可获得它的所有的数字特征,从而使“统计平均”化为“时间平均”,使实际测量和计算的问题大为简化。第21页,本讲稿共85页 注注意意:具具有有各各态态历历经经性性的的随随机机过过程程必必定定是是平平稳稳随随机机过过程程,但但平平稳稳随随机机过过程程不不一一定定是是各各态态历历经经的的。在在通通信信系系统统中中所所遇遇到到的随机信号和噪
15、声,的随机信号和噪声,一般均能满足各态历经条件。一般均能满足各态历经条件。第22页,本讲稿共85页 2.2.3平稳随机过程自相关函数的性质平稳随机过程自相关函数的性质 对于平稳随机过程而言,它的自相关函数是特别重要的一个函数。其一,平稳随机过程的统计特性,如数字特征等,可通过自相关函数来描述;其二,自相关函数与平稳随机过程的谱特性有着内在的联系。因此,我们有必要了解平稳随机过程自相关函数的性质。设(t)为实平稳随机过程,则它的自相关函数 R()=E(t)(t+)(2.2-8)具有下列主要性质:(1)R(0)=E2(t)=S(t)的平均功率(2.2-9)(2)R()=E2(t)(t)的直流功率
16、(2.2-10)第23页,本讲稿共85页 这里利用了当时,(t)与(t+)没有依赖关系,即统计独立,且认为(t)中不含周期分量。(3)R()=R(-)的偶函数 (2.2-11)这一点可由定义式(2.2-8)得证。(4)|R()|R(0)R()的上界 (2.2-12)考虑一个非负式即可得证。(5)R(0)-R()=2方差,(t)的交流功率(2.2-13)当均值为0时,有R(0)=2。第24页,本讲稿共85页 2.2.4平稳随机过程的功率谱密度平稳随机过程的功率谱密度 随机过程的频谱特性是用它的功率谱密度来表述的。我们知道,随机过程中的任一实现是一个确定的功率型信号。而对于任意的确定功率信号f(t
17、),它的功率谱密度为 式中,FT()是f(t)的截短函数fT(t)(见图 2-2)所对应的频谱函数。我们可以把f(t)看成是平稳随机过程(t)中的任一实现,因而每一实现的功率谱密度也可用式(2.2-14)来表示。由于(t)是无穷多个实现的集合,哪一个实现出现是不能预知的,因此,某一实现的功率谱密度不能作为过程的功率谱密度。过程的功率谱密度应看做是任一实现的功率谱的统计平均,即 第25页,本讲稿共85页图 2-2 功率信号f(t)及其截短函数第26页,本讲稿共85页(t)的平均功率S则可表示成 虽然式(2.2-15)给出了平稳随机过程(t)的功率谱密度P(),但我们很难直接用它来计算功率谱。那么
18、,如何方便地求功率谱P()呢?我们知道,确知的非周期功率信号的自相关函数与其谱密度是一对傅氏变换关系。对于平稳随机过程,也有类似的关系,即 P()=其傅里叶反变换为第27页,本讲稿共85页于是 R0 因为R(0)表示随机过程的平均功率,它应等于功率谱密度曲线下的面积。因此,P()必然是平稳随机过程的功率谱密度函数。所以,平稳随机过程的功率谱密度P()与其自相关函数R()是一对傅里叶变换关系,即 第28页,本讲稿共85页于是 因为R(0)表示随机过程的平均功率,它应等于功率谱密度曲线下的面积。因此,P()必然是平稳随机过程的功率谱密度函数。所以,平稳随机过程的功率谱密度P()与其自相关函数R()
19、是一对傅里叶变换关系,即第29页,本讲稿共85页或简记为 R()P()关系式(2.2-18)称为维纳-辛钦关系,在平稳随机过程的理论和应用中是一个非常重要的工具。它是联系频域和时域两种分析方法的基本关系式。根据上述关系式及自相关函数R()的性质,不难推演功率谱密度P()有如下性质:第30页,本讲稿共85页 (1)P()0,非负性;(2.2-20)(2)P(-)=P(),偶函数。(2.2-21)因此,可定义单边谱密度P()为 P1()=0W 0W 0 例 2-1某随机相位余弦波(t)=Acos(ct+),其中A和c均为常数,是在(0,2)内均匀分布的随机变量。(1)求(t)的自相关函数与功率谱密
20、度;(2)讨论(t)是否具有各态历经性。第31页,本讲稿共85页 解 (1)先考察(t)是否广义平稳。(t)的数学期望为(t)的自相关函数为第32页,本讲稿共85页见(t)的数学期望为常数,而自相关函数只与时间间隔有关,所以(t)为广义平稳随机过程。根据平稳随机过程的相关函数与功率谱密度是一对傅里叶变换,即R()P(),则因为第33页,本讲稿共85页 cosc (-c)+(+c)所以,功率谱密度为 P()=(-c)+(+c)平均功率为 S=R(0)=(2)现在来求(t)的时间平均。根据式(2.2-6)可得第34页,本讲稿共85页 比较统计平均与时间平均,得a=,R()=,因此,随机相位余弦波是
21、各态历经的。第35页,本讲稿共85页2.3高斯随机过程高斯随机过程 2.3.1定义定义 若随机过程(t)的任意n维(n=1,2,)分布都是正态分布,则称它为高斯随机过程或正态过程。其n维正态概率密度函数表示如下:fn(x1,x2,xn;t1,t2,tn)式中,ak=E(tk),2k=E(tk)-ak2,|B|为归一化协方差矩阵的行列式,即第36页,本讲稿共85页1 b12 b1nB21 1 b2nBn1 bn2 1|B|jk为行列式|B|中元素bjk的代数余因子,bjk为归一化协方差函数,且第37页,本讲稿共85页 2.3.2重要性质重要性质 (1)由式(2.3-1)可以看出,高斯过程的n维分
22、布完全由n个随机变量的数学期望、方差和两两之间的归一化协方差函数所决定。因此,对于高斯过程,只要研究它的数字特征就可以了。(2)如果高斯过程是广义平稳的,则它的均值与时间无关,协方差函数只与时间间隔有关,而与时间起点无关,由性质(1)知,它的n维分布与时间起点无关。所以,广义平稳的高斯过程也是狭义平稳的。(3)如果高斯过程在不同时刻的取值是不相关的,即对所有jk有bjk=0,这时式(2.3-1)变为第38页,本讲稿共85页fn(x1,x2,xn;t1,t2,tn)=(2.3-2)也就是说,如果高斯过程在不同时刻的取值是不相关的,那么它们也是统计独立的以后分析问题时,会经常用到高斯过程中的一维分
23、布。例如,高斯过程在任一时刻上的样值是一个一维高斯随机变量,其一维概率密度函数可表示为=f(x1,t1)f(x2,t2)f(xn,tn)第39页,本讲稿共85页 式中,a为高斯随机变量的数学期望,2为方差。f(x)曲线如图 2-3所示。由式(2.3-3)和图2-3可知f(x)具有如下特性:(1)f(x)对称于x=a这条直线。(2)且有第40页,本讲稿共85页图2-3 正态分布的概率第41页,本讲稿共85页 3)a表示分布中心,表示集中程度,f(x)图形将随着的减小而变高和变窄。当a=0,=1时,称f(x)为标准正态分布的密度函数。当我们需要求高斯随机变量小于或等于任意取值x的概率P(x)时,还
24、要用到正态分布函数。正态分布函数是概率密度函数的积分,即 这个积分无法用闭合形式计算,我们要设法把这个积分式和可以在数学手册上查出积分值的特殊函数联系起来,一般常用以下几种特殊函数:第42页,本讲稿共85页(1)误差函数和互补误差函数。误差函数的定义式为 它是自变量的递增函数,erf(0)=0,erf()=1,且erf(-x)=-erf(x)。我们称1-erf(x)为互补误差函数,记为erfc(x),即 erfc(x)=1-erf(x)=它是自变量的递减函数,erfc(0)=1,erfc()=0,且erfc(-x)=2-erfc(x)。当x1时(实际应用中只要x2)即可近似有第43页,本讲稿共
25、85页 (2)概率积分函数和Q函数。概率积分函数定义为(x)=(2.3-10)这是另一个在数学手册上有数值和曲线的特殊函数,有()=1。Q函数是一种经常用于表示高斯尾部曲线下的面积的函数,其定义为 比较式(2.3-8)与式(2.3-10)和式(2.3-11),可得 第44页,本讲稿共85页 现在让我们把以上特殊函数与式(2.3-6)进行联系,以表示正态分布函数F(x)。若对式(2.3-6)进行变量代换,令新积分变量t=(z-a)/,就有dz=dt,再与式(2.3-10)联系,则有 F(x)=(2.3-15)若对式(2.3-6)进行变量代换,令新积分变量t=(z-a)/,就有dz=dt,再利用式
26、(2.3-5),则不难得到 第45页,本讲稿共85页 用误差函数或互补误差函数表示F(x)的好处是,它简明的特性有助于今后分析通信系统的抗噪声性能。F(X)=第46页,本讲稿共85页 2.3.3高斯白噪声高斯白噪声 信号在信道中传输时,常会遇到这样一类噪声,它的功率谱密度均匀分布在整个频率范围内,即 P()=(2.3-17)这种噪声被称为白噪声,它是一个理想的宽带随机过程。式中n0为一常数,单位是瓦/赫。显然,白噪声的自相关函数可借助于下式求得,即 R()=第47页,本讲稿共85页 这说明,白噪声只有在=0时才相关,而它在任意两个时刻上的随机变量都是互不相关的。图 2-4 画出了白噪声的功率谱
27、和自相关函数的图形。如果白噪声又是高斯分布的,我们就称之为高斯白噪声。由式(2.3-18)可以看出,高斯白噪声在任意两个不同时刻上的取值之间,不仅是互不相关的,而且还是统计独立的。应当指出,我们所定义的这种理想化的白噪声在实际中是不存在的。但是,如果噪声的功率谱均匀分布的频率范围远远大于通信系统的工作频带,我们就可以把它视为白噪声。第 3章将要讨论的热噪声和散弹噪声就是近似白噪声的例子。第48页,本讲稿共85页2.4随机过程通过线性系统随机过程通过线性系统 通信的目的在于传输信号,信号和系统总是联系在一起的。通信系统中的信号或噪声一般都是随机的,因此在以后的讨论中我们必然会遇到这样的问题:随机
28、过程通过系统(或网络)后,输出过程将是什么样的过程?这里,我们只考虑平稳过程通过线性时不变系统的情况。随机信号通过线性系统的分析,完全是建立在确知信号通过线性系统的分析原理的基础之上的。我们知道,线性系统的响应vo(t)等于输入信号vi(t)与系统的单位冲激响应h(t)的卷积,即vo(t)=vi(t)*h(t)=第49页,本讲稿共85页 若 vo(t)Vo(),vi(t)Vi(),h(t)H(),则有 Vo()=H()Vi()(2.4-2)若线性系统是物理可实现的,则 vo(t)=或 如果把vi(t)看作是输入随机过程的一个样本,则vo(t)可看作是输出随机过程的一个样本。显然,输入过程i(t
29、)的每个样本与输出过程o(t)的相应样本之间都满足式(2.4-4)的关系。这样,就整个过程而言,便有第50页,本讲稿共85页 o(t)=(2.4-5)假定输入i(t)是平稳随机过程,现在来分析系统的输出过程o(t)的统计特性。我们先确定输出过程的数学期望、自相关函数及功率谱密度,然后讨论输出过程的概率分布问题。1.输出过程o(t)的数学期望对式(2.4-5)两边取统计平均,有 Eo(t)=eh()i(t-)d=式中利用了平稳性假设Ei(t-)=Ei(t)=a(常数)。又因为第51页,本讲稿共85页 H(W)=求得H(0)=所以Eo(t)=aH(0)由此可见,输出过程的数学期望等于输入过程的数学
30、期望与直流传递函数H(0)的乘积,且Eo(t)与t无关。2.输出过程输出过程o(t)的自相关函数的自相关函数 第52页,本讲稿共85页Ro(t1,t1+)=Eo(t1)o(t1+)=E 根据平稳性 Ei(t1-)i(t1+-)=Ri(+-)有Ro(t1,t1+)=h()h()Ri(+-)dd=Ro 可见,o(t)的自相关函数只依赖时间间隔而与时间起点t1无关。由以上输出过程的数学期望和自相关函数证明,若线性系统的输入过程是平稳的,那么输出过程也是平稳的。第53页,本讲稿共85页 3.输出过程输出过程o(t)的功率谱密度的功率谱密度 对式(2.4-7)进行傅里叶变换,有令则有Po()=即 Po(
31、)=第54页,本讲稿共85页 可见,系统输出功率谱密度是输入功率谱密度Pi()与系统功率传输函数|H()|2的乘积。这是十分有用的一个重要公式。当我们想得到输出过程的自相关函数Ro()时,比较简单的方法是先计算出功率谱密度Po(),然后求其反变换,这比直接计算Ro()要简便得多。例 2 带带限限白白噪噪声声。试求功率谱密度为n0/2的白噪声通过理想矩形的低通滤波器后的功率谱密度、自相关函数和噪声平均功率。理想低通的传输特性为H()=K0e-jwt 0 其他第55页,本讲稿共85页 解 由上式得|H()|2=,|H。输出功率谱密度为Po()=|H()|2Pi()=,|H 可见,输出噪声的功率谱密
32、度在|H内是均匀的,在此范围外则为零,如图 2-5(a)所示,通常把这样的噪声称为带限白噪声。其自相关函数为第56页,本讲稿共85页图2-5 带限白噪声的功率谱和自相关函数第57页,本讲稿共85页 式 中,H=2fH。由 此 可 见,带 限 白 噪 声 只 有 在=k/2fH(k=1,2,3,)上得到的随机变量才不相关。它告诉我们,如果对带限白噪声按抽样定理抽样的话,则各抽样值是互不相关的随机变量。这是一个很重要的概念。如图 2-5(b)所示,带限白噪声的自相关函数Ro()在=0 处有最大值,这就是带限白噪声的平均功率:Ro(0)=n0fH 4.输出过程输出过程o(t)的概率分布的概率分布 从
33、原理上看,在已知输入过程分布的情况下,通过式(2.4-5),即第58页,本讲稿共85页 总可以确定输出过程的分布。其中一个十分有用的情形是:如果线性系统的输入过程是高斯型的,则系统的输出过程也是高斯型的。因为从积分原理来看,上式可表示为一个和式的极限,即第59页,本讲稿共85页 由于i(t)已假设是高斯型的,所以,在任一时刻的每项i(t-k)h(k)k都是一个高斯随机变量。因此,输出过程在任一时刻得到的每一随机变量,都是无限多个高斯随机变量之和。由概率论得知,这个“和”的随机变量也是高斯随机变量。这就证明,高斯过程经过线性系统后其输出过程 仍为高斯过程。更一般地说,高斯过程经线性变换后的过程仍
34、为高斯过程。但要注意,由于线性系统的介入,与输入高斯过程相比,输出过程的数字特征已经改变了。第60页,本讲稿共85页2.5窄带随机过程窄带随机过程 随机过程通过以fc为中心频率的窄带系统的输出,即是窄带过程。所谓窄带系统,是指其通带宽度ffc,且fc远离零频率的系统。实际中,大多数通信系统都是窄带型的,通过窄带系统的信号或噪声必是窄带的,如果这时的信号或噪声又是随机的,则称它们为窄带随机过程。如用示波器观察一个实现的波形,则如图2-6(b)所示,它是一个频率近似为fc,包络和相位随机缓变的正弦波。第61页,本讲稿共85页图2-6 窄带过程的频谱和波形示意第62页,本讲稿共85页 因此,窄带随机
35、过程(t)可用下式表示:(t)=a(t)cosct+(t),a(t)0(2.5-1)等价式为(t)=c(t)cosct-s(t)sinct (2.5-2)其中c(t)=a(t)cos(t)(2.5-3)s(t)=a(t)sin(t)(2.5-4)式中,a(t)及(t)分别是(t)的随机包络和随机相位,c(t)及s(t)分别称为(t)的同相分第63页,本讲稿共85页 量和正交分量,它们也是随机过程,显然它们的变化相对于载波cosct的变化要缓慢得多。由式(2.5-1)至(2.5-4)看出,(t)的统计特性可由a(t),(t)或c(t),s(t)的统计特性确定。反之,如果已知(t)的统计特性则可确
36、定a(t),(t)以及c(t),s(t)的统计特性。第64页,本讲稿共85页 2.5.1同相和正交分量的统计特性同相和正交分量的统计特性 设窄带过程(t)是平稳高斯窄带过程,且均值为零,方差为2。下面将证明它的同相分量c(t)和正交分量s(t)也是零均值的平稳高斯过程,而且与(t)具有相同的方差。1.数学期望数学期望 对式(2.5-2)求数学期望:E(t)=Ec(t)cosct-Es(t)sinct (2.5-5)可得第65页,本讲稿共85页 Ec(t)=0 Es(t)=0 (2.5-6)2.自相关函数自相关函数 R(t,t+)=E(t)(t+)=Ec(t)cosct-s(t)sinct c(
37、t+)cosc(t+)-s(t+)sinc(t+)=Rc(t,t+)cosct cosc(t+)-Rcs(t,t+)cosctsinc(t+)-Rsc(t,t+)sinctcosc(t+)+Rs(t,t+)sinctsinc(t+)第66页,本讲稿共85页 式中Rc(t,t+)=Ec(t)c(t+)Rcs(t,t+)=Ec(t)s(t+)Rsc(t,t+)=Es(t)c(t+)Rs(t,t+)=Es(t)s(t+)因为(t)是平稳的,故有 R(t,t+)=R()这就要求式(2.5-7)的右边也应该与t无关,而仅与时间间隔有关。若取使sinct=0 的所有t值,则式(2.5-7)应变为 第67页
38、,本讲稿共85页 R()=Rc(t,t+)cosc-Rcs(t,t+)sinc(2.5-8)这时,显然应有 Rc(t,t+)=Rc()Rcs(t,t+)=Rcs()所以,式(2.5-8)变为 R()=Rc()cosc-Rcs()sinc (2.5-9)再取使cosct=0的所有t值,同理有 R()=Rs()cosc+Rsc()sinc (2.5-10)第68页,本讲稿共85页 其中应有 Rs(t,t+)=Rs()Rsc(t,t+)=Rsc()由以上的数学期望和自相关函数分析可知,如果窄带过程(t)是平稳的,则c(t)与s(t)也必将是平稳的。进一步分析,式(2.5-9)和式(2.5-10)应同
39、时成立,故有 Rc()=Rs()(2.5-11)Rcs()=-Rsc()(2.5-12)可见,同相分量c(t)和正交分量s(t)具有相同的自相关函数,而且根据互相关函数的性质,应有第69页,本讲稿共85页 Rcs()=Rsc(-)将上式代入式(2.5-12),可得 Rsc()=-Rsc(-)(2.5-13)同理可推得Rcs()=-Rcs(-)(2.5-14)式(2.5-13)、(2.5-14)说明,c(t)、s(t)的互相关函数Rsc()、Rcs()都是的奇函数,在=0时 Rsc(0)=Rcs(0)=0 (2.5-15)于是,由式(2.5-9)及式(2.5-10)得到第70页,本讲稿共85页
40、Rsc(0)=Rcs(0)=0 (2.5-15)于是,由式(2.5-9)及式(2.5-10)得到 R(0)=Rc(0)=Rs(0)(2.5-16)即2=2c=2s (2.5-17)这表明(t)、c(t)和s(t)具有相同的平均功率或方差(因为均值为0)。另外,因为(t)是平稳的,所以(t)在任意时刻的取值都是服从高斯分布的随机变量,故在式(2.5-2)中有第71页,本讲稿共85页 取t=t1=0 时,(t1)=c(t1)取t=t2=32c时,(t2)=s(t2)所以c(t1),s(t2)也是高斯随机变量,从而c(t)、s(t)也是高斯随机过程。又根据式(2.5-15)可知,c(t)、s(t)在
41、同一时刻的取值是互不相关的随机变量,因而它们还是统计独立的。上所述,我们得到一个重要结论:一个均值为零的窄带平稳高斯过程(t),它的同相分量c(t)和正交分量s(t)也是平稳高斯过程,而且均值都为零,方差也相同。此外,在同一时刻上得到的c和s是互不相关的或统计独立的。第72页,本讲稿共85页 2.5.2包络和相位的统计特性包络和相位的统计特性 由上面的分析可知,c和s的联合概率密度函数为 f(c,s)=f(c)f(s)=设a,的联合概率密度函数为f(a,),则利用概率论知识,有 f(a,)=f(c,s)根据式(2.5-3)和式(2.5-4)在t时刻随机变量之间的关系 c=acos s=asin
42、第73页,本讲稿共85页 得到 Cos sin-asin acos=于是f(a,)=af(c,s)=注意,这里a0,而在(0,2)内取值。再利用概率论中边际分布知识将f(a,)对积分,可求得包络a的一维概率密度函数为第74页,本讲稿共85页 可见,a服从瑞利分布。同理,f(a,)对a积分可求得相位的一维概率密度函数为f()=可见,服从均匀分布。第75页,本讲稿共85页 综上所述,我们又得到一个重要结论:一个均值为零,方差为2的窄带平稳高斯过程(t),其包络a(t)的一维分布是瑞利分布,相位(t)的一维分布是均匀分布,并且就一维分布而言,a(t)与(t)是统计独立的,即有下式成立:f(a,)=f
43、(a)f()(2.5-23)第76页,本讲稿共85页2.6正弦波加窄带高斯噪声正弦波加窄带高斯噪声 信号经过信道传输后总会受到噪声的干扰,为了减少噪声的影响,通常在接收机前端设置一个带通滤波器,以滤除信号频带以外的噪声。因此,带通滤波器的输出是信号与窄带噪声的混合波形。最常见的是正弦波加窄带高斯噪声的合成波,这是通信系统中常会遇到的一种情况,所以有必要了解合成信号的包络和相位的统计特性。设合成信号为 r(t)=A cos(ct+)+n(t)(2.6-1)第77页,本讲稿共85页 式中,n(t)=nc(t)cosct-ns(t)sinct为窄带高斯噪声,其均值为零,方差为2n;正弦信号的A,c均
44、为常数,是在(0,2)上均匀分布的随机变量。于是 r(t)=Acos+nc(t)cosct-Asin+ns(t)sinct =zc(t)cosct-zs(t)sinct =z(t)cosct+(t)(2.6-2)式中zc(t)=A cos+nc(t)(2.6-3)zs(t)=Asin+ns(t)(2.6-4)第78页,本讲稿共85页合成信号r(t)的包络和相位为z(t)=利用上一节的结果,如果值已给定,则zc、zs是相互独立的高斯随机变量,且有 Ezc=Acos Ezs=Asin第79页,本讲稿共85页 所以,在给定相位的条件下的zc和zs的联合概率密度函数为 f(zc,zs/)=利用上一节相
45、似的方法,根据式(2.6-3)、(2.6-4)可以求得在给定相位的条件下的z和的联合概率密度函数为f(z,/)=f(zc,zs/)=zf(zc,zs/)第80页,本讲稿共85页求条件边际分布,有由于故有 第81页,本讲稿共85页 式中,I0(x)为零阶修正贝塞尔函数。当x0时,I0(x)是单调上升函数,且有I0(0)=1。因此 f(z/)=由上式可见,f(z/)与无关,故正弦波加窄带高斯过程的包络概率密度函数为 这个概率密度函数称为广义瑞利分布,也称莱斯(Rice)密度函数。上式存在两种极限情况:第82页,本讲稿共85页 (1)当信号很小,A0,即信号功率与噪声功率之比 =r0时,x值很小,有
46、I0(x)=1,这时合成波r(t)中只存在窄带高斯噪声,式(2.6-8)近似为式(2.5-21),即由莱斯分布退化为瑞利分布。(2)当信噪比r很大时,有I0(x),这时在zA附近,f(z)近似于高斯分布,即 f(z)由此可见,信号加噪声的合成波包络分布与信噪比有关。小信噪比时,它接近于瑞利分布;大信噪比时,它接近于高斯分布;在一般情况下它是莱斯分布。图 2-7(a)给出了不同的r值时f(z)的曲线。第83页,本讲稿共85页 关于信号加噪声的合成波相位分布f(),由于比较复杂,这里就不再演算了。不难推想,f()也与信噪比有关。小信噪比时,f()接近于均匀分布,它反映这时窄带高斯噪声为主的情况;大信噪比时,f()主要集中在有用信号相位附近。图 2-7(b)给出了不同的r值时f()的曲线。第84页,本讲稿共85页图 2 7 正弦波加窄带高斯过程的包络与相位分布 第85页,本讲稿共85页