基于裁剪DFT扩展FBMC系统的PAPR抑制仿真平台

资 源 简 介

本项目是一个基于MATLAB环境开发的物理层仿真工具，旨在通过裁剪离散傅里叶变换扩展（Pruned DFT-spreading）技术解决滤波器组多载波（FBMC）系统固有的高峰均功率比（PAPR）难题。在传统FBMC系统中，多子载波的相干叠加会导致信号在时域产生巨大的峰值，从而降低射频功放的线性效率并增加非线性失真。本项目实现的Pruned DFT-s技术通过对DFT扩展过程进行逻辑优化，根据子载波映射规律去除多余的计算分支，不仅在数学上等效于一种特殊的预编码结构以降低PAPR，还显著减少了发射端的运算开

详 情 说 明

项目介绍

本项目是一个基于 MATLAB 环境开发的物理层仿真工具，专注于研究和评估裁剪离散傅里叶变换扩展（Pruned DFT-spreading）技术在滤波器组多载波（FBMC）系统中的应用。

在现代无线通信中，FBMC 系统因其极低的带外辐射和较高的频谱效率被视为 5G 及未来通信的关键候选波形。然而，与 OFDM 类似，FBMC 同样面临高峰均功率比（PAPR）的挑战，这会导致射频功放效率下降。本项目通过引入 Pruned DFT-s 技术，在发射端对数据进行基于子块的 DFT 预编码扩展。这种方法通过逻辑优化去除了冗余计算，在降低信号峰值功率的同时，保持了较低的计算复杂度。

本仿真平台提供了一个标准化的对比环境，能够量化分析标准 FBMC、经典 DFT-s-FBMC 以及 Pruned-DFT-s-FBMC 在 PAPR、功率谱密度（PSD）以及误码率（BER）等多个维度的性能表现。

功能特性

1. 全链路仿真：实现了从随机比特生成、16-QAM 调制、OQAM 映射到 FBMC 调制与解调的完整物理层链路。
2. 多方案对比：在一套代码中集成了三种波形方案（Standard/Classic DFT-s/Pruned DFT-s），方便直接进行性能基准测试。
3. PAPR 统计分析：通过大量的蒙特卡洛实验，自动生成互补累积分布函数（CCDF）曲线，直观展示不同技术对峰值功率的抑制效果。
4. 频谱特性评估：利用周期图法计算并绘制发射信号的功率谱密度，用于分析不同波形的带外泄露情况。
5. 鲁棒性验证：在加性高斯白噪声（AWGN）信道下进行位错误率（BER）仿真，验证引入预编码技术后对系统检测性能的影响。
6. 高度参数化：支持自定义子载波数量（M）、重叠因子（K）以及裁剪粒度（Prune Factor）等关键系统参数。

系统要求

1. 软件环境：MATLAB R2018b 或更高版本。
2. 必备工具箱：

- Communications Toolbox - Signal Processing Toolbox

使用方法

1. 准备环境：确保 MATLAB 已安装上述工具箱。
2. 运行仿真：在 MATLAB 指令窗口中直接运行仿真主程序脚本。
3. 结果查看：程序运行完成后，将自动弹出四张分析图表：

- PAPR CCDF 性能对比图 - 发射信号功率谱密度（PSD）对比图 - AWGN 信道下的 BER 曲线图 - 时域信号瞬时功率包络图

1. 日志输出：实时查看 MATLAB 控制台输出的仿真进度，包括 PAPR、PSD 和 BER 阶段的执行提示。

核心功能与逻辑实现

程序逻辑严格遵循数字信号处理流程，分为参数初始化、PAPR 统计循环、频谱计算、BER 遍历仿真以及多维度绘图五个主要部分。

1. 信号调制逻辑实现

系统采用 PHYDYAS 原型滤波器，基于重叠因子 K=4 构建滤波器系数。调制过程遵循 OQAM 逻辑，将输入的 QAM 复数符号拆分为实部和虚部，并交替在时域进行平移。每个子载波上的符号通过多相网络（PPN）进行滤波。

1. 三种预编码方案的实现细节

• Standard FBMC：数据直接映射到活跃子载波进行调制。
• Classic DFT-s-FBMC：在调制前，对整个活跃子载波块执行 N 点 FFT 变换，实现时域信号在频域的均匀扩展，以降低峰值。
• Pruned-DFT-s-FBMC：将活跃子载波进一步划分为更小的子块（Sub-blocks），并在每个子块内执行裁剪规模的 DFT 变换。这种方式通过降低子载波间的相关性，进一步优化了 PAPR 表现。

1. PAPR 计算逻辑

程序通过对调制后的时域信号序列进行统计，计算信号的瞬时峰值功率与平均功率的比值。实验通过设置 500 次迭代来实现蒙特卡洛统计，最终计算出 PAPR 超过特定阈值的概率。

1. 误码率（BER）仿真流程

在不同的信噪比（SNR）条件下，程序将生成的比特流通过完整的调制、信道（AWGN）及解调链路。在接收端，针对 Pruned 方案执行对应的逆裁剪 DFT（IDFT）还原数据，最后通过判定比特错误数来计算 BER。

算法与关键函数分析

1. PHYDYAS 滤波器设计

代码中使用了经典的 PHYDYAS 频率响应系数（H1=0.9719, H2=0.7071, H3=0.2351）构建原型滤波器。通过公式计算时域脉冲响应，并进行能量归一化处理，确保了子载波间的正交性和良好的频域封闭性。

1. FBMC 调制函数 (多相网络 PPN 结构)

该函数通过 IFFT 变换和周期扩展实现多载波合成。关键步骤包括：

• 子载波映射：将活跃载波放置在频谱中心以模拟护带。
• 相位旋转：实现 OQAM 所需的 pi/2 相位偏移。
• 重叠相加：通过滑动窗口机制，将滤波后的符号在时域进行 K 倍重叠叠加。

1. FBMC 解调函数 (分析滤波器组 AFB)

作为调制的逆过程，该函数通过匹配滤波和 PPN 叠加汇总，恢复各个子载波的分量，并通过取实部的方式完成 OQAM 信号的判决。

1. 裁剪变换 (Pruning Logic)

在主程序逻辑中直接嵌入了裁剪逻辑。通过设置 prune_factor 参数，将数据矩阵切片，分别对每个切片进行 FFT 运算。这一逻辑实现了预编码阶段的计算优化。

1. 辅助功能函数

包含了计算比特错误和二进制转十进制的逻辑，用于支持 BER 仿真的数据转换。