第二章 平稳随机过程的谱分析
- 格式:doc
- 大小:1003.34 KB
- 文档页数:37
第二章 随机过程分析1.1 学习指导 1.1.1 要点随机过程分析的要点主要包括随机过程的概念、分布函数、概率密度函数、数字特征、通信系统中常见的几种重要随机过程的统计特性。
1. 随机过程的概念随机过程是一类随时间作随机变化的过程,它不能用确切的时间函数描述。
可从两种不同角度理解:对应不同随机试验结果的时间过程的集合,随机过程是随机变量概念的延伸。
2. 随机过程的分布函数和概率密度函数如果ξ(t )是一个随机过程,则其在时刻t 1取值ξ(t 1)是一个随机变量。
ξ(t 1)小于或等于某一数值x 1的概率为P [ ξ(t 1) ≤x 1],随机过程ξ(t )的一维分布函数为F 1(x 1, t 1) = P [ξ(t 1) ≤x 1] (2-1)如果F 1(x 1, t 1)的偏导数存在,则ξ(t )的一维概率密度函数为1111111(,)(, ) (2 - 2)∂=∂F x t f x t x对于任意时刻t 1和t 2,把ξ(t 1) ≤x 1和ξ(t 2) ≤x 2同时成立的概率{}212121122(, ; , )(), () (2 - 3)F x x t t P t x t x ξξ=≤≤称为随机过程ξ(t )的二维分布函数。
如果2212122121212(,;,)(,;,) (2 - 4)F x x t t f x x t t x x ∂=∂⋅∂存在,则称f 2(x 1, x 2; t 1, t 2)为随机过程ξ(t )的二维概率密度函数。
对于任意时刻t 1,t 2,…,t n ,把{}n 12n 12n 1122n n ()(),(),,() (2 - 5)=≤≤≤F x x x t t t P t x t x t x ξξξ,,,;,,,称为随机过程ξ(t )的n 维分布函数。
如果n n 12n 12n n 12n 12n 12n(x )() (2 - 6)∂=∂∂∂F x x t t t f x x x t t t x x x ,,,;,,,,,,;,,,存在,则称f n (x 1, x 2, …, x n ; t 1, t 2, …, t n )为随机过程ξ(t )的n 维概率密度函数。
随机信号分析实验报告实验一:平稳随机过程的数字特征实验二:平稳随机过程的谱分析实验三:随机信号通过线性系统的分析实验四:平稳时间序列模型预测班级:姓名:学号:一、实验目的1、加深理解平稳随机过程数字特征的概念2、掌握平稳随机序列期望、自相关序列的求解3、分析平稳随机过程数字特征的特点二、实验原理平稳随机过程数字特征求解的相关原理三、实验过程function y = experimentnumber = 49; %学号49I = 8; %幅值为8u = 1/number;Ex = I*0.5 + (-I)*0.5;N = 64;C0 = 1; %计数p(1) = exp(-u);for m = 2:Nk = 1:m/2;p(m) = exp(-u*m) + sum((u*m).^(2*k)./factorial(2*k)*exp(-u*m));2222()[()()]{()()}{()()}X R m E X n X n m I P X n X n m I I P X n X n m I =+=+=-+=-E[X(n)]= I P{X(n)=+I}+(-I)P{X(n)=-I}=0⨯⨯0m >当时,/2220(){()()}(2)!m k m k m P X n X n m I e P k λλ⎢⎥⎣⎦-=+===∑222()(1)(21)X R m I P I P I P =--=-2()()X X X C m R m m =-me I m n X n X E m R λ22)]()([)(-=+=end;pp = [fliplr(p) C0 p];Rx = (2*pp - 1)*I^2;m = -N:N;Kx = Rx - Ex*Ex;rx = Kx/25;subplot(211), plot(m,Rx); axis([-N N 0 I*I]); title('自相关序列');subplot(212), plot(m,rx); axis([-N N 0 1]); title('自相关序数');四、实验结果及分析自相关序列的特点分析:m>0时Rx(m)随着m的增大而减小,m<0时Rx(m)随着m的增大而增大。
第二章平稳随机过程的谱分析平稳随机过程第二章平稳随机过程的谱分析本章要解决的问题:●随机信号是否也可以应用频域分析方法?●傅里叶变换能否应用于随机信号?● 相关函数与功率谱的关系● 功率谱的应用● 采样定理● 白噪声的定义2.1 随机过程的谱分析2.1.1 预备知识1、付氏变换:对于一个确定性时间脉冲x(t),设x(t)是时间t 的非周期实函数,且x(t) 满足狄利赫利条件(有限个极值,有限个断点,断点为有限值)且绝对可积,能量有限,则x(t)傅里叶变换存在。
即:满足上述三个条件的x(t)的傅里叶变换为:其反变换为:2、帕赛瓦等式由上面式子可以重新得到:——称为非周期性三十天拉热函数的帕塞瓦(Parseval)等式。
物理意义:若x(t)表示的是电压(或电流) ,则上式左边代表x(t)在时间(-∞, ∞) 区间的总能量(单位阻抗)。
因此,等式右边的被积函数X X (ω)2表示了信号x(t)能量按频率分布的情况,故称X X (ω)2为能量谱密度。
2.1.2、随机过程的功率谱密度变换一个信号的惟教变换是否存在,可能需要满足三个条件,那么随机信号是否满足这三个条件从而存在付氏呢?随机信号持续时间无限长,因此,对于非0的样本函数,它的能量一般也是无限的,因此,其付氏变换不牵涉到。
但是注意到它的平均功率是有限的,在特定的条件下,仍然洪可以利用博里叶变换这一工具。
为了将傅里叶变换方法常量应用于随机过程,必须对过程的待测函数做某些限制,最简单的一种方法是应用截取函数。
截取函数x T (t):图2.1 x (t)及其截取函数当x(t)为有限值时,裁取函数x T (t)满足绝对可积条件。
因此,x T (t)的傅里叶变换存在,有很明显,式的变化)x T (t)也应满足帕塞瓦等式,即:(注意积分区间和表达用2T 除上式等号用的两端,可以得到等号于两边取集合平均,可以得到:令T→∞,再取极限,便可得到随机过程的平均功率。
实验二平稳随机过程的谱分析谱分析是对平稳随机过程的频率特性进行研究的一种方法。
它通过分析随机过程在不同频率下的能量分布,可以揭示出随机过程的主要频率成分和其相应的能量。
在实验二中,我们将以一个平稳随机过程为例,详细介绍谱分析的方法和步骤,并通过具体的实例来说明如何进行谱分析。
首先,我们需要明确谱密度函数的概念。
谱密度函数描述了随机过程在各个频率上的能量分布,其定义为随机过程在单位频率范围内的功率谱与单位频率之比。
一般地,谱密度函数可以通过傅里叶变换和自相关函数计算得到。
接下来,我们需要计算随机过程的自相关函数。
自相关函数反映了随机过程在不同时刻之间的相关性,其定义为随机过程在不同时刻的取值之积的期望。
通过计算自相关函数,我们可以得到随机过程的自相关系数和自相关函数的性质。
然后,我们可以通过自相关函数计算随机过程的功率谱密度函数。
功率谱密度函数描述了随机过程在各个频率上的能量分布,其定义为自相关函数的傅里叶变换。
通过计算功率谱密度函数,我们可以得到随机过程的频谱特性。
在进行谱分析时,我们需要选择适当的算法和工具进行计算。
常见的算法包括周期图法、Welch法和傅里叶变换法。
周期图法是一种通过周期图对随机过程进行频谱分析的方法,其步骤包括选择窗函数、计算周期图和计算功率谱密度函数。
Welch法是一种通过分段计算随机过程的频谱的方法,其步骤包括选择窗函数、选择段数、计算每一段的频谱并对它们求平均。
傅里叶变换法是一种通过对随机过程进行傅里叶变换得到频谱的方法,其步骤包括对随机过程进行傅里叶变换和计算功率谱密度函数。
最后,我们可以通过绘制频谱图来直观地表示随机过程的频谱特性。
频谱图是将频率作为横坐标、功率谱密度函数的取值作为纵坐标,以直方图或曲线的形式展示出来。
通过观察频谱图,我们可以得到随机过程的主要频率成分和其相应的能量。
综上所述,谱分析是一种揭示平稳随机过程频率特性的重要方法。
通过计算自相关函数和功率谱密度函数,并绘制频谱图,可以得到平稳随机过程的主要频率成分和其相应的能量,进而对随机过程进行频域分析。
第二章 平稳过程的谱分析§1谱理论简介我们知道,由Wold 分解定理,一个平稳过程t Y 可以找到一个平稳的(,)ARM A p q 来近似。
且已知1,,T y y ,当T →∞,我们可以一致的估计(,)ARM A p q 模型中的未知参数,并由此来把握平稳过程t Y 。
现在,我们换一个角度看t Y ,把所有二阶矩平稳过程看成为一个Hilbert 空间,那么,由Hilbert 空间的谱表示定理,任何一个二阶矩平稳过程t Y 都可以表示成为一个右连续的正交增量过程的R —S 积分,即,()i tt Y edz πωπω-=⎰,()()()z A iB ωωω=+。
满足:[()()]0i j E dA dA ωω=, [()()]0i j E dB dB ωω=,i j ∀≠。
(正交增量性)[()][()]0E dA E dB ωω==, [()][()]Var dA Var dB ωω=,且右连续是指均方收敛,即,2[()()]0E A A ωδω+-→,0δ↓。
( 参见MIT 教本)将t Y 改写成,0cos()()sin()()t Y t dA t dB ππωωωω=+⎰⎰。
定义[()][()]Var dA Var dB ωω==2()dF ω,[0,]ωπ∀∈。
那么由(),()A B ωω的正交增量性和右连续性,知()F ω是一个[0,]π上的非减右连续的函数。
称()F ω为t Y 的谱分布函数。
又将()dF ω写成,()()dF f d ωωω=,则()f ω就称为t Y 的谱密度函数。
注意,()F ω或()f ω是由(),()A B ωω唯一决定的,也就是由t Y 唯一决定的。
这里唯一性指的是几乎处处唯一。
反过来也正确。
任给一个谱密度函数()f ω或谱分布函数()F ω,可以决定一个唯一的右连续的正交增量过程,()()()z A iB ωωω=+,并由()z ω决定一个唯一的平稳过程t Y 。
实验二平稳随机过程的谱分析一、实验目的1、复习信号处理的采样定理2、理解功率谱密度函数与自相关函数的关系3、掌握对功率谱密度函数的求解和分析二、实验设备计算机、Matlab软件三、实验内容与步骤已知平稳随机过程的相关函数为:RX(τ)=1-|τ|/T |τ|<T=0 |τ|>=TT=学号*3设计程序求:1.利用采样定理求R1(m)2.利用RX(τ)求SX(w),3.利用功率谱密度采样定理求S(w)(离散时间序列的功率谱密度)4.利用IFFT求R(m)5.利用求出的R1(m),用FFT求S1(w)6.比较上述结果。
四、实验原理平稳随机过程的谱分析和付立叶变换1、TTTTSadjTRFTSTXX222)(sin4)2()exp()/1(2)}({)(ωωωτωτττω==--==⎰2、如果时间信号的采样间隔为T0,那么在频谱上的采样间隔1/(N*T0),保持时域和频域的采样点一致N3、注意实际信号以原点对称,画图时是以中心对称,注意坐标的变换五、实验报告要求1、打印所求出的R1(m)、R(m)、S1(w)、S(w)序列,并绘图。
采样点数根据采样定理求出,并在程序中设置为可任意键盘输入的值,以便了解采样点数变化和由采样所得序列能否正确恢复原始信号的关系。
2、附上程序和必要的注解。
六、实验过程function y = experiment2close all;clc;number = 41;T = number*3;T0 = 0.1 %input('采样间隔T0=');t = -T: T0: T;t1 = -2*T: T0: 2*T;n = T/T0;Rx1 = 1 - abs(t)/T;Rx = [zeros(1, n) Rx1 zeros(1, n)];figure(1),subplot(211), plot(t1, Rx);title('自相关函数') ; %自相关函数F = 1/(2*T0);F0 = 1/(4*T);f = -F: F0: F;w = 2* pi* f;a = w*T/2;Sx = T*sin(a).*sin(a)./(a.*a);Sx(2*n + 1) = T;subplot(212), plot(f, Sx); title('功率谱密度函数') ; %功率谱密度函数figure(2),R1 = Rx;subplot(211),plot(R1); title('自相关序列') ; %自相关序列S1 = T0*abs(fft(R1));S1 = fftshift(S1);subplot(212), plot(S1); title('自相关序列FFT得到功率谱密度函数') ; %自相关序列FFT得到功率谱密度函数figure(3),S = Sx;subplot(211), plot(S); title('功率谱密度函数采样序列') % 功率谱密度函数采样序列R = 1/T0*abs(ifft(S));R = ifftshift(R);subplot(212), plot(R); title('功率谱密度序列IFFT得到自相关序列') %功率谱密度序列IFFT得到自相关序列七、实验结果及分析-150-100-5005010015000.51自相关函数-5-4-3-2-1012345020406080功率谱密度函数05001000150020002500300000.51自相关序列050010001500200025003000020406080自相关序列FFT 得到功率谱密度函数050010001500200025003000020406080功率谱密度函数采样序列05001000150020002500300000.51功率谱密度序列IFFT 得到自相关序列八、实验心得体会通过本次对平稳随机过程的谱分析的实验,进一步熟悉了Matlab 软件的使用操作,加深了书本上的理论知识,如信号处理的采样定理的理解,掌握了功率谱密度函数与自相关函数的关系,以及对功率谱密度函数的求解和分析方法。
第二章平稳随机过程的谱分析本章要解决的问题:●随机信号是否也可以应用频域分析方法?●傅里叶变换能否应用于随机信号?●相关函数与功率谱的关系●功率谱的应用●采样定理●白噪声的定义2.1 随机过程的谱分析2.1.1 预备知识1、付氏变换:对于一个确定性时间信号x(t),设x(t)是时间t的非周期实函数,且x(t) 满足狄利赫利条件(有限个极值,有限个断点,断点为有限值)且绝对可积,能量有限,则x(t)傅里叶变换存在。
即:满足上述三个条件的x(t)的傅里叶变换为:其反变换为:2、帕赛瓦等式由上面式子可以得到:——称为非周期性时间函数的帕塞瓦(Parseval)等式。
物理意义:若x(t)表示的是电压(或电流),则上式左边代表x(t)在时间(-∞,∞)区间的总能量(单位阻抗)。
因此,等式右边的被积函数2)(ωXX表示了信号x(t)能量按频率分布的情况,故称2)(ωXX为能量谱密度。
2.1.2、随机过程的功率谱密度一个信号的付氏变换是否存在,需要满足三个条件,那么随机信号是否满足这三个条件从而存在付氏变换呢?随机信号持续时间无限长,因此,对于非0的样本函数,它的能量一般也是无限的,因此,其付氏变换不存在。
但是注意到它的平均功率是有限的,在特定的条件下,仍然可以利用博里叶变换这一工具。
为了将傅里叶变换方法应用于随机过程,必须对过程的样本函数做某些限制,最简单的一种方法是应用截取函数。
x(t):截取函数T图2.1 x(t)及其截取函数x(t)满足绝对可积条件。
因此,当x(t)为有限值时,裁取函数Tx(t)的傅里叶变换存在,有Tx(t)也应满足帕塞瓦等式,即:(注意积分区间和表达很明显,T式的变化)用2T 除上式等号的两端,可以得到等号两边取集合平均,可以得到:令∞→T ,再取极限,便可得到随机过程的平均功率。
交换求数学期望和积分的次序,可以得到:(注意这里由一条样本函数推广到更一般的随机过程,即下面式子对所有的样本函数均适用)ωωπd TT X E dt t XE TX T TT T 2]),([lim 21)]([21lim22⎰⎰∞∞-∞→-∞→=上式等号的左边表示的正是随机过程消耗在单位电阻上的平均功率(包含时间平均和统计平均),以后我们将简称它为随机过程的功率并记为Q 。
再看等式的右边,它当然也存在,并且等于Q 。
又因为2),(ωT XX非负,所以极限TT XE XT 2]),([lim2ω∞→必定存在,记为)(ωX S :ωωπd S dt t XE TQ X TT T ⎰⎰∞∞--∞→==)(21)]([21lim2注意:(1)Q 为确定性值,不是随机变量(2))(ωX S 为确定性实函数。
(见式)● 两个结论: 1.><=)]([2t X E A Q 式中,><>=<∞→.21lim.TA T表示时间平均。
它说明,随机过程的平均功率可以通过对过程的均方值求时间平均来得到,即对于一般的随机过程(例如,非平稳随机过程)求平均功率,需要既求时间平均,又求统计平均。
显然, Q 不是随机变量。
若随机过程为平稳的,则)0()]([)]([22X R t X E t X E A Q =>=<=这是因为均方值与时间t 无关,其时间平均为它自身。
由于已经对2),(ωT XX求了数学期望,所以)(ωX S 不再具有随机性,它是ω的确定性函数。
● 功率谱密度:)(ωX S 描述了随机过程X(t)的功率在各个不同频率上的分布——称)(ωX S 为随机过程X(t)的功率谱密度。
● 对)(ωX S 在X(t)的整个频率范围内积分,便可得到X(t)的功率。
● 对于平稳随机过程,则有:⎰∞∞-=ωωπd S t X E X )(21)]([22.1.3、功率谱密度的性质证明:证明:因为2),(ωT XX进行了取模运算,这是ω的实函数,所以)(ωX S 也是ω的实函数,且为确定性实函数。
证明:因此:即:得:证明:对于平稳随机过程,有:⎰∞∞-=ωωπd S t X E X )(21)]([22.2 联合平稳随机过程的互功率谱密度2.2.1、互谱密度可由单个随机过程的功率谱密度的概念,以及相应的分析方法推广而来。
考虑两个平稳实随机过程X(t)、Y(t), 它们的样本函数分别为)(t x 和)(t y ,定义两个截取函数()t x T 、()t y T 为:因为()t x T 、()t y T 都满足绝对可积的条件,所以它们的傅里叶变换存在。
在时间范围(-T ,T)内,两个随机过程的互功率)(T Q XY 为:(注意()t x T 、()t y T 为确定性函数,所以求平均功率只需取时间平均)由于()t x T 、()t y T 的傅里叶变换存在,故帕塞瓦定理对它们也适用,即:注意到上式中,)(t x 和)(t y 是任一样本函数,因此,具有随机性,取数学期望,并令∞→T ,得:])()(21[lim )]([lim dt t y t x TE Q T Q E TT T XY XY T ⎰-∞→∞→===]),(21[lim dt t t R TTT XYT ⎰-∞→=ωωωπd TT X T XE Y XT 2)],(),([lim 21*⎰∞∞-∞→定义互功率谱密度为:得:同理,有:又知以上定义了互功率和互功率谱密度,并且导出了它们之间的关系。
2.2.2、互谱密度和互相关函数的关系平稳随机过程的自相关函数与其功率谱密度之间互为傅里叶变换,互相关函数与互谱密度之间也存在着类似关系。
定义:对于两个实随机过程X(t)、Y(t),其互谱密度)(ωXY S 与互相关函数),(τ+t t R XY 之间的关系为若X(t)、Y(t)各自平稳且联合平稳,则有即:式中,><.A 表示时间平均。
显然:证明:略,参见自相关函数和功率谱密度关系的证明。
结论:对于两个联合平稳(至少是广义联合平稳)的实随机过程,它们的互谱密度与其互相关函数互为傅里叶变换。
2.3.3、互谱密度的性质互功率谱密度和功率谱密度不同,它不再是频率ω的正的、实的和偶函数。
性质1:)()()(*ωωωYX YX XY S S S =-= 证明:⎰∞∞--=ττωωτd eR S j XY XY )()(=⎰∞∞---ττωτd e R j YX )( 令ττ-==⎰∞∞-ττωτd eR j YX )(=)(*ωYXS=⎰∞∞---τττωd eR j YX )()(=)(ω-YX S性质2:)(Re[)](Re[ωω-=XY XY S S ;)(Re[)](Re[ωω-=YX YX S S证明:式中Re[·]表示实部。
亦即互谱密度的实部为ω的偶函数。
ττωωτd eR S j XY XY ⎰∞∞--=)()(=τωτωττd j R XY )]sin()[cos (⎰∞∞--+所以:τωττωd R S XYXY ⎰∞∞-=cos )()](Re[ 令ττ-==τωττd R XY⎰∞∞--cos )(=)](Re[ω-XY S其它同理可证。
性质3:证明:类似性质2证明。
性质4:若X(t)与Y(t)正交,则有证明:若X(t)与Y(t)正交,则0),(),(2121==t t R t t R YX XY 所以,0)()(==ωωYX XY S S性质5:若X(t)与Y(t)不相关,X(t)、Y(t)分别具有常数均值Xm 和Y m ,则证明:因为X(t)与Y(t)不相关,所以Y X m m t Y t X E =)]()([21ττωωτd eR S j XY XY ⎰∞∞--=)()(=τωτd em m j YX ⎰∞∞--=)(2ωδπY X m m (注意)(21ωπδ↔) 性质6:式中,A<∙>表示时间平均。
这给出了一般的随机过程(包含平稳)的互谱密度与互相关函数的关系式。
[例2.2] 设两个随机过程X(t)和Y(t)联合平稳,其互相关函数)(τXY R 为:求互谱密度)(ωXY S ,)(ωYX S 解:先求)(ωXY S :再求)(ωYX S2.3 功率谱密度与自相关函数之间的关系确定信号:x(t)↔ )(ωj X 。
随机信号:平稳随机过程的自相关函数↔功率谱密度。
1.定义:若随机过程X(t)是平稳的,自相关函数绝对可积,则自相关函数与功率谱密度构成一对付氏变换,即:这一关系就是著名的维纳—辛钦定理、或称为维纳—辛钦公式。
2. 证明:下面就来推导这一关系式。
证明方法类似式的证明。
因为:由(3.1.14)式 2[(,)]()lim2XX T E XT S Tωω→∞= *1lim[(,)(,)]2XXT E XT XT Tωω→∞==1lim2T T→∞121122[()()]T T j t j t TTE X t edt X t e dt ωω---⎰⎰交换积分和数学期望顺序=21()12121lim[()()]2T T j t t TTT E X t X t edt dt Tω----→∞⎰⎰=⎰⎰----∞→-T T TT t t j X T dt dt et t R T21)(1212)(21limω设12t t -=τ,12t t u +=,则22ut τ+=,21τ-=u t所以:2121212121),(),(21=-=∂∂=u t t J τ t1t2-TT2T2Tu -2T τ-=T u 2τ+-=T u 2τ+=T u 2τ--=T u 2τ图2.2则du e R d TS j X TT T T X ωτττττω--+-∞→⎰⎰=)(21{21lim)(2022 })(210222du eR d j X TT T ωτττττ--+--⎰⎰+=})(2121{lim 2222du eR d T j X TT T T T ωτττττ---+-∞→⎰⎰=τττωτd eR T Tj X TTT --∞→⎰-)()2(21lim22=τττωτd eR Tj X TTT --∞→⎰-)()21(lim22 (1)=⎰∞∞--ττωτd eR j X )(-22lim)()2T j X TT R ed Tωττττ--→∞⎰(注意T ∞→,02→Tτ;且∞→τ时,0)(→τX R 。
因此,通常情况下,第二项为0=⎰∞∞--ττωτd eR j X )(证毕。
推论:对于一般的随机过程X(t),有:则平均功率为:ωωπd S dt t XE TX TT T ⎰⎰∞∞--∞→=)(21)]([21lim2(0=τ)——时间平均加统计平均。
利用自相关函数和功率谱密度皆为偶函数的性质,又可将维纳—辛钦定理表示成:3.单边功率谱由于实平稳过程x(t)的自相关函数)(τX R 是实偶函数,功率谱密度也一定是实偶函数。