【MATLAB源码-第452期】基于MATLAB的CPFSK与GFSK调制的多解调算法误码率、频偏与相偏鲁棒性对比仿真

【MATLAB源码-第452期】基于MATLAB的CPFSK与GFSK调制的多解调算法误码率、频偏与相偏鲁棒性对比仿真 操作环境MATLAB 2024a1、算法描述摘要连续相位频移键控是一类 具有相位连续特性的数字调制方式。与普通相位突变型调制相比它在包络稳定性、频谱控制和功率放大器适应性方面具有一定优势。本文围绕CPFSK与GFSK两类连续相位调制信号构建了一个基于MATLAB的调制、信道、解调和性能评估仿真系统。系统主要完成CPFSK与GFSK调制波形生成、频谱估计、AWGN信道与瑞利衰落信道下的误码率统计以及频偏和相位偏差条件下的鲁棒性分析。解调算法方面本文对鉴频解调、差分解调、相干解调和Viterbi解调四类方法进行对比。系统通过多次Monte Carlo随机仿真输出不同信噪比、不同偏差强度下的BER曲线并保存时域波形、频谱、误码率和单次解调统计图像。仿真结果能够直观反映不同解调算法在复杂信道和同步偏差条件下的性能差异。整体来看该系统结构清晰参数集中图像输出完整适合作为连续相位调制解调算法性能分析和教学演示的仿真平台。CPM信号具有连续相位和记忆特性这一点也与MATLAB通信工具箱对连续相位调制的定义一致。关键词 CPFSKGFSK连续相位调制鉴频解调差分解调相干解调Viterbi解调误码率频偏鲁棒性相位偏差1 引言数字通信 系统的核心任务是将离散比特信息转换为适合信道传输的波形并在接收端尽可能准确地恢复原始比特。调制方式直接决定信号带宽、抗噪声能力、接收机复杂度和硬件实现难度。在低功耗无线通信、短距离通信和物联网通信场景中频移键控及其连续相位形式仍然具有重要价值。其主要原因是这类信号通常具有恒包络特征对功率放大器线性度要求相对较低适合功耗受限和成本敏感的终端设备。CPFSK是连续相位频移键控的基础形式。它通过控制瞬时频率的正负偏移来表示二进制信息同时保持相位连续。GFSK则在CPFSK基础上引入高斯滤波器对基带符号进行平滑处理使频率变化更加缓和从而改善频谱旁瓣扩展。Bluetooth Core Specification中也采用GFSK作为其物理层调制相关内容之一说明该类调制方式在实际短距离通信系统中仍具有工程意义。单纯讨论调制方式并不完整。接收端解调算法会显著影响系统误码性能。对于连续相位调制信号常见的解调方法包括鉴频解调、差分解调、相干解调和基于最大似然思想的Viterbi解调。鉴频解调结构简单适合快速实现但抗噪声能力有限。差分解调不依赖绝对相位具有一定抗相位偏差能力。相干解调利用参考波形相关判决理论性能较好但对同步条件要求更高。Viterbi解调利用状态路径搜索思想能够考虑连续相位调制的记忆 特性但计算复杂度高于普通判决方法。本文所研究的MATLAB系统正是围绕上述问题展开。系统不是只画一条BER曲线而是同时覆盖调制波形、频谱特征、AWGN信道、瑞利信道、频偏、相偏和多算法解调对比。这样的仿真结构比单一调制演示更有实际价值也更适合作为通信系统课程设计或算法性能分析项目。2 系统模型系统采用二进制随机比特作为输入。比特序列首先经过符号映射将0和1分别映射为-1和1。对于CPFSK调制系统按照调制指数对符号序列进行过采样并通过相位累加生成连续相位信号。其关键思想是当前输出相位不是孤立决定的而是由历史符号共同累积形成。因此CPFSK信号具有记忆性。对于GFSK调制系统在符号过采样之后增加高斯频率脉冲成形环节。高斯滤波器由BT积和脉冲跨度控制。BT积越小滤波越强频率变化越平滑频谱旁瓣通常会被压低但符号间相关性也会增强。系统中设置BT为0.3调制指数为0.5每比特采样点数为8。这组参数能够体现GFSK相对CPFSK的频谱平滑特点也便于在时域和频域图像中观察差异。信道方面系统包含AWGN信道和瑞利衰落信道。AWGN信道主要用于模拟理想加性白噪声环境。MATLAB通信文档中也将AWGN信道描述为向通过信道的信号添加白高斯噪声的模型。 瑞利信道用于模拟多径衰落下幅度和相位随机变化的影响。项目代码中对瑞利信道采用每比特一个平坦衰落系数并在接收端设置理想均衡。这样做的目的不是完整模拟实际接收机的信道估计误差而是突出衰落环境对不同解调算法误码率的影响。偏差模型方面系统分别设置归一化频偏和固定相位偏差。频偏通过复指数旋转叠加到基带信号上相位偏差则通过固定相位旋转实现。之后再叠加AWGN噪声并比较四种解调算法在偏差条件下的BER变化。这一部分对工程应用有意义因为实际通信系统中本振误差、晶振偏差、移动环境和同步残差都会导致载波频率或相位偏差。3 调制算法设计CPFSK调制模块的核心是相位连续累加。系统先将二进制比特转换为双极性符号再按照过采样倍数扩展为采样序列。每个采样点对相位的贡献由调制指数和符号极性决定。最后通过复指数函数生成单位幅度复基带信号。由于输出信号幅度恒定信息主要承载在相位变化和瞬时频率变化中。GFSK调制模块在CPFSK基础上增加高斯滤波。系统先生成高斯脉冲再将过采样符号序列与该脉冲卷积。卷积后的信号不再是突变的矩形频率脉冲而是经过平滑处理的频率控制序列。随后同样通过相位累加生成复基带信号。这样可以降低频谱扩展使信号在频域上更加集中。从工程角度看CPFSK实现简单瞬时频率变化直接适合用于理解连续相位调制基本原理。GFSK通过高斯滤波改善频谱表现更接近实际低功耗无线通信中的波形选择。研究中同时保留CPFSK和GFSK有利于学生或工程人员理解“频率平滑”和“频谱控制”之间的关系。4 解调算法设计系统实现了四类解调算法。第一类是鉴频解调。该方法计算相邻采样点之间的相位差并将相位差累加作为判决度量。由于频移键控的信息体现在瞬时频率正负变化上因此相位差可以近似反映频率变化方向。鉴频解调实现简单不需要复杂状态搜索适合低复杂度接收机。但它对噪声较敏感特别是在低信噪比或衰落环境下瞬时相位差容易出现较大扰动。第二类是差分解调。该方法利用一个比特周期内首尾采样点之间的相位差进行判决。相较于普通鉴频方法差分解调不强依赖绝对载波相位因此对固定相位偏差具有一定适应性。但是如果存在较强噪声、频偏或相位旋转累积差分判决仍可能出现明显误码。第三类是相干解调。系统为当前相位状态下的比特0和比特1分别构造参考波形并计算接收片段与参考波形的相关值。相关值较大的一支被判为发送比特。相干解调充分利用波形匹配信息理论上可获得更好的性能。但它需要较准确的相位状态跟踪。如果同步误差较大相干参考波形与接收信号失配性能会下降。第四类是Viterbi解调。连续相位调制信号存在相位记忆当前符号不仅影响当前波形也影响后续相位状态。Viterbi算法将相位状态离散化通过路径度量递推寻找最优状态序列。系统中每个分支度量由接收片段与候选参考波形的相关输出得到再通过动态规划更新路径。最终从最优终止状态反向回溯得到比特序列。Viterbi方法通常性能较强但需要保存前驱状态和分支判决计算量和存储量高于鉴频、差分和普通相干解调。这四类算法形成了较清晰的复杂度梯度。鉴频和差分属于低复杂度方法相干解调属于中等复杂度方法Viterbi解调属于更高性能但复杂度更高的方法。系统将它们放在同一仿真框架中对比是合理的。5 仿真设计系统主程序采用集中式流程。首先设置随机种子保证仿真结果可复现。然后读取系统参数创建图像和结果输出文件夹 。接着生成短比特序列分别通过CPFSK和GFSK调制生成波形并绘制时域实部和归一化瞬时频率曲线。该图主要用于观察两种调制方式在时域上的差异。第二步系统生成较长比特序列用于频谱估计。频谱估计采用分段加窗和FFT平均思想。系统使用类汉宁窗降低频谱泄漏并通过fftshift将频谱中心移动到零频附近。输出图像用于比较CPFSK和GFSK的功率谱分布。第三步系统分别在AWGN信道下测试CPFSK和GFSK并对四类解调算法统计BER。信噪比设置为0 dB到12 dB步进为2 dB。每个信噪比点进行多次Monte Carlo随机仿真每次生成3000个随机比特。这样能够降低单次随机结果带来的偶然性。第四步系统在瑞利信道下测试CPFSK的四类解调算法性能。由于瑞利衰落会引入更明显的幅度和相位扰动所以该部分通常比AWGN场景更能体现算法抗信道变化能力。第五步系统固定Eb/N0为10 dB分别测试频偏和相位偏差对CPFSK解调性能的影响。频偏范围为0到0.2步进为0.025。相位偏差范围为0度到120度步进为15度。该设计能够较直观地展示算法对同步偏差的敏感程度。第六步系统生成单次解调统计图。该图包括接收波形、发送比特、四类算法判决序列和误码柱状图。它不是单纯给出BER曲线而是进一步展示具体判决差异有利于理解误码产生的位置和算法输出差别。6 结果分析从时域波形角度看CPFSK和GFSK都保持连续相位特征但二者瞬时频率变化方式不同。CPFSK的频率变化更加直接符号切换位置的频率变化较明显。GFSK由于经过高斯滤波瞬时频率曲线更加平滑。这种平滑会降低频谱突变带来的高频分量因此GFSK通常具有更好的频谱集中性。从频谱图看GFSK相较CPFSK应表现出更明显的旁瓣抑制趋势。这个结论与GFSK引入高斯滤波的基本原理一致也与关于GFSK在AWGN信道中进行性能分析的研究方向相符。已有研究也常从高斯滤波、功率谱密度和BT积角度分析GFSK性能。从AWGN信道下的BER曲线看随着Eb/N0升高各类解调算法的误码率一般会下降。这是数字通信仿真的基本规律信噪比越高噪声扰动越小接收端判决越可靠。四类算法之间的差异主要来自判决信息利用程度不同。鉴频解调主要依靠相邻采样点相位差信息利用较简单差分解调依靠符号首尾相位差对绝对相位不敏感相干解调利用参考波形相关能够更充分利用波形结构Viterbi解调进一步利用状态记忆适合连续相位调制的路径搜索。从瑞利信道结果看误码率通常会明显劣于AWGN场景。即使代码中加入了理想均衡衰落信道下的噪声增强和相位扰动仍会影响判决可靠性。低复杂度算法在这种环境下更容易出现性能退化。Viterbi解调由于考虑状态路径理论上更适合处理具有记忆特性的信号但它也不是无条件最优。如果信道估计、同步和噪声模型与算法假设不一致性能仍会受到影响。频偏仿真用于反映载波同步误差对系统的影响。频偏会使接收信号相位随时间持续旋转。对于依赖相位轨迹或参考波形匹配的算法频偏会造成累计失配。相干解调和Viterbi解调虽然理论性能强但如果没有频偏补偿其参考波形会逐渐偏离接收信号。鉴频和差分方法对部分固定相位误差不那么敏感但对持续频率偏移同样会受到影响。相位偏差仿真用于反映固定载波相位误差的影响。固定相位偏差不会像频偏那样随时间持续累积但会改变接收信号整体相位位置。差分解调因为主要使用相对相位变化对固定相位偏差有一定天然适应性。相干解调如果参考相位没有同步则会出现相关度下降。Viterbi解调同样需要参考状态与接收信号相位关系基本一致否则路径度量会被破坏。综合来看本系统的仿真结果可以支持三个判断。第一GFSK相较CPFSK更适合强调频谱平滑和带外抑制的场景。第二低复杂度解调算法适合快速实现但在低信噪比、衰落和同步偏差条件下性能有限。第三Viterbi解调能够利用连续相位调制的状态记忆适合用于性能上限分析但实际工程中需要权衡计算复杂度、同步精度和实现成本。7 系统特点该系统的第一个特点是结构清晰。主程序只负责组织整体流程具体功能被拆分到function文件夹中。调制、信道、解调、绘图、频谱估计和结果输出分别由独立函数完成。这种结构便于阅读、维护和二次修改。第二个特点是参数集中。采样点数、调制指数、BT积、比特数、Monte Carlo次数、信噪比范围、频偏范围、相位偏差范围和输出路径都集中在cfg_sys函数中。后续如果需要调整仿真规模或测试条件只需要修改配置函数不需要在多个脚本中反复查找参数。第三个特点是对比维度较完整。系统不仅比较CPFSK和GFSK还比较四类解调算法不仅分析AWGN信道还加入瑞利信道不仅分析信噪比变化还分析频偏和相位偏差。这种多维对比比单一BER曲线更有说服力。第四个特点是输出结果直观。系统生成8张PNG图像包括波形对比、频谱对比、AWGN误码率、瑞利误码率、频偏影响、相位偏差影响和单次解调统计。图像命名清晰便于直接用于报告、答辩或技术交流。第五个特点是结果可复现。主程序固定随机种子并采用统一参数配置。虽然Monte Carlo仿真本身具有随机性但固定种子能够保证同一环境下重复运行时结果一致方便排查问题和比较算法变化。第六个特点是算法层次清楚。系统从简单鉴频、差分判决到相干相关判决再到Viterbi路径搜索形成了由低复杂度到高复杂度的算法体系。该体系适合用来说明“复杂度提升是否换来性能收益”这一通信系统设计中的核心问题。第七个特点是工程表达明确。系统没有停留在公式层面而是将调制、噪声、衰落、频偏、相偏和误码统计都落实到可运行代码中。对于通信类课程设计或算法演示来说这种完整链路比单独讲某一个函数更有价值。8 结论本文围绕CPFSK与GFSK连续相位调制系统完成了一个多算法解调性能仿真平台的分析。系统以MATLAB为实现工具构建了从比特生成、调制、信道叠加、偏差建模、解调判决到结果绘图和数据保存的完整流程。解调算法覆盖鉴频、差分、相干和Viterbi四种方法能够从不同复杂度层面对连续相位调制信号进行性能比较。从系统功能看该项目能够完成三类关键任务第一比较CPFSK和GFSK在时域波形和频谱上的差异第二比较不同解调算法在AWGN和瑞利信道下的BER表现第三分析频偏和相位偏差对解调性能的影响。从工程价值看该系统适合用于通信系统仿真教学、调制解调算法对比、同步误差影响分析和低功耗无线通信波形研究。最终可以得出结论GFSK通过高斯滤波改善频谱平滑性适合频谱受限场景鉴频和差分解调实现简单适合低复杂度系统相干解调在同步较好时具有较强判决能力Viterbi解调能够利用连续相位调制的记忆特性适合进行高性能接收机分析。实际系统设计中不能只看某一个信噪比点的BER而应同时考虑信道类型、同步误差、复杂度、参数可调性和实现成本。本文系统正好提供了一个较完整的仿真框架用于支撑上述分析。参考文献Bluetooth SIG. Bluetooth Core Specification Version 5.4. Bluetooth Special Interest Group, 31 Jan. 2023.MathWorks. “Continuous-Phase Modulation.” Communications Toolbox Documentation, MathWorks, 2024.MathWorks. “comm.CPMModulator.” Communications Toolbox Documentation, MathWorks, 2024.MathWorks. “AWGN Channel.” Communications Toolbox Documentation, MathWorks, 2024.Zhang, Xiaoyue. “Performance Analysis of GFSK Modulation over AWGN Channel.” Applied and Computational Engineering, 2024.Aldana-López, Rodrigo, Omar Longoria-Gandara, Jose Valencia-Velasco, and Javier Vázquez Castillo. “Orthogonal Space-Time Bluetooth System for IoT Communications.” IoT, vol. 7, no. 1, 2025.“Digital Linear GFSK Demodulator for IoT Devices.” IET Communications, 2018.Forkan, Md., Amit Halder, Mst. Nazmun Nahar Ripa, and Md. Riyad Tanshen. “Performance Analysis of BPSK 8-FSK Modulation Technique through AWGN Channel in Wireless Communication System.” Engineering Research and Applications, AIJR Publisher, 2022.ETSI. Digital Video Broadcasting; Second Generation DVB Interactive Satellite System; Part 2: Lower Layers for Satellite Standard. ETSI EN 301 545-2, version 1.4.1, 2020.Shafi, M. M. Sadman. “MATLAB-Simulated Dataset for Automatic Modulation Classification in Wireless Fading Channels.” Data in Brief, 2026.2、仿真结果演示3、关键代码展示略4、MATLAB 源码获取V点击下方名片关注公众号获取