基于OQAM调制的FBMC与OFDM系统仿真及性能对比分析
项目介绍
本项目设计并实现了一个完整的通信系统物理层仿真链路,重点在于对比分析FBMC-OQAM(基于偏移正交幅度调制的滤波器组多载波)与传统OFDM(正交频分复用)技术的性能。作为5G及未来6G移动通信的候选波形,FBMC通过移除循环前缀(CP)并利用具有良好时频聚焦特性的原型滤波器,旨在提高频谱效率并降低带外辐射。
本仿真平台完全基于MATLAB构建,通过高度模块化的代码实现了从比特生成、调制、波形合成、信道传输到接收解调的完整过程,并提供了多维度的可视化分析工具,用于评估两种波形在频谱特性、误码率(BER)及峰值平均功率比(PAPR)方面的表现。
功能特性
- 双系统并行仿真:同时构建了标准OFDM系统和FBMC-OQAM系统,确保在相同参数配置(如子载波数、调制阶数)下进行公平对比。
- 先进的波形设计:
*
OFDM:支持可配置的循环前缀(CP)插入与去除,采用标准IFFT/FFT处理流程。
*
FBMC:实现了基于PHYDYAS原型滤波器的多载波调制。包含OQAM预处理,严格遵循实部/虚部的时间交错(Staggering)机制,保证复数域的正交性。
- 信道环境:目前主要针对加性高斯白噪声(AWGN)信道环境进行性能评估。
- 全面的性能分析指标:
*
频谱分析:计算并对比功率谱密度(PSD),展示FBMC在旁瓣抑制方面的优势。
*
可靠性分析:通过蒙特卡洛模拟计算误码率(BER),在0-20dB信噪比范围内绘制性能曲线。
*
PAPR分析:统计信号的峰值平均功率比,生成互补累积分布函数(CCDF)曲线。
- 可视化展示:自动生成包含PSD对比、BER曲线、接收星座图及PAPR分布的综合图表。
算法与实现细节
本节详细说明代码中的核心逻辑与算法实现,基于实际程序内容。
1. 系统参数配置
程序首先初始化关键通信参数,定义了子载波数量(N=64)、重叠因子(K=4,用于FBMC滤波器设计)、调制阶数(16-QAM)以及仿真符号数。特别地,对于FBMC系统,代码加载了PHYDYAS原型滤波器的频率响应系数(k=4的标准系数:P0至P3),这是实现低带外辐射的关键。
2. 误码率(BER)仿真链路
主程序包含一个SNR循环(0dB至20dB),在每个信噪比点执行以下操作:
- 比特生成:生成足够长度的随机比特流。
- OFDM链路:执行QAM映射、IFFT变换、添加循环前缀、AWGN信道加噪、去除CP、FFT变换及QAM解调。
- FBMC链路:执行QAM映射、OQAM交错调制(将复数符号分解为时间偏移的实数值)、通过综合滤波器组、信道加噪、通过分析滤波器组及OQAM解调。
- 误差统计:比较发射与接收比特,由于FBMC滤波器组存在固有延迟,代码中包含对齐截断处理逻辑,计算有效数据段的误码率。
3. FBMC-OQAM 核心处理逻辑
代码中详细实现了OQAM的调制预处理(Staggering):
- QAM分解:将标准的QAM复数符号分解为实部和虚部。
- 时间交错:为了满足正交性,对实部和虚部在时间轴上进行偏移。
- 相位规则:根据子载波索引(奇/偶)应用特定的相位规则。偶数子载波按照[实部, 虚部]排列,奇数子载波按照[虚部, 实部]排列,并结合相位旋转因子(j的幂次)以维持正交性。
- 多相相干:利用重叠因子K和PHYDYAS系数构建滤波器组结构。
4. 信号特性分析与可视化
为了直观展示性能差异,代码在仿真结束后运行一次高信噪比(30dB)的单次传输,并进行特殊的统计分析:
- PSD分析:使用Welch方法计算发射信号的功率谱密度,并归一化频率轴。通过绘图清晰展示OFDM较高的带外旁瓣与FBMC极低的带外泄漏。
- PAPR统计:通过多次随机试验运行仅发射机的子函数,计算信号瞬时功率与平均功率的比值,统计其超过特定阈值的概率(CCDF)。
- 星座图:绘制接收端的复数信号散点图,直观反映信号在经过信道和解调后的聚焦程度。
系统要求
- MATLAB R2018b 或更高版本。
- DSP System Toolbox(用于
pwelch等频谱分析函数)。 - Communications Toolbox(用于
qammod、qamdemod、biterr、awgn等基础通信函数)。
使用方法
- 确保MATLAB环境已安装上述工具箱。
- 直接运行主程序脚本。
- 程序将输出初始化信息、实时打印不同SNR下的BER统计结果。
- 运行结束后,会自动弹出一个包含四个子图的综合性能对比窗口:
* 左上:功率谱密度(Power Spectral Density)对比。
* 右上:误码率(BER)随SNR变化的曲线。
* 左下:接收信号星座图。
* 右下:PAPR的互补累积分布函数(CCDF)。