ESC
输入关键词搜索文章
目录

随机信号分析

通信原理 · 第二章
平稳随机过程、自相关与功率谱密度
7核心节
4核心分布
2关键定理
1实用公式
导论
学习目标与课程定位

第一章建立了通信系统的整体框架(信源→信道→信宿),引入了信息量、误码率等性能指标。但通信系统中传输的信号不可避免地受到噪声干扰,而噪声的本质是随机的——你无法预知下一个时刻噪声的精确取值,只能描述其统计规律。

本章回答的核心问题:如何用数学工具描述和分析随机信号?这是后续所有章节(信道分析、基带传输、调制解调、最佳接收)的数学基础。

PDF第二章 随机信号分析 开篇p.1

pdf/通信原理/第二章.pdf · p.1

打开原文

本章学习目标
  1. 理解随机过程的定义及其与随机变量的区别
  2. 掌握平稳随机过程(严平稳/广义平稳)与各态历经性
  3. 熟练运用自相关函数的性质与维纳-辛钦定理
  4. 掌握高斯过程、窄带过程、白噪声的统计特性
  5. 理解正弦波加窄带高斯过程的莱斯分布
  6. 能分析随机过程通过线性系统的输入输出关系

前置知识回顾

本章需要概率论中以下概念:

概念公式说明
期望$E[X] = \int_{-\infty}^{\infty} x f(x) dx$随机变量的"平均值"
方差$D[X] = E[(X - E[X])^2]$偏离均值的程度
协方差$\text{Cov}(X,Y) = E[XY] - E[X]E[Y]$两个随机变量的线性关联
正态分布$f(x) = \frac{1}{\sqrt{2\pi}\sigma} e^{-\frac{(x-\mu)^2}{2\sigma^2}}$最重要的连续分布

关键直觉:概率论研究的是"静态"的随机变量,而本章要研究的是随时间演化的随机变量族——也就是随机过程

§2.1
随机过程的基本概念

为什么需要随机过程?

如果你测量一个电阻的热噪声电压,每次测量得到一条随时间变化的波形。重复实验,每次得到的波形都不同。你无法写出一条确定的函数 $v(t)$ 来描述它,但你可以描述所有可能波形的统计规律。这种"一族时间函数"就是随机过程。

PDF2.1 随机过程的定义p.2

pdf/通信原理/第二章.pdf · p.2

打开原文

定义

设随机实验 $E$ 的样本空间为 $S$,每个实验结果 $\xi \in S$ 对应一个时间函数 $X(\xi, t)$,则称 $\{X(\xi, t), t \in T\}$随机过程,简记为 $X(t)$

两个基本属性(对理解非常关键):
视角固定量变化量结果
固定 $\xi$(选定一次实验)样本时间 $t$ 变化得到一条确知的样本函数
固定 $t = t_1$(选定一个时刻)时间$\xi$ 变化得到一个随机变量 $X(t_1)$

类比:想象拍电影。固定一部电影(固定 $\xi$),随着时间推进得到一帧帧画面(样本函数);固定某一时刻(固定 $t_1$),所有电影在这一帧的画面构成一个随机变量。

分布函数与概率密度

既然固定时刻 $t_1$$X(t_1)$ 是随机变量,它就有概率分布:

$$F(x, t_1) = P\{X(t_1) \leq x\}$$
$$f(x, t_1) = \frac{\partial F(x, t_1)}{\partial x}$$

要描述两个时刻之间的关联,需要二维分布:

$$f(x_1, x_2; t_1, t_2) = \frac{\partial^2 F(x_1, x_2; t_1, t_2)}{\partial x_1 \partial x_2}$$

以此类推,有 $n$ 维分布。维度越高,描述越完整,但计算也越复杂。实际工程中,通常只用一二阶数字特征(均值、相关函数)就够了。

数字特征

(1)均值(数学期望)

$$a(t) = E[X(t)] = \int_{-\infty}^{\infty} x \cdot f(x, t) dx$$

物理意义:随机过程在时刻 $t$ 的所有样本的"中心位置"。它本身是 $t$ 的确定函数,所有样本在它附近起伏。

(2)方差

$$\sigma^2(t) = D[X(t)] = E[(X(t) - a(t))^2]$$

