随机过程¶
基本概念¶
随机过程是在时间进程中处于不同时刻的随机变量的集合。
随机过程其实是随机变量在连续时间上的一次“升维”
分布函数/概率密度函数¶
设\(\xi(t)\)表示一个随机过程,则它在任意时刻\(t_{1}\)的值\(\xi\left(t_{1}\right)\)是一个随机变量
定义一维分布函数 $$ F_{1}\left(x_{1},t_{1}\right)=P[\xi(t_{1})\leqslant x_{1}] $$ 一维概率密度函数 $$ f_1(x_1,t_1)=\frac{\partial F_1(x_1,t_1)}{\partial x_1} $$
定义二维分布函数 $$ F_{2}(x_{1},x_{2};t_{1},t_{2})=P|\xi(t_{1})\leqslant x_{1},\xi(t_{2})\leqslant x_{2}| $$ 二维概率密度函数 $$ f_2(x_1,x_2;t_1,t_2)=\frac{\partial^2F_2(x_1,x_2;t_1,t_2)}{\partial x_1\cdot\partial x_2} $$
数字特征¶
数学期望¶
方差¶
可以证明:
自相关函数¶
自相关函数反映不同时刻随机过程取值的相关性
自协方差函数¶
自协方差函数反映不同时刻随机过程的统计相关性
且有 $$ B(t_1,t_2)=R(t_1,t_2)-a(t_1)a(t_2) $$
互相关函数¶
平稳随机过程¶
定义¶
若一个随机过程\(\xi(t)\)的统计特性与时间起点无关,即时间平移不影响其任何统计特性,则称该随机过程是在严格意义下的平稳随机过程,简称严平稳随机过程。 $$ f_{n}(x_{1},x_{2},\cdots,x_{n};t_{1},t_{2},\cdots,t_{n})=\f_{n}(x_{1},x_{2},\cdots,x_{n};t_{1}+\Delta,t_{2}+\Delta,\cdots,t_{n}+\Delta) $$
定义广义平稳过程,如果随机过程\(\xi(t)\)满足
期望与\(t\)无关 $$ E\left[\xi(t)\right]=\int_{-\infty}^{\infty}x_{1}f_{1}\left(x_{1}\right)\mathrm{d}x_{1}=a $$ 自相关函数只与时间间隔\(\tau\)有关 $$ \begin{aligned} R(t_{1},t_{2})&=E[\xi(t_{1})\xi(t_{1}+\tau)]\\&=\int_{-\infty}{\infty}\int_{-\infty}=R(\tau) \end{aligned} $$}x_{1}x_{2}f_{2}\left(x_{1},x_{2};\tau\right)\mathrm{d}x_{1}\mathrm{d}x_{2
各态历经性¶
假设\(x(t)\)是平稳过程\(\xi\left(t\right)\)的任意一次实现
定义 $$ \overline{a}=x(t)=\lim_{T\to\infty}\frac{1}{T}\int_{-\frac{2}{T}}^{\frac{2}{T}} x(t)dt \ \overline{R(\tau)}=x(t)x(t+\tau)=\lim_{T\to\infty}\frac{1}{T}\int_{-\frac{2}{T}}^{\frac{2}{T}} x(t)x(t+\tau)dt $$ 若 $$ \begin{cases}a=\overline{a}\\R(\tau)=\overline{R(\tau)}&\end{cases} $$ 则称该平稳过程具有各态历经性
各态历经性的含义:平稳过程中的任一次实现都经历了随机过程的所有可能状态.
充分条件:具有各态历经的随机过程一定是平稳过程.
判断一个随机过程是否广义平稳只需验证下式成立与否
$$ \begin{cases}a(t)=a\\R(t_1,t_1+\tau)=R(\tau)&\end{cases} $$ 判断一个平稳过程是否各态历经只需验证下式成立与否
平稳过程的自相关函数¶
定义平稳过程的自相关函数 $$ R(\tau)=E[\xi(t)\xi(t+\tau)] $$ 1. 在时域中计算平均功率\(R(0)=E[\xi^2(t)]\)且为最大值 1. 偶函数\(R(\tau)=R(-\tau)\) 1. 直流功率\(R\left(\infty\right)=E^{2}\left[\xi\left(t\right)\right]=a^{2}\)
$$ \lim_{\tau\to\infty}R(\tau)=\lim_{\tau\to\infty}E\left[\xi(t)\xi(t+\tau)\right]=E\left[\xi(t)\right]\cdot E\left[\xi(t+\tau)\right]=E^{2}\left[\xi(t)\right] $$ 上式利用了当\(\tau\to\infty\)时,\(\xi(t)\)与\(\xi(t+\tau)\)没有任何依赖关系,即统计独立
- 交流功率\(R(0)-R(\infty)=\sigma^{2}\)
平稳过程的功率谱函数¶
随机过程中的任一样本是一个确定的功率型信号
定义功率谱密度 $$ P_{x}(f)=\lim_{T\to\infty}\frac{\mid X_{T}(f)\mid^{2}}{T} $$
维纳-辛钦定理:平稳过程的功率谱密度与其自相关函数是一对傅里叶变换关系 $$ \begin{cases} P_\xi(\omega)=\int_{-\infty}^\infty R(\tau)\mathrm{e}^{-\mathrm{j}\omega\tau}\mathrm{d}\tau\\ R(\tau)=\frac{1}{2\pi}\int_{-\infty}^\infty P_\xi(\omega)\mathrm{e}^{\mathrm{j}\omega\tau}\mathrm{d}\omega& \end{cases} $$
或 $$ \begin{cases} P_\xi(f)&=\int_{-\infty}^\infty R(\tau)\mathrm{e}^{-\mathrm{j}\omega\tau}\mathrm{d}\tau\\ R(\tau)&=\int_{-\infty}^\infty P_\xi(f)\mathrm{e}^{\mathrm{j}\omega\tau}\mathrm{d}f&\end{cases} $$ 简记为 $$ R(\tau)\Leftrightarrow P_{\xi}(f) $$
-
在频域中计算平稳过程的平均功率: $$ R(0)=\int_{-\infty}^{\infty}P_{\xi}(f)\mathrm{d}f $$
-
各态历经过程的任一样本函数的功率谱密度等于过程的功率谱密度
-
功率谱密度具有非负性\(P_{\xi}(f)\geq0\)和实偶性\(P_{\xi}(-f)=P_{\xi}(f)\)
高斯随机过程¶
定义高斯过程 $$ f_{n}(x_{1},x_{2},\cdots,x_{n};t_{1},t_{2},\cdots,t_{n})=\ \frac{1}{(2\pi)^\frac{n}{2}\sigma_{1}\sigma_{2}\cdots\sigma_{n}\mid B\mid^\frac{1}{2}}\mathrm{exp}\left[-\frac{1}{2\mid B\mid}\sum_{j=1}{n}\sum_{k=1}\right)\right] $$}\mid B\mid_{jk}\left(\frac{x_{j}-a_{j}}{\sigma_{j}}\right)\left(\frac{x_{k}-a_{k}}{\sigma_{k}
- \(n\)维分布只依赖各个随机变量的均值、方差和归一化协方差,我们只需要研究这些数字特征
- 广义平稳的高斯过程也是严平稳的
- 如果高斯过程在不同时刻的取值不相关,那么它们是统计独立的
- 高斯过程经线性变换后(输入线性系统后的输出)仍是高斯过程
高斯随机变量¶
一维概率密度函数¶
高斯过程在任一时刻上的取值是一个正态分布的随机变量高斯随机变量,其一维概率密度函数为 $$ f(x)=\frac{1}{\sqrt{2\pi}\sigma}\exp\left(-\frac{\left(x-a\right){2}}{2\sigma\right) $$ 1. }\(f(a+x)=f(a-x)\)
-
\(\int_{-\infty}^{\infty}f(x)\mathrm{d}x=1\)
-
当\(a=0\),\(\sigma=1\)时,称为标准化的正态分布 $$ f(x)=\frac{1}{\sqrt{2\pi}}\exp\left(-\frac{x^{2}}{2}\right) $$
正态分布函数¶
定义正态分布的概率密度函数的积分为正态分布函数 $$ F(x)=P(\xi\leqslant x)=\int_{-\infty}{x}\frac{1}{\sqrt{2\pi}\sigma}\mathrm{exp}\left[-\frac{(z-a)z $$}}{2\sigma^{2}}\right]\mathrm{d
我们可以用以下几种函数来进行变量代换来方便计算
- 使用误差函数\(\operatorname{erf}(x)\)
令\(t=\frac{(z-a)}{\sqrt{2}\sigma}\),则\(\mathrm{d}z=\sqrt{2}\sigma\mathrm{d}t\),有: $$ F(x)=\frac{1}{2}\cdot\frac{2}{\sqrt{\pi}}\int_{-\infty}{(x-a)/\sqrt{2}\sigma}\mathrm{e}\right) $$ 其中 $$ \mathrm{erf}(x)=\frac{2}{\sqrt{\pi}}\int_{0}}}\mathrm{d}t=\frac{1}{2}+\frac{1}{2}\mathrm{erf}\left(\frac{x-a}{\sqrt{2}\sigma{x}\mathrm{e}t $$}}\mathrm{d
- 使用互补误差函数\(erfc(x)\)
此时 $$ F(x)=1-\frac{1}{2}\mathrm{erfc}\left(\frac{x-a}{\sqrt{2}\sigma}\right) $$
- 使用\(Q(X)\)
那么 $$ P(\xi>x)=Q\left(\frac{x-a}{\sigma}\right) $$
平稳随机过程通过线性系统¶
关系式¶
\(v_{\mathrm{i}}(t)\)是输入随机过程的一个样本,\(v_{o}(t)\)是输出随机过程的一个样本,随机过程\(\xi_{\mathrm{i}}\left(t\right)\)的样本包括[\(v_{\mathrm{i},n}\left(t\right)\) ,\(n=1,2,\cdots\)]
确知信号通过线性时不变系统时,可由系统的单位冲激响应表征 $$ v_{\mathrm{o}}(t)=h(t)*v_{\mathrm{i}}(t)=\int_{-\infty}^{\infty}h_{\mathrm{i}}(\tau)v_{\mathrm{i}}(t-\tau)\mathrm{d}\tau $$ 于是输入与输出随机过程也应满足 $$ \xi_{o}(t)=\int_{-\infty}^{\infty}h(\tau)\xi_{i}(t-\tau)\mathrm{d}\tau $$
输出过程的统计特性¶
- 输出过程\(\xi_{o}(t)\)的均值
输入过程是平稳的,有 $$ E\left[\xi_{\mathrm{i}}\left(t-\tau\right)\right]=E\left[\xi_{\mathrm{i}}\left(t\right)\right]=a $$ 所以 $$ \label{输出过程均值} E\left[\xi_{\mathrm{o}}(t)\right]=a\cdot\int_{-\infty}^{\infty}h(\tau)\mathrm{d}\tau=a\cdot H(0) $$
上式表明输出过程的均值是一个常数,其中\(H(0)\)是直流增益
- 输出过程\(\xi_{o}(t)\)的自相关函数
输入过程是平稳的,有 $$ E[\xi_{\mathrm{i}}(t_{1}-\alpha)\xi_{\mathrm{i}}(t_{1}+\tau-\beta)]=R_{\mathrm{i}}(\tau+\alpha-\beta) $$ 所以 $$ \label{输出过程自相关函数} R_{\mathrm{o}}(t_{1},t_{1}+\tau)=\int_{-\infty}{\infty}\int_{-\infty}(\tau) $$ 上面}h(\alpha)h(\beta)R_{\mathrm{i}}(\tau+\alpha-\beta)\mathrm{d}\alpha\mathrm{d}\beta=R_{\mathrm{o}\(\ref{输出过程均值}\)和\(\ref{输出过程自相关函数}\)两式表明,若线性系统的输入过程是平稳的,那么输出过程也是平稳的
- 功率谱密度
对上式进行傅里叶变换 $$ P_{o}(f)=\int{\infty}R_{_o}(\tau)\mathrm{e}\tau $$}\omega\tau}\mathrm{d
令\(\tau^{\prime}=\tau+\alpha-\beta\) $$ P_{o}\left(f\right)=\int_{-\infty}{\infty}h\left(\alpha\right)\mathrm{e}}\omega\alpha}\mathrm{d}\alpha\int_{-\infty{\infty}h\left(\beta\right)\mathrm{e}}\omega\beta}\mathrm{d}\beta\int_{-\infty{\infty}R_{i}\left(\tau}\right)\mathrm{e{-\mathrm{j}\omega\tau\left(f\right)\ =H^{*}\left(f\right)\cdot H\left(f\right)\cdot P_{i}\left(f\right)=\mid H\left(f\right)\mid^{2}P_{i}\left(f\right) $$ 上式表明,输出过程的功率谱密度是输入过程的功率谱密度乘系统频率响应模值的平方}}\mathrm{d}\tau^{\prime}P_{o
- 输出过程\(\xi_{o}(t)\)的概率分布
上式表明,高斯过程经线性变换后的过程仍为高斯过程
窄带随机过程¶
若随机过程\(\xi(t)\)的谱密度集中在中心频率\(f_\mathrm{c}\)附近相对窄的频带范围\(\Delta f\)内,即满足 $$ \begin{cases} \Delta f\ll f_{\mathrm{c}}\\ f_\mathrm{c}远离零频率 \end{cases} $$ 则称该\(\xi(t)\)为窄带随机过程
一个典型的窄带随机过程的频谱密度和波形
窄带随机过程\(\xi(t)\)可表示为 $$ \xi(t)=a_{\xi}(t)\cos\left[\omega_{\mathrm{c}}t+\varphi_{\xi}(t)\right]\quad a_{\xi}(t)\geqslant0 $$
其中\(a_{\xi}(t)\)称为窄带随机过程\(\xi(t)\)的随机包络 \(\varphi_{\xi}(t)\)称为随机相位 \(\omega_{\mathrm{c}}\)称为中心频率
上式展开得 $$ \xi(t)=\xi_{\mathrm{c}}(t)\cos\omega_{\mathrm{c}}t-\xi_{\mathrm{s}}(t)\sin\omega_{\mathrm{c}}t $$ 其中 $$ \xi_{\mathrm{c}}(t)=a_{\xi}(t)\cos\varphi_{\xi}(t) $$
分别称为\(\xi\left(t\right)\)的同相分量和正交分量
统计特性¶
现研究典型窄带随机过程\(\xi(t)\)的统计特性,其中\(\xi(t)\)是一个均值为0,方差为\(\sigma_{\xi}^{2}\)的平稳高斯窄带随机过程 $$ \xi(t)=\xi_{\mathrm{c}}(t)\cos\omega_{\mathrm{c}}t-\xi_{\mathrm{s}}(t)\sin\omega_{\mathrm{c}}t $$
- \(\xi_c(t)\)和\(\xi_s(t)\)的统计特性
若窄带过程\(\xi(t)\)是平稳的,则\(\xi_c(t)\)和\(\xi_s(t)\)也是平稳的
\(\xi_{c}(t)\)和\(\xi_{s}(t)\)为零均值,等方差,不相关,独立的平稳高斯过程。
- \(a_\xi\left(t\right)\)和\(\varphi_\xi\left(t\right)\)的统计特性
此\(\xi(t)\)的窄带平稳高斯过程的包络\(a_\xi\left(t\right)\)的一维分布是瑞利分布,相位\(\varphi_\xi\left(t\right)\)的一维分布是均匀分布,且就一维分布而言二者统计独立 $$ \begin{cases} f(a_\xi)=\frac{a_\xi}{\sigma_\xi2}\exp\left(-\frac{a_\xi2}{2\sigma_\xi^2}\right),a_\xi\geq0\\ f(\varphi_\xi)=\frac{1}{2\pi},0\leq\varphi_\xi\leq2\pi\\ f(a_\xi,\varphi_\xi)=f(a_\xi)\cdot f(\varphi_\xi) \end{cases} $$
正弦波加窄带高斯噪声¶
数学模型¶
带通滤波器的输出是正弦波已调信号与窄带高斯噪声的混合波形
设正弦波加窄带高斯噪声的混合信号为 $$ r(t)=A\cos(\omega_{\mathrm{c}}t+\theta)+n(t) $$ 其中\(n(t)\)为窄带高斯噪声,其均值为\(0\) 方差为\(\sigma_{n}^{2}\) $$ n(t)=n_{\mathrm{c}}(t)\cos\omega_{\mathrm{c}}t-n_{\mathrm{s}}(t)\sin\omega_{\mathrm{c}}t $$ 于是 $$ \begin{aligned} r(t)&=\left[A\mathrm{cos}\theta+n_\mathrm{c}(t)\right]\mathrm{cos}\omega_\mathrm{c}t-\left[A\mathrm{sin}\theta+n_\mathrm{s}(t)\right]\mathrm{sin}\omega_\mathrm{c}t\\ &=z_{\mathrm{c}}(t)\cos\omega_{\mathrm{c}}t-z_{\mathrm{s}}(t)\sin\omega_{\mathrm{c}}t\\ &=z(t)\cos\left[\omega_{\mathrm{c}}t+\varphi(t)\right]\ \end{aligned} $$ 其中 $$ z_{\mathrm{c}}(t)=A\mathrm{cos}\theta+n_{\mathrm{c}}(t) $$
\(r(t)\)的包络和相位分别为 $$ z(t)=\sqrt{z_{\mathrm{c}}{2}(t)+z_{\mathrm{s}}\quad z\geqslant0 $$}(t)
统计特性¶
随机包络\(z(t)\)的统计特性:莱斯分布/广义瑞利分布 $$ f(z)=\frac{z}{\sigma_n2}\exp\left(-\frac{z2+A2}{2\sigma_n2}\right)I_0\left(\frac{Az}{\sigma_n^2}\right),z\geq0 $$ 大信噪比环境:\(r= \frac {A^2}{2\sigma _n^2}\gg 1\) 近似为高斯分布
小信噪比环境:\(A\to0\)退化为瑞利分布
随机相位\(\varphi(t)\)的统计特性:
小信噪比环境:\(f(\varphi)\)服从均匀分布
大信噪比环境:集中在有用信号相位附近,在\(0\)相位附近取值集中。
高斯白噪声和带限白噪声¶
白噪声¶
定义白噪声:噪声的功率谱密度在所有频率上均为一常数,即其双边功率谱密度为: $$ P_n(f)=\frac{n_0}{2}\quad(-\infty<f<+\infty) $$ 或单边功率谱密度为: $$ P_n(f)=n_0\quad(0<f<+\infty) $$ 取傅里叶反变换得白噪声的自相关函数 $$ R(\tau)=\frac{n_0}{2}\delta(\tau) $$
- 白噪声仅在\(\tau=0\)时相关,而在任意两个时刻的随机变量都不相关
-
白噪声的平均功率无穷大,即\(R(0)=\frac{n_0}{2}\delta(0)=\infty\)
-
如果白噪声取值的概率服从高斯分布,则称为高斯白噪声,其通常用来作为通信信道中的噪声模型
- 实际中不存在理想白噪声,因为其功率为无穷大,若噪声功率谱均匀分布的频率范围远大于通信系统的工作频带,则视为白噪声。
低通白噪声¶
定义低通白噪声:白噪声通过理想矩形的低通滤波器或理想低通信道,则输出的噪声称为低通白噪声
假设理想低通滤波器传输特性:模为 \(1\) ,截止频率为\(|f|\leqslant f_{H}\)
则低通白噪声的功率谱密度
$$
P_n(f)=
\begin{cases}\frac{n_0}{2}&|f|\leqslant f_H\\0&其他&
\end{cases}
$$
自相关函数为
$$
R(\tau)=n_0f_H\frac{\sin2\pi f_H\tau}{2\pi f_H\tau}
$$

此低通白噪声也通常称作带限白噪声,其在\(\tau=\frac{k}{2f_{H}}(k=1,2,3)\)上得到的随机变量不相关
带通白噪声¶
定义带通白噪声:白噪声通过理想矩形的带通滤波器或理想带通信道,则输出的噪声称为带通白噪声
假设理想带通滤波器的传输特性为:
$$
H(f)=\begin{cases}1&f_c-\frac{B}{2}\leqslant\mid f\mid\leqslant f_c+\frac{B}{2}\\0&\text{其他}&\end{cases}
$$
则带通白噪声的功率谱密度
$$
P_n(f)=
\begin{cases}\frac{n_0}{2}&f_\mathrm{c}-\frac{B}{2}\leqslant\mid f\mid\leqslant f_\mathrm{c}+\frac{B}{2}\\0&\text{其他}&
\end{cases}
$$
自相关函数
$$
\begin{aligned}
R\left(\tau\right)&=\int_{-\infty}{\infty}P_{n}\left(f\right)\mathrm{e}f\\
&=\int_{-f_{c}-\frac{B}{2}}}2\pi f\tau}\mathrm{d{-f_{c}+\frac{B}{2}}\frac{n_{0}}{2}\mathrm{e}}2\pi f\tau}\mathrm{d}f+\int_{f_{c}-\frac{B}{2}{f_{c}+\frac{B}{2}}\frac{n_{0}}{2}\mathrm{e}f\\
&=n_0B\frac{\sin\pi B\tau}{\pi B\tau}\mathrm{cos}2\pi f_\mathrm{c}\tau
\end{aligned}
$$
}2\pi f\tau}\mathrm{d
通常,带通滤波器\(B\ll f_{c}\) 因此称为窄带滤波器,相应地把带通白噪声称为窄带高斯白噪声,其统计特性与前面描述的窄带随机过程相同,因此有 $$ n(t)=n_{\mathrm{c}}(t)\cos\omega_{\mathrm{c}}t-n_{\mathrm{s}}(t)\sin\omega_{\mathrm{c}}t $$
参考上式,平均功率为 $$ N=n_0B $$
创建日期: 2025-05-26
作者:
