BOC信号码跟踪抖动仿真：理论模型与Matlab实现

简介：本资源面向卫星导航信号处理方向的研究生、工程师及MATLAB仿真初学者，聚焦BOC调制解调系统中码跟踪性能的关键指标——码跟踪抖动标准差，系统仿真其随载噪比（C/N₀）变化的规律，并对比Unambiguous与KFPELTA两类典型跟踪环路的抗噪能力差异。压缩包共21个文件（2.3MB），含14个核心MATLAB脚本（实现主跟踪、对比分析与函数封装）、4个预存数据.mat文件（支撑不同BOC阶数如BOC10_5/BOC14_2的仿真复用）、2张结果图（直观展示抖动曲线）及1段AVI操作录像（完整演示MATLAB环境配置、路径设置与一键运行流程）。已有253人学习下载，配套录像特别强调当前文件夹路径设置要点，有效规避常见运行报错；所有代码模块清晰、注释完备，支持快速复现、参数调整与算法对比验证。

1. BOC信号为什么不能直接套用BPSK的抖动结论

1.1 BOC调制到底改了什么

做接收机基带算法这一行，绕不开BOC调制。BOC的全称是Binary Offset Carrier，二进制偏移载波，它和传统BPSK-R最直观的区别，就是在PRN码上额外乘了一个方波子载波。用BOC(m,n)表示时，m和n分别是相对于1.023 MHz的倍数，比如BOC(1,1)代表1.023 MHz子载波和1.023 Mcps码率，BOC(14,2)则对应14×1.023 MHz子载波和2×1.023 Mcps码率。这个额外的子载波把本来集中在主瓣的能量推到两侧，形成典型的双峰谱，所以BOC又被称为分裂谱调制。

这样做的直接收益有两个：一是频谱和原有BPSK信号的重叠程度降低，能够改善不同系统、不同信号之间的互干扰；二是等效带宽变大，信号理论上具备更好的码相位鉴别能力。但代价也相当明显，BOC信号的相关函数不再像BPSK那样是一个干净的单峰，而是出现多个幅度相近的波峰，整个相关曲线的包络变窄，但主峰周围有大量副峰。这会让捕获和跟踪阶段面临副峰锁定的风险，也让码环行为变得更加微妙。

1.2 从S曲线到鉴别器增益：抖动仿真的底层逻辑

码跟踪环本质上是一个负反馈系统。本地生成的BOC码与接收信号做相关，超前和滞后相关结果经过鉴别器生成误差信号，再由环路滤波器驱动码NCO，让本地码相位持续逼近输入码相位。系统跟踪精度的高低，并不完全取决于瞬时相关峰高度，而是取决于鉴别器过零点附近的增益和环路带宽内噪声的共同作用。

问题在于，鉴别器的增益来自S曲线，而S曲线的形状由相关函数、相关间距、前端滤波带宽共同决定。BOC信号的相关函数有过零点密集、多峰共存的特点，它的S曲线在零附近虽然通常比BPSK更陡，但也更容易受到相关间距和带宽变化的影响。同一个环路带宽、同一个载噪比，BPSK和BOC的码跟踪抖动标准差会有明显差别；即使同样是BOC，BOC(1,1)和BOC(6,1)结果也可能相差很大。因此工程上不能拿BPSK的结论直接外推，必须对具体调制方式做完整的仿真统计。

2. 仿真链路搭建：把码环跟踪的每个环节拆开来看

2.1 仿真链路总览

整个仿真链路可以拆成三大块：信号源模块、接收端跟踪环路模块、数据统计模块。信号源先生成PRN码，完成BOC调制，再加入高斯白噪声；接收端完成载波剥离、相关积分、码鉴别、环路滤波和码NCO更新；数据统计模块在环路进入稳态之后记录码相位误差序列，并计算标准差。

为了防止载波相位对码环统计造成干扰，我采用零中频复基带模型，本地载波用给定频率的复指数直接解调。实际工程中载波环和码环会存在耦合，但在只评估码跟踪抖动的前提下，这种简化是常见做法，也能让代码结构更清晰。参数上选择BOC(1,1)信号，码率1.023 Mcps，子载波1.023 MHz，采样率20.46 MHz，PRN码长1023，相关间距0.1 chip，相干积分时间1 ms，码环带宽1 Hz。

