1.算法描述
随着全球导航卫星系统的高速发展,导航系统的数量也越来越多,比如使用最广泛的GPS导航系统,以及越来越备受关注的中国北斗导航系统等.因此导航频段变得越来越拥挤,且各个频段内的信号相互干扰,在如此情况下,一种二进制偏移载波(Binary Offset Carrier,BOC)调制信号被成功提出,用来解决该问题.而伴随着对BOC信号的深入研究,一系列BOC衍生信号也被提出来应用于各个导航系统中.此外,实际环境中存在着许多干扰因素,给BOC及其衍生信号的捕获带来了困难,因此也成为了国内外学者研究的热点与难点.
BOC 调制就是在原有 BPSK 调制的基础上，再加上一个二进制副载波对 BPSK 信号进行二次扩频。由于方波易于生产，因此采用方波来替代正弦波，这样可以节约硬件资源。一般用 BOC（sf ,cf）来表示sf 代表副载波频率，cf 代表伪码速率。因为sf 和cf都是 1.023MHz 的整数倍，所以文献中用 BOC(m,n)的表示形式，其中 m 表示的是副载波频率，n 表示的是扩频码速率，它们分别表示 1.023MHz 的 m 倍和 n 倍。
BOC 调制的频谱可以分为两个对称独立的边带，它的自相关函数是一种多峰结构，其中主峰较窄，所以有可能对接收机处理方式做一些改变，以便在实现时获得性能的改善。另外一种可行的方法是使鉴相器工作在一个带宽信号的相对函数上，进而保持对主峰的粗跟踪，之后再转入用双边带信号工作的鉴相器，来作精跟踪。由此可见，BOC 调制信号的频谱中上边带和下边带的冗余，自相关函数中的多峰，如果加以开发利用有可能在接收机的信号截获、码跟踪以及数据解调等处理中带来利益。
BOC 调制信号是在载波调制前对扩频信号用副载波进行二次扩频。其流程如下图所示：
正、余弦BOC调制信号的码片波形可以分别表示为
正弦BOC调制信号BOCs(m,n)和余弦BOC调制信号BOCc(m,n)的自相关函数分别为
boc和bpsk区别很多:
1 功率谱不同，boc将基带信号能量搬移到±fs处，bpsk基带信号能量集中在0频。boc的这种特性使得一个频点可以发射多个信号分量。
2 相关函数不同，boc有多个峰，会使得接收机产生误锁，而bpsk不会。
总结，boc是的导航信号调制方式设计更具灵活性，但接收处理更有挑战。bpsk接收处理简单。
2.仿真效果预览
matlab2022a仿真结果如下：
3.MATLAB部分代码预览
%输出即将赋值的可见星号码
PRN = Visible_sv(ii);
fprintf('%02d ', PRN);
%产生随机的多普勒偏移频率
frqshift(ii) = settings.acqSearchBand * 1e3 * rand(1,1) - ...
(settings.acqSearchBand/2) * 1e3;
%此时载波频率
frqBins(ii) = settings.IF + frqshift(ii);
%产生随机导航数据为1或-1
NavData = randi([0,1],1,ceil(navdatLength+1))*2 - 1;
%每个采样点所对应的导航数据号
dataValueIndex = floor((ts * (1:settings.msToProcess*1e-3*settings.samplingFreq)) / ...
(1/50));
%每个采样点所对应的导航数据
NavDataS = NavData(dataValueIndex + 1);
clear dataValueIndex
%利用generateCAcode函数产生每个可见星对应的C/A码
caCode1023 = generateCAcode(PRN);
%每个采样点所对应的C/A码片序号
codeValueIndex = floor((ts * (1:settings.msToProcess*1e-3*settings.samplingFreq)) / ...
(1/settings.codeFreq));
%产生每个采样点所对应的随机相位
codePhase(ii) = floor(1023*rand(1,1));
%每个随机相位下采样点对应的C/A码片数据
codePhase_o=codePhase(ii);
longCaCode = caCode1023((rem(codeValueIndex + ...
settings.CAcodeLength - codePhase(ii), settings.CAcodeLength) + 1));
clear caCode1023
%产生副载波并与伪码相乘得到BOC信号码片数据
subCarr = square(settings.subcarrFreq * phasePoints); %sign(sin(settings.subcarrFreq * phasePoints));
longBocCode = longCaCode.*subCarr;
clear longCaCode subCarr
%产生正弦余弦载波信号
sinCarr = sin(frqBins(ii) * phasePoints);
cosCarr = cos(frqBins(ii) * phasePoints);
%产生正交以及同向的BOC调制信号
I1 = cosCarr .* longBocCode .* NavDataS;
Q1 = sinCarr .* longBocCode .* NavDataS;
clear NavDataS longBocCode sinCarr cosCarr
%多个可见星正交以及同向的BOC调制信号叠加，即模拟实际情况
signalI = signalI + I1;
signalQ = signalQ + Q1;
clear I1 Q1
end
A48
