Costas环载波同步仿真：BPSK/QPSK/MSK/GMSK的Simulink实现

简介：这是一套面向通信与信号处理方向学习者的 MATLAB/Simulink 仿真资源，重点围绕 MSK、GMSK、QPSK、BPSK 四种调制方式下的 Costas 环载波同步问题，提供可直接运行的仿真模型，适合本科、硕士阶段的课程作业、科研入门以及教师备课参考。压缩包内共 14 个文件，以 8 个 mdl 仿真模型为主体，覆盖上述四种调制方式的多种环路实现，另有 m 脚本用于基带信号生成，txt 与 doc 文档补充参数说明和使用注意，png 图片保存了运行结果波形，压缩包整体大小约 1013KB，文件组织清晰、针对性明确。该资源已有 355 人学习/下载。代码基于 MATLAB 2014/2019a 编写，内含运行结果，打开模型即可查看不同调制方式下 Costas 环的同步效果；同时文档中对模型结构、信号流程等做了配套说明，便于理解载波同步环路中鉴相器、环路滤波器、压控振荡器等模块的连接与参数设置。对于正在研究数字调制解调、锁相环、Costas 环，或需要在 Simulink 中完成相关仿真的读者，这套文件能提供直观的模型参考、调试思路和结果对比，实用价值较强。

1. 突发通信里为什么离不开 Costas 环

做无线通信物理层仿真的人，迟早会撞上载波同步这道坎。发射端把基带信号搬到中频或射频，接收端下变频时本振不可能和发射端完全同频同相，残留频偏和相偏会让星座图转圈。BPSK、QPSK 这类恒包络或准恒包络调制，信息全部承载在相位上，相位一转，解调就废。MSK 和 GMSK 虽然相位连续，但对频偏同样敏感，只是表现方式不同——星座图不转圈，而是眼图张开度变差、误码率抬出一条平台。

Costas 环的价值在于它不用额外导频，直接从调制信号里把载波误差提出来。这份资源里就是一套基于 Simulink 的搭好的 Costas 环解调模型，覆盖 BPSK、QPSK、MSK、GMSK 四种调制方式，配了配套的 m 脚本和说明文档。适合三类人：刚接触载波同步、想看看环路在 Simulink 里怎么落地的新手；做课程设计或毕业设计、需要可跑的仿真框架的本科生和研究生；以及做软件无线电或通信算法验证、想快速对比不同调制方式下环路性能差异的工程师。模型从 MATLAB 2014a 到 2019a 都能打开，拿到手改改参数就能跑。

2. Costas 环的结构拆解与环路滤波器参数设计

2.1 三种鉴相器的工作机制

Costas 环核心就三块：鉴相器、环路滤波器、压控振荡器（Simulink 里通常用 VCO 模块或自己搭的 NCO）。鉴相器看起来是简单的乘加运算，但不同调制方式的鉴相表达式差异很大，直接决定环路能不能锁定。

对 BPSK，接收信号与本地载波相乘后分成 I/Q 两路，鉴相误差为sign(I) * Q。这个表达式的含义是：I 路先做硬判决，判决结果再和 Q 路相乘。当环路锁定且无噪声时，Q 路趋于零，误差信号也趋于零；频偏存在时，误差信号会随相位差线性变化，经过环路滤波器后控制 VCO 拉回。对 QPSK，鉴相误差为sign(I) * Q - sign(Q) * I，相当于把四象限符号判决纳入误差计算。

MSK 和 GMSK 的鉴相器要复杂得多。MSK 本质上是偏移 QPSK（OQPSK），I/Q 两路在时间上错开一个符号周期 Tb，叠加了半正弦加权。直接套用 QPSK 的鉴相结构当然也能锁，但误差信号里会混入 ISI 分量，导致稳态相位抖动偏大。通常的做法是在鉴相前加匹配滤波器，或者在误差输出端加低通，让符号间干扰被平均掉。GMSK 因为高斯滤波进一步压缩了相位路径，鉴相器输出的误差信号更弱，需要环路滤波器提供更高的直流增益。

2.2 环路滤波器系数的手算与脚本化

环路滤波器是 Costas 环稳定性的命门。一阶环只能跟踪固定频偏，二阶环可以跟踪频率斜升。Simulink 里最常见的结构是比例积分滤波器，传递函数为：

F(s) = (tau2 * s + 1) / (tau1 * s)

离散化后就是用两个乘法器加一个累加器搭成的 IIR 结构。环路带宽 BL 和阻尼系数 ξ 直接决定捕获速度与稳态抖动，经典设计中取 ξ = 0.707，环路带宽取符号速率的 1% 到 3%。