物理意义:样本在时刻 $t$ 相对均值的偏离程度,即"起伏有多剧烈"。

(3)自相关函数

$$R(t_1, t_2) = E[X(t_1) X(t_2)]$$

物理意义:衡量 $X(t)$$t_1$$t_2$ 两个时刻取值之间的线性关联程度。若 $R$ 很大,说明两个时刻的值"同步变化"的趋势强。

(4)协方差函数

$$C(t_1, t_2) = E[(X(t_1) - a(t_1))(X(t_2) - a(t_2))] = R(t_1, t_2) - a(t_1) \cdot a(t_2)$$

当均值 $a(t) = 0$ 时,协方差函数 $=$ 自相关函数。通信中常假设均值为零,所以两者等价。

记忆技巧$R$ 是"原点矩"(直接相乘),$C$ 是"中心矩"(减去均值再相乘)。两者关系:$C = R - a_1 \cdot a_2$
随机过程的样本函数族JSXGraph
随机过程的样本函数族
§2.2
平稳随机过程

为什么需要"平稳"?

一般随机过程的均值 $a(t)$ 和自相关 $R(t_1, t_2)$ 都是时间的复杂函数,分析起来极其困难。如果随机过程的统计特性不随时间平移而改变,即"平稳",那均值就退化为常数,自相关就只取决于时间差——大大简化分析。

PDF2.2 平稳随机过程p.5

pdf/通信原理/第二章.pdf · p.5

打开原文

严平稳 vs 广义平稳

严平稳(狭义平稳):随机过程的 $n$ 维概率密度函数与时间起点无关:

$$f(x_1, x_2, \ldots, x_n; t_1, t_2, \ldots, t_n) = f(x_1, x_2, \ldots, x_n; t_1 + \tau, t_2 + \tau, \ldots, t_n + \tau)$$

这要求所有阶的统计特性都不随时间变化,条件非常苛刻。

广义平稳(宽平稳):只需要满足两个条件:

  1. 均值为常数$E[X(t)] = a$(与 $t$ 无关)
  2. 自相关函数只与时间间隔有关$R(t_1, t_2) = R(t_2 - t_1) = R(\tau)$
关系:严平稳 $\Rightarrow$ 广义平稳(条件更强推出更弱);但广义平稳 $\nRightarrow$ 严平稳。例外:高斯过程中两者等价。工程中说的"平稳"默认指广义平稳。

各态历经性(遍历性)

问题:即使知道平稳过程,要算 $E[X(t)]$ 仍然需要"无数条样本函数"做统计平均,实际中不可能。

核心思想:如果随机过程有遍历性,那么时间平均 = 统计平均

$$\lim_{T \to \infty} \frac{1}{2T} \int_{-T}^{T} X(t) dt = E[X(t)]$$

左端是时间平均(一条样本就够),右端是统计平均(需要所有样本)。这意味着:用一条足够长的样本就能推断整个随机过程的统计特性。

"两个骰子"类比(来自 [W05]):有两个骰子,一个只能掷出 1-3,另一个只能掷出 4-6。等概率随机选一个来掷。统计均值为 $E = (2 + 5)/2 = 3.5$。但如果你选定了一个骰子(一条样本),你永远只能看到 1-3 或 4-6——遍历不到全部状态。所以这个平稳过程不具遍历性。遍历性要求任意一条样本都能"经历"所有可能状态。

自相关函数的性质(平稳过程)

对平稳过程,自相关函数 $R(\tau)$ 有以下重要性质:

性质公式直觉
$\tau=0$ 给出平均功率$R(0) = E[X^2(t)]$信号自身的均方值=功率
偶函数$R(\tau) = R(-\tau)$过去和未来的相关性对称
$\tau=0$ 处取最大值$R(0) \geq |R(\tau)|$任何时刻与自己最相关
非负定其傅里叶变换 $\geq 0$功率谱密度不可能为负
平稳过程自相关函数 R(τ) 典型形状JSXGraph
平稳过程自相关函数 R(τ) 典型形状

频谱特性——功率谱密度与维纳-辛钦定理

随机信号无法直接做傅里叶变换(能量无限,不满足绝对可积),但我们关心的是功率的频域分布。