这里尤其要说明采样率的选择。BOC(1,1)的子载波是1.023 MHz方波，如果采样率只取2.046 MHz，码片内只有两个采样点，子载波的正负半周刚好能表示，但相关峰形状会受到严重影响，S曲线畸变，抖动统计结果不可信。取20倍码率以上是比较稳妥的做法，既能刻画方波子载波的跳变沿，又不会让单次仿真数据量失控。

2.2 载噪比控制：最容易出错的一步

仿真里要扫的变量是载噪比，单位是dB-Hz。载噪比和信噪比最根本的区别在于，信噪比的参考带宽是采样率或前端带宽，而载噪比的参考带宽是1 Hz。也就是说C/N0 = P_signal / N0，其中P_signal是信号功率，N0是噪声功率谱密度，单位是W/Hz，与采样带宽无关。

在做离散复基带仿真时，信号样本功率为Ps，噪声样本方差为Pn，噪声功率谱密度N0 = Pn / fs，因此C/N0的线性值等于Ps * fs / Pn。反过来，给定目标载噪比C/N0_lin，噪声方差应设为：

Pn = Ps * fs / C/N0_lin

举个例子，目标C/N0是35 dB-Hz，线性值就是10^3.5 ≈ 3162。如果采样率是20.46 MHz，信号功率归一化为1，那么噪声方差应该是20.46e6 / 3162 ≈ 6470。这个数值比信号功率大几个数量级，很多人在这一步直接把噪声方差设为1，结果实际载噪比完全不对，扫出来的曲线整体偏移好几个dB。这个坑我后面还会再提，因为它的错误表现并不直观，曲线形状看起来依然正常，但数值完全不可信。

2.3 环路的离散化实现要点

码环的离散化实现主要涉及码NCO和鉴别器。码NCO负责产生本地码相位，每个采样点累加码相位增量，当相位超过一个码片时，就切换PRN码序号并完成BOC子载波符号翻转。相关器分为早、即时、迟三路，早期和迟期相对即时码相差±d/2 chip。

鉴别器推荐用非相干超前减滞后功率鉴别器：

D = (IE^2 + QE^2 - IL^2 - QL^2) / (IE^2 + QE^2 + IL^2 + QL^2)

这个鉴别器在BOC信号上稳定性较好，不依赖载波相位的绝对跟踪状态，也不会因为信号幅度波动而出现增益漂移。分母做归一化之后，鉴别器输出的量纲仍然是码相位误差，环路滤波器可以直接复用。为了防止副峰锁定，每次试验的初始码相位误差都设置在±0.3 chip以内，确保环路从相关主峰附近开始收敛。这个细节在正式项目中非常关键，因为BOC信号副峰吸引力很大，初始误差超过半码片，环路很可能锁到副峰上，抖动统计值会变成错误结果。

3. 码跟踪抖动标准差的理论模型：公式、数值积分和对照逻辑

3.1 热噪声抖动公式的来源与每一项的含义

为了验证蒙特卡洛统计结果是否合理，仿真里通常叠加一条理论曲线。码环热噪声抖动标准差可以写成：

sigma_DLL = sqrt( Bn * (1 - 0.5*Bn*Ti) * ∫ Gs(f) * sin^2(π*f*d*Tc) df / ( (2π)^2 * (C/N0) * [ ∫ f*Gs(f)*sin(π*f*d*Tc) df ]^2 ) )

积分范围是[-beta_r/2, beta_r/2]，beta_r是接收机前端双端带宽。公式里每一项的含义都不容忽视：Bn是码环单边等效噪声带宽，单位为Hz；Ti是相干积分时间，单位为s；d是早迟间距，单位是chip；Tc是码片持续时间，单位为s；Gs(f)是归一化信号功率谱密度；C/N0是载噪比线性值。