% costas_loop_param.m % 输入：符号速率 Rs，环路带宽因子 BL_factor，采样率 Fs % 输出：比例系数 Kp，积分系数 Ki Rs = 1e6; % 1 Msps BL_factor = 0.02; % 环路带宽 = 2% * Rs zeta = 0.707; % 阻尼系数 Kd = 1; % 鉴相器增益（BPSK时为1） K0 = 1; % VCO增益归一化 BL = BL_factor * Rs; wn = BL / (zeta + 1/(4*zeta)); % 自然角频率 K = Kd * K0; % 双线性变换后的数字滤波器系数 T = 1 / Fs; c1 = 2*zeta*wn*T / K; c2 = (wn*T)^2 / K;

这段脚本算出的c1和c2就是 Simulink 中比例支路和积分支路的增益。需要注意，Fs是环路更新速率，不是系统采样率。如果环路滤波器放在 1 MHz 符号率下，但模型步长是 1e-8 秒，那么Fs应该设为实际环路更新频率，否则系数差几个数量级，环路要么不收敛要么振荡。实际调参时我会先在 MATLAB 命令行跑一遍这段脚本，把c1、c2打印出来再填进 Simulink 常量模块，而不是在模型里现算。

2.3 VCO/NCO 的建模差异

Simulink 自带的Voltage-Controlled Oscillator模块直接把控制电压映射到输出频率，但在数字通信仿真里更推荐自己搭 NCO。原因是 VCO 模块的输出是模拟连续相位，而数字接收机里通常用相位累加器 + 正余弦查找表的结构，方便观测相位误差和做时钟同步。

NCO 参数据实如下：phase_increment对应中心频率 f_c，控制字delta_phase的更新公式为delta_phase_new = f_c/Fs + c1*error + c2*integral。VCO 模块的Quiescent frequency要设为下变频后的残余频偏（通常是 0 或一个很小的值），Output amplitude设为 1。如果模型中 I/Q 两路共用同一个 NCO 输出，注意 Sine 和 Cosine 的相位差必须是 90 度，不要用两个独立 VCO 模块分别生成，否则两路载波相位会漂移。

3. 从 BPSK 到 GMSK：Simulink 模型里改了什么

3.1 BPSK 与 QPSK 的 Costas 环模型结构

压缩包里的costas_pYll_BPSK.mdl是最底的参考模型。结构上，接收信号进Product模块和本地载波相乘，乘出来的 I 路和 Q 路分别过低通滤波器，I 路进Sign模块做硬判决，判决输出再和 Q 路相乘得到误差信号。误差信号经环路滤波器后送到 VCO/NCO，形成闭合回路。

QPSK 的模型costas_pYll_qpsk_1.mdl、costas_pYll_qpsk_2.mdl、costas_pYll_qpsk_3.mdl是三版不同的实现。qpsk_1用的是经典四相乘子结构，qpsk_2改成复数乘法加angle提取相位差，qpsk_3是加了自动增益控制（AGC）的版本。实际用下来，qpsk_2在低信噪比下抖动最小，因为atan2型鉴相器不会因符号判决错误产生 180 度翻转。

% 设置 BPSK 模型的仿真参数（MATLAB 命令行执行） % 用 sim() 函数批量跑不同 Eb/N0 下的误码率 EbN0_dB = 0:2:10; ber = zeros(size(EbN0_dB)); for k = 1:length(EbN0_dB) assignin('base', 'EbN0_current', EbN0_dB(k)); simOut = sim('costas_pYll_BPSK.mdl', 'StopTime', '1e-4'); % 对比发送 bit 与解调 bit，统计误码 rxBits = simOut.rxBits(100:end); % 丢掉前100个符号，避开捕获段 txBits = simOut.txBits(100:end); ber(k) = sum(rxBits ~= txBits) / length(rxBits); end

这段脚本说明了一点：环路捕获期间的错判不能计入误码统计，必须丢弃前导符号。BPSK 在正常参数下几百个符号就能锁定，但 QPSK 在低信噪比下可能需要上千符号。私信或文档里提到的「内含运行结果」，指的就是simOut里保存的误码率曲线和星座图数据。

3.2 MSK 的 I/Q 时序错位处理

costas_pYll_msk1.mdl和costas_pYll_msk2.mdl对应两种 MSK 实现。第一种是正交法，I 路和 Q 路分别用半正弦脉冲成形，两路错开 Tb 输入；第二种是相位累加法，直接生成连续相位信号再正交分解。对于 Costas 环来说，两种方式下鉴相器的输入信号等效，但要注意 I/Q 采样时序。