推导思路(三步走,来自 [W04] 的完整推导链):

  1. 截断:取 $X(t)$ 的一段 $x_T(t)$,在 $[-T, T]$ 内等于 $x(t)$,之外为零
  2. 傅里叶变换:截断信号满足可积条件,可做 FT 得 $F_x(\omega, T)$
  3. 取极限:功率谱密度定义为
    $$S_X(\omega) = \lim_{T \to \infty} \frac{1}{2T} E[|F_x(\omega, T)|^2]$$

维纳-辛钦定理:平稳随机过程的自相关函数与功率谱密度是一对傅里叶变换:

$$\boxed{R(\tau) \;\stackrel{\mathscr{F}}{\longleftrightarrow}\; S(\omega)}$$

$S(\omega) = \int_{-\infty}^{\infty} R(\tau) e^{-j\omega\tau} d\tau$$R(\tau) = \frac{1}{2\pi} \int_{-\infty}^{\infty} S(\omega) e^{j\omega\tau} d\omega$

意义:这个定理把时域的相关分析和频域的功率分析桥接起来了。已知 $R(\tau)$ 可以求 $S(\omega)$,反之亦然。在通信工程中,设计滤波器时经常在频域操作,而测量信号时用时域操作,维纳-辛钦定理让两者可以自由切换。
§2.3
高斯随机过程

为什么高斯过程这么重要?

通信中的热噪声、散弹噪声等,根据中心极限定理,都近似服从高斯分布。高斯过程是最常用的噪声模型。

PDF2.3 高斯随机过程p.15

pdf/通信原理/第二章.pdf · p.15

打开原文

定义

若随机过程 $X(t)$ 的任意 $n$ 维分布都服从正态分布,则称为高斯随机过程

四条重要性质

  1. 广义平稳 $\Leftrightarrow$ 严平稳:高斯过程中两种平稳等价(因为高斯分布完全由均值和协方差确定)
  2. 不相关 $\Leftrightarrow$ 独立:对高斯随机变量,不相关就意味着统计独立
  3. 高斯过程之和仍为高斯:多个高斯噪声叠加还是高斯噪声
  4. 通过线性系统后仍为高斯:这对分析噪声经过滤波器后的特性至关重要

一维概率密度

$$f(x) = \frac{1}{\sqrt{2\pi}\sigma} \exp\left[-\frac{(x - a)^2}{2\sigma^2}\right]$$

其中 $a$ 为均值,$\sigma^2$ 为方差。曲线关于 $x = a$ 对称,呈钟形。

误差函数

正态分布的累积分布函数 $\Phi(x)$ 没有闭合表达式,引入误差函数来简化表示:

$$\text{erf}(x) = \frac{2}{\sqrt{\pi}} \int_0^x e^{-t^2} dt$$
$$\text{erfc}(x) = 1 - \text{erf}(x) = \frac{2}{\sqrt{\pi}} \int_x^{\infty} e^{-t^2} dt$$

关键性质:

  • $\text{erf}(0) = 0$$\text{erf}(\infty) = 1$(递增函数)
  • $\text{erfc}(0) = 1$$\text{erfc}(\infty) = 0$(递减函数)
  • $\text{erfc}(-x) = 2 - \text{erfc}(x)$
为什么要用误差函数? 在计算通信系统误码率时,最终都会化为 $\text{erfc}$ 的形式。它简洁的特性有助于分析系统的抗噪性能。
高斯概率密度与误差函数JSXGraph
高斯概率密度与误差函数
§2.4
窄带随机过程

什么是窄带?

若随机过程的功率谱集中在某个中心频率 $\omega_c$ 附近,且带宽 $\Delta\omega \ll \omega_c$,则称为窄带随机过程。通信中的接收信号经过带通滤波器后就是典型的窄带过程。

PDF2.4 窄带随机过程p.20

pdf/通信原理/第二章.pdf · p.20

打开原文

两种表示方式

窄带过程 $\xi(t)$ 可以写成:

方式一:包络 + 相位表示

$$\xi(t) = a_\xi(t) \cos[\omega_c t + \varphi_\xi(t)]$$

其中 $a_\xi(t) \geq 0$ 是包络,$\varphi_\xi(t) \in [0, 2\pi)$ 是相位,它们都是随时间慢变化的(相对于 $\omega_c$)。

方式二:同相 + 正交分量表示

展开上式可得:

