用MATLAB和Pluto SDR复现通信原理实验：正弦波、方波收发实测与波形畸变分析

基于Pluto SDR的通信原理实验：波形传输畸变分析与MATLAB实战指南

实验背景与硬件准备

在通信工程的教学体系中，理论知识与实践操作往往存在一道难以跨越的鸿沟。当学生在课本上学习奈奎斯特采样定理时，很难直观理解"为什么方波传输会出现吉布斯现象"；当讨论带限信道特性时，也无法亲眼目睹高频分量被滤除后的波形变化。这正是我们引入Pluto SDR硬件平台的意义所在——它将抽象的理论公式转化为可视化的波形对比，让通信原理变得触手可及。

Pluto SDR（Software Defined Radio）是一款由Analog Devices推出的低成本软件定义无线电设备，其核心优势在于：

• 硬件参数：支持70MHz至6GHz频段，20MHz瞬时带宽，12位ADC/DAC
• 开发接口：通过USB 3.0接口与MATLAB/Simulink、Python等环境直接交互
• 教学适用性：体积仅信用卡大小，适合实验室多台配置

实验前需要准备的环境包括：

1. Pluto SDR设备（建议使用原厂固件v0.31以上版本）
2. MATLAB R2020b及以上（需安装Communications Toolbox和Support Package for Pluto SDR）
3. 50Ω阻抗匹配的天线或同轴电缆直连方案

提示：为避免2.4GHz WiFi干扰，实验频率建议选择2.5GHz以上频段，如2.6GHz或3.5GHz

正弦波传输实验与信道响应分析

基础正弦信号生成

我们从最简单的单频正弦波开始，建立对SDR系统的直观认识。在MATLAB中生成1MHz正弦波的典型代码如下：

fs = 1e6; % 采样率1MHz f0 = 100e3; % 信号频率100kHz t = 0:1/fs:1e-3; % 1ms时间序列 tx_signal = 0.9*sin(2*pi*f0*t); % 生成正弦波（幅度略小于1防止DAC饱和）

参数配置对传输质量的影响可以通过下表对比：

接收信号处理与畸变检测

接收端需要特别关注自动增益控制(AGC)带来的影响。建议采用以下处理流程：

1. 初始帧丢弃：前3-5帧数据因AGC收敛不稳定，应舍弃
2. 幅度归一化：根据接收信号峰值动态调整显示范围
3. 相位对齐：通过互相关算法补偿传输延迟

rx = sdrrx('Pluto','SamplesPerFrame',4000); for i=1:5 % 丢弃前5帧用于AGC稳定 rx(); end valid_data = rx(); % 获取有效数据 % 时域对齐处理 [corr,lags] = xcorr(valid_data, tx_signal); [~,idx] = max(abs(corr)); delay = lags(idx); aligned_rx = valid_data(delay+1:delay+length(tx_signal));

典型正弦波传输结果可能呈现以下特征：

• 包络波动：反映信道多径效应或天线阻抗失配
• 相位抖动：本振频率偏移或时钟不同步导致
• 谐波失真：发射机非线性引入的二次/三次谐波

方波传输与吉布斯现象实证

方波谐波特性分析

方波在数学上可表示为无限奇次谐波的叠加：

方波 = 4/π * (sin(ωt) + 1/3 sin(3ωt) + 1/5 sin(5ωt) + ...)

当通过有限带宽信道传输时，高频谐波被滤除，导致波形出现以下典型畸变：

1. 上升沿/下降沿变缓：缺失高频分量使跳变不再陡峭
2. 过冲振荡：吉布斯现象在跳变处产生约9%的幅度震荡
3. 平台波动：剩余谐波的相位干扰导致平顶部分不平坦

MATLAB实现与参数优化

生成5周期方波的代码需要特别注意占空比控制：

f_square = 50e3; % 方波频率50kHz duty_cycle = 50; % 占空比百分比 t_square = 0:1/fs:100e-6; % 覆盖5个周期 tx_square = 0.9*square(2*pi*f_square*t_square, duty_cycle);

接收端观测到的畸变程度与系统带宽直接相关，可通过调整采样率验证：

注意：实际可用带宽受Pluto硬件限制，最大瞬时带宽约20MHz

多波形对比与教学演示方案

综合实验设计

为全面展示不同波形的传输特性，建议按以下顺序进行实验演示：

1. 单频正弦波：验证系统基本功能，测量信噪比
2. 双音测试：使用两个相近频率正弦波检测非线性失真
3. 三角波：观察中等谐波含量信号的传输特性
4. 锯齿波：分析非对称波形的畸变模式
5. 方波：极端案例展示吉布斯现象

课堂互动环节设计

可引导学生思考以下问题并通过实验验证：

• 为什么方波接收波形会出现振铃效应？
• 如何通过调整采样率改善锯齿波的阶梯状畸变？
• 对比实际测量结果与傅里叶级数理论预测的差异

以下MATLAB代码可同时生成多种波形供对比：

waveforms = struct(); waveforms.sine = sin(2*pi*50e3*t); waveforms.triangle = sawtooth(2*pi*50e3*t, 0.5); waveforms.sawtooth = sawtooth(2*pi*50e3*t); waveforms.square = square(2*pi*50e3*t); figure; for i = 1:4 subplot(2,2,i); plot(t(1:200), waveforms.(fields(i))(1:200)); title(fieldnames(waveforms){i}); end

高级实验拓展方向

信道模拟与失真补偿

在基础实验之上，可以引入更复杂的信道条件：

1. 多径环境模拟：

   delayed_wave = [zeros(1,50), tx_signal(1:end-50)]; multipath_signal = tx_signal + 0.3*delayed_wave;

2. 加性高斯白噪声：

   SNR = 15; % dB noisy_signal = awgn(tx_signal, SNR, 'measured');

实时频谱监测方案

利用MATLAB的DSP System Toolbox实现实时频谱分析：

spectrumAnalyzer = dsp.SpectrumAnalyzer(... 'SampleRate', fs, ... 'FrequencySpan', 'Span and center frequency', ... 'Span', 500e3, ... 'CenterFrequency', 0); while true rx_data = rx(); spectrumAnalyzer(rx_data); end

自动测量指标计算

编写函数量化波形失真程度：

function [thd, enob] = analyze_waveform(signal, fundamental_freq, fs) % 计算总谐波失真(THD) thd = thd(signal, fs, 5); % 估算有效位数(ENOB) snr_val = snr(signal, fs); enob = (snr_val - 1.76)/6.02; end

实验中发现Pluto SDR在2.5GHz频段典型性能为：

• THD：-45dBc（正弦波，0dB发射增益）
• ENOB：约9.5位（受12位ADC限制）