这个公式的物理含义可以这样理解：分子里的sin²项表示早迟相关器对噪声的采样作用，它和信号功率谱密度共同决定有多少热噪声进入码环；分母里那个带sin的单重积分，本质上是环路S曲线在零附近的斜率。它越大，说明码相位误差转换成鉴别器输出的增益越高，同样的噪声对应到码相位域就越小。两者一比，得到的就是码相位误差的方差。

3.2 用Matlab数值积分算理论曲线

理论上Gs(f)可以用BOC方波功率谱的解析表达式直接算，但工程上更推荐直接用仿真的BOC码序列做周期图估计，再对频谱做归一化。这样做的好处是，如果信号模型里加了滤波、限带或非理想子载波，理论曲线能自动跟随实际信号模型，不容易出现“公式是公式、仿真是仿真”的两张皮。

归一化规则是sum(Gs)*df = 1，也就是离散功率谱所有频率点之和乘以频率分辨率等于1。之后用trapz做数值积分，注意积分范围要和前端带宽保持一致。我工程里常用的理论曲线计算函数大致是：

function sigmaDll = dllJitterTheoryBn(bn, cn0Lin, d, br, psdFreq, psdVal, tc, ti) % psdFreq: 频率轴，单位Hz % psdVal: 归一化功率谱密度，单位1/Hz % tc: 码片持续时间，单位s % ti: 相干积分时间，单位s df = psdFreq(2) - psdFreq(1); idx = abs(psdFreq) <= br/2; f = psdFreq(idx); G = psdVal(idx); num = trapz(f, G .* (sin(pi*f*d*tc)).^2); den = trapz(f, f .* G .* sin(pi*f*d*tc)); sigmaDll = sqrt(bn * (1 - 0.5*bn*ti) * num / ((2*pi)^2 * cn0Lin * den^2)); end

要特别提醒的是，psdFreq必须包含负频率。如果只取正频率，分子和分母都会少了一半信号能量，理论值会明显偏大。这是我实际调试中踩过的问题，当时理论曲线始终比仿真高一倍，查了半天才发现积分区间少了一半。

3.3 为什么理论曲线和蒙特卡洛统计不能要求完全重合

理论公式建立在若干理想假设上：噪声是高斯白噪声、环路处于线性工作区、鉴别器S曲线可以被一阶泰勒展开近似、早迟相关器输出完全对称。实际蒙特卡洛仿真中，环路滤波器的离散化、积分时间与环路带宽的比值、鉴别器归一化方式、初始相位误差残留，都会带来额外偏差。因此，理论曲线和仿真曲线的合理关系是趋势一致、量级接近，而不是强行要求逐点重合。

通常我会把理论曲线和仿真曲线画在同一张图上，如果偏差在30%以内，就认为仿真链路基本正确；如果偏差超过一倍，优先检查噪声功率换算、PSD归一化和鉴别器增益。这个判断标准是我实际项目中的经验值，不一定适用于所有场景，但至少能帮你快速定位问题。

4. Matlab工程实现：主程序、关键函数与参数配置

4.1 工程文件组织

为了让代码能复现，我会把工程整理成五个文件，职责划分清楚。initParams.m集中定义所有系统参数；generateBOC.m生成PRN码和BOC基带波形；runTracking.m执行一次给定载噪比下的DLL闭环跟踪并返回码相位误差序列；computeDllJitter.m计算理论曲线；plotCNoCurve.m负责扫描载噪比并画图。这样每个文件都能单独调试，也方便换参数验证。

我自己的习惯是先用脚本形式跑通单点，再封装成函数。因为码环仿真涉及的变量非常多，一开始全部塞进一个脚本会很难排查问题。先固定一个载噪比，比如35 dB-Hz，跑一次跟踪，画一下S曲线和码相位误差收敛曲线，确认环路稳定了，再写载噪比扫描循环。这个顺序能省下大量调代码的时间。

4.2 BOC信号生成与PSD估计

BOC波形生成的代码不复杂，关键是码NCO和子载波的相位要严格对齐。一种常见的实现是先生成PRN码序列，再按采样点索引计算对应的码片序号和子载波符号：