$$\xi(t) = \xi_c(t) \cos\omega_c t - \xi_s(t) \sin\omega_c t$$

其中 $\xi_c(t) = a_\xi(t)\cos\varphi_\xi(t)$同相分量$\xi_s(t) = a_\xi(t)\sin\varphi_\xi(t)$正交分量

两条核心结论

结论一(同相分量与正交分量的统计特性):

$\xi(t)$ 是均值为零、方差为 $\sigma_\xi^2$ 的窄带平稳高斯过程,则:

  • $\xi_c(t)$$\xi_s(t)$ 各自也是平稳高斯过程
  • 均值都为零:$E[\xi_c] = E[\xi_s] = 0$
  • 方差相同:$\sigma_c^2 = \sigma_s^2 = \sigma_\xi^2$
  • 同一时刻上 $\xi_c$$\xi_s$ 互不相关(也是统计独立的,因为高斯)
直觉:窄带高斯噪声被"解调"到基带后,I路和Q路是两个独立的零均值高斯噪声,方差与原噪声相同。

结论二(包络与相位的分布):

包络 $a_\xi(t)$ 服从瑞利分布

$$f(a) = \frac{a}{\sigma_\xi^2} \exp\left(-\frac{a^2}{2\sigma_\xi^2}\right), \quad a \geq 0$$

相位 $\varphi_\xi(t)$ 服从 $[0, 2\pi)$ 上的均匀分布

$$f(\varphi) = \frac{1}{2\pi}, \quad 0 \leq \varphi < 2\pi$$

包络与相位在同一时刻统计独立

瑞利分布的直觉:如果把 $\xi_c$$\xi_s$ 看成二维平面上互相垂直的高斯随机变量,包络 $a = \sqrt{\xi_c^2 + \xi_s^2}$ 就是从原点到该点的距离。距离的分布在原点为零(概率为零),在某个值处达到峰值,然后拖尾——这就是瑞利分布的形状。

窄带随机过程:时域波形与包络JSXGraph
窄带随机过程:时域波形与包络
窄带高斯过程包络的瑞利分布JSXGraph
窄带高斯过程包络的瑞利分布
§2.5
白噪声

理想白噪声

定义:功率谱密度在整个频域内均匀分布的噪声:

$$P(\omega) = \frac{n_0}{2} \quad (-\infty < \omega < \infty)$$

其中 $n_0$ 为正常数(单位:W/Hz)。称"白"是因为类比白光——所有频率分量强度相同。

自相关函数(由维纳-辛钦定理):

$$R(\tau) = \frac{1}{2\pi} \int_{-\infty}^{\infty} \frac{n_0}{2} e^{j\omega\tau} d\omega = \frac{n_0}{2} \delta(\tau)$$
关键含义:白噪声只在 $\tau = 0$ 时相关,任意两个不同时刻都完全不相关。这就像理想的"无记忆"噪声——这一刻的值对下一刻没有任何暗示。

物理现实:理想白噪声功率无限大,现实中不存在。但实际噪声带宽远大于信号带宽时,可近似为白噪声。

PDF2.5 白噪声p.25

pdf/通信原理/第二章.pdf · p.25

打开原文

带限白噪声

白噪声通过理想低通滤波器(截止频率 $f_0$)后变为带限白噪声:

$$P(f) = \begin{cases} \frac{n_0}{2}, & |f| \leq f_0 \\ 0, & |f| > f_0 \end{cases}$$

自相关函数:

$$R(\tau) = n_0 f_0 \cdot \frac{\sin(2\pi f_0 \tau)}{2\pi f_0 \tau} = n_0 f_0 \cdot \text{sinc}(2f_0 \tau)$$

$\tau = k/(2f_0)$$k \neq 0$ 的整数)时,$R(\tau) = 0$,即这些时刻上互不相关

与抽样定理的联系:如果以 $f_s = 2f_0$(奈奎斯特速率)对带限白噪声采样,各采样值互不相关。这在后续 PCM 系统分析中很重要。
白噪声功率谱密度与带限白噪声自相关JSXGraph
白噪声功率谱密度与带限白噪声自相关
§2.6
正弦波加窄带高斯过程

为什么研究这个?

通信接收端收到的信号 = 有用的正弦信号 + 信道中的窄带高斯噪声。分析这个混合信号的统计特性,直接关系到接收机能正确检测信号的概率(误码率)。

