随机信号处理实验报告：从原理到实践，避坑指南

在华东理工大学的随机信号处理实验中，经常遇到理论与实践脱节的问题。很多同学仅仅停留在公式推导层面，缺乏对随机信号本质的理解，导致实验结果与预期偏差较大。本文将从实验原理、代码实现、以及常见问题三个方面进行深入剖析，帮助大家更好地完成随机信号处理实验。

随机信号基础理论回顾

随机过程的概念

随机过程可以看作是一系列随机变量的集合，这些随机变量随着时间或其他参数变化而变化。它不同于确定性信号，其每一个时刻的值都具有不确定性，需要用概率统计的方法进行描述。理解随机过程的关键在于掌握其统计特性，例如均值、方差、自相关函数等。在实验中，我们需要通过采样获取离散的随机信号，并利用这些样本估计其统计特性。

功率谱密度 (PSD)

功率谱密度描述了随机信号功率在不同频率上的分布情况。它是随机信号分析的重要工具，可以用于识别信号中的主要频率成分，并进行信号的滤波和降噪。计算 PSD 的常用方法是 Welch 估计法，该方法通过对信号进行分段、加窗、求周期图，然后进行平均，从而降低估计的方差。在实际实验中，窗函数的选择会影响 PSD 的分辨率和方差，需要根据具体情况进行权衡。常用的窗函数包括矩形窗、汉明窗、汉宁窗等。

实验中常见的随机信号

• 高斯白噪声： 具有均匀功率谱密度的随机信号，在各个频率上的功率相等。它是一种理想化的噪声模型，常用于模拟实际环境中的背景噪声。
• 正弦信号加噪声： 模拟实际通信系统中接收到的信号，其中正弦信号代表有用信号，噪声代表干扰。
• 语音信号： 语音信号具有时变性和非平稳性，是一种典型的随机信号。对其进行分析可以提取出语音的特征，用于语音识别、语音合成等应用。

基于 Python 的随机信号处理实验实现

实验环境搭建

推荐使用 Anaconda 管理 Python 环境，并安装以下常用库：

conda create -n random_signal python=3.8
conda activate random_signal
pip install numpy scipy matplotlib

代码示例：生成高斯白噪声并计算 PSD

import numpy as np
import matplotlib.pyplot as plt
from scipy import signal

# 设置参数
fs = 1000  # 采样率
t = np.arange(0, 1, 1/fs)  # 时间向量
n = len(t)  # 样本数

# 生成高斯白噪声
noise = np.random.normal(0, 1, n)  # 均值为 0，方差为 1 的高斯白噪声

# 计算 PSD (Welch 估计)
f, Pxx = signal.welch(noise, fs, nperseg=256)  # nperseg：每个段的长度

# 绘制 PSD
plt.figure()
plt.semilogy(f, Pxx)  # 使用对数坐标显示功率谱
plt.xlabel('Frequency [Hz]')
plt.ylabel('PSD [V**2/Hz]')
plt.title('Power Spectral Density of Gaussian White Noise')
plt.grid()
plt.show()

代码示例：正弦信号加噪声的滤波

import numpy as np
import matplotlib.pyplot as plt
from scipy import signal

# 设置参数
fs = 1000  # 采样率
t = np.arange(0, 1, 1/fs)  # 时间向量
f0 = 50  # 正弦信号频率
A = 1  # 正弦信号幅度

# 生成正弦信号
signal_pure = A * np.sin(2 * np.pi * f0 * t)

# 添加高斯白噪声
noise = 0.5 * np.random.normal(0, 1, len(t))
signal_noisy = signal_pure + noise

# 设计 Butterworth 带通滤波器
cutoff = [40, 60]  # 截止频率
order = 4  # 滤波器阶数
b, a = signal.butter(order, cutoff, btype='bandpass', fs=fs)  # 计算滤波器系数

# 应用滤波器
signal_filtered = signal.lfilter(b, a, signal_noisy)  # 使用 lfilter 进行滤波

# 绘制结果
plt.figure(figsize=(12, 8))
plt.subplot(3, 1, 1)
plt.plot(t, signal_pure)
plt.title('Pure Sine Wave')
plt.subplot(3, 1, 2)
plt.plot(t, signal_noisy)
plt.title('Noisy Sine Wave')
plt.subplot(3, 1, 3)
plt.plot(t, signal_filtered)
plt.title('Filtered Sine Wave')
plt.tight_layout()
plt.show()

实验避坑经验总结

• 数据预处理： 在进行信号处理之前，务必对数据进行预处理，例如去除直流分量、归一化等。避免因数据量纲不同导致计算错误。
• 参数选择： 不同的实验需要选择合适的参数，例如窗函数的类型和长度、滤波器的阶数和截止频率等。参数选择不当可能导致实验结果不理想。
• 结果分析： 对实验结果进行仔细分析，理解其物理意义，并与理论预期进行对比。如果结果与预期偏差较大，需要仔细检查实验步骤和代码实现。
• 版本控制: 使用 Git 进行版本控制，避免代码丢失或混乱。同时，注释要清晰明了，方便自己和他人理解。

通过以上分析，相信大家可以更好地完成华东理工大学的随机信号处理实验，并对随机信号处理的原理和应用有更深入的理解。随机信号处理实验报告的撰写也会更加得心应手。