前几天跟几个学生调Simulink里的调制系统仿真他们一开始盯着MSK和GMSK的数学公式头疼不已一个问“相位连续”到底是怎么个连续法另一个搞不明白“高斯滤波为什么能让频谱变窄”。我说别急着翻书先把仿真模型搭起来看到相位轨迹和频谱图之后这些问题自己就通了。今天这篇就把我在Simulink里做MSK与GMSK调制系统仿真的完整过程写出来从信号模型、模块搭建、参数设置到结果解读最后再聊几个我踩过的坑。这篇文章适合通信专业学生、做物理层算法验证的工程师以及准备通信岗位面试需要快速建立直观理解的朋友。1. 为什么MSK和GMSK是移动通信里绕不开的一对“恒定包络兄弟”1.1 移动通信里调制方式选型的两个硬约束蜂窝移动通信系统从诞生那天起就面对两个很现实的问题第一频率是稀缺资源运营商花大价钱拍下来的频谱必须榨干用尽所以调制信号的带宽要窄、带外泄漏要低第二手机终端是电池供电的射频功放的效率直接决定待机时间和发热所以调制信号最好能让功放工作在高效率状态。这两个约束放在一起会逼着系统工程师做出一个重要取舍。像QAM、QPSK这类线性调制频率效率确实高但它们的包络不是恒定的信号幅度会随着数据变化起伏。带外泄漏和包络起伏这两件事就像跷跷板的两头你想压住带外辐射就得把信号通过带限滤波器结果滤波后的包络会剧烈波动你要让功放高效工作就得让功放工作在接近饱和的区段可是包络起伏的信号一旦进饱和区就会被削顶产生严重失真和邻道干扰。所以工程上需要一类信号包络恒定相位连续频谱相对集中。MSK和GMSK就是这类信号里最有代表性的两个。GSM当年选择GMSK不是没有道理的它能让手机里的功放工作在效率较高的非线性区同时把频谱控制在一个可接受的范围内。正是因为这些实际工程需求MSK和GMSK成了移动通信教科书里绕不开的内容。1.2 它们在调制家族里的位置我把常见数字调制方式的关系整理成一张表方便对照调制方式恒定包络相位连续性频谱集中度典型应用QPSK / QAM否有跳变可能较高Wi-Fi、4G/5G数据信道普通FSK是通常不连续一般低速率遥控、简单数传MSK是连续较好卫星测控、AIS等GMSK是连续好GSM、蓝牙BR、卫星通信从数学上看MSK是连续相位频移键控CPFSK在调制指数( h 0.5 )时的特例也可以看成是OQPSK的一种特殊实现区别在于MSK的基带脉冲不是矩形而是半正弦波形。GMSK则是在MSK的基础上把要发送的基带信号先经过一个高斯低通滤波器再把滤波后的信号送去做连续相位调频。用生活类比可能更好理解MSK像是一辆每次换挡都瞬间完成的车动力传递很快但顿挫感明显GMSK相当于给换挡机构加了一套缓震装置加速过程更平顺油耗频谱占用下降了但动力响应会有点“拖泥带水”这就是码间串扰的来源。1.3 Simulink仿真对理解这些问题有多大帮助公式和文字只能告诉你“GMSK的频谱比MSK更集中”但为什么高斯滤波之后频谱就集中了集中在什么程度BT积从0.5改到0.3又会牺牲多少误码性能这些问题在公式里很难形成直观感受但放到Simulink里跑一遍几分钟就能看到结果。Simulink的价值不在于把模块连起来而在于它是一个可以反复试验的“数字实验台”。你可以把MSK的调制指数改大改小可以切换高斯滤波器BT积可以加噪声观察误码率甚至可以中途换一根连接线从理想信道改成衰落信道所有中间信号都能用各类示波器模块直接查看。这种“动手改参数、立刻看结果”的反馈路径是理解调制系统最有效率的方式也是我在实际项目里验证算法想法时依赖它的原因。2. 建模仿真前先把MSK的相位连续和GMSK的高斯整形吃透2.1 MSK频率差、相位连续和最小移频到底怎么算很多教材把MSK定义为“调制指数为0.5的连续相位FSK”这个定义很简洁但初看容易懵。拆开讲就清楚了。对于二进制频移键控发送频率会在两个频率之间切换一个代表“1”一个代表“0”。一般FSK的频率切换是瞬间完成的相位在切换点可能发生跳变频谱会被拉得很宽。MSK则要求相位在切换点保持连续也就是说信号的时间波形不能有突变的相位拐点。MSK的“最小频率间隔”是一个关键数字。设符号周期为( T_b )调制指数( h 0.5 )则峰值频偏为[ \Delta f \frac{h}{2T_b} \frac{1}{4T_b} ]所以两个发送频率分别为( f_c \frac{1}{4T_b} )和( f_c - \frac{1}{4T_b} )两者之差正好是( \frac{1}{2T_b} )。这个间隔之所以叫“最小”是因为在这个间隔下两个频率的信号在一个符号周期内保持正交而再小的话接收端就很难区分同时相位连续这个条件也会破坏。另外一个理解MSK的重要视角是正交分解。MSK信号可以写成[ s(t) d_I(t)\cos\left(\frac{\pi t}{2T_b}\right)\cos(2\pi f_c t) d_Q(t)\sin\left(\frac{\pi t}{2T_b}\right)\sin(2\pi f_c t) ]其中( d_I(t) )和( d_Q(t) )都是±1的二进制数据序列且Q路数据在时间上比I路延迟一个符号周期( T_b )。这个式子说明MSK本质上和OQPSK是一家人只是基带脉冲从矩形变成了半正弦。这个等价关系对仿真特别重要因为后面自己搭MSK调制器时就是按这个结构来的。每个符号周期内MSK信号的相位会线性增加或减少( \frac{\pi}{2} )所以把相位轨迹画出来是一条斜率不变但方向可以切换的折线。折线在符号边界处没有跳变只是斜率改变这就叫相位连续。2.2 GMSK高斯滤波器、BT积与码间干扰的代价GMSK和MSK的区别只差一个前置的高斯低通滤波器。这个滤波器的频率响应是[ H(f) e^{-\alpha^2 f^2} ]其中( \alpha )与滤波器的3dB带宽( B )有关。工程上更喜欢用归一化的“BT积”来描述它也就是滤波器的3dB带宽( B )乘以符号周期( T_b )。BT积是一个无量纲的数它决定了滤波强度BT积很大比如接近无穷大时高斯滤波器几乎不工作GMSK退化成MSKBT减小滤波器带宽变窄对信号的高频成分压得更多功率谱更集中但滤波必然带来时域扩展一个符号的能量会拖到相邻符号产生码间干扰。GSM系统用的BT积是0.3蓝牙经典BR标准在教材里通常用0.5做演示。BT0.3时高斯滤波器的冲激响应在时间上会扩展到三到四个符号周期所以接收端必须用均衡或最大似然序列估计来弥补这个ISI否则误码性能会下降。这也是为什么GSM接收机里普遍存在均衡器不是发射机做得不好而是工程上选择了“用带宽换频谱效率再用接收机算法把性能捞回来”的策略。2.3 接收端为什么常见差分检测但仿真BER常用相干解调调制系统仿真不只是发射机的事接收机怎么选会直接影响结果的解读。实际硬件中GMSK/MSK接收机有两大类相干解调和非相干差分检测。相干解调性能好但需要恢复载波的频率和相位在快速移动和衰落环境下做起来很费劲差分检测对相位同步要求低实现简单但噪声会被相乘放大性能有一定损失。做Simulink仿真时我一般分成两个层次来对待如果目的是对比MSK和GMSK在AWGN信道下的理论误码性能就直接用Communications Toolbox里的CPM Demodulator Baseband做相干解调这样结果稳定也方便和理论曲线对照如果目的是验证某个低功耗蓝牙或卫星接收机的算法架构那就得自己搭差分检测链路。本文后面主要按前者来展开因为这是理解调制系统原理最高效的路径。3. 在Simulink里搭出完整的MSK/GMSK收发链路3.1 两条技术路线选择进入实操环节。在Simulink里做MSK/GMSK调制系统仿真有两条路线我建议都要了解路线实现方式优点缺点适用场景路线ACommunications Toolbox的CPM调制/解调模块搭建快、参数直观、支持MSK/GMSK一键切换、BER统计稳定模块封装较深底层波形细节看不全链路级性能验证、参数扫描、课程设计路线B自己搭正交调制器半正弦加权、高斯滤波、相位积分信号每一根线都可见能彻底搞懂原理模块多、数值问题多、调试费时间入门教学、算法原型验证、FPGA参考模型实际项目里我先用路线A把链路跑通建立整体概念再用路线B只搭MSK发射机来验证“半正弦加权Q路延迟”的等价关系。这样做既省时间又不会让细节成为黑盒。3.2 路线ACPM调制解调模块的搭建步骤新建Simulink模型后按下面步骤搭发射机到接收机的完整链路在模块库中搜索并添加Random Integer GeneratorSet size设为2Sample time设为1秒。这里把符号周期设成1秒是为了方便观察采样率后续再提高。添加CPM Modulator Baseband按表格设置参数。添加AWGN Channel噪声模式选择“Eb/N0 (dB)”Eb/N0先设15相当于一个较高信噪比。添加CPM Demodulator Baseband参数与调制器保持一致输出类型选Bit。添加Error Rate Calculation接收端输出接入Rx端口把发射端原始比特延迟后接Tx端口再用Display显示误码数。在信号线上接Time Scope和Spectrum Scope分别查看时域波形和功率谱密度。MSK和GMSK在CPM模块里的参数设置如下参数MSKGMSKM-ary number22Modulation index0.50.5Frequency pulse shapeRectangularGaussianPulse length (symbol intervals)13或4Bandwidth-time product不适用0.3Samples per symbol1616需要提醒一句不同MATLAB版本里模块名称和参数名可能略有差异。比如有的版本叫“Pulse length”有的叫“Pulse length in symbol intervals”BT积的选项有的叫“Bandwidth time product”。如果找不到直接在模块参数对话框下方看帮助说明即可。3.3 同步和BER统计的接法新手经常在这里卡住明明发射机和接收机参数都对BER却怎么都不归零。原因几乎都是解调延迟没对齐。CPM解调器内部有滤波器、Viterbi解码等处理输出数据会比发射端延迟若干个符号周期。最稳妥的做法是先把AWGN Channel的Eb/N0设得很大比如100相当于无噪声跑一次仿真通过调整Error Rate Calculation模块里的Delay参数让误码数变为0。记下这个Delay值再加噪声做正式统计。如果你懒得手调可以用一个Transport Delay模块在参考比特线上补延迟配合示波器观察两个序列对齐的情况。这里还有一个很实用的技巧把发射机的数据源和接收机的解调输出都接到同一个Data Type Conversion再进示波器用“Logic Analyzer”或“Time Scope”同时观察原始比特和解调比特。看到两者完全重合后再去设置延迟参数就非常直观。3.4 路线B自己搭MSK正交调制器来验证原理如果你像我一样总觉得CPM模块是个黑盒那建议抽半小时搭一个MSK正交调制器。结构就是前面提到的正交分解式。在Simulink里可以这样搭用Bernoulli Binary Generator生成随机比特Sample time设为1秒。用Unipolar to Bipolar Converter把0/1变成±1或者用Gain配合Bias来做。把数据分成I、Q两路。I路直接输出Q路用Delay模块延迟1个符号周期正常情况还可以做符号翻转具体取决于数据映射关系。生成半正弦加权因子用一个斜坡信号源计算时间( t )分别给I路乘( \cos(\pi t / (2T_b)) )给Q路乘( \sin(\pi t / (2T_b)) )。如果模型里已经是离散时间可以直接用Repeating Sequence产生半正弦包络。复基带信号的实部为I路加权结果虚部为Q路加权结果。把这个复信号直接给示波器看就能看到MSK特有的圆形相位轨迹。这里要特别提醒自己搭的时候所有模块的采样时间必须统一。比如符号周期设为1秒、每符号16个采样点那么总体采样时间应该是1/16秒。如果哪个模块Sample Time不一致Simulink会报错或者更隐蔽地出现信号错位导致波形看起来很怪。至于GMSK我不建议在Simulink里用“高斯滤波器调频器”的方式从零搭发射机因为模拟滤波器和连续时间调频的数值误差很容易让结果失真。直接用CPM模块设置成Gaussian脉冲是最省心也最可靠的做法。4. 仿真参数怎么定才不容易出“假结果”4.1 每符号采样点数和数据长度仿真结果可信不可信首先看两个数字每符号采样点数Samples per symbol和统计符号数。Samples per symbol决定波形离散化的精细程度。一般来说至少8个点推荐16个点严谨一点可以设32。如果这个值太小高斯滤波器的输出会被严重“锯齿化”频谱上会出现不该有的周期波纹。我在实际仿真里观察过从8提高到16功率谱旁瓣会有明显变化从16再提高到32变化就很小了。所以16是一个性价比很高的经验值。统计符号数取决于你要看什么。如果只是定性观察波形几百个符号就够如果要统计BER到( 10^{-4} )这种量级至少需要几万个符号否则误码数太少统计误差很大。仿真停止时间可以这样算符号周期×符号数。比如符号周期1秒、统计10000个符号仿真停止时间就是10000秒。这个时间看起来大但由于是离散系统Simulink跑起来并不慢。4.2 高斯滤波器BT积和脉冲长度L在CPM模块里GMSK除了BT积还有一个容易被忽略的参数Pulse length也就是高斯频率脉冲的持续时间用符号数表示。这个参数本质上是把理论上无限长的高斯脉冲截断到有限长度。BT0.3时高斯脉冲的主要能量集中在3到4个符号以内所以Pulse length设3或4就够。设太小等于把高斯脉冲硬切掉一部分发射频谱会出现鼓起的不正常旁瓣设太大计算量增加但改善有限。如果你在频谱上看到GMSK旁瓣有异常波纹可以先把Pulse length从3改到4或5看看是否改善。BT积这个参数我建议你在仿真中一定要做一次扫描。从0.5、0.4、0.3到0.2每改一次就看一下功率谱和BER的变化。这个过程做完你对“频带利用率和误码性能的权衡”会有非常具体的认识。4.3 Eb/N0和SNR为什么必须用统一参考对比MSK和GMSK的误码性能时一个最常见的错误是用SNR信号噪声功率比而不是Eb/N0每比特能量与噪声功率谱密度之比。原因是两者占用的带宽不同同样的SNR对应的比特能量可能完全不同比较起来没有意义。在AWGN Channel模块里我建议直接选“Eb/N0 (dB)”模式。对二进制MSK/GMSK每个符号正好一个比特所以“Bits per symbol”设为1。如果使用模块默认的SNR模式需要自己换算很容易出错。还有一种更隐蔽的情况你在Simulink里给数据加了增益或者改变了信号功率AWGN Channel内部会基于输入信号功率计算噪声如果功率归一化没做对仿真出来的BER曲线整体偏移但看起来又很“合理”这种假结果最坑人。我的做法是先把发射机输出信号功率归一到1然后在AWGN Channel里用Eb/N0模式并且每换一组参数比如从MSK切到GMSK、改变BT积都保持同样的参考设置。4.4 仿真步长和数值稳定性Simulink的CPM调制解调器内部是离散时间系统建议把求解器设成“Fixed-step”并手动指定步长。步长就是采样周期等于符号周期除以Samples per symbol。比如符号周期1秒、每符号16个点固定步长就是1/16秒。很多人在默认“Variable-step”下也能跑通但变步长会在信号斜率变化剧烈的地方自动加密计算点不仅慢而且可能在相位轨迹上制造难以解释的细节。设成定步长后结果可复现性大大提高也方便你把模型自动生成代码或者做批量参数扫描。5. 从波形、频谱到误码率看懂两种调制的本质差异5.1 看频率轨迹和相位轨迹最直观的区别把MSK和GMSK的调制输出接到Time Scope上重点看频率轨迹和相位轨迹。MSK的基带频率脉冲是矩形意味着在每个符号周期内频率值恒定符号切换时频率是瞬间跳变的。虽然相位没有跳变但看频率轨迹就会看到明显的方波形状。对应的相位轨迹是一条斜率恒定的折线方向可以在符号边界快速切换看起来有棱角。GMSK经过高斯滤波后频率脉冲变成了钟形曲线频率切换是平滑过渡的。相位轨迹在符号边界处会显现出圆弧过渡整个轨迹比MSK圆润很多。这个区别用一句话总结GMSK把MSK的所有“棱角”都磨平了代价是符号之间互相重叠。在Simulink里看相位时记得先把信号转成角度并做去卷绕unwrap否则你会发现相位一直在0到2π之间跳变没法看出连续的效果。可以用一个MATLAB Function模块写unwrap(angle(u))再送给Scope。5.2 看功率谱密度带外泄漏差距有多大调制方式的频谱特性最直观的就是看功率谱密度。用Spectrum Scope同时显示MSK和GMSK的PSD可以看到两个关键差异。第一主瓣宽度不同。MSK的主瓣相对宽GMSK在BT0.3时主瓣明显收窄。第二旁瓣衰落速度不同。MSK虽然相位连续旁瓣衰减比普通FSK好但和GMSK相比还是有明显差距。GMSK的旁瓣下降得非常快带外辐射被大幅压低。GSM之所以能用200kHz的信道间隔承载270.833kbps的码流靠的就是GMSK这种优异的频谱集中能力。看频谱时要用足够长的数据做平均否则谱线毛糙难以对比。建议在Spectrum Scope里开启“Averaging”至少累计几千个符号再做观察。这个数字不是随便拍的。5.3 看误码率曲线恒定包络的代价误码率是最终裁判。在AWGN信道下MSK采用相干解调时误码性能和BPSK完全一致理论曲线就是[ P_e Q\left(\sqrt{\frac{2E_b}{N_0}}\right) ]我经常用这条理论曲线来检验自己的Simulink链路有没有搭对如果MSK的仿真BER点和理论曲线叠不上那问题基本出在延迟对齐、噪声设置或参数一致性上。GMSK因为高斯滤波引入了可控的码间干扰误码性能会有损失。BT0.3时在BER等于( 10^{-4} )这个水平上GMSK比MSK要损失1到1.5dB左右。BT0.5时损失小一些大约零点几dB。如果接收端用Viterbi均衡或者最大似然序列估计这种损失可以被部分弥补这也是为什么GSM系统敢用GMSK。把这两张图放在脑子里你就能理解调制选型从来不是“哪个好”的问题而是“在什么约束下选哪个”的问题MSK解调简单、性能好但频谱不够紧凑GMSK把频谱压下来了但需要更复杂的接收机来换回性能。6. 仿真过程中我踩过的几个坑以及对应解法6.1 坑一CPM调制端和解调端参数对不上BER直接卡在0.5我最早做这个仿真时犯过一个低级错误调制端选了Gaussian脉冲、BT0.3解调端却忘了改默认还在Rectangular脉冲上。结果BER显示0.5看起来就像没信号一样。