fs = 20.46e6; fc = 1.023e6; fsc = 1.023e6; codeLen = 1023; numSamples = round(fs / 1000); % 1ms信号 prn = generatePrn(length(codeLen), 1); t = (0:numSamples-1) / fs; codeIdx = mod(floor(t * fc), codeLen) + 1; codeSeq = prn(codeIdx); subcarrier = sign(sin(2*pi*fsc*t)); bocSignal = codeSeq .* subcarrier;

PSD估计可以用periodogram或者pwelch。因为理论计算需要的是信号本身的功率谱密度而不是加窗后的谱分辨率，我通常直接对整段BOC序列做FFT，然后归一化。注意BOC是方波子载波，如果直接用periodogram默认的汉明窗，谱峰会变宽，理论曲线和仿真曲线会系统性偏离，所以我在这个场景里不用加窗。

4.3 跟踪环路主循环

码环主循环是整个工程的核心，逻辑上分为相关运算、鉴别器、环路滤波和NCO更新四步。下面给出一阶码环滤波器的简化版本，目的是把热噪声抖动的统计逻辑讲清楚：

function codeErrorSeq = runTracking(bocSignal, cn0Lin, params) fs = params.fs; fc = params.fc; fsc = params.fsc; ti = params.ti; d = params.d; bn = params.bn; numSteps = length(bocSignal) / round(fs*ti); codeNco = 0; codePhase = 0; codeErrorSeq = zeros(numSteps, 1); loopGain = bn * 4 * ti; % 一阶码环的简化增益 for k = 1:numSteps % 取一段相干积分数据 idx = (k-1)*round(fs*ti)+1 : k*round(fs*ti); blk = bocSignal(idx); % 生成本地早、即时、迟码 [early, prompt, late] = localCode(codePhase-d/2, codePhase, codePhase+d/2, params); % 相关积分 ie = real(blk .* conj(early)); % 简化为复相关 ... % 鉴别器 disc = (ie^2 + qe^2 - il^2 - ql^2) / (ie^2 + qe^2 + il^2 + ql^2); % 环路滤波与NCO更新 codeNco = codeNco + loopGain * disc; codePhase = codePhase + codeNco * round(fs*ti) / fc; codeErrorSeq(k) = codePhase; % 记录相对初始相位 end end

这段代码省略了本地码生成和正交下变频的具体写法，工程里需要补充载波NCO和正交分量累加。完整实现时，伪码同步部分要保证本地码率和接收信号码率一致，否则码相位会线性漂移，抖动统计值里混入多普勒频移导致的趋势项，曲线在低载噪比段会明显抬高。

4.4 参数设置与循环逻辑

载噪比扫描范围我通常取20到50 dB-Hz，步进2 dB。每个载噪比点做20次蒙特卡洛试验，每次试验跑1000个相干积分周期，也就是1秒仿真数据，然后丢弃前面100个周期，取稳态后900个周期的码相位误差序列计算标准差。这样可以消除初始收敛段对统计结果的影响。

参数集中定义的好处是，换BOC阶数、换相关间距、换环路带宽都只需要改一行。下面是我常用的参数表。

5. 仿真结果解读：曲线形态、参数敏感性和工程判断

5.1 抖动标准差随载噪比变化的总体趋势

拿BOC(1,1)信号跑完扫描，得到的曲线和我预期一致：载噪比越低，码跟踪抖动标准差越大，而且大致遵循平方根反比关系。也就是说，载噪比下降10 dB，抖动标准差约增大为原来的3.16倍。这是因为热噪声抖动主项和C/N0是平方根关系，在中等和高载噪比区间，曲线在双对数坐标下近似一条斜率为-0.5的直线。

下面是一组示例结果，对应BOC(1,1)、相关间距0.1 chip、码环带宽1 Hz、前端带宽20 MHz。数值是我在类似配置下得到的量级参考，不同实现可能略有出入。

从表格看，蒙特卡洛统计值和理论值很接近，差值基本在10%以内。在25 dB-Hz这种较低载噪比段，仿真值比理论值略高，是因为环路逐渐接近非线性工作区，鉴别器增益下降，理论公式的线性化假设不再完全成立。这是正常现象，不需要过度担心。

