简介这是一份基于MATLAB的GMSK调制解调实验课程设计文档面向通信工程专业本科生及需要完成同类课题的开发者系统讲解了高斯滤波最小频移键控的调制解调原理、Simulink模型搭建与瑞利信道下的仿真验证。文档共1个docx文件压缩包大小约1.81MB篇幅完整包含摘要、设计原理、功能模块设计、波形/眼图/误码率结果分析和结论等结构可直接作为课程设计说明书模板或仿真实战参考。内容从GMSK基本概念出发覆盖信号发生、调制解调、误码率统计和波形观察四个核心模块并重点分析BT参数对频谱效率和抗干扰能力的影响帮助读者在Simulink中复现完整链路并理解关键性能指标。已有464人学习下载适合希望系统掌握GMSK仿真流程、撰写实验报告或搭建通信系统模型的读者。1. GMSK调制解调实验到底要你交出什么GMSK是GSM时代蜂窝物理层的事实标准蓝牙里的GFSK也是它的近亲理解GMSK等于同时拿下连续相位调制 频谱成形两个通信基本功。用MATLAB做这个实验目的不是造一个能发射的信号而是看清三件事高斯滤波器如何用带宽换频谱纯度信息如何藏在每比特±π/2的相位增量里以及差分检测为什么在无线调制和解调电路原理课程里被反复强调。适合通信原理课设、准备面试的本科生或想把手写调制器方法捡起来的工程师。下面给的代码不挑MATLAB版本装好即可运行只有一处用到通信工具箱会给出无工具箱替代写法。2. GMSK调制把频率轨迹从折线磨成平滑曲线的建模与代码2.1 MSK到GMSK高斯滤波器到底改在哪一处MSK是连续相位FSK里最容易被接受的基准调制指数 h0.5瞬时频率在每个比特周期内恒定比特边界只会发生相位连续但频率跳变。相位轨迹是折线频率轨迹是方波。方波意味着高频分量丰富频谱旁瓣按 f⁻⁴ 衰减对相邻信道干扰大。GMSK在频率成形之前插入一个高斯型低通滤波器把方波形的瞬时频率磨成高斯脉冲形状的过渡曲线。更准确地说GMSK的基带等效模型是NRZ数据先经过预调制高斯滤波器输出作为瞬时频率再积分得到相位。高斯滤波器传递函数是H(f) exp( -ln2 · f² / (2B²) )其中 B 是该滤波器的 3dB 带宽BT乘积就是 B 和比特周期 Tb 的乘积。BT 是GMSK唯一真正需要手工调的参数它同时决定频谱和误码率。BT 越大滤波越弱频谱越接近MSKBT 越小频谱越紧但单个比特的能量被拉到相邻比特上眼图闭合得越厉害。这张取舍表通常直接写进实验报告BT 取值典型场景频谱特征代价∞即MSK理论基准主瓣最宽旁瓣按 f⁻⁴ 衰减无ISI眼图全开0.5蓝牙GFSK同族主瓣收窄旁瓣较低ISI轻微一比特差分可用0.3GSM主瓣窄旁瓣极低ISI明显解调需差分或均衡实验里最常见的做法是把同一套代码跑三遍BT 取 0.3、0.5 和一个很大的值近似MSK用 pwelch 量频谱、用误码率曲线看代价这组对照本身就有讨论价值。2.2 两种生成高斯成形脉冲的MATLAB写法调制器里唯一需要制造的元件就是高斯预调制滤波器。手头有通信工具箱时直接用gaussdesign它返回一个单位DC增益的对称FIRsum(g) 1这个性质后面归一化相位要用。bwT 0.3; % BT乘积 L 4; % 截断长度单位比特周期 ovs 16; % 每比特采样点数 g gaussdesign(bwT, L, ovs); % 单位DC增益sum(g)1没有通信工具箱时可以按高斯滤波器的解析式自己采样生成。3dB带宽B对应的高斯标准差是σ sqrt(ln2) / (2·π·B)把它作为滤波器时域宽度直接写函数function h gmsk_gauss_lpf(bwT, L, ovs) % 手写高斯预调制滤波器与gaussdesign等效 Tb 1; % 归一化比特周期 sigma sqrt(log(2)) / (2*pi*bwT); % 3dB带宽→标准差 t ((0 : L*ovs) - L*ovs/2) / ovs; % 时间轴单位比特周期 h exp(-t.^2 / (2*sigma^2)); % 高斯冲激响应 h h / sum(h); % 强制单位DC增益 end两个版本的共同点是sum(h)1。这个性质不是形式主义它保证一个持续整个比特周期的1输入滤波后在该比特周期内的积分恰好是1于是相位累加器里每个比特贡献的相位增量恰为 ±π/2。节外生枝说一句gaussdesign的BT是3dB带宽不是噪声等效带宽别拿0.5/Tb这种经验值去替代。2.3 NRZ保持 相位累加最小可运行调制器完整的复基带调制器分四步生成 ±1 数据 → 把每个比特保持 ovs 个采样得到NRZ波形 → 过高斯滤波器得到频率轨迹 → 累加积分并映射到复指数。注意一个新手高发错误不要用upsample造冲激串再接高斯滤波器那等于拿滤波器的冲激响应当频率脉冲在BT较大时相位轨迹和频谱都会失真必须先做NRZ hold。Nbits 2000; % 比特数 rng(42); % 固定种子结果可复现 bits 2*(rand(1, Nbits) 0.5) - 1; % 等概率 ±1 % —— 1) NRZ保持每个比特重复ovs次 —— nrz reshape(repmat(bits, ovs, 1), 1, []); % —— 2) 高斯滤波得到频率轨迹 —— g gaussdesign(bwT, L, ovs); freq filter(g, 1, nrz); % 因果滤波存在群延迟 % —— 3) 补偿FIR对称延迟把判决点拉回比特网格 —— delay (length(g) - 1) / 2; % g 对称延迟为半窗长 freq [freq(delay1:end), zeros(1, delay)]; % —— 4) 相位累加 复指数 —— phase cumsum(freq) / ovs * (pi/2); % 每比特相位增量归一化到±π/2 s exp(1j * phase); % 复基带GMSK信号逐项说明第1步把每个比特扩成 ovs 个相同采样这是NRZ而不是冲激序列物理意义是持频一个周期第2步filter是因果卷积输出头部会晚delay个采样不补偿的话后续所有判决点整体平移误码率直接变0.5第3步左移delay个采样末尾补零只影响最后一点点相位统计误码时本来就要裁边。cumsum是离散积分/ ovs把每比特内的样本数折算成比特周期乘以π/2是因为调制指数 h0.5即每个比特相位变化恰为 ±90°。exp(1j*phase)得到单位功率的复基带信号mean(abs(s).^2)应等于1这是后面加噪时能量归一化的前提。验证调制做对没有最简单的是看相位轨迹的平滑度。画plot(unwrap(angle(s)))BT0.3 时应该是连续S形曲线而不是折线换成bwT100重跑一遍折线特征立刻出现这就是向MSK退化。3. GMSK解调差分与相干两条路线的实现和取舍3.1 一比特差分为什么它和GMSK是天生一对GMSK的信息不在绝对相位里而在每个比特周期结束时的相对相位增量1 比特让相位增加 π/2-1 比特让相位减小 π/2。这决定了差分解调是成本最低的天然匹配方案——接收机把当前采样与延迟一个比特的采样共轭相乘相位差就直接暴露了数据。与DPSK不同MSK/GMSK用一比特差分检测时不需要在发送端做差分编码因为相位增量本身就是信息这是个常见的概念混淆点。% 加AWGN信道Eb/N010 dB 时任取一个信噪比验证 EbN0_dB 10; rho 10^(EbN0_dB/10); % 线性Eb/N0 sigma sqrt(ovs / rho); % 每采样复噪声标准差 nz sigma/sqrt(2) * (randn(1, length(s)) 1j*randn(1, length(s))); r s nz; % 接收信号单位功率 % —— 一比特差分解调 —— rd r .* conj([zeros(1, ovs), r(1:end-ovs)]); % 延迟一个比特后共轭相乘 idx 2*ovs : ovs : length(r) - ovs; % 取每个比特结束时刻的差值 d_hat imag(rd(idx)) 0; % 相位增量为正→1 tx bits(2 : length(d_hat)1); % 判决点从第2比特开始对齐 ber mean(d_hat ~ (tx 0)); % 误码统计代码逻辑rd(n)是 r(n) 与 r(n-ovs) 的共轭乘积它的辐角就是跨越一个比特周期的相位差。判决点取2*ovs : ovs : end即每个比特周期的末端此时相位增量已经完整积分完比取比特中心更严格。噪声方差sigma sqrt(ovs/rho)的推导见第4章这里先记住复基带每采样噪声功率是sigma^2信号功率是1。若误码率算出来在0.5附近先查对齐idx起点、tx的下标、以及发送端freq的延迟补偿三处只要错一处结果就是随机的。3.2 两比特差分BT0.3时更稳的退路当 BT 压到 0.3高斯脉冲的能量拖到前后各2~3个比特一比特差分的判决点被ISI污染误码率曲线在高信噪比段出现平台。常见做法是把积分窗加长到两个比特跨两个比特的相位差只有三种结果两个比特同号则总增量为 ±π异号则互相抵消为0。判决看实部而非虚部得到的是相邻比特是否同号的翻转标志还需要做一次递推恢复原始数据。r2e r .* conj([zeros(1, 2*ovs), r(1:end-2*ovs)]); % 两比特延迟共轭相乘 idx2 3*ovs : ovs : length(r) - ovs; u_hat real(r2e(idx2)) 0; % 同号→1异号→0 % 差分译码由翻转标志递推原比特±1域 u_hat 2*u_hat - 1; % 转±11同号-1异号 d_hat2 ones(1, length(u_hat)); for k 2 : length(u_hat) d_hat2(k) d_hat2(k-1) * u_hat(k-1); % 当前符号 前一符号 × 翻转标志 end ber2 mean(d_hat2(3:end) ~ (bits(1:length(d_hat2)-2) 0));递推里有个相位模糊初始符号是任取的恢复出的序列可能整体取反所以统计前要裁掉首尾并对齐。两比特差分比一比特多付出一个比特的观测延迟、多一重递推换来的是BT0.3时显著更低的地板这是GSM时代接收机设计里真实用过的取舍写进报告比只报一条曲线更有说服力。3.3 相干解调的简版实现与性能对照仿真环境里可以免费拥有理想载波相位这是实测给不了的便宜。简版相干接收按观察窗 候选相位距离实现已知上一比特结束时的相位当前比特只有 ±π/2 两种走向把接收信号在比特中心的相位与两个候选相位分别求距离选近的那个。这套走查式判决丢掉积分增益但足够演示相干至少不比差分差的趋势严格最优要上Laurent分解加Viterbi均衡那是课设加分项。phi_ref 0; % 理想初始相位 bits_co zeros(1, Nbits); for k 3 : Nbits r_center r((k-1)*ovs round(ovs/2)); % 比特中心采样 phi_m angle(r_center); % 两种候选相位在上一状态上 ±π/2用复平面距离而非直接角度差 d_plus abs(angle(exp(1j*(phi_m - phi_ref - pi/2)))); d_minus abs(angle(exp(1j*(phi_m - phi_ref pi/2)))); bits_co(k) d_plus d_minus; phi_ref phi_ref (2*bits_co(k)-1)*pi/2; % 用判决结果更新参考相位 end ber_co mean(bits_co(3:end) ~ (bits(3:end) 0));exp(1j*(...))再取辐角是为了避免角度差跨越 ±π 时出现±2π的跳变比直接写abs(phi_m - phi_ref - pi/2)稳健。真实系统里理想相位并不存在需要Costas环或数据辅助估计所以这个数值只作性能上界。三条路线的对照关系如下解调方式复杂度需要载波同步BT0.3时表现一比特差分最低纯复数乘加不需要10dB附近约10⁻³量级两比特差分低多一重递推不需要高SNR优于一比特地板更低相干(理想相位)中需要(仿真免费)比差分再低1~2dB4. 从跑通到成报告参数表和误码率测量脚本4.1 先把参数写死结果才有说服力实验报告最怕我调了个参数跑出来了但说不清出处。下面这张表是脚本里所有应当显式声明的量Mente系数和rng种子一起写进报告复现性和说服力同时到位。参数推荐取值说明BT0.3 / 0.5 / 100(近似∞)三档对照报告核心变量Nbits≥ 1e5扫描BER时统计到1e-3需要约10倍错误样本ovs16低于8会直接看到眼图闭合、BER爬升L4截断小于3频谱出现截断纹波Eb/N0范围0~14 dB步进2 dB覆盖差分检测的可测区间蒙特卡洛每组独立3次取平均消除种子偶然性画errorbar更好Nbits 这里有个矛盾调制器演示跑2000比特足够画BER曲线至少1e5比特。建议脚本里把Nbits做成变量演示和扫描分开别用同一个参数跑到底。4.2 Eb/N0到复基带噪声方差最容易翻车的换算许多同学直接在复基带信号上加sigmasqrt(1/(2*ratio))之类的高斯白噪声导致BER曲线的横轴其实是信噪比而不是 Eb/N0和别人结果无法对比。这里给出自洽的换算信号s已归一化到单位功率每个比特能量 Eb P·Tb 1Tb 在代码里是归一化的1。离散复基带里每采样噪声方差与 Eb/N0 的关系是σ² ovs / ρρ 10^(Eb/N0/10)推导很朴素匹配滤波对 ovs 个样本取平均后信号幅度保持1噪声方差变成 σ²/ovs于是滤波后信噪比 1/(σ²/ovs) ρ恰为目标的 Eb/N0。反过来验算一遍Eb/N010dB即ρ10ovs16时 σ²1.6匹配滤波后SNR16/1.610 ✓。代码里按实部虚部分量各取sigma/sqrt(2)即本章前一节nz那行。4.3 BER扫描脚本与曲线判读把第3章的一比特差分包进循环就是完整的扫描实验EbN0_dB 0 : 2 : 14; ber zeros(size(EbN0_dB)); for m 1 : length(EbN0_dB) rho 10^(EbN0_dB(m)/10); sigma sqrt(ovs / rho); nz sigma/sqrt(2) * (randn(1, length(s)) 1j*randn(1, length(s))); r s nz; rd r .* conj([zeros(1, ovs), r(1:end-ovs)]); idx 2*ovs : ovs : length(r) - ovs; d_hat imag(rd(idx)) 0; tx bits(2 : length(d_hat)1); ber(m) mean(d_hat ~ (tx 0)); end semilogy(EbN0_dB, ber, o-); grid on; xlabel(E_b/N_0 (dB)); ylabel(BER);判读时有三个经验值可以自检。其一MSK一比特差分的理论误码率约 0.5·exp(-Eb/N0)8dB时约7e-4GMSK在BT0.3时还会恶化1~2dB所以一比特差分曲线应在10dB附近进入10⁻³量级。其二若曲线在某个高信噪比处掉不下去说明存在错误平台先怀疑ISI或对齐而不是噪声。其三把三档BT的曲线画在一起MSK应最低、BT0.3应最高但对应频谱恰恰相反——这组对比就是报告结论的素材。改bwT时记得重新生成g和delay否则曲线全部错位。5. 交报告前自检眼图、频谱与三个高频翻车点5.1 三张图证明调制解调链路自洽误码率能说明能用但有说服力的报告还需要波形证据三张图各管一件事。眼图证明解调判决点选对了把差分解调输出imag(rd)喂给eyediagram横轴一个比特周期BT0.3时眼图应明显收窄但仍张开张开度随Eb/N0提高而改善。频谱证明高斯滤波器在起作用pwelch(s, hamming(1024), 512, 1024, ovs)横轴以比特率为单位对照BT0.3和100两条曲线主瓣收窄和旁瓣滚降的差异一目了然。相位轨迹证明调制模型没建错plot(unwrap(angle(s)))平滑S形是GMSK折线是MSK。这三张图各自只对应一个plot系列函数冲突起来时优先查图而不是查代码。5.2 三个高频翻车点和它们的样子第一个是滤波器延迟没补偿。症状是误码率恒在0.5附近、眼图相位整体平移半个码元周期。检查方法是把发送bits和d_hat画在同一张图错位对齐看错误集中在固定偏移处就是这个问题。第二个是把复基带GMSK当实信号用sign(real(r))直接判决。GMSK相位每比特走±90°实部在0附近徘徊这种判决等价于把圆上的信号压到一条线上任何信噪比下都没正确率。第三个是噪声方差没按第4.2节的公式换算或者加噪前忘了验证mean(abs(s).^2)1表现为BER曲线横向漂移2~3dB和理论对不上。最后提供一个快速定位工具给接收信号人为加一个小的载波频偏r_f r .* exp(1j*2*pi*df*Tb*(0:length(r)-1)/ovs)用一比特差分检测扫df*Tb从0到0.1观察BER随频偏的上升斜率。一比特差分的判决相位里会多出一项 2π·Δf·Tb 的固定偏置当它接近 π/2 时判决直接失效这个实验可以1分钟讲清楚频偏容限为什么是比特率的几十分之一比背公式直观得多。把这段结果放在报告最后就算参数选得中庸这份工作量的含金量也够了。本文还有配套的精品资源点击获取 GMSK调制解调MATLAB实现：从原理到误码率分析
简介这是一份基于MATLAB的GMSK调制解调实验课程设计文档面向通信工程专业本科生及需要完成同类课题的开发者系统讲解了高斯滤波最小频移键控的调制解调原理、Simulink模型搭建与瑞利信道下的仿真验证。文档共1个docx文件压缩包大小约1.81MB篇幅完整包含摘要、设计原理、功能模块设计、波形/眼图/误码率结果分析和结论等结构可直接作为课程设计说明书模板或仿真实战参考。内容从GMSK基本概念出发覆盖信号发生、调制解调、误码率统计和波形观察四个核心模块并重点分析BT参数对频谱效率和抗干扰能力的影响帮助读者在Simulink中复现完整链路并理解关键性能指标。已有464人学习下载适合希望系统掌握GMSK仿真流程、撰写实验报告或搭建通信系统模型的读者。1. GMSK调制解调实验到底要你交出什么GMSK是GSM时代蜂窝物理层的事实标准蓝牙里的GFSK也是它的近亲理解GMSK等于同时拿下连续相位调制 频谱成形两个通信基本功。用MATLAB做这个实验目的不是造一个能发射的信号而是看清三件事高斯滤波器如何用带宽换频谱纯度信息如何藏在每比特±π/2的相位增量里以及差分检测为什么在无线调制和解调电路原理课程里被反复强调。适合通信原理课设、准备面试的本科生或想把手写调制器方法捡起来的工程师。下面给的代码不挑MATLAB版本装好即可运行只有一处用到通信工具箱会给出无工具箱替代写法。2. GMSK调制把频率轨迹从折线磨成平滑曲线的建模与代码2.1 MSK到GMSK高斯滤波器到底改在哪一处MSK是连续相位FSK里最容易被接受的基准调制指数 h0.5瞬时频率在每个比特周期内恒定比特边界只会发生相位连续但频率跳变。相位轨迹是折线频率轨迹是方波。方波意味着高频分量丰富频谱旁瓣按 f⁻⁴ 衰减对相邻信道干扰大。GMSK在频率成形之前插入一个高斯型低通滤波器把方波形的瞬时频率磨成高斯脉冲形状的过渡曲线。更准确地说GMSK的基带等效模型是NRZ数据先经过预调制高斯滤波器输出作为瞬时频率再积分得到相位。高斯滤波器传递函数是H(f) exp( -ln2 · f² / (2B²) )其中 B 是该滤波器的 3dB 带宽BT乘积就是 B 和比特周期 Tb 的乘积。BT 是GMSK唯一真正需要手工调的参数它同时决定频谱和误码率。BT 越大滤波越弱频谱越接近MSKBT 越小频谱越紧但单个比特的能量被拉到相邻比特上眼图闭合得越厉害。这张取舍表通常直接写进实验报告BT 取值典型场景频谱特征代价∞即MSK理论基准主瓣最宽旁瓣按 f⁻⁴ 衰减无ISI眼图全开0.5蓝牙GFSK同族主瓣收窄旁瓣较低ISI轻微一比特差分可用0.3GSM主瓣窄旁瓣极低ISI明显解调需差分或均衡实验里最常见的做法是把同一套代码跑三遍BT 取 0.3、0.5 和一个很大的值近似MSK用 pwelch 量频谱、用误码率曲线看代价这组对照本身就有讨论价值。2.2 两种生成高斯成形脉冲的MATLAB写法调制器里唯一需要制造的元件就是高斯预调制滤波器。手头有通信工具箱时直接用gaussdesign它返回一个单位DC增益的对称FIRsum(g) 1这个性质后面归一化相位要用。bwT 0.3; % BT乘积 L 4; % 截断长度单位比特周期 ovs 16; % 每比特采样点数 g gaussdesign(bwT, L, ovs); % 单位DC增益sum(g)1没有通信工具箱时可以按高斯滤波器的解析式自己采样生成。3dB带宽B对应的高斯标准差是σ sqrt(ln2) / (2·π·B)把它作为滤波器时域宽度直接写函数function h gmsk_gauss_lpf(bwT, L, ovs) % 手写高斯预调制滤波器与gaussdesign等效 Tb 1; % 归一化比特周期 sigma sqrt(log(2)) / (2*pi*bwT); % 3dB带宽→标准差 t ((0 : L*ovs) - L*ovs/2) / ovs; % 时间轴单位比特周期 h exp(-t.^2 / (2*sigma^2)); % 高斯冲激响应 h h / sum(h); % 强制单位DC增益 end两个版本的共同点是sum(h)1。这个性质不是形式主义它保证一个持续整个比特周期的1输入滤波后在该比特周期内的积分恰好是1于是相位累加器里每个比特贡献的相位增量恰为 ±π/2。节外生枝说一句gaussdesign的BT是3dB带宽不是噪声等效带宽别拿0.5/Tb这种经验值去替代。2.3 NRZ保持 相位累加最小可运行调制器完整的复基带调制器分四步生成 ±1 数据 → 把每个比特保持 ovs 个采样得到NRZ波形 → 过高斯滤波器得到频率轨迹 → 累加积分并映射到复指数。注意一个新手高发错误不要用upsample造冲激串再接高斯滤波器那等于拿滤波器的冲激响应当频率脉冲在BT较大时相位轨迹和频谱都会失真必须先做NRZ hold。Nbits 2000; % 比特数 rng(42); % 固定种子结果可复现 bits 2*(rand(1, Nbits) 0.5) - 1; % 等概率 ±1 % —— 1) NRZ保持每个比特重复ovs次 —— nrz reshape(repmat(bits, ovs, 1), 1, []); % —— 2) 高斯滤波得到频率轨迹 —— g gaussdesign(bwT, L, ovs); freq filter(g, 1, nrz); % 因果滤波存在群延迟 % —— 3) 补偿FIR对称延迟把判决点拉回比特网格 —— delay (length(g) - 1) / 2; % g 对称延迟为半窗长 freq [freq(delay1:end), zeros(1, delay)]; % —— 4) 相位累加 复指数 —— phase cumsum(freq) / ovs * (pi/2); % 每比特相位增量归一化到±π/2 s exp(1j * phase); % 复基带GMSK信号逐项说明第1步把每个比特扩成 ovs 个相同采样这是NRZ而不是冲激序列物理意义是持频一个周期第2步filter是因果卷积输出头部会晚delay个采样不补偿的话后续所有判决点整体平移误码率直接变0.5第3步左移delay个采样末尾补零只影响最后一点点相位统计误码时本来就要裁边。cumsum是离散积分/ ovs把每比特内的样本数折算成比特周期乘以π/2是因为调制指数 h0.5即每个比特相位变化恰为 ±90°。exp(1j*phase)得到单位功率的复基带信号mean(abs(s).^2)应等于1这是后面加噪时能量归一化的前提。验证调制做对没有最简单的是看相位轨迹的平滑度。画plot(unwrap(angle(s)))BT0.3 时应该是连续S形曲线而不是折线换成bwT100重跑一遍折线特征立刻出现这就是向MSK退化。3. GMSK解调差分与相干两条路线的实现和取舍3.1 一比特差分为什么它和GMSK是天生一对GMSK的信息不在绝对相位里而在每个比特周期结束时的相对相位增量1 比特让相位增加 π/2-1 比特让相位减小 π/2。这决定了差分解调是成本最低的天然匹配方案——接收机把当前采样与延迟一个比特的采样共轭相乘相位差就直接暴露了数据。与DPSK不同MSK/GMSK用一比特差分检测时不需要在发送端做差分编码因为相位增量本身就是信息这是个常见的概念混淆点。% 加AWGN信道Eb/N010 dB 时任取一个信噪比验证 EbN0_dB 10; rho 10^(EbN0_dB/10); % 线性Eb/N0 sigma sqrt(ovs / rho); % 每采样复噪声标准差 nz sigma/sqrt(2) * (randn(1, length(s)) 1j*randn(1, length(s))); r s nz; % 接收信号单位功率 % —— 一比特差分解调 —— rd r .* conj([zeros(1, ovs), r(1:end-ovs)]); % 延迟一个比特后共轭相乘 idx 2*ovs : ovs : length(r) - ovs; % 取每个比特结束时刻的差值 d_hat imag(rd(idx)) 0; % 相位增量为正→1 tx bits(2 : length(d_hat)1); % 判决点从第2比特开始对齐 ber mean(d_hat ~ (tx 0)); % 误码统计代码逻辑rd(n)是 r(n) 与 r(n-ovs) 的共轭乘积它的辐角就是跨越一个比特周期的相位差。判决点取2*ovs : ovs : end即每个比特周期的末端此时相位增量已经完整积分完比取比特中心更严格。噪声方差sigma sqrt(ovs/rho)的推导见第4章这里先记住复基带每采样噪声功率是sigma^2信号功率是1。若误码率算出来在0.5附近先查对齐idx起点、tx的下标、以及发送端freq的延迟补偿三处只要错一处结果就是随机的。3.2 两比特差分BT0.3时更稳的退路当 BT 压到 0.3高斯脉冲的能量拖到前后各2~3个比特一比特差分的判决点被ISI污染误码率曲线在高信噪比段出现平台。常见做法是把积分窗加长到两个比特跨两个比特的相位差只有三种结果两个比特同号则总增量为 ±π异号则互相抵消为0。判决看实部而非虚部得到的是相邻比特是否同号的翻转标志还需要做一次递推恢复原始数据。r2e r .* conj([zeros(1, 2*ovs), r(1:end-2*ovs)]); % 两比特延迟共轭相乘 idx2 3*ovs : ovs : length(r) - ovs; u_hat real(r2e(idx2)) 0; % 同号→1异号→0 % 差分译码由翻转标志递推原比特±1域 u_hat 2*u_hat - 1; % 转±11同号-1异号 d_hat2 ones(1, length(u_hat)); for k 2 : length(u_hat) d_hat2(k) d_hat2(k-1) * u_hat(k-1); % 当前符号 前一符号 × 翻转标志 end ber2 mean(d_hat2(3:end) ~ (bits(1:length(d_hat2)-2) 0));递推里有个相位模糊初始符号是任取的恢复出的序列可能整体取反所以统计前要裁掉首尾并对齐。两比特差分比一比特多付出一个比特的观测延迟、多一重递推换来的是BT0.3时显著更低的地板这是GSM时代接收机设计里真实用过的取舍写进报告比只报一条曲线更有说服力。3.3 相干解调的简版实现与性能对照仿真环境里可以免费拥有理想载波相位这是实测给不了的便宜。简版相干接收按观察窗 候选相位距离实现已知上一比特结束时的相位当前比特只有 ±π/2 两种走向把接收信号在比特中心的相位与两个候选相位分别求距离选近的那个。这套走查式判决丢掉积分增益但足够演示相干至少不比差分差的趋势严格最优要上Laurent分解加Viterbi均衡那是课设加分项。phi_ref 0; % 理想初始相位 bits_co zeros(1, Nbits); for k 3 : Nbits r_center r((k-1)*ovs round(ovs/2)); % 比特中心采样 phi_m angle(r_center); % 两种候选相位在上一状态上 ±π/2用复平面距离而非直接角度差 d_plus abs(angle(exp(1j*(phi_m - phi_ref - pi/2)))); d_minus abs(angle(exp(1j*(phi_m - phi_ref pi/2)))); bits_co(k) d_plus d_minus; phi_ref phi_ref (2*bits_co(k)-1)*pi/2; % 用判决结果更新参考相位 end ber_co mean(bits_co(3:end) ~ (bits(3:end) 0));exp(1j*(...))再取辐角是为了避免角度差跨越 ±π 时出现±2π的跳变比直接写abs(phi_m - phi_ref - pi/2)稳健。真实系统里理想相位并不存在需要Costas环或数据辅助估计所以这个数值只作性能上界。三条路线的对照关系如下解调方式复杂度需要载波同步BT0.3时表现一比特差分最低纯复数乘加不需要10dB附近约10⁻³量级两比特差分低多一重递推不需要高SNR优于一比特地板更低相干(理想相位)中需要(仿真免费)比差分再低1~2dB4. 从跑通到成报告参数表和误码率测量脚本4.1 先把参数写死结果才有说服力实验报告最怕我调了个参数跑出来了但说不清出处。下面这张表是脚本里所有应当显式声明的量Mente系数和rng种子一起写进报告复现性和说服力同时到位。参数推荐取值说明BT0.3 / 0.5 / 100(近似∞)三档对照报告核心变量Nbits≥ 1e5扫描BER时统计到1e-3需要约10倍错误样本ovs16低于8会直接看到眼图闭合、BER爬升L4截断小于3频谱出现截断纹波Eb/N0范围0~14 dB步进2 dB覆盖差分检测的可测区间蒙特卡洛每组独立3次取平均消除种子偶然性画errorbar更好Nbits 这里有个矛盾调制器演示跑2000比特足够画BER曲线至少1e5比特。建议脚本里把Nbits做成变量演示和扫描分开别用同一个参数跑到底。4.2 Eb/N0到复基带噪声方差最容易翻车的换算许多同学直接在复基带信号上加sigmasqrt(1/(2*ratio))之类的高斯白噪声导致BER曲线的横轴其实是信噪比而不是 Eb/N0和别人结果无法对比。这里给出自洽的换算信号s已归一化到单位功率每个比特能量 Eb P·Tb 1Tb 在代码里是归一化的1。离散复基带里每采样噪声方差与 Eb/N0 的关系是σ² ovs / ρρ 10^(Eb/N0/10)推导很朴素匹配滤波对 ovs 个样本取平均后信号幅度保持1噪声方差变成 σ²/ovs于是滤波后信噪比 1/(σ²/ovs) ρ恰为目标的 Eb/N0。反过来验算一遍Eb/N010dB即ρ10ovs16时 σ²1.6匹配滤波后SNR16/1.610 ✓。代码里按实部虚部分量各取sigma/sqrt(2)即本章前一节nz那行。4.3 BER扫描脚本与曲线判读把第3章的一比特差分包进循环就是完整的扫描实验EbN0_dB 0 : 2 : 14; ber zeros(size(EbN0_dB)); for m 1 : length(EbN0_dB) rho 10^(EbN0_dB(m)/10); sigma sqrt(ovs / rho); nz sigma/sqrt(2) * (randn(1, length(s)) 1j*randn(1, length(s))); r s nz; rd r .* conj([zeros(1, ovs), r(1:end-ovs)]); idx 2*ovs : ovs : length(r) - ovs; d_hat imag(rd(idx)) 0; tx bits(2 : length(d_hat)1); ber(m) mean(d_hat ~ (tx 0)); end semilogy(EbN0_dB, ber, o-); grid on; xlabel(E_b/N_0 (dB)); ylabel(BER);判读时有三个经验值可以自检。其一MSK一比特差分的理论误码率约 0.5·exp(-Eb/N0)8dB时约7e-4GMSK在BT0.3时还会恶化1~2dB所以一比特差分曲线应在10dB附近进入10⁻³量级。其二若曲线在某个高信噪比处掉不下去说明存在错误平台先怀疑ISI或对齐而不是噪声。其三把三档BT的曲线画在一起MSK应最低、BT0.3应最高但对应频谱恰恰相反——这组对比就是报告结论的素材。改bwT时记得重新生成g和delay否则曲线全部错位。5. 交报告前自检眼图、频谱与三个高频翻车点5.1 三张图证明调制解调链路自洽误码率能说明能用但有说服力的报告还需要波形证据三张图各管一件事。眼图证明解调判决点选对了把差分解调输出imag(rd)喂给eyediagram横轴一个比特周期BT0.3时眼图应明显收窄但仍张开张开度随Eb/N0提高而改善。频谱证明高斯滤波器在起作用pwelch(s, hamming(1024), 512, 1024, ovs)横轴以比特率为单位对照BT0.3和100两条曲线主瓣收窄和旁瓣滚降的差异一目了然。相位轨迹证明调制模型没建错plot(unwrap(angle(s)))平滑S形是GMSK折线是MSK。这三张图各自只对应一个plot系列函数冲突起来时优先查图而不是查代码。5.2 三个高频翻车点和它们的样子第一个是滤波器延迟没补偿。症状是误码率恒在0.5附近、眼图相位整体平移半个码元周期。检查方法是把发送bits和d_hat画在同一张图错位对齐看错误集中在固定偏移处就是这个问题。第二个是把复基带GMSK当实信号用sign(real(r))直接判决。GMSK相位每比特走±90°实部在0附近徘徊这种判决等价于把圆上的信号压到一条线上任何信噪比下都没正确率。第三个是噪声方差没按第4.2节的公式换算或者加噪前忘了验证mean(abs(s).^2)1表现为BER曲线横向漂移2~3dB和理论对不上。最后提供一个快速定位工具给接收信号人为加一个小的载波频偏r_f r .* exp(1j*2*pi*df*Tb*(0:length(r)-1)/ovs)用一比特差分检测扫df*Tb从0到0.1观察BER随频偏的上升斜率。一比特差分的判决相位里会多出一项 2π·Δf·Tb 的固定偏置当它接近 π/2 时判决直接失效这个实验可以1分钟讲清楚频偏容限为什么是比特率的几十分之一比背公式直观得多。把这段结果放在报告最后就算参数选得中庸这份工作量的含金量也够了。本文还有配套的精品资源点击获取