专利名称：基于码长特征的频谱分离系数检测方法
技术领域：
本发明涉及的是一种卫星导航技术领域的检测方法，用于全球导航卫星系统 (Global Navigation Satellite Systems, GNSS)的基于码长特征的频谱分离系数检测方法。
背景技术：
在全世界的卫星导航系统频谱资源有限的背景下，由于有众多卫星导航系统的存 在，使得简单要求各卫星导航系统的授权信号频谱完全分离是不可能的，因此，部分的频谱 重叠是不可避免和必需的；另一方面，在ITU框架下，已经证明并且达成国际共识的是，在 频谱重叠的条件下，现有卫星导航信号可以互存而不影响彼此的正常使用。因此，对于现有 的卫星导航信号和设计中的卫星导航信号，评估该信号和其它卫星导航信号的兼容性(兼 容性是指两个或多个系统同时工作时，不会引起冲突，相对于单系统工作的情况不会产生 显著的性能下降)是一个重要的研究课题。经对现有技术文献的检索发现，F. SoualIe，Th. B. Burger.撰文Impact of Galileo Spreading Code Selection and Data Rate onto Navigation Signal Interference [C] //ION GPS/GNSS 2003. Portland JON, 2003 1035-1043. “伽利略扩 频码的选择和数据速率对导航信号干扰的影响”，文中提出以频谱分离系数(Spectral Separation Coefficient (SSC))来做为全球导航卫星系统兼容性评估的主要参数之一，并 指出扩频码和数据速率对频谱分离系数的影响，但文章给出的基于频谱分离系数定义的检 测方法过程复杂，且没有定量的衡量分类标准，不具有明确的可操作性，因此寻找有效的频 谱分离系数检测方法对于全球导航卫星系统兼容性评估问题具有重要的意义。在进一步的检索中，尚未发现任何与本文类似的卫星导航信号基于码长特征的频 谱分离系数检测方法的文献和报道。