MSK 的符号周期 Tb 等于比特周期，I/Q 两路各传 0.5 个比特率，实际符号率是比特率的一半。错误做法是把 QPSK 模型直接照搬，没有把 Q 路延迟 Tb。这样环路能锁，但解调出来的数据码型是错的，星座图上四个点不落在 ±π/4 而是落在坐标轴上。检查方法：解调后把 I/Q 两路重新合成复数信号，画散点图，正确的 MSK 星座图应该呈现四个旋转点轨迹，而不是静止的四个点。

3.3 GMSK 的额外麻烦

costas_pYll_gmsk.mdl和costas_pYll_gmsk1.mdl的差异在环路滤波器带宽。GMSK 经过高斯滤波后相位路径变平滑，高频分量减少，鉴相器输出的误差信号频谱更集中。一个常见问题是：把 BPSK 的环路参数直接套到 GMSK 上，环路带宽偏大，噪声通过环路滤波器进入 VCO，星座图发散但误码率看着还行——错误的稳定。

GMSK 模型里通常要加一个phase unwrap处理，因为高斯滤波让相位变化率降低，atan2输出不会跨过 ±π 边界，但如果前端有残留频偏，相位累加器会溢出。我一般会在环路滤波器后加一个Modulo 2π模块，把相位误差限幅在 -π 到 π 之间，防止 VCO 控制字跳变。

表：四种调制的 Costas 环关键差异

GSMK 的 BT 值（通常 0.3 或 0.5）也影响环路设计。BT 越小，相位路径越平滑，需要的环路带宽越窄，捕获时间越长。如果仿真时间不够长，环路还没锁住仿真就结束了，输出误码率会虚高——先加长StopTime验证锁定了再统计。

4. 仿真收敛判据与四个高频排查点

4.1 怎么判断环路真锁定了

很多人把「环路滤波器输出不再变化」当作锁定判据，这对 BPSK 还行，对 QPSK/MSK 是不准确的。正确的做法是看 VCO 控制电压和发射机中心频率的差。锁定时这个差值应该等于初始频偏的估计值，而不是零。理想情况下，环路锁定后 VCO 输出频率 = 发射频率 - 残留频偏补偿量。

% 从 simOut 里提取 VCO 控制电压 t = simOut.tout; vco_ctrl = simOut.vcoCtrl; figure(1); plot(t, vco_ctrl); xlabel('时间 (s)'); ylabel('VCO 控制电压'); % 锁定后控制电压应趋于一条水平直线 % 直线的斜率应该约为0，否则表示频率还有斜升未消除

BPSK 在理想信道下锁定后vco_ctrl是常值，QPSK 会有小的纹波，MSK/GMSK 的纹波更大。如果vco_ctrl持续单调上升或下降而不收敛，说明环路带宽太小或初始频偏超出 VCO 的牵引范围。

4.2 高频排查点

排查一：初始频偏过大时环路失锁。Costas 环的捕获范围大约在 ±BL/2 以内，超过这个范围 VCO 拉不回来。解决方法是先扫频粗估频偏，把残余频偏压到捕获范围内再启动环路。Simulink 里可以加一个Freq Offset Estimator模块做粗同步，但更直接的做法是在发射端加一个已知导频段，用导频估频偏。

排查二：环路滤波器输出极性接反。误差信号经过滤波器进 VCO 时，如果极性反了，环路变成正反馈，VCO 只会加速偏离。症状是星座图越转越快。BPSK/QPSK 的sign(I)*Q鉴相器，极性方向要和 VCO 的频率控制方向匹配。检查方法：把初始相位误差设一个小值，观察误差信号是收敛到 0 还是发到 ±π。

排查三：同样一套模型，改采样率后不收敛。Simulink 里Sine Wave或VCO模块的输出采样率决定了 NCO 的相位步进。采样率从 1e6 改成 2e6 后，环路滤波器的Fs参数也要同步更新。手动改容易漏。建议把采样率设成基带工作区的全局变量，NCO、环路滤波器、示波器全部引用同一个变量。

排查四：离散模型步长太大导致相位在符号周期内跳变。固定步长仿真中，步长大于符号周期的 1/16 时，NCO 的相位量化噪声会直接进入鉴相器，环路稳定点出现随机偏移。通常把最大步长设为符号周期的 1/32 或更小。

4.3 捕获失败的快速实验方法

