Alex_McAvoy

想要成为渔夫的猎手

平稳过程的谱密度

【相关函数的谱密度】

平稳过程的相关函数可视为一表示位移的时间函数,在时域上描述了随机过程的统计特征,因此,对于平稳过程的相关函数,利用 Fourier 分析的方法进行研究,便可在频域上描述平稳过程的统计特征,进而得到平稳过程谱密度这一概念

{X(t),<t<+} 是均方连续的平稳过程,则其相关函数

RX(τ)=12π+ejωπdFX(ω),<τ<+

被称为 {X(t),<t<+} 相关函数的谱展开式,或谱分解式

其中,FX(ω) 被称为 {X(t),<t<+}谱函数,其在 (,+) 上非负、有界、单调不减、右连续,且有

FX()=0FX(+)=2πRX(0)

若存在函数 SX(ω),使得

FX(ω)=ωSX(λ)dλ,<ω<+

则称 SX(ω){X(t),<t<+}谱密度

【功率谱密度】

谱密度的概念来自于无线电技术,在物理学中表示功率谱密度

x(t),<t<+ 为一确定性信号,如果 x(t) 满足 Dirichlet 条件且绝对可积,那么 x(t) 具有频谱

F(ω)=+ejωtx(t)dt

则在 x(t)F(ω) 中有 Parseval 等式成立

+x2(t)dt=12π+|Fx(ω)|2dω

等式左边表示 x(t)(,+) 上的总能量,相应地,|Fx(ω)|2 被称为 x(t) 的能谱密度

在实际应用中,很多信号的总能量是无限的,不满足绝对可积的条件,此时通常转为研究 x(t)(,+) 上的平均功率,即

limT+12TTTx2(t)dt

x(t) 的截尾函数

xT(t)={x(t),|t|T0,|t|>T

其在 (,+) 上绝对可积,此时 xT(t) 的 Fourier 变换为

Fx(ω,T)=+ejωtxT(t)dt=TTejωtx(t)dt

那么,Parseval 等式就变为

TTx2(t)dt=12π+|Fx(ω,T)|2dω

将 Parseval 等式两边同时除以 2T,再令 T+,即可得 x(t)(,+) 上的平均功率

limT+12TTTx2(t)dt=12π+limT+12T|Fx(ω,T)|2dω

相应地,能谱密度为

Sx(ω)=limT+12T|Fx(ω,T)|2=limT+12T|TTejωtx(t)dt|2

称为确定性信号 x(t)ω 处的功率谱密度


对于平稳过程,类似地,有如下定义

{X(t),<t<+} 平稳过程,称

limT+E[12TTTX2(t)dt]

{X(t),<t<+}平均功率,称

limT+12TE[|FX(ω,T)|2]

{X(t),<t<+}功率谱密度,其中

FX(ω,T)=TTejωtX(t)dt

对于平稳过程 {X(t),<t<+},若 RX(τ) 绝对可积,那么 {X(t),<t<+}谱密度是功率谱密度,即

SX(ω)=limT+12TE[|TTejωtX(t)dt|2]=limT+12TE[|FX(ω,T)|2]

【Wiener-Khintchine 公式】

{X(t),<t<+} 是均方连续的平稳过程,且 RX(τ) 绝对可积,即

+|RX(τ)|dτ<+

FX(ω),且有 Wiener-Khintchine 公式

RX(ω)=12π+ejωπSX(ω)dω,<τ<+SX(ω)=+ejωπRX(τ)dτ,<ω<+

当取 ω=0 时,有

RX(0)=12π+SX(ω)dωSX(0)=+RX(τ)dτ

其中,RX(0) 说明了功率谱密度曲线下的总面积(平均功率)等于平稳过程的均方值,SX(0) 说明了功率谱密度的零频率分量等于相关函数曲线下的总面积

可以发现,Wiener-Khintchine 公式给出了平稳过程相关函数和谱密度的转换关系,揭示了从时间角度描述平稳过程的统计规律和从频率角度描述平稳过程的统计规律间的联系

在实际应用中,根据实际情况去选择时域方法或等价的频域方法来解决问题

【谱密度的计算】

平稳过程谱密度的计算,包括由相关函数计算谱密度和由谱密度计算相关函数两方面

由 Wiener-Khintchine 公式可知,实际上是计算 Fourier 变换和逆 Fourier 变换的问题,因此,计算方法有两种:

  1. 直接计算积分
  2. 利用 Fourier 变换性质与常用的相关函数和谱密度变换结果进行计算

下图给出了常用的相关函数 RX(τ) 和谱密度 SX(ω) 的变换


例如:已知平稳过程的相关函数为 RX(τ)=5+4e3|τ|cos22τ,求其谱密度 SX(ω)

由于

RX(τ)=5+4e3|τ|cos22τ=5+2e3|τ|+2e3|τ|cos4τ

利用 Fourier 变换性质与上表,有:

SX(ω)=F(RX(τ))=5F(1)+2F(e3|τ|)+2F(e3|τ|cos4τ)=10πδ(ω)+269+ω2+2(39+(ω4)2+39+(ω+4)2)=10πδ(ω)+129+ω2+69+(ω4)2+69+(ω+4)2

其中,F 代表 Fourier 变换

【互谱密度】

定义

{X(t),<t<+}{Y(t),<t<+} 是联合平稳的平稳过程,若互相关函数 RXY(τ) 绝对可积,即

+|RXY(τ)|dτ<+

则称

SXY(ω)=+ejωtRXY(τ)dτ,<ω<+

{X(t),<t<+}{Y(t),<t<+}互谱密度

互谱密度一般是 ω 的复函数,不像谱密度那样具有物理意义,其作用在于:将时域上描述 {X(t),<t<+}{Y(t),<t<+} 相互关系的互相关函数,转换到频域上来研究他们的相关关系

性质

互谱密度具有如下简单的性质:

1)SXY(ω)=SYX(ω)

2)RXY(τ)SXY(ω) 是一对 Fourier 变换,即

RXY(τ)=12π+ejωτSXY(ω)dω,<τ<+SXY(ω)=+ejωτRXY(τ)dτ,<ω<+

3)若 {X(t),<t<+}{Y(t),<t<+} 是实联合平稳的平稳过程,则 SXY(ω) 的实部 Re(SXY(ω)) 是偶函数,虚部 Im(SXY(ω)) 是奇函数

4)|SXY(ω)|2SX(ω)SY(ω),|SYX(ω)|2SX(ω)SY(ω)

感谢您对我的支持,让我继续努力分享有用的技术与知识点!
0 comments
Anonymous
Markdown is supported

Be the first person to leave a comment!