简介GSM通信系统MATLAB仿真程序面向通信工程专业学生、研究人员与移动通信入门者提供了一套可视化复现GSM物理层关键技术的完整仿真链路涵盖卷积编码、交织、GMSK调制、频分多址与时分多址接入、衰落信道模拟、维特比译码与误码率评估等环节。资源压缩包共31个文件其中28个m脚本/函数文件对应发射、接收、信道模型、数据生成与性能分析等独立模块另含3个sim模型文件用于系统级仿真整体体积仅30KB结构清晰便于逐模块运行与拆解。目前已有470人浏览学习适合课程设计、毕业设计或自学移动通信原理时参考。通过运行该程序可直观观察调制星座图、频谱图及交织前后差异灵活修改参数计算误码率与符号错误率深入理解信道编码与交织如何对抗衰落和干扰既巩固GSM系统知识也提升MATLAB编程与通信仿真能力。1. GSM仿真程序是什么一张误码率曲线背后的完整物理层链路如果你正在做通信方向的课程设计或者刚接手物理层算法岗大概率会遇到这样一个场景网上散落着各种GSM仿真程序 matlab 工程打开却分不清哪个能一次跑通、哪个只是挂了个名字。我这次拆的这份工程属于少有的把GSM物理层端到端做全的——从卷积编码、交织、突发成型、GMSK调制到多径信道、定时同步、信道估计和MLSE均衡全链路可执行。它不是带界面的大屏Demo跑完能直接拿到误码率BER曲线、星座图和眼图能用于课程实验也能作为真实GSM物理层算法验证的骨架。适合通信工程学生、刚入门的算法工程师以及要在MATLAB里对比其它调制方式的人。2. 架构设计GSM协议栈中的物理层模块化拆解与目录规划2.1 物理层仿真覆盖哪些协议点GSM系统结构分MS、BTS、BSC、MSC等多个网元但这份仿真程序聚焦的是空中接口的物理层链路也就是从发射端信源比特到接收端判决比特这段。整个链路包含四个大环节信源编码之前的信息比特这里用随机数模拟、信道编码卷积码、交织与突发成型、GMSK调制经过信道后接收端做匹配滤波、定时同步、信道估计和Viterbi均衡译码。协议站点的选择很典型GSM全速率语音信道采用的是(2,1,4)卷积码码率1/2约束长度5生成多项式写作八进制[27 31]交织器是按块交织设计的发射端把卷积码输出按块重排接收端逆序复原调制部分采用GMSK带宽时间积BT固定为0.3符号速率约270.833kbit/s。这些特征决定了仿真程序的结构也决定了代码里很多常量的来源——比如训练序列长度是26比特、一个正常突发是148比特3尾58数据26训练58数据3尾。2.2 源码目录结构与职责边界工程代码按“发射机 / 信道 / 接收机 / 工具”四大块组织每一块都是独立的函数文件主脚本只负责把它们串起来。这种分层方式和实际GSM协议栈的物理层处理流程一一对应拿到源码后看目录就能猜到每个文件的作用二次开发时替换某一个模块也不会影响其它部分。gsm_sim/ ├─ main.m 主脚本参数配置与BER扫描 ├─ tx/ │ ├─ conv_encoder.m 卷积编码含生成多项式定义 │ ├─ interleaver.m 块交织与解交织 │ ├─ burst_formatter.m 突发组装加尾比特、插训练序列 │ └─ gmsk_modulator.m GMSK调制 ├─ channel/ │ ├─ multipath_channel.m 典型城市多径信道 │ └─ awgn_channel.m 加性白高斯噪声 ├─ rx/ │ ├─ freq_offset_comp.m 频偏补偿 │ ├─ symbol_sync.m 训练序列互相关定时同步 │ ├─ channel_estimator.m 最小二乘信道估计 │ └─ mlse_equalizer.m MLSE均衡 └─ utils/ ├─ ber_calc.m 误码率统计 └─ plot_routine.m 绘制星座图、眼图、误码率曲线比较关键的设计决策是接收机单独拆出symbol_sync.m和channel_estimator.m而不是把它们揉进均衡器里。原因在于GSM突发结构里训练序列正好位于突发中间定时同步、信道估计、均衡三个阶段共用一个训练序列拆开后每段逻辑都能单独调试——信道估计不对时可以单独打印est_h出来看不必等整条链路跑完。2.3 模块间的数据流与关键中间信号主脚本里传的是复基带波形不是比特只有收发两端的最外层才处理比特。发射端burst_formatter.m输出的是按符号速率排列的比特序列含训练序列和尾比特gmsk_modulator.m把它变成复基带采样点tx_signal信道模块对tx_signal做卷积加噪得到rx_signal接收端的symbol_sync.m找到突发起始位置channel_estimator.m用训练序列估计出信道冲激响应est_hmlse_equalizer.m在est_h上做网格搜索输出软/硬判决比特。% main.m 核心流程节选 conf.BT 0.3; % GSM标准带宽时间积 conf.fs 270833 * 4; % 采样率每符号4个采样点 conf.training_seq load_training_seq(1); % 载入26位训练序列 tx_bits randi([0 1], 10000, 1); encoded conv_encoder(tx_bits, conf); interleaved interleaver(encoded, conf, encode); burst burst_formatter(interleaved, conf); tx_signal gmsk_modulator(burst, conf);这段代码传达的信息很清楚每一层之间传递的是完整数据单元不是零散比特。conf结构体贯穿所有模块里面放的都是仿真级参数——BT、采样率、训练序列、突发长度把conf单独抽出来后续做参数扫描时就不用改函数签名。3. 核心实现GMSK调制、突发成型与Viterbi均衡的MATLAB代码3.1 GMSK调制高斯滤波与相位累积的落地写法GMSK的本体是CPFSK连续相位频移键控它先把0/1比特映射成±1通过高斯低通滤波器做脉冲成型再对成型脉冲做积分得到相位轨迹最后用复指数表达基带信号。这份工程里的gmsk_modulator.m跳过了方便但黑匣子化的工具箱函数直接手写高斯滤波器方便你在实验中改BT观察频谱变化。function ph gmsk_modulator(data, conf) % data: 差分编码后的 0/1 比特序列列向量 % conf: 包含 BT、Tb、fs、Nsyn 等仿真参数 Tb 1 / 270833; % GSM符号周期 Ls round(conf.fs * Tb); % 每个符号的采样点数 B conf.BT / Tb; % 高斯滤波器3dB带宽 K sqrt(log(2)) / (sqrt(2) * pi * B); data double(data(:)) * 2 - 1; % 0 - -11 - 1 span conf.Nsyn; % 滤波器截断符号数常见取4 n (-span*Ls : span*Ls) * Tb / Ls; hg exp(-(n / K).^2); % 高斯冲激响应 hg hg / sum(hg); % 归一化直流增益 up zeros(1, length(data) * Ls); up(1 : Ls : end) data; % 上采样其余点补零 y conv(up, hg); % 高斯成型 y y(1 : length(up)); % 相位增量 (pi/(2*Tb)) * 高斯脉冲积分 ph cumsum(y) * (pi / 2) * (Tb / Ls); end逻辑说明cumsum(y)是相位累积的核心每次累加相当于对高斯脉冲做时间积分得到的ph就是连续变化的调制相位最后把ph写成exp(1j*ph)就是在基带做正交调制。Ls决定一个符号的采样点数采样率越高波形重建越精确仿真时间成比例增长工程里通常取4到8。参数说明conf.BT 0.3是GSM标准参数BT越大频谱滚降越快、邻道泄漏越小但符号间干扰越大conf.Nsyn控制滤波器截断长度太短会看到明显的频谱拖尾太长相位轨迹变缓。常见做法是固定BT用这组代码画相位轨迹图验证——正常GSM相位轨迹在最差路径上的最大相位偏差应该在0.6左右。3.2 突发成型与训练序列帧结构里的确定性成分GSM正常突发不是裸数据它把148比特排成“尾比特3 数据57 训练序列26 数据57 尾比特3”的结构接收端的定时同步和信道估计全靠中间这26位训练序列。工程里的实现直接用掩码拼接把卷积编码、交织后的数据切两段中间插入训练序列。function burst burst_formatter(interleaved_bits, conf) tail zeros(1, 3); % 尾比特固定为0 n length(interleaved_bits); half floor(n / 2); burst [tail, ... interleaved_bits(1 : half - 3)., ... conf.training_seq, ... interleaved_bits(half - 2 : end)., ... tail]; end逻辑说明tail放在突发两端作为均衡器训练和路径收尾的已知状态接收端Viterbi译码时知道起点是0能有效避免起始状态不确定训练序列放在正中间这样离两端数据距离最近信道的时变影响最小。参数说明conf.training_seq是从GSM标准的8组训练序列中任取一组长度26但实际起作用的是中间16位两端各5位作为保护带。这段话值得记住做信道估计时卷积矩阵的行数用的是训练序列的总长度但有效估计长度是26 - LchLch为信道记忆深度所以不要把训练序列全部当成有效数据。3.3 接收机侧信道估计与MLSE均衡接收端最关键的链路是“粗同步→信道估计→均衡”。粗同步用训练序列与接收符号做互相关峰值位置就是突发头部信道估计用最小二乘法解卷积方程均衡则用Viterbi网格搜索。工程里这三步对应三个函数核心代码集中在均衡器里。function bits mlse_equalizer(rx_sym, est_h, Lch) % rx_sym: 定时同步后的符号序列按符号速率采样 % est_h : 信道冲激响应长度 Lch1 N_APP 2^Lch; % 网格状态数 n_sym length(rx_sym); metric zeros(N_APP, 1); survivor zeros(N_APP, n_sym); for n 1 : n_sym new_metric inf(N_APP, 2); new_surv zeros(N_APP, 2); for state 1 : N_APP % 状态 最近Lch个发送比特的编码 for bit [0 1] prev_state floor((state-1)/2) 1; prev_bits de2bi(prev_state - 1, Lch, left-msb); tx_pattern [prev_bits, bit] * 2 - 1; expected sum(est_h .* tx_pattern); cand metric(prev_state) abs(rx_sym(n) - expected)^2; if cand new_metric(state, bit1) new_metric(state, bit1) cand; new_surv(state, bit1) prev_state; end end end metric min(new_metric, [], 2); % survivor 路径回溯略 end bits traceback(survivor); end逻辑说明Viterbi均衡把多径信道看成有限状态机state表示最近的Lch个发送比特expected是假设信道冲激响应为est_h时的理想接收值cand是累计路径度量取最小度量对应的路径最后回溯得到发送比特序列。这段代码对应的是BPSK调制的MLSE简化模型GMSK实际应用时会把调制符号先化成正交基展开但状态搜索框架完全一致。参数说明Lch是信道记忆深度典型城市信道6条多径时取5到6Lch越大均衡性能越好但状态数按2^Lch指数增长。工程实践里Lch6时状态数64MATLAB单Symbol循环还能接受不要盲目把Lch调到8以上仿真时间会指数增长收益却非常有限。4. 运行指南从零跑通仿真脚本并复现典型城市信道曲线4.1 环境准备与文件清单运行这份工程只需要MATLAB基础环境最好装上Communications Toolbox因为polytrellis、convenc、vitdec这些函数在工具箱里是现成的。如果没装工具箱代码里的手写均衡器仍然能跑但卷积码和Viterbi译码要换成自实现版本——这也是很多同学翻车的地方装了低版本MATLAB工具箱函数名冲突导致报错。文件清单如下文件名职责改动频率main.m参数配置、循环扫描、画图调用高tx/gmsk_modulator.mGMSK调制低rx/mlse_equalizer.m均衡核心低channel/multipath_channel.m多径信道抽头中utils/plot_routine.m绘图高4.2 运行主脚本与观察指标在MATLAB里切换到工程根目录运行main。脚本会先打印当前配置然后逐个Eb/N0点进行Monte Carlo仿真最后弹出三张图误码率曲线、接收星座图、眼图。 cd /path/to/gsm_sim main跑通一次需要等多久取决于nbits和 Eb/N0点数。拿默认的nbits10000来说Eb/N0从0到10dB共6个点全链路跑完大约两三分钟。如果直接把nbits改成100000单点误码率统计会更平滑但每点耗时翻十倍做好跑一晚上的准备。误码率曲线怎么看理想GSM在AWGN信道下Eb/N0大于6dB时BER应该低于1e-3多径信道下因为引入了符号间干扰同Eb/N0下BER会抬高半个到一个量级这是正常现象。如果你的曲线在多径信道下仍然和AWGN几乎重合反而说明信道模块根本没生效——常见原因就是信道抽头幅度没归一化。4.3 参数调整实验Eb/N0扫描与多径时延改变跑通一次只是开始这份工程最有意思的地方在于改参数看变化。第一个建议实验是改main.m里的EbN0dB数组范围比如EbN0dB 0:1:12把间隔从2dB改成1dB曲线更密能更精确看到香农极限附近的行为。第二个建议实验是修改multipath_channel.m的时延与功率参数% 典型城市TU信道6条多径最大时延约1.2us tap_power_dB [0, -3, -6, -8, -10, -12]; % 抽头平均功率 tap_delay_us [0, 0.1, 0.3, 0.5, 0.8, 1.2]; % 时延单位us把时延全部缩小到原来的三分之一相当于信道从“频率选择性”变成“平衰落”BER会显著下降反过来把时延拉到2us以上均衡器跟不上BER曲线会出现错误底。通过这种改动你能直观体会为什么GSM实际部署中要限制小区半径——保护周期和均衡能力都有限时延扩展太大会直接把链路打崩。5. 常见问题与排查相位模糊、定时偏移与信噪比换算5.1 误码率曲线出现“肩膀”不是继续下降而是停在1e-3现象BER曲线在Eb/N0达到8dB后不再下降形成一条平线看起来像有系统误差。原因这是GMSK相位模糊的典型症状。GMSK接收端解调得到的符号会有一个不确定的90度或180度相位偏转如果没有做差分编码判决器会直接判错一部分比特。解决在发射端gmsk_modulator之前加一步差分编码接收端在mlse_equalizer输出后做差分判决。差分编码的公式很直白d(n) b(n) xor b(n-1)判决时反过来b_est(n) d_est(n) xor d_est(n-1)。5.2 星座图整体旋转固定角度BER却不低现象接收星座图的四个散点不是理想位置而是绕原点旋转了比如45度但误码率统计出来居然还能接受。原因帧边界没有对齐。训练序列位置偏移会导致信道估计矩阵错位估计出的est_h代入了错误的相位旋转。解决先用xcorr做一次粗同步把接收符号序列和本地训练序列调制符号做互相关峰值位置就是训练序列起点。注意这里一定要用调制后的复数符号做相关不能用0/1比特——用比特相关只会得到0查不出任何东西。5.3 多径信道下BER居高不下均衡器似乎没起作用现象AWGN信道下仿真结果正常换成multipath_channel后BER从1e-3跳到1e-1加了均衡器也没改善。原因均衡器是接了但信道估计没做或者est_h传了一个全1向量进mlse_equalizer。这是代码结构问题——est_h是channel_estimator的输出必须显式传给均衡器。解决在main.m里打上断点逐个检查est_h的复数幅度。典型城市信道下est_h应该有6个非零抽头幅度依次递减。如果看到est_h只有一个非零值说明信道估计矩阵构造有误检查卷积矩阵的行列数是否对应length(train_seq) - Lch。5.4 曲线整体向右移了约3dB现象仿真结果和理论曲线形状完全一致但整条BER曲线比理论值右移了3dB。原因Eb/N0与SNR换算错误。很多人直接把awgn函数的SNR参数填成Eb/N0忘了码率已经不完全是1以及采样率和符号率的差别。GSM的卷积码码率是1/2也就是说每1个信息比特对应2个信道比特换算时必须计入码率因子。解决统一用SNR EbN0dB 10*log10(Rb/Bn)换算其中Rb是信息比特率Bn是噪声等效带宽。工程里更省事的做法是先固定发送功率用awgn(rx_signal, snr, measured)再在注释里写明换算关系避免每次扫参时来回改。5.5 Viterbi均衡器回溯时比特顺序颠倒前段对后段错现象均衡后前58个数据比特译码正确后58个全错而且错得规律。原因网格状态映射用到的de2bi默认是右最高位right-msb而发射端突发成型时数据是按左最高位排列的。两者不一致导致状态编号和实际发送比特串的映射错位。解决所有de2bi和bi2de调用统一加选项left-msb保证发射端、均衡器、译码器三处比特顺序基准一致。这个坑在MATLAB老版本里尤其致命——有些版本默认是left-msb有些是right-msb同一段代码换版本就翻车。6. 二次开发方向从GSM向5G NR物理层演进的基本动作这份工程最大的价值在于它把物理层收发链路的接口定义清楚了发射端输入比特输出基带波形接收端输入基带波形输出判决比特。基于这个接口改造方向非常明确。第一个方向是替换调制器。把gmsk_modulator换成QPSK或OFDM映射接收端同步均衡对应调整。常见做法是保留突发结构、训练序列位置和均衡器框架只把相位累积替换成星座点映射这样做的工程量大约只有两个函数但能快速对比不同调制方式在同一信道下的性能差距。验证方法是画BER曲线——QPSK的理论曲线在AWGN下大约比GMSK好2dB左右如果你的结果偏离太大说明替换时忘了做幅度归一化。第二个方向是替换信道编码。把 (2,1,4) 卷积码换成Turbo码或LDPC码需要在conv_encoder.m和译码链路里同时改。经验是先用MATLAB自带的comm.LDPCEncoder跑通再用手写实现替换避免一上来就陷入编码矩阵构造的泥潭。改动后一定要验证编码增益——好一点的编码在BER1e-4处能比卷积码多带来3dB以上的增益。第三个方向是把仿真波形送进SDR硬件。用comm.SDRTransmitter把tx_signal重采样到设备采样率再用接收端程序对接真实射频信号。这里最大的坑是频偏——仿真里可以不考虑真实传播必然有几十到几百赫兹的频偏需要把freq_offset_comp.m的补偿循环次数调大否则星座图会持续旋转。还有一个我自己踩过的教训每改动一个模块先跑一个全零比特序列的回环测试——发射机输入全0接收机输出必须全0。任何一处相位翻转或比特方向错误都会让这个测试立刻暴露。从那以后我每次拿到新的通信仿真项目都强制走一遍全零回环再开始改参数。希望帮到你。本文还有配套的精品资源点击获取
简介GSM通信系统MATLAB仿真程序面向通信工程专业学生、研究人员与移动通信入门者提供了一套可视化复现GSM物理层关键技术的完整仿真链路涵盖卷积编码、交织、GMSK调制、频分多址与时分多址接入、衰落信道模拟、维特比译码与误码率评估等环节。资源压缩包共31个文件其中28个m脚本/函数文件对应发射、接收、信道模型、数据生成与性能分析等独立模块另含3个sim模型文件用于系统级仿真整体体积仅30KB结构清晰便于逐模块运行与拆解。目前已有470人浏览学习适合课程设计、毕业设计或自学移动通信原理时参考。通过运行该程序可直观观察调制星座图、频谱图及交织前后差异灵活修改参数计算误码率与符号错误率深入理解信道编码与交织如何对抗衰落和干扰既巩固GSM系统知识也提升MATLAB编程与通信仿真能力。1. GSM仿真程序是什么一张误码率曲线背后的完整物理层链路如果你正在做通信方向的课程设计或者刚接手物理层算法岗大概率会遇到这样一个场景网上散落着各种GSM仿真程序 matlab 工程打开却分不清哪个能一次跑通、哪个只是挂了个名字。我这次拆的这份工程属于少有的把GSM物理层端到端做全的——从卷积编码、交织、突发成型、GMSK调制到多径信道、定时同步、信道估计和MLSE均衡全链路可执行。它不是带界面的大屏Demo跑完能直接拿到误码率BER曲线、星座图和眼图能用于课程实验也能作为真实GSM物理层算法验证的骨架。适合通信工程学生、刚入门的算法工程师以及要在MATLAB里对比其它调制方式的人。2. 架构设计GSM协议栈中的物理层模块化拆解与目录规划2.1 物理层仿真覆盖哪些协议点GSM系统结构分MS、BTS、BSC、MSC等多个网元但这份仿真程序聚焦的是空中接口的物理层链路也就是从发射端信源比特到接收端判决比特这段。整个链路包含四个大环节信源编码之前的信息比特这里用随机数模拟、信道编码卷积码、交织与突发成型、GMSK调制经过信道后接收端做匹配滤波、定时同步、信道估计和Viterbi均衡译码。协议站点的选择很典型GSM全速率语音信道采用的是(2,1,4)卷积码码率1/2约束长度5生成多项式写作八进制[27 31]交织器是按块交织设计的发射端把卷积码输出按块重排接收端逆序复原调制部分采用GMSK带宽时间积BT固定为0.3符号速率约270.833kbit/s。这些特征决定了仿真程序的结构也决定了代码里很多常量的来源——比如训练序列长度是26比特、一个正常突发是148比特3尾58数据26训练58数据3尾。2.2 源码目录结构与职责边界工程代码按“发射机 / 信道 / 接收机 / 工具”四大块组织每一块都是独立的函数文件主脚本只负责把它们串起来。这种分层方式和实际GSM协议栈的物理层处理流程一一对应拿到源码后看目录就能猜到每个文件的作用二次开发时替换某一个模块也不会影响其它部分。gsm_sim/ ├─ main.m 主脚本参数配置与BER扫描 ├─ tx/ │ ├─ conv_encoder.m 卷积编码含生成多项式定义 │ ├─ interleaver.m 块交织与解交织 │ ├─ burst_formatter.m 突发组装加尾比特、插训练序列 │ └─ gmsk_modulator.m GMSK调制 ├─ channel/ │ ├─ multipath_channel.m 典型城市多径信道 │ └─ awgn_channel.m 加性白高斯噪声 ├─ rx/ │ ├─ freq_offset_comp.m 频偏补偿 │ ├─ symbol_sync.m 训练序列互相关定时同步 │ ├─ channel_estimator.m 最小二乘信道估计 │ └─ mlse_equalizer.m MLSE均衡 └─ utils/ ├─ ber_calc.m 误码率统计 └─ plot_routine.m 绘制星座图、眼图、误码率曲线比较关键的设计决策是接收机单独拆出symbol_sync.m和channel_estimator.m而不是把它们揉进均衡器里。原因在于GSM突发结构里训练序列正好位于突发中间定时同步、信道估计、均衡三个阶段共用一个训练序列拆开后每段逻辑都能单独调试——信道估计不对时可以单独打印est_h出来看不必等整条链路跑完。2.3 模块间的数据流与关键中间信号主脚本里传的是复基带波形不是比特只有收发两端的最外层才处理比特。发射端burst_formatter.m输出的是按符号速率排列的比特序列含训练序列和尾比特gmsk_modulator.m把它变成复基带采样点tx_signal信道模块对tx_signal做卷积加噪得到rx_signal接收端的symbol_sync.m找到突发起始位置channel_estimator.m用训练序列估计出信道冲激响应est_hmlse_equalizer.m在est_h上做网格搜索输出软/硬判决比特。% main.m 核心流程节选 conf.BT 0.3; % GSM标准带宽时间积 conf.fs 270833 * 4; % 采样率每符号4个采样点 conf.training_seq load_training_seq(1); % 载入26位训练序列 tx_bits randi([0 1], 10000, 1); encoded conv_encoder(tx_bits, conf); interleaved interleaver(encoded, conf, encode); burst burst_formatter(interleaved, conf); tx_signal gmsk_modulator(burst, conf);这段代码传达的信息很清楚每一层之间传递的是完整数据单元不是零散比特。conf结构体贯穿所有模块里面放的都是仿真级参数——BT、采样率、训练序列、突发长度把conf单独抽出来后续做参数扫描时就不用改函数签名。3. 核心实现GMSK调制、突发成型与Viterbi均衡的MATLAB代码3.1 GMSK调制高斯滤波与相位累积的落地写法GMSK的本体是CPFSK连续相位频移键控它先把0/1比特映射成±1通过高斯低通滤波器做脉冲成型再对成型脉冲做积分得到相位轨迹最后用复指数表达基带信号。这份工程里的gmsk_modulator.m跳过了方便但黑匣子化的工具箱函数直接手写高斯滤波器方便你在实验中改BT观察频谱变化。function ph gmsk_modulator(data, conf) % data: 差分编码后的 0/1 比特序列列向量 % conf: 包含 BT、Tb、fs、Nsyn 等仿真参数 Tb 1 / 270833; % GSM符号周期 Ls round(conf.fs * Tb); % 每个符号的采样点数 B conf.BT / Tb; % 高斯滤波器3dB带宽 K sqrt(log(2)) / (sqrt(2) * pi * B); data double(data(:)) * 2 - 1; % 0 - -11 - 1 span conf.Nsyn; % 滤波器截断符号数常见取4 n (-span*Ls : span*Ls) * Tb / Ls; hg exp(-(n / K).^2); % 高斯冲激响应 hg hg / sum(hg); % 归一化直流增益 up zeros(1, length(data) * Ls); up(1 : Ls : end) data; % 上采样其余点补零 y conv(up, hg); % 高斯成型 y y(1 : length(up)); % 相位增量 (pi/(2*Tb)) * 高斯脉冲积分 ph cumsum(y) * (pi / 2) * (Tb / Ls); end逻辑说明cumsum(y)是相位累积的核心每次累加相当于对高斯脉冲做时间积分得到的ph就是连续变化的调制相位最后把ph写成exp(1j*ph)就是在基带做正交调制。Ls决定一个符号的采样点数采样率越高波形重建越精确仿真时间成比例增长工程里通常取4到8。参数说明conf.BT 0.3是GSM标准参数BT越大频谱滚降越快、邻道泄漏越小但符号间干扰越大conf.Nsyn控制滤波器截断长度太短会看到明显的频谱拖尾太长相位轨迹变缓。常见做法是固定BT用这组代码画相位轨迹图验证——正常GSM相位轨迹在最差路径上的最大相位偏差应该在0.6左右。3.2 突发成型与训练序列帧结构里的确定性成分GSM正常突发不是裸数据它把148比特排成“尾比特3 数据57 训练序列26 数据57 尾比特3”的结构接收端的定时同步和信道估计全靠中间这26位训练序列。工程里的实现直接用掩码拼接把卷积编码、交织后的数据切两段中间插入训练序列。function burst burst_formatter(interleaved_bits, conf) tail zeros(1, 3); % 尾比特固定为0 n length(interleaved_bits); half floor(n / 2); burst [tail, ... interleaved_bits(1 : half - 3)., ... conf.training_seq, ... interleaved_bits(half - 2 : end)., ... tail]; end逻辑说明tail放在突发两端作为均衡器训练和路径收尾的已知状态接收端Viterbi译码时知道起点是0能有效避免起始状态不确定训练序列放在正中间这样离两端数据距离最近信道的时变影响最小。参数说明conf.training_seq是从GSM标准的8组训练序列中任取一组长度26但实际起作用的是中间16位两端各5位作为保护带。这段话值得记住做信道估计时卷积矩阵的行数用的是训练序列的总长度但有效估计长度是26 - LchLch为信道记忆深度所以不要把训练序列全部当成有效数据。3.3 接收机侧信道估计与MLSE均衡接收端最关键的链路是“粗同步→信道估计→均衡”。粗同步用训练序列与接收符号做互相关峰值位置就是突发头部信道估计用最小二乘法解卷积方程均衡则用Viterbi网格搜索。工程里这三步对应三个函数核心代码集中在均衡器里。function bits mlse_equalizer(rx_sym, est_h, Lch) % rx_sym: 定时同步后的符号序列按符号速率采样 % est_h : 信道冲激响应长度 Lch1 N_APP 2^Lch; % 网格状态数 n_sym length(rx_sym); metric zeros(N_APP, 1); survivor zeros(N_APP, n_sym); for n 1 : n_sym new_metric inf(N_APP, 2); new_surv zeros(N_APP, 2); for state 1 : N_APP % 状态 最近Lch个发送比特的编码 for bit [0 1] prev_state floor((state-1)/2) 1; prev_bits de2bi(prev_state - 1, Lch, left-msb); tx_pattern [prev_bits, bit] * 2 - 1; expected sum(est_h .* tx_pattern); cand metric(prev_state) abs(rx_sym(n) - expected)^2; if cand new_metric(state, bit1) new_metric(state, bit1) cand; new_surv(state, bit1) prev_state; end end end metric min(new_metric, [], 2); % survivor 路径回溯略 end bits traceback(survivor); end逻辑说明Viterbi均衡把多径信道看成有限状态机state表示最近的Lch个发送比特expected是假设信道冲激响应为est_h时的理想接收值cand是累计路径度量取最小度量对应的路径最后回溯得到发送比特序列。这段代码对应的是BPSK调制的MLSE简化模型GMSK实际应用时会把调制符号先化成正交基展开但状态搜索框架完全一致。参数说明Lch是信道记忆深度典型城市信道6条多径时取5到6Lch越大均衡性能越好但状态数按2^Lch指数增长。工程实践里Lch6时状态数64MATLAB单Symbol循环还能接受不要盲目把Lch调到8以上仿真时间会指数增长收益却非常有限。4. 运行指南从零跑通仿真脚本并复现典型城市信道曲线4.1 环境准备与文件清单运行这份工程只需要MATLAB基础环境最好装上Communications Toolbox因为polytrellis、convenc、vitdec这些函数在工具箱里是现成的。如果没装工具箱代码里的手写均衡器仍然能跑但卷积码和Viterbi译码要换成自实现版本——这也是很多同学翻车的地方装了低版本MATLAB工具箱函数名冲突导致报错。文件清单如下文件名职责改动频率main.m参数配置、循环扫描、画图调用高tx/gmsk_modulator.mGMSK调制低rx/mlse_equalizer.m均衡核心低channel/multipath_channel.m多径信道抽头中utils/plot_routine.m绘图高4.2 运行主脚本与观察指标在MATLAB里切换到工程根目录运行main。脚本会先打印当前配置然后逐个Eb/N0点进行Monte Carlo仿真最后弹出三张图误码率曲线、接收星座图、眼图。 cd /path/to/gsm_sim main跑通一次需要等多久取决于nbits和 Eb/N0点数。拿默认的nbits10000来说Eb/N0从0到10dB共6个点全链路跑完大约两三分钟。如果直接把nbits改成100000单点误码率统计会更平滑但每点耗时翻十倍做好跑一晚上的准备。误码率曲线怎么看理想GSM在AWGN信道下Eb/N0大于6dB时BER应该低于1e-3多径信道下因为引入了符号间干扰同Eb/N0下BER会抬高半个到一个量级这是正常现象。如果你的曲线在多径信道下仍然和AWGN几乎重合反而说明信道模块根本没生效——常见原因就是信道抽头幅度没归一化。4.3 参数调整实验Eb/N0扫描与多径时延改变跑通一次只是开始这份工程最有意思的地方在于改参数看变化。第一个建议实验是改main.m里的EbN0dB数组范围比如EbN0dB 0:1:12把间隔从2dB改成1dB曲线更密能更精确看到香农极限附近的行为。第二个建议实验是修改multipath_channel.m的时延与功率参数% 典型城市TU信道6条多径最大时延约1.2us tap_power_dB [0, -3, -6, -8, -10, -12]; % 抽头平均功率 tap_delay_us [0, 0.1, 0.3, 0.5, 0.8, 1.2]; % 时延单位us把时延全部缩小到原来的三分之一相当于信道从“频率选择性”变成“平衰落”BER会显著下降反过来把时延拉到2us以上均衡器跟不上BER曲线会出现错误底。通过这种改动你能直观体会为什么GSM实际部署中要限制小区半径——保护周期和均衡能力都有限时延扩展太大会直接把链路打崩。5. 常见问题与排查相位模糊、定时偏移与信噪比换算5.1 误码率曲线出现“肩膀”不是继续下降而是停在1e-3现象BER曲线在Eb/N0达到8dB后不再下降形成一条平线看起来像有系统误差。原因这是GMSK相位模糊的典型症状。GMSK接收端解调得到的符号会有一个不确定的90度或180度相位偏转如果没有做差分编码判决器会直接判错一部分比特。解决在发射端gmsk_modulator之前加一步差分编码接收端在mlse_equalizer输出后做差分判决。差分编码的公式很直白d(n) b(n) xor b(n-1)判决时反过来b_est(n) d_est(n) xor d_est(n-1)。5.2 星座图整体旋转固定角度BER却不低现象接收星座图的四个散点不是理想位置而是绕原点旋转了比如45度但误码率统计出来居然还能接受。原因帧边界没有对齐。训练序列位置偏移会导致信道估计矩阵错位估计出的est_h代入了错误的相位旋转。解决先用xcorr做一次粗同步把接收符号序列和本地训练序列调制符号做互相关峰值位置就是训练序列起点。注意这里一定要用调制后的复数符号做相关不能用0/1比特——用比特相关只会得到0查不出任何东西。5.3 多径信道下BER居高不下均衡器似乎没起作用现象AWGN信道下仿真结果正常换成multipath_channel后BER从1e-3跳到1e-1加了均衡器也没改善。原因均衡器是接了但信道估计没做或者est_h传了一个全1向量进mlse_equalizer。这是代码结构问题——est_h是channel_estimator的输出必须显式传给均衡器。解决在main.m里打上断点逐个检查est_h的复数幅度。典型城市信道下est_h应该有6个非零抽头幅度依次递减。如果看到est_h只有一个非零值说明信道估计矩阵构造有误检查卷积矩阵的行列数是否对应length(train_seq) - Lch。5.4 曲线整体向右移了约3dB现象仿真结果和理论曲线形状完全一致但整条BER曲线比理论值右移了3dB。原因Eb/N0与SNR换算错误。很多人直接把awgn函数的SNR参数填成Eb/N0忘了码率已经不完全是1以及采样率和符号率的差别。GSM的卷积码码率是1/2也就是说每1个信息比特对应2个信道比特换算时必须计入码率因子。解决统一用SNR EbN0dB 10*log10(Rb/Bn)换算其中Rb是信息比特率Bn是噪声等效带宽。工程里更省事的做法是先固定发送功率用awgn(rx_signal, snr, measured)再在注释里写明换算关系避免每次扫参时来回改。5.5 Viterbi均衡器回溯时比特顺序颠倒前段对后段错现象均衡后前58个数据比特译码正确后58个全错而且错得规律。原因网格状态映射用到的de2bi默认是右最高位right-msb而发射端突发成型时数据是按左最高位排列的。两者不一致导致状态编号和实际发送比特串的映射错位。解决所有de2bi和bi2de调用统一加选项left-msb保证发射端、均衡器、译码器三处比特顺序基准一致。这个坑在MATLAB老版本里尤其致命——有些版本默认是left-msb有些是right-msb同一段代码换版本就翻车。6. 二次开发方向从GSM向5G NR物理层演进的基本动作这份工程最大的价值在于它把物理层收发链路的接口定义清楚了发射端输入比特输出基带波形接收端输入基带波形输出判决比特。基于这个接口改造方向非常明确。第一个方向是替换调制器。把gmsk_modulator换成QPSK或OFDM映射接收端同步均衡对应调整。常见做法是保留突发结构、训练序列位置和均衡器框架只把相位累积替换成星座点映射这样做的工程量大约只有两个函数但能快速对比不同调制方式在同一信道下的性能差距。验证方法是画BER曲线——QPSK的理论曲线在AWGN下大约比GMSK好2dB左右如果你的结果偏离太大说明替换时忘了做幅度归一化。第二个方向是替换信道编码。把 (2,1,4) 卷积码换成Turbo码或LDPC码需要在conv_encoder.m和译码链路里同时改。经验是先用MATLAB自带的comm.LDPCEncoder跑通再用手写实现替换避免一上来就陷入编码矩阵构造的泥潭。改动后一定要验证编码增益——好一点的编码在BER1e-4处能比卷积码多带来3dB以上的增益。第三个方向是把仿真波形送进SDR硬件。用comm.SDRTransmitter把tx_signal重采样到设备采样率再用接收端程序对接真实射频信号。这里最大的坑是频偏——仿真里可以不考虑真实传播必然有几十到几百赫兹的频偏需要把freq_offset_comp.m的补偿循环次数调大否则星座图会持续旋转。还有一个我自己踩过的教训每改动一个模块先跑一个全零比特序列的回环测试——发射机输入全0接收机输出必须全0。任何一处相位翻转或比特方向错误都会让这个测试立刻暴露。从那以后我每次拿到新的通信仿真项目都强制走一遍全零回环再开始改参数。希望帮到你。本文还有配套的精品资源点击获取
平方复杂度无法被常规优化消除｜TaoToken 统一 Key 通道实测)