如果环路完全锁不住，先用排除法。把发射端的调制信号换成等幅单音，也就是把信息序列全置 1，让 BPSK 模型退化成纯载波恢复。单音下 Costas 环的鉴相器输出应该是一个光滑的误差曲线，VCO 控制电压快速收敛。单音能锁、调制不能锁，问题在鉴相器；单音都锁不住，问题在环路滤波器或 NCO 本身。这一步能省至少半天排查时间。

5. 环路带宽缩放与参数批跑技巧

5.1 按符号速率缩放环路参数

一套调试好的参数在 1 Mbps 下工作正常，换到 10 Mbps 不能直接用。正确的缩放关系：环路带宽 BL 和符号速率 Rs 的比值保持恒定，c1、c2按 Rs 的比例缩放。推导上，环路滤波器的积分时间常数随符号率缩短，所以积分系数要乘以缩放倍数。

% bpsk_loop_scale.m % 从 1Mbps 缩放到 10Mbps 的参数重算 Rs_old = 1e6; Rs_new = 10e6; scale = Rs_new / Rs_old; % 旧参数（在 1Mbps 下调好的） c1_old = 0.001; c2_old = 1e-5; % 新参数 c1_new = c1_old * scale; c2_new = c2_old * scale^2; fprintf('c1_new = %g, c2_new = %g\n', c1_new, c2_new);

c2按平方缩放的原因是积分支路的时间常数涉及两次积分运算，频率从 1M 变 10M，相位累加的积分时间缩短 10 倍，而误差积累的权重还要再乘一个时间因子。实际应用中，我会先按这个公式算出初始值，再在小范围内微调c2（正负 30% 是正常搜索区间）。

表：不同符号率下的经验参数（BPSK 模型，ξ=0.707）

表中的捕获时间与符号速率无关，因为环路带宽和符号率等比缩放。想让捕获更快就增大 BL_factor，代价是稳态相位抖动变大，星座图点弥散更明显。工程上一般先用大带宽快捕，锁定后再切小带宽精跟踪，这就是双模式环路。Simulink 里实现不难：用Switch模块切换两套滤波器系数，控制信号来自一个简单的定时器。

5.2 用脚本驱动 sim 做蒙特卡洛扫描

手工改 Simulink 参数再点运行，效率太低了。把需要扫描的变量全部assignin到 base workspace，用sim()在循环里跑，最后把结果汇总画图。注意每次sim()调用后要清理 workpace 里的临时变量，不然变量残留会影响下一次仿真。

% monte_carlo_costas.m % 扫描不同初始频偏下 QPSK Costas 环的捕获时间 freq_offsets = 100:100:1000; % 频偏值，单位 Hz capture_time = zeros(size(freq_offsets)); for k = 1:length(freq_offsets) assignin('base', 'freq_offset', freq_offsets(k)); simOut = sim('costas_pYll_qpsk_3.mdl', 'StopTime', '1e-3'); vco_ctrl = simOut.vcoCtrl; t = simOut.tout; % 捕获时间定义为控制电压第一次进入 ±5% 误差带 % 且之后不再超出 idx = find(abs(vco_ctrl - mean(vco_ctrl(end-100:end))) > 0.05*abs(mean(vco_ctrl(end-100:end))), 1, 'last'); if isempty(idx) capture_time(k) = t(end); else capture_time(k) = t(idx); end end figure(2); plot(freq_offsets, capture_time); xlabel('初始频偏 (Hz)'); ylabel('捕获时间 (s)');

这段脚本里的捕获判据用了「最后一段的均值」作为锁定参考，比直接用理论频偏值可靠，因为环路锁定点的 VCO 控制电压不为零。批跑时如果模型里有随机噪声源，每次仿真的结果会有抖动，同一组参数跑 10 次取均值更有参考意义。

5.3 波形记录与回放调试

Simulink 的To Workspace模块把中间信号存到工作区后，可以离线做频谱分析。锁定后 VCO 控制电压的功率谱应该没有明显的离散谱线——如果有，说明还有周期性干扰进入环路。一个实用技巧是把error_signal导出来做 FFT，看它的频谱峰值是否在符号速率附近。如果在，说明鉴相器前面的低通滤波器带宽太宽，符号间干扰泄漏到环路里了，把滤波器带宽从 0.5Rs 压到 0.3Rs 再试。这一招对 GMSK 尤其有效，高斯滤波后的信号带宽窄，环路前端滤波器的带宽匹配比 QPSK 更敏感。

本文还有配套的精品资源，点击获取