PDF2.6 正弦波加窄带高斯过程p.30

pdf/通信原理/第二章.pdf · p.30

打开原文

模型

设合成信号:

$$r(t) = A\cos(\omega_c t + \theta) + n(t)$$

其中 $A$ 为信号振幅(常数),$\theta$$(0, 2\pi)$ 上均匀分布,$n(t)$ 是均值为零、方差为 $\sigma_n^2$ 的窄带平稳高斯过程。

包络服从莱斯分布

合成信号的包络 $z(t)$ 的概率密度函数为:

$$f(z) = \frac{z}{\sigma_n^2} \exp\left(-\frac{z^2 + A^2}{2\sigma_n^2}\right) I_0\left(\frac{Az}{\sigma_n^2}\right), \quad z \geq 0$$

其中 $I_0(\cdot)$ 是第一类零阶修正贝塞尔函数。这个分布称为广义瑞利分布莱斯分布

两种极端情况

条件退化分布物理含义
$A = 0$(无信号)瑞利分布纯噪声
$A \gg \sigma_n$(强信号)近似高斯分布信号主导

直觉:想象在二维平面上,噪声把信号点"模糊"成一个云团。没有信号时($A=0$),云团中心在原点,距离服从瑞利分布;有信号时,云团中心偏移到 $(A, 0)$,距离服从莱斯分布。信号越强,云团越集中,接收端越容易区分信号和噪声。

§2.7
随机过程通过线性系统

问题

随机信号 $X_i(t)$ 通过冲激响应为 $h(t)$ 的线性时不变系统,输出 $X_o(t)$ 的统计特性如何?

PDF2.7 随机过程通过线性系统p.35

pdf/通信原理/第二章.pdf · p.35

打开原文

输入输出关系

$$X_o(t) = \int_{-\infty}^{\infty} h(\tau) X_i(t - \tau) d\tau$$

三条核心结论

(1)输出仍为平稳过程

若输入 $X_i(t)$ 是平稳的,则输出 $X_o(t)$ 也是平稳的,且:

  • 输出均值:$E[X_o(t)] = E[X_i(t)] \cdot H(0)$,其中 $H(0) = \int h(\tau) d\tau$ 是系统的直流增益
  • 输出自相关:$R_o(\tau) = h(\tau) * h(-\tau) * R_i(\tau)$

(2)功率谱密度的关系

$$\boxed{P_o(\omega) = |H(\omega)|^2 \cdot P_i(\omega)}$$

输出功率谱 = 输入功率谱 × 系统功率传输函数

这是整章最实用的公式之一。它在频域中直观地表明:系统像一个"频率选择器",对每个频率分量按 $|H(\omega)|^2$ 进行加权。

(3)高斯过程经线性系统仍为高斯

这与 2.3 节中高斯过程的第四条性质呼应:高斯性在通过线性系统后保持不变。这意味着白噪声通过滤波器后,虽然功率谱变了,但分布类型仍是高斯的。

速查
复习速查表
概念核心公式/性质关键词
随机过程固定 $t$ → 随机变量;固定 $\xi$ → 样本函数双重性
均值$a(t) = E[X(t)]$一阶矩
自相关$R(t_1, t_2) = E[X(t_1)X(t_2)]$二阶矩
广义平稳$E[X(t)] = a$$R = R(\tau)$均值常数 + 相关只与 $\tau$ 有关
各态历经时间平均 = 统计平均一条样本即可
维纳-辛钦$R(\tau) \leftrightarrow S(\omega)$时域频域桥梁
高斯过程不相关⇔独立;线性变换后仍高斯最重要的噪声模型
窄带过程包络 → 瑞利分布;相位 → 均匀分布I/Q 两路独立高斯
白噪声$P(\omega) = n_0/2$$R(\tau) = \frac{n_0}{2}\delta(\tau)$频域平坦 + 时域冲击
莱斯分布正弦波 + 窄带高斯 → 包络服从莱斯$A=0$ 退化为瑞利
线性系统$P_o(\omega) = |H(\omega)|^2 \cdot P_i(\omega)$频域加权
下一步:第三章将利用本章建立的功率谱密度、噪声模型等工具,分析信道对信号的影响,并引出香农公式。