发明内容
本发明的目的在于克服现有技术中的不足，提供一种基于码长特征的频谱分离系 数检测方法。本发明解决了现有卫星导航信号频谱分离系数的检测过程复杂且检测效率不 高的技术问题，在不引入可见的计算误差的前提下，大大提高频谱分离系数的检测效率，使 得全球导航卫星系统的兼容性评估有快速、准确的实现途径。本发明是通过以下技术方案实现的本发明包括步骤如下步骤一、系统启动，在检测设备上设定需要检测的目标频段，得到位于该频段上的 所有卫星信号的集合M= {xi(t),x2(t),x3(t),…xn (t)}，并且得到初始集合的大小N;步骤二、根据目标频段的系统信息，获得GPS(Global Positioning System全球定 位系统)和Galileo (伽利略)系统的相应频段的卫星导航信号的信息列表；步骤三、将集合M中属于步骤二中所得信息表中的卫星信号的成员移入集合Mk，Mk的集合大小为Nk = k ；集合M的大小变为N = N-k,此时集合M为待检测集合，Mk为已知信 号集合；步骤四、将已知信号集合Mk中的元素导航信号xk(t) (k = 1，. . . Nk)的码长特征参 数《 k初始化为，其中Lk为第k个信号伪随机码码长，fkc为第k个信号码片速率，
Jkc
fkd为第k个信号导航数据速率；将α k < 1的信号放入短码集合A，假设集合元素为n，将 1的信号放入长码集合B，假设集合元素为m ；步骤五、对待检测集合M中的元素Xi (t)(i = 1，...，N)根据码长特征参数、逐 一处理，计算频谱分离系数κ i
Lf.,1)将待检测导航信号Xi (t)的码长特征参数a i初始化为f
Jic2)根据待检测导航信号的码长特征参数a i，计算导航信号\(0与集合A中每个 元素的频谱分离系数kAi (i = 1，. . .，η)；计算导航信号Xi (t)与集合B中每个元素的频谱 分离系数 kBi(i = 1，. . .，m)；3)取频谱分离系数kAi(i = 1，...，η)和kBi(i = 1，...，m)中的最大值，做为待 检测导航信号Xi (t)的频谱分离系数κ i ；步骤六、设置系统最终频谱分离系数阈值Lim，默认阈值为-50dB/Hz，如果Ici > Lim,则检测结果为导航信号Xi (t)不能达到频谱分离系数检测标准，如果Ici ( Lim,则检测 结果为导航信号Xi (t)频谱分离系数检测达标，符合卫星导航信号兼容性要求，基于码长特 征的频谱分离系数检测完成。所述的频谱分离系数κ i，待检测导航信号的码长特征参数a i < 1并且干扰信号 属于集合A时，、为
/2K = J Gs (f公式(1)
-A/2其中，I为接收机前端带宽，Gs(f)为期望信号的功率谱密度，Gi (f)为干扰信号 的功率谱密度，Δ f为多普勒频移。所述的频谱分离系数κ i，待检测导航信号的码长特征参数a i > 1时，无论干扰 信号属于集合A或者B，Ki为:
A/2K= f Gj/XU/H/公式(2)
-A/2其中，β,为接收机前端带宽，Gs。(f)为期望信号的码片功率谱密度，Gic;(f)为干扰 信号的码片功率谱密度。计算频谱分离系数时，待检测导航信号为期望信号，已知信号集合A、B中的信号 均为干扰信号。所述的功率谱密度是指用密度的概念表示信号功率在各频率点的分布情况，对 功率谱在频域上积分就得到信号的功率，使用周期性伪随机码且包含数据调制的卫星导航 信号的功率谱密度Gx (f)
其中L为信号伪随机码码长，fc为信号码片速率，fd为导航数据速率，DFT(Ck)为 伪随机码的离散傅立叶变换，Φ为调制波形的傅立叶变换。本发明首先在检测设备上设定需要检测的目标频段，得到位于该频段上的所有卫 星信号，以码长特征对卫星信号进行分类，将待检测卫星信号分别与已分类的已知卫星信 号进行频谱分离系数的检测，最终得到频谱分离系数的检测结果。本发明对卫星导航信号进行基于码长特征的频谱分离系数的检测，以便快速判定 不同卫星导航信号间的兼容性。故首先在检测设备上设定需要检测的目标频段，得到位于 该频段上的所有卫星信号，以码长特征对卫星信号进行分类，将待检测卫星信号分别与已 分类的已知卫星信号进行频谱分离系数的检测，最终得到频谱分离系数的检测结果。该检 测方法以码长特征参数代入频谱分离系数检测，且实施的每一步骤都有定量的衡量分类标 准，在不引入可见的计算误差的前提下，大大提高频谱分离系数的检测效率，使得全球导航 卫星系统的兼容性评估有快速、准确的实现途径。


图1实施例中待检测导航信号的功率谱密度曲线示意图；图2实施例中待检测信号的码片功率谱密度曲线示意图；图3实施例中待检测信号与GPS CA信号间的频谱分离系数曲线示意图。
具体实施例方式下面结合附图对本发明的实施例作详细说明本实施例在以本发明方案为前提下 进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的 实施例。实施例设计待检测卫星信号X1 (t)工作在Ll频段上，X1 (t)信号的伪随机码码长L为 2046，信号码片速率f。为1. 023MHz、导航数据速率fd为50Hz。步骤一系统启动，在检测设备上设定需要检测的目标频段为Li，得到位于该频 段上的所有卫星信号的集合 M = {Xl (t)，GPS CA, GPS LlP (Y)，GPS LlC, Galieo El OS}，并 且得到初始集合的大小N = 5 ；步骤二 根据目标频段的系统信息，获得GPS(Global Positioning System全球定 位系统)和Galileo (伽利略)系统的Ll频段的卫星导航信号的信息列表；步骤三将集合M中属于步骤二中所得信息表中的卫星信号的成员移入集合Mk，Mk ={GPSCA, GPS LlP(Y)，GPS LlC, Galileo El OS}的集合大小为 Nk = 4 ；集合 M 的大小变 为N = 1，此时集合M为待检测集合，Mk为已知信号集合；步骤四将集合#中的元素导航信号&(0&=1，...炉)的码长特征参数Cik初 始化为^t (Lk为第k个信号伪随机码码长，fkc为第k个信号码片速率，fkd为第k个信号导航数据速率)，将ak< 1的信号放入集合A (短码集合)，将Cik彡1的信号放入集合 B (长码集合)，具体计算数值如表(1)所示； 表 1步骤五对待检测集合M中的元素X1 (t)进行如下处理1)将待检测导航信号X1 (t)的码长特征参数α ！初始化为= 2046 * 50 /1.023 x IO6 = 0.1 ；
fic2)根据待检测导航信号的码长特征参数α i，计算导航信号\(0与集合A中每个 元素的频谱分离系数kAi (i = 1，. . .，η)；计算导航信号Xi (t)与集合B中每个元素的频谱 分离系数 kBi(i = 1，...，m)。因为待检测导航信号的码长特征参数α ！ < 1，则以公式⑴计算导航信号X1 (t) 与集合A中元素GPS CA的频谱分离系数kA1，结果为-50dB/Hz到-75dB/Hz的周期波动值， 如图(3)所示；以公式(2)计算导航信号Xi(t)与集合B中每个元素的频谱分离系数kB1 =-71. 16dB/Hz, kB2 = -65. 23dB/Hz, kB3 = -68. 04dB/Hz ；3)取频谱分离系数 kAi (i = 1，· · ·，η)和 kBi (i = 1，· · ·，m)中的最大值 _50dB/ Hz，做为导航信号X1 (t)的Ic1频谱分离系数；步骤六设置系统最终频谱分离系数阈值Lim为-50dB/Hz，因为Ic1彡Lim，所以检 测结果为导航信号Xi (t)频谱分离系数检测达标，符合卫星导航信号兼容性要求，基于码长 特征的频谱分离系数检测完成。本实施例的频谱分离系数的检测方法，相比现有技术检测方法，得到频谱分离系 数误差仅为0. 232dB/Hz，而运算速度提高了 1000倍以上。
权利要求
一种基于码长特征的频谱分离系数检测方法，其特征在于，包括步骤如下步骤一、系统启动，在检测设备上设定需要检测的目标频段，得到位于该频段上的所有卫星信号的集合M＝{x1(t)，x2(t)，x3(t)，…xn(t)}，并且得到初始集合的大小N；步骤二、根据目标频段的系统信息，获得GPS和Galileo系统的相应频段的卫星导航信号的信息列表；步骤三、将集合M中属于步骤二中所得信息表中的卫星信号的成员移入集合Mk，Mk的集合大小为Nk＝k；集合M的大小变为N＝N k，此时集合M为待检测集合，Mk为已知信号集合；步骤四、将集合Mk中的元素导航信号xk(t)(k＝1，...Nk)的码长特征参数αk初始化为其中Lk为第k个信号伪随机码码长，fkc为第k个信号码片速率，fkd为第k个信号导航数据速率；将αk＜1的信号放入短码集合A，假设集合元素为n，将αk≥1的信号放入长码集合B，假设集合元素为m；步骤五、对所述的待检测集合M中的元素xi(t)(i＝,1，...，N)根据码长特征参数αi逐一处理，计算频谱分离系数κi；步骤六、设置系统最终频谱分离系数阈值Lim，默认阈值为 50dB/Hz，如果ki＞Lim，则检测结果为导航信号xi(t)不能达到频谱分离系数检测标准，如果ki≤Lim，则检测结果为导航信号xi(t)频谱分离系数检测达标，符合卫星导航信号兼容性要求，基于码长特征的频谱分离系数检测完成。FDA0000024096860000011.tif
2.根据权利要求1所述的基于码长特征的频谱分离系数检测方法，其特征是，所述的 计算频谱分离系数k,为1)将导航信号Xi(t)的码长特征参数a ,初始化为f ； 2)根据待检测导航信号的码长特征参数a,，计算导航信号Xi(t)与集合A中每个元素 的频谱分离系数kAi(i = 1，. . .，n)；计算导航信号Xi(t)与集合B中每个元素的频谱分离 系数 kBi(i = 1, . . . , m)；取所述的频谱分离系数kAi(i = 1，. . .，n)和kBi(i = 1，. . .，m)中的最大值，做为导航 信号Xi(t)的频谱分离系数Ki。
3.根据权利要求2所述的基于码长特征的频谱分离系数检测方法，其特征是，所述的 频谱分离系数k ,，待检测导航信号的码长特征参数a,< 1并且干扰信号属于集合A时， K i为 其中，为接收机前端带宽，Gs(f)为期望信号的功率谱密度，Gi(f)为干扰信号的功 率谱密度，Af为多普勒频移；所述的频谱分离系数k”待检测导航信号的码长特征参数、> 1时，无论干扰信号 属于集合A或者B，、为: 其中，为接收机前端带宽，Gs。(f)为期望信号的码片功率谱密度，GiJf)为干扰信号 的码片功率谱密度；计算频谱分离系数时，待检测导航信号为期望信号，已知信号集合A、B中的信号均为 干扰信号。
4.根据权利要求3所述的基于码长特征的频谱分离系数检测方法，其特征是，所述的 功率谱密度是指用密度的概念表示信号功率在各频率点的分布情况，对功率谱在频域上 积分就得到信号的功率，使用周期性伪随机码且包含数据调制的卫星导航信号的功率谱密 SGx(f)其中L为信号伪随机码码长，fc为信号码片速率，fd为导航数据速率，DFT(Ck)为伪随 机码的离散傅立叶变换，Φ为调制波形的傅立叶变换。
全文摘要
一种卫星导航技术领域的基于码长特征的频谱分离系数检测方法。包括步骤如下得到所有卫星信号初始集合的大小N；获得信号的信息列表；集合M为待检测集合，Mk为已知信号集合；码长特征参数αk初始化为将αk＜1的信号放入短码集合A，假设集合元素为n，将αk≥1的信号放入长码集合B，假设集合元素为m；如果ki＞Lim，则检测结果为导航信号xi(t)不能达到频谱分离系数检测标准，如果ki≤Lim，则检测结果为导航信号xi(t)频谱分离系数检测达标，符合卫星导航信号兼容性要求，基于码长特征的频谱分离系数检测完成。本发明在不引入可见的计算误差的前提下，大大提高频谱分离系数的检测效率，使得全球导航卫星系统的兼容性评估有快速、准确的实现途径。
文档编号G01S19/30GK101900820SQ20101024560
公开日2010年12月1日 申请日期2010年8月5日 优先权日2010年8月5日
发明者刘卫, 刘莉, 张炎华, 战兴群, 杜刚, 牛满仓, 翟传润 申请人:上海交通大学