随机过程笔记：2.谱分析 - 无推导版

谱分析Spectral Analysis of Stochastic Process

谱，从某种角度出发，进行分解，以把握特征。

1.引入：确定性信号的分解

对于确定性周期信号：$X(t),deterministic,Periodic:X(t+T)=X(t)$

，有：$x(t)={\sum }_k{\alpha }_k{e}^{j{\omega }_{k}t},{\omega }_k=\frac{2k\pi }{T},\frac{2\pi }{T}$为基频Base Frequency，其中${\alpha }_k=\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}x(t)e^{-j{\omega }_{k}t}dt$

1. 展开仅仅成立在$t\in [-\frac{T}{2},\frac{T}{2}]$
2. 对于非周期函数也可在区间做傅里叶级数展开，区间之外傅里叶级数是其周期延拓
3. 基函数$\frac{1}{\sqrt{T}}{e}^{j{\omega }_{k}t},k\in (-\infty ,+\infty )$是规范正交基。此时与对应函数做内积就可以直接得到系数，相当于在对应方向上的投影。

当$T\to \infty$，则有傅里叶变换：
$$\begin{array}{r}\underset{T\to \infty }{\lim }x(t)=\underset{T\to \infty }{\lim }\frac{1}{T}\sum _{k=-\infty }^{+\infty }[{\int }_{-\frac{T}{2}}^{\frac{T}{2}}x(s)e^{-j{\omega }_{k}s}ds]e^{j{\omega }_{k}t}=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }X(\omega )e^{j\omega t}d\omega \end{array}$$
傅里叶变换对：
$$\{\begin{array}{rl}x(t)& =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }X(\omega )e^{j\omega t}d\omega \\ X(\omega )& ={\int }_{-\infty }^{+\infty }x(t)e^{-j\omega t}dt\end{array}$$

2.随机信号的分解——功率谱密度

2.1.功率谱密度与相关函数

相比确定信号，随机信号可能存在一个问题，积分是否收敛？

对于傅里叶变换积分收敛，存在一个条件：$x(t)\in L^1(\mathbb{R})\iff {\int }_{-\infty }^{+\infty }|x(t)|dt<\infty$。该条件对于确定性信号，是普遍满足的。对于不满足的情况，通常也会做一些处理。

例如cos(t)，引入广义函数$\frac{1}{2}[\delta (\omega -1)+\delta (\omega +1)]$

面对这样的问题，可以提供两种解决办法。物理的角度：可以想到的是，希望做傅里叶的是有衰减趋势的函数，可以考虑二阶的函数。大部分相关函数是衰减的（也有周期振荡的）。下面以复随机信号(宽平稳)为例，根据推导可得：
$$\begin{array}{rl}& \underset{T\to \infty }{\lim }\frac{1}{T}E|{\int }_{-\frac{T}{2}}^{\frac{T}{2}}x(t)e^{-j\omega t}dt{|}^2={\int }_{-\infty }^{+\infty }{R}_x(u)e^{-j\omega u}du=S_x(\omega )\end{array}$$
上述结果即为功率谱密度Power Spectral Density，简称PSD。得到随机信号的傅里叶变换对（由相关函数的傅里叶变换得到功率谱密度）：
$$\{\begin{array}{rl}{S}_x(\omega )=& {\int }_{-\infty }^{+\infty }{R}_x(u)e^{-j\omega u}du\\ R_x(u)=& \frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{S}_x(\omega )e^{j\omega u}d\omega \end{array}$$

2.2.功率谱密度的性质

• 功率：量纲是$I^{2}T$焦耳$J$，即能量。换个角度，可得到$R_x(0)=\frac{1}{2\pi }{\int }_{-\infty }^{\infty }{S}_x(\omega )d\omega =E|x(t){|}^2$，为功率，则$S_x(\omega )$单位为功率除以频率，也就是功率乘以时间，故量纲是$J$，即能量。

• 谱：功率谱密度反映的是随机过程在每个频点上功率的大小，是一个二阶量。

  • $S_{\alpha x}(\omega )=|\alpha {|}^2{S}_x(\omega )$

  • $R_x(0)-R_x(\tau )\ge \frac{1}{4}[R_x(0)-R_x(2\tau )]$

  • $R_x(0)-R_x(\tau )\ge \frac{1}{4^n}[R_x(0)-R_x(2^n\tau )]$

• $S_{x+y}\ne S_x(\omega )+S_y(\omega )$

• 密度：体现在是常数。

• 关于$S_x(\omega )\ge 0$，从另外一个角度分析。相关函数$R_x(t)$是正定的，根据Bochner的结果，其傅里叶变换也是正的。

• 如果是实变量，功率谱是偶函数：$S_x(\omega )=S_x(-\omega )$。实信号没有负频率的说法，其频率负轴是频率正轴的镜像。

• 随机过程通过线性系统：有$Y(t)=(h\ast x)(t)={\int }_{-\infty }^{\infty }h(t-\tau )x(\tau )d\tau$，其中$h(t)$是系统的冲激响应，$H(\omega )$是传递函数。

X(t)

H(LTI：线性时不变系统)

Y(t)

宽平稳的随机过程通过线性系统仍然是宽平稳。
$$\begin{array}{r}{S}_y(\omega )=S_x(\omega )|H(\omega ){|}^2\end{array}$$
Wiener-Khinchine Relation.

Wiener：Cybernetics控制论，数学，美

Khinchine：排队论之父，前苏联