基于码元序列估计的北斗卫星导航系统b1频带弱d1信号捕获方法
【专利摘要】本发明公开了一种基于码元序列估计的北斗卫星导航系统B1频带弱D1信号捕获方法，首先，将每1ms中频卫星信号与本地信号作相关运算，然后对信号中纽曼霍夫曼(NH)码元序列以及导航信息位进行假设，得到不同的扩展NH码元序列，根据得到的扩展NH码元序列，对每1ms的相关运算结果进行相干和非相干累加，将累加结果的最大值作为检测量和阈值进行比较，最后，判断是否捕获到可视卫星信号。本发明解决了信号中NH码元以及导航信息位跳变引起的相邻伪随机噪声(PRN)码周期间比特符号翻转问题，同时可以直接获得位同步信息，减小卫星定位时间。
【专利说明】
基于码元序列估计的北斗卫星导航系统B1频带弱D1信号捕获 方法
技术领域
[0001] 本发明涉及卫星导航基带信号处理技术领域，具体涉及基于码元序列估计的北斗 卫星导航系统B1频带弱D1信号捕获方法。
【背景技术】
[0002] 北斗卫星导航系统B1频段播发的D1导航电文信号进行了 NH(纽曼霍夫曼）码的二 次编码调制，将数据速率由50bps提高至lkbps，提高了信号抗干扰和抗多径的性能。同时， NH码的调制使得信号中PRN(伪随机噪声）码的周期长度扩大了 20倍，减小了不同卫星的信 号间互相关值，提高了信号的互相关性能。对于高强度卫星信号来说，只需要对一个PRN码 周期长度的信号进行捕获操作处理即可得到所检测的卫星是否可视及其多普勒频移和码 延时等信息。但是对于弱卫星信号来说，往往需要对多个PRN码周期长度的信号进行处理才 能够得到所需的参数。由于北斗卫星导航系统B1频段D1信号中采用的NH码的码片宽度与 PRN码周期均为lms，使得单个PRN码周期长度信号中可能存在比特符号翻转。比特符号翻转 有NH码元跳变和导航电文信息位翻转两种来源。虽然对于单个PRN码周期长度信号中的比 特符号翻转问题可以用DBZP(二倍分组块补零）的方法来解决，但对多个PRN码周期信号相 关运算结果进行积分累加时，由于相邻PRN码周期调制的NH码元以及导航信息数据位可能 是不同的，进行直接累加会导致累加结果变小，从而影响最终捕获结果的判断。
【发明内容】
[0003] 要解决的技术问题是现有技术对多个PRN码周期信号相关运算结果进行积分累加 时，最终捕获结果的判断不准确。
[0004] 为解决上述技术问题，本发明采用的技术方案是：基于码元序列估计的北斗卫星 导航系统B1频带弱D1信号捕获方法，包括以下步骤:步骤1:将每lms中频卫星信号与本地信 号作相关运算，得到每lms的相关运算结果;步骤2:对中频卫星信号中纽曼霍夫曼(NH)码元 序列以及导航信息位进行假设，得到不同的扩展NH码元序列;步骤3:根据步骤2得到的扩展 NH码元序列，对步骤1得到的每lms的相关运算结果进行相干累加和非相干累加;步骤4:将 累加结果的最大值作为检测量和阈值进行比较，判断是否捕获到可视卫星信号。
[0005] 本方法解决了北斗卫星信号中相邻PRN码间高频率比特符号翻转而影响接收机的 捕获结果判断的技术问题。
[0006] 进一步，步骤1的具体步骤如下：1)设定预检测积分时间为Mms;设定多普勒频移搜 索步长为△ f，多普勒频移搜索范围为(fmin，fmax);检测阈值设置为Vt;将多普勒频移搜索值 为$的本地载波与延时了？个采样点的本地PRN码信号相乘，得到lms本地信号l(n);2)取 Mms连续的中频采样信号r(n)，将每lms的中频信号分别与本地信号做相关运算，得到M个相 关结果为价万,？)；步骤2的具体方法如下:根据北斗卫星导航系统B1频带D1信号中采用的NH 码序列，设置t/= 个的扩展NH码元序列ENHj，每个ENHj序列由原始NH码序列与Mms信号 中导航信息位的不同取值假设相结合得到，且均包含(M+19)个码片;假设Mms信号中的起始 NH码元是原始NH码序列中第ini个码元，ini = {1，2，…，20}，对于每一种起始NH码元序号 ini的假设，由J个ENH序列可以得到J个ENHj子序列每个子序列是由相应ENHj序列 中第ini到第ini+M-1个码元组成。
[0007] 进一步，步骤3的具体步骤如下：1)将EW/r序列中的码元分别与步骤1中得到的M 个值对应相乘，得到新的相关结果序列；2)在信号中起始NH码元序号为ini的假设 下，共得到J个新的相关结果序列；3)将得到的每个相关结果序列分成K段，对每段的N个相 关结果值进行相干累加并取模;4)在信号中起始NH码元序号为ini的假设下，对于每一段的 相干累加，由J组相关结果序列可以得到J个相干累加模值，取其中最大值作为本段相干累 加模值结果;5)对得到的K个模值结果进行非相干累加，作为信号中起始NH码元序号为ini 的假设下的非相干累加结果;得到的20个非相干累加结果分别对应于20种信号中起始NH码 元序号的假设;6)取最大值作为估计参数为(U)时的最终检测量。
[0008] 进一步，步骤4的具体步骤如下：1)由步骤3得到了参数估计值为(IJ)时检测量； 2)将检测量与阈值进行比较，检测量大于阈值时，表示捕获成功，跳转步骤3);检测量小于 阈值时，表示未检测到信号，跳转步骤1);3)捕获信号中的多普勒频移和码延时(万，？)，同时 得到了信号中的起始NH码元序号，即导航信息位同步信息;若在多普勒频移和码延时搜索 范围之内没有捕获到卫星信号，则宣布捕获失败。
[0009] 进一步，所述lms本地信号通过将延时7的本地PRN码采样信号与多普勒频移为万 的本地复载波采样信号相乘得到，表示为：
[0010] = -f)e~J2：：(/" "77,，，；; (1)
[0011 ]其中c (n)为本地PRN码采样信号，f IF为中频频率，Ts为采样间隔。
[0012] 进一步，所述中频采样信号的模型为：
[0013] /(") = -sj2P^c{nT^ - r)d{nTs - r)cos(2/T(//; + fd)nTs + (p) + \v{n) (2)
[0014]其中Pr为接收信号功率，c(n)为PRN码采样信号，d(n)为导航信息位的采样，t为接 收信号中的码延时，fd为接收信号中的多普勒频移，f为信号初始相位，w(n)为叠加在信号 上的高斯白噪声；
[0015] 所述每lms的中频采样信号分别与本地信号做相关运算，得到相关结果为：
[0016] /?(./；,.,r) = ^r(/〇/(?) (3) n
[0017] 为了便于表述，将估计参数为时得到的M个相关结果表示为h，其中 i = l，2,…，M〇
[0018] 进一步，将NH码序列表示为NH={c1C2...c2()}，其中 CiG{-l，l}，设置个 的扩展NH码元序列ENHj，每个ENHj序列由原始NH码序列与Mms信号中导航信息位的不同取值 假设相结合得到，每个扩展码元序列中包含(M+19)个码片，表示为：
[0019] ENH： = [c, ... c20 d2c{ ... djC^. ? ?? dg^c^ ... dg^c^ d-^cx ... ^cn]Dd(M_12())] (4)
[0020] 其中cF' = [1右右…心… 为第j种假设下信号中包含的导航信息位，
为Mms信号中导航信息位最大比特数；
[0021 ]所述相干累加并取模的模值为： N
[0022] S；'LJ - =|R^ENHfJ| (5) i-1
[0023] 其中，k=l，2, ? ? ?，K, j = l，2, ? ? ?，J，ini = l，2, ? ? ?，20;ENHfJ表示为已顺化^序 列中第k段子序列的转置且ENK，= [<# ... <# ...:<c=v]7' ;
[0024]所述20个非相干累加结果的计算公式表示为：
[0025] .5 .= max .{S? | /'" = 1，…，20;; = max j 各max (S':"J ;;| (6)--
[0026] 本发明的优点是:本发明采用了先猜后检的方法，通过对信号中NH码元序列和导 航信息位进行假设和估计，扩展了信号积分累加时间，解决了信号中NH码元和导航信息比 特未知导致的捕获性能削弱问题，优点是能够较好的捕获到弱信号北斗卫星，使得接收机 可以在不同环境下正常运行。本发明中通过对每一段相关结果进行相干累加取最大模值的 方法，使累加结果数目不会随着非相干累加运算的进行而增加，节约了存储和运算资源。本 发明在正确检测到卫星多普勒频移和码延时信息的同时，获得了信号中导航信息位同步信 息，可以直接传递到跟踪过程中进行确认，减小接收机定位所需时间。
【附图说明】：
[0027]图1是本发明中基于码元序列估计的弱信号捕获原理图。
【具体实施方式】：
[0028]如图1所示，本发明描述的是一种北斗卫星导航系统B1频带D1导航电文信号捕获 方法，该发明通过对信号中NH码元序列和导航信息位的假设，有效的增加信号累加积分时 间，使接收机可以在弱信号环境下成功捕获到可视卫星。该方法的创新性在于对信号中NH 码元序列以及导航信息位进行假设，得到不同的扩展NH码元序列，根据得到的扩展NH码元 序列，对每lms的相关运算结果进行相干和非相干累加，通过对每一段相关结果进行相干累 加取最大模值的方法，使累加结果的数目不会随着非相干累加运算的进行而增加，节约了 存储和运算资源，同时在正确检测到卫星多普勒频移和码延时信息的同时，获得了信号的 位同步信息，可以直接传递到跟踪过程中进行确认，减小接收机定位所需时间。
[0029]本发明方法具体步骤如下：
[0030]步骤1:设定预检测积分时间为Mms;设定多普勒频移搜索步长为A f，多普勒频移 搜索范围为(fmin，fmax);检测阈值设置为Vt;将多普勒频移搜索值为I的本地载波与延时7 的本地PRN码采样信号相乘，得到lms本地信号l(n);
[0031] 具体过程如下：
[0032] 将延时f的本地PRN码采样信号与多普勒频移为I的本地复载波采样信号相乘， 得到lms本地信号表示为：
[0033] l(n) = c(nT - T)e j2z{J" l)"r- (1)
[0034] 其中c (n)为本地PRN码采样信号，f IF为中频频率，Ts为采样间隔。
[0035] 步骤2:取Mms连续的中频采样信号r(n)，将得到的Mms信号中每lms信号分别与本 地信号做相关运算，得到M个相关结果为
[0036] 具体过程如下；
[0037] 首先得到中频采样信号模型为
[0038] r(n) = yj2Pltc(nT - T)d{nl\ - r)cos(2;r(/,v + fr!)nT^ + (p) + \v(n) ( 2 )
[0039] 其中PR为接收信号功率，c(n)为PRN码采样信号，d(n)为导航信息位的采样，t为接 收信号中的码延时，fd为接收信号中的多普勒频移，为信号初始相位，w(n)为叠加在信号 上的尚斯白噪声。
[0040] 然后将每lms的中频采样信号分别与本地信号做相关运算，得到相关结果为：
[0041] R(f：rT) ^ (3) n
[0042] 为了便于表述，将估计参数为时得到的M个相关结果i?(U)表示为心，其中 i = l，2,…，M〇
[0043] 步骤3:根据北斗卫星导航系统B1频带D1信号中采用的NH码序列，设置j 个 的扩展NH码元序列ENHj，每个ENHj序列由原始NH码序列与Mms信号中导航信息位的不同取值 假设相结合得到，且均包含(M+19)个码片;假设Mms信号中的起始NH码元是原始NH码序列中 第ini个码元，ini = {1，2，…，20}，对于每一种起始NH码元序号ini的假设，由J个ENH序列可 以得到J个ENHj子序列，每个子序列是由相应ENHj序列中第ini到第ini +M-1个码元组 成;将i?皿f序列中的码元分别与步骤2中得到的M个值对应相乘，得到新的相关结 果序列;在信号中起始NH码元序号为ini的假设下，共得到J个新的相关结果序列。将得到的 每个相关结果序列分成K段，对每段的N个相关结果值进行相干累加并取模;在信号中起始 NH码元序号为ini的假设下，对于每一段的相干累加，由J组相关结果序列可以得到J个相干 累加模值，取其中最大值作为本段相干累加模值结果，然后对得到的K个模值结果进行非相 干累加，作为信号中起始NH码元序号为ini的假设下的非相干累加结果。最终得到的20个非 相干累加结果分别对应于20种信号中起始NH码元序号的假设，取最大值作为估计参数为 (万,？)时的最终检测量。
[0044] 具体过程如下；
[0045] 尽管Mms信号中包含的NH码元是由固定的原始NH码序列构成，导航信息位的未知 特性使得信号中NH码元的取值存在不确定性，由于最终的累加量是以绝对值形式进行检测 的，因而信号中导航信息位的取值共有^/刮种可能性。将NH码序列表示为NH = {C1C2. ? .C2〇}，其中CiE {-I，1}，考虑到彳目号中导航彳目息位的不同取值，可以得到j = 种 扩展NH码元序列EN出，每个扩展码元序列中包含(M+19)个码片，表示为：
[0046] 五身丑-[('1 ……" ? C'p-Cj …fip:Crn0.卿-⑶.)] (4)
[0047] 其中心…为第j种假设下信号中包含的导航信息位
为Mms信号中导航信息位最大比特数，如上所述，由于最终的累加量是以绝对值形式进行检 测的，将信号中第一个导航信息位比特设为1，则对信号中导航信息位的取值只需要考虑 2H种可能组合。在本专利中侧重于考虑在增加信号积分累加时间时信号中相邻PRN码周期 间的比特翻转问题，假设了步骤1中计算相关结果的每lms信号中有一个完整的PRN码周期， 因而Mms信号中包含了M个NH码元。
[0048] 假设Mms信号中的起始NH码元为原始NH码序列中第ini个码元且ini = {1，2，…， 20}，则在每一种起始NH码元序号假设下由J个EN氏序列可以得到J个子序列￡胃「+，每个子 序列是由相应EN氏序列中第ini到第ini+M-1个码元组成;将^^/广序列中的码元分别与步 骤2中得到的M个及(1,9值对应相乘，得到新的相关结果序列。由此可知，对于每种信号中 起始NH码元序号为ini的假设，可以得到J个新的相关结果序列。
[0049] 由于对相关结果进行累加时有相干累加和非相干累加两种方式，假设将得到的每 个相关结果序列分成K段，对每段的N个相关结果进行相干累加，然后对相干结果进行非相 干累加，可以得到第k个相干累加结果的模值为：
[0050] S^J = =|R^.ENHfJ| (5)
[0051 ]其中，k=l，2, ? ? ?，K, j = l，2, ? ? ?，J，ini = l，2, ? ? ? JCLENHfW表示为已顺…^序 列中第k段子序列的转置且…。
[0052] 在信号中起始NH码元序号为ini的假设下，对于同一序号段相关结果的相干累加， 由J组相关结果序列可以得到J个相干累加模值，取其中最大值作为本序号段相干累加模值 结果，然后对得到的K个模值结果进行非相干累加，作为信号中起始NH码元序号为ini的假 设下的非相干累加结果。由上述可知在估计参数为(万时，可以得到20个非相干累加结 果，分别对应于Mms信号中20种起始NH码元的假设，取最大值作为估计参数为(I，？)时的最 终检测量。上述操作可以表示为：
[0053] S - maxjs? \ini - 1,",,2〇|= max }| (6)
[0054] 步骤4:由步骤3得到了参数估计值为(万,？)时检测量，将检测量与阈值进行比较， 若检测量大于阈值，则表示捕获成功，得到信号中的多普勒频移和码延时(U)，同时得到 了信号中的起始NH码元序号，即导航信息位同步信息;如果检测量小于阈值，表示参数估计 值为(U)时，没有检测到信号，改变参数估计值(万，？)，重复步骤1~3进行搜索检测直至捕 获成功;若在多普勒频移和码延时搜索范围之内没有捕获到卫星信号，则宣布捕获失败。
[0055] 以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应 视为本发明的保护范围。
【主权项】
1. 基于码元序列估计的北斗卫星导航系统BI频带弱Dl信号捕获方法，其特征是，包括 W下步骤： 步骤1:将每Ims中频卫星信号与本地信号作相关运算，得到每Ims的相关运算结果； 步骤2:对中频卫星信号中纽曼霍夫曼(NH)码元序列W及导航信息位进行假设，得到不 同的扩展NH码元序列； 步骤3:根据步骤2得到的扩展NH码元序列，对步骤1得到的每Ims的相关运算结果进行 相干累加和非相干累加； 步骤4:将累加结果的最大值作为检测量和阔值进行比较，判断是否捕获到可视卫星信 号。2. 根据权利要求1所述基于码元序列估计的北斗卫星导航系统Bl频带弱Dl信号捕获方 法，其特征是，步骤1的具体步骤如下： 1) 设定预检测积分时间为Mms;设定多普勒频移捜索步长为Af,多普勒频移捜索范围 为(f"in，f"ax);检测阔值设置为Vt;将多普勒频移捜索值为方的本地载波与延时了？个采样 点的本地PRN码信号相乘，得到Ims本地信号Kn); 2) 取Mms连续的中频采样信号Hn),将每Ims的中频信号分别与本地信号做相关运算， 得到M个相关结果为/?(7六；)； 步骤2的具体方法如下： 根据北斗卫星导航系统Bl频带Dl信号中采用的NH码序列，设置苦1个的扩展NH码 元序列ENHj,每个ENHj序列由原始NH码序列与Mms信号中导航信息位的不同取值假设相结合 得到，且均包含(M+19)个码片;假设Mms信号中的起始NH码元是原始NH码序列中第ini个码 元，ini = {1，20}，对于每一种起始NH码元序号ini的假设，由J个ENH序列可W得到J个 ENH牙序列心VHf，每个子序列是由相应ENHj序列中第in巧Ij第ini+M-1个码元组成。3. 根据权利要求2所述基于码元序列估计的北斗卫星导航系统Bl频带弱Dl信号捕获方 法，其特征是，步骤3的具体步骤如下： 1) 将序列中的码元分别与步骤1中得到的1个巧7^，訂值对应相乘，得到新的相 关结果序列； 2) 在信号中起始N曲马元序号为ini的假设下，共得到J个新的相关结果序列； 3) 将得到的每个相关结果序列分成K段，对每段的N个相关结果值进行相干累加并取 模； 4) 在信号中起始NH码元序号为ini的假设下，对于每一段的相干累加，由J组相关结果 序列可W得到J个相干累加模值，取其中最大值作为本段相干累加模值结果； 5) 对得到的K个模值结果进行非相干累加，作为信号中起始NH码元序号为ini的假设下 的非相干累加结果;得到的20个非相干累加结果分别对应于20种信号中起始NH码元序号的 假设； 6) 取最大值作为估计参数为(方，?)时的最终检测量。4. 根据权利要求3所述基于码元序列估计的北斗卫星导航系统Bl频带弱Dl信号捕获方 法，其特征是，步骤4的具体步骤如下： 1) 由步骤3得到了参数估计值为(万，r肘检测量； 2) 将检测量与阔值进行比较，检测量大于阔值时，表示捕获成功，跳转步骤3);检测量 小于阔值时，表示未检测到信号，跳转步骤1); 3) 捕获信号中的多普勒频移和码延时(}；；，〇,同时得到了信号中的起始NH码元序号，即 导航信息位同步信息；若在多普勒频移和码延时捜索范围之内没有捕获到卫星信号，则宣 布捕获失败。5. 根据权利要求1或2所述的基于码元序列估计的北斗卫星导航系统Bl频带弱Dl信号 捕获方法，其特征是:所述Ims本地信号通过将延时？的本地PR师马采样信号与多普勒频移为 充的本地复载波采样信号相乘得到，表示为：(1) 其中C (n)为本地PRN码采样信号，f IF为中频频率，Ts为采样间隔。6. 根据权利要求1或2所述的基于码元序列估计的北斗卫星导航系统Bl频带弱Dl信号捕获方法-甘赔县.日K未*搞?巧拴倍皂? 首刑* ? (2) 其中扣为接收信号功率，c(n)为PR师马采样信号，d(n)为导航信息位的采样，T为接收信 号中的码延时，fd为接收信号中的多普勒频移，與为信号初始相位，w(n)为叠加在信号上的 局斯白噪声；所述每Ims的中师吉皂 A功I 与号做相关运算，得到相关结果为： (3) 为了便于表述，将估计参数为(万，。时得到的M个相关结果/?(方，。表示为Ri，其中i = 1， 2,... ,M。7. 根据权利要求3所述的基于码元序列估计的北斗卫星导航系统Bl频带弱Dl信号捕获 方法，其特征是:将畑码序列表示为NH={：cic2...C2〇}，其中(：1￡{-1，^，设置^ = 2^等1个的 扩展NH码元序列ENHj，每个ENHj序列由原始NH码序列与Mms信号中导航信息位的不同取值假 设相结合得到，每个扩展码元序列中包含(M+19)个码片，表示为：其中4' =[:^;'與…^^'…^^^为第巧中假设下信号中包含的导航信息位，为Mms信号中导航信息位最大比特数； 寅巾nit而照的照化击.樹 其中，k=l，2，...，KJ = I，2，...，J，ini = l，2，...，20;ENHr'表示为 ENHini'j 序列中 第k段子序列的转置且ENKfj' =战端…式说… <(鮮f ; 所述20个非相干累加结果的计算公式表示为：
【文档编号】G01S19/30GK105911569SQ201610208635
【公开日】2016年8月31日
【申请日】2016年4月6日
【发明人】樊祥宁
【申请人】东南大学