5.2 环路带宽、相关间距和BOC阶数的影响

扫完主曲线之后，我通常还会做几组参数敏感性测试。第一组是改码环带宽，把Bn从1 Hz改到2 Hz和0.5 Hz。结果是码环带宽越大，抖动标准差越大，因为环路对噪声的等效积分带宽变宽了。但实际系统里不能一味减小带宽，因为带宽太小会导致环路拉不住动态应力，高动态场景下会出现常值偏置甚至失锁。所以码跟踪抖动标准差曲线和动态应力曲线是一对矛盾，最终选择是折中结果。

第二组是改相关间距d。d从0.1 chip改成0.5 chip之后，热噪声抖动明显增大，这是因为增大早迟间距会降低S曲线在零附近的斜率，也就是鉴别器增益下降。但如果d取得太小，信号经过前端限带后相关峰变圆，S曲线在零附近反而可能变平坦。过小的相关间距对多径误差也未必有利。通常在BOC(1,1)这类信号上，d取0.1到0.2 chip是比较常见的范围。

第三组是换BOC阶数，比如BOC(6,1)。由于子载波频率更高，信号的Gabor带宽更大，理论上相同载噪比下热噪声抖动会更低。但随之而来的是前端带宽要求更高，采样率也要相应提高，副峰数量增多，跟踪环路实现难度加大。这种曲线对比可以直接支撑信号体制选型的讨论，做星座设计和频率兼容分析时特别有用。

5.3 结果报表与工程用途

仿真完成后，我习惯把结果整理成一张双对数坐标图，横轴是载噪比C/N0，纵轴是码跟踪抖动标准差，单位用chips或换算成米。BOC(1,1)的1 chip对应约293米，换算后能直观给出测距误差水平。

这张曲线图在工程上的用途非常明确：第一，接收机链路预算时，给定前端载噪比指标，可以查出码环能达到的稳态跟踪精度；第二，环路带宽设计时，可以根据载噪比和动态指标折中选带宽；第三，做抗干扰分析时，可以评估干扰导致载噪比恶化后，码环性能退化了多少。仿真脚本一旦搭好，这些分析都是几分钟内改参数重跑的事。

6. 配套操作录像的复现过程与避坑清单

6.1 录像里操作的完整步骤

这份配套仿真操作录像的核心目的，是让你能照着屏幕操作把整条曲线跑出来。录像从打开Matlab开始，第一步是运行initParams.m初始化参数，第二步运行generateBOC.m生成BOC(1,1)基带信号，并短暂显示时域波形和频谱，确认子载波和频谱双峰没有异常。第三步是单点测试，先把载噪比固定在35 dB-Hz，运行runTracking.m观察码相位误差收敛曲线，确认环路稳定。第四步运行plotCNoCurve.m，扫描20到50 dB-Hz的全部载噪比点，等待进度条走完，最终生成结果图和表格。

录像里特意展示了如何修改相关间距和环路带宽后重新出图，因为这是使用这套仿真工具最频繁的操作。如果你手上只有标题里提到的这套仿真，没有完整工程，按上面的步骤从零把脚本组织起来，也能得到同样的结构。

6.2 我踩过的几个坑

先说噪声功率换算。前面提到过，复基带仿真里如果直接把噪声方差设成1，而信号功率也做了归一化，实际载噪比会被采样率严重拉低，导致曲线整体右移，非常隐蔽。我的解决方法是把载噪比换算写成独立函数，每次设置参数时打一次日志，确认理论C/N0和实际产生C/N0的误差在0.1 dB以内。

第二个坑是统计时长不够。早期我每点只跑200 ms，蒙特卡洛抖动统计值很不稳定，曲线毛刺很多。后来改为每点跑1秒、重复20次后，曲线才变得平滑。原因很简单，码环本身是低通系统，等效噪声带宽越窄，需要的统计时间越长。如果环路带宽只有1 Hz，相关时间常数为0.5秒，只统计

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