高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

利用部分可信信号的导航终端欺骗干扰检测方法

王环宇 林红磊 欧钢 唐小妹

高健, 张耀宗, 孟祥蕊, 马芳卉. 几类指标为2的不可约拟循环码的重量分布[J]. 电子与信息学报, 2022, 44(12): 4312-4318. doi: 10.11999/JEIT211104
引用本文: 王环宇, 林红磊, 欧钢, 唐小妹. 利用部分可信信号的导航终端欺骗干扰检测方法[J]. 电子与信息学报, 2024, 46(10): 4053-4061. doi: 10.11999/JEIT240067
GAO Jian, ZHANG Yaozong, MENG Xiangrui, MA Fanghui. Weight Distributions of Some Classes of Irreducible Quasi-cyclic Codes of Index 2[J]. Journal of Electronics & Information Technology, 2022, 44(12): 4312-4318. doi: 10.11999/JEIT211104
Citation: WANG Huanyu, LIN Honglei, OU Gang, TANG Xiaomei. The Spoofing Detection Method of Navigation Terminal Using Partial Authenticated Signals[J]. Journal of Electronics & Information Technology, 2024, 46(10): 4053-4061. doi: 10.11999/JEIT240067

利用部分可信信号的导航终端欺骗干扰检测方法

doi: 10.11999/JEIT240067
详细信息
    作者简介:

    王环宇:女,博士生,研究方向为卫星导航信号认证

    林红磊:男,副研究员,研究方向为卫星导航系统与信号处理

    欧钢:男,教授,研究方向为卫星导航定位、综合导航定位授时

    唐小妹:女,研究员,研究方向为导航信号体制设计、导航安全对抗等

    通讯作者:

    林红磊 linhonglei@nudt.edu.cn

  • 中图分类号: TN967.1

The Spoofing Detection Method of Navigation Terminal Using Partial Authenticated Signals

  • 摘要: 导航信号认证服务处于初步部署阶段,认证信号对地覆盖重数无法满足独立定位授时需求,现有研究对这一阶段利用部分通过认证的信号,即可信信号,实现欺骗检测的方法关注度较低。针对这一现状,该文根据欺骗攻击原理,提出以可信信号为基准,基于可信信号伪距残差的欺骗检测方法,建立该场景下的欺骗检测模型,并分析影响所提方法检测性能的因素。经过仿真,在可信卫星数目为3颗、用户定位精度约10 m条件下,当欺骗导致的定位偏差为100 m时,该方法的平均欺骗检测概率可达0.96。此外,该文对算法欺骗检测盲区进行了分析,证明所提算法对于绝大部分欺骗导致的定位结果均有效。
  • 线性码的重量分布是编码理论中的一个重要研究问题。近几年,线性码的重量分布,尤其是循环码的重量分布,被国内外的编码学者广泛关注与研究。Ding等人[1]研究了不可约循环码的重量分布,利用不可约循环码构造了一些参数较好的线性码。基于指数和理论,文献[2,3]在有限域上某些可约循环码的重量分布研究中也取得了很大进展。

    少重量线性码,如常重量线性码[4]、2-重量线性码[5]、3-重量线性码[6]以及其他少重量线性码[7,8]等,是重要的线性码类,可用于认证码[9]、结合方案以及秘密共享方案的构造。Ding[10]基于不可约循环码的重量分布构造了几类3-重量线性码。Schmidt等人[11]基于离散傅里叶变换和高斯和理论得到了不可约循环码至多有两个重量的充要条件。Zhou等人[12]构造了7类3-重量循环码并分析了由这些3-重量循环码得到的秘密共享方案的结构。

    拟循环码是一类重要的线性码,它与卷积码和低密度校验码密切相关。指标为2的拟循环码也称为分块长度相等的双循环码。Borges等人[13]给出了2元域上双循环码的显式生成元并确定了双循环码及其对偶码生成元之间的关系。Gao等人[14,15]给出了4元双循环码的生成元以及与对偶码生成元之间的关系,并且证明了4元双循环码是渐进优的。Patanker等人[16]利用高斯和确定了几类2元双循环码的重量分布。

    本文主要利用有限域上指标为2的不可约拟循环码构造少重量的线性码。首先,基于有限域上的高斯周期,本文给出了几类指标为2的不可约拟循环码的重量分布;其次,基于不可约拟循环码的重量分布,本文构造了几类2-重量线性码和3-重量线性码,其中包括3类最优的2-重量线性码。

    p是一个素数,q=pmr=qt,其中mt是正整数。令Fq表示q元有限域。 Fql的线性子空间称为有限域Fq上码长为l的线性码。Ai表示线性码C中Hamming重量为i的码字个数。 定义码长为l的线性码C的重量分布多项式为1+A1x++Alxl

    l=m0+m1。如果C中的任意码字c=(c1,0,c1,1,,c1,m01|c2,0,c2,1,,c2,m11)经循环移位T作用后有T(c)=(c1,m01,c1,0,,c1,m02|c2,m11,c2,0,, c2,m12)C则称线性码C是分块长度为(m0,m1)的双循环码。若m0=m1,则称线性码C是指标为2的拟循环码。

    Fr=Fqt表示qt元有限域且ζ为其本原元。设θFr=Fr{0}ord(θ)=nh(x)θ1Fq上的极小多项式,deg(h(x))=s,则Fq[x]/Fq[x]h(x)h(x)=Fqs是一个qs元有限域且θFqs=Fqs{0}。因为deg(h(x))=s,所以s|t(qs1)|(qt1)。令Fqs=w,其中w=ζqt1qs1,则wFqs的本原元。

    假设θ=wN,其中nN=qs1。定义指标为2的不可约拟循环码为

    C={(Tr(αβ0),Tr(αβ0θ),,Tr(αβ0θn1)|Tr(αβ1),Tr(αβ1θ),,Tr(αβ1θn1))|αFqs}
    (1)

    其中,β0,β1FqsθTr()表示从FqsFq的迹映射。此外,CFq上的校验多项式为θ1的极小多项式h(x)且维数为dim(C)=s

    Trq/qpp表示从FqFp的迹映射,χ是从Fq到模长为1的复数组成的乘法群的映射,对任意的x,yFqχ(x+y)=χ(x)χ(y)。设bFq,定义

    χb(c)=e2π1Trq/p(bc)/p,cFq
    (2)

    Fq上的加法特征。如果b=0,则χ0(c)=1χ0称为Fq上的平凡加法特征;如果b=1,则χ1称为Fq上的标准加法特征。

    定义Ci(N,qs)=wiwNi=0,1,,N,其中wN表示Fqs的一个子群,陪集Ci(N,qs)称为Fqs的阶为N的分圆类。根据加法特征和分圆类,定义Fqs上的高斯周期为

    ηi(N,qs)=xCi(N,qs)χ(x),i=0,1,,N

    其中,χFqs上的标准加法特征。一般情况下,高斯周期的值很难计算,但是可以通过高斯多项式ψ(N,qs)(X)[17]得到一些特殊的值,其中高斯多项式定义为

    ψ(N,qs)(X)=N1i=0(Xηi(N,qs))

    引理1[1] 当N=2时,高斯周期的值为

    η0(2,qs)={1+(1)sm1qs22,p1(mod4)1+(1)sm1(1)smqs22,p3(mod4)

    η1(2,qs)=1η0(2,qs)

    引理2[1] 当N=3时,高斯多项式ψ(3,qs)(X)的分解如下:

    (1)如果p2(mod3),则sm为偶数并且

    ψ(3,qs)(X)={33(3X+1+2qs2)(3X+1qs2)2sm233(3X+12qs2)(3X+1+qs2)2sm2

    (2)如果p1(mod3),sm0(mod3),则

    ψ(3,qs)(X)=127(3X+1c1qs3)(3X+1+12(c1+9d1)qs3)(3X+1+12(c19d1)qs3)

    其中,c1d1满足4qs3=c21+27d21,c11(mod3),gcd(c1,p)=1

    引理3[1] 当N=4时,高斯多项式ψ(4,qs)(X)的分解如下:

    (1)如果p3(mod4),则sm为偶数并且

    ψ(3,qs)(X)={44(4X+1+3qs2)(4X+1qs2)3sm244(4X+13qs2)(4X+1+qs2)3sm2

    (2)如果p1(mod4),sm0(mod4),则

    ψ(3,qs)(X)=44((4X+1)+qs2+2u1qs4)((4X+1)+qs22u1qs4)3×44((4X+1)qs2+4v1qs4)((4X+1)qs24v1qs4)

    其中,u1,v1满足qs2=u21+4v21,u11(mod4),gcd(u1,p)=1

    为了确定指标为2的不可约拟循环码的重量分布,本文还需要以下引理。

    引理 4[1] 设e1qs1的正因数并且0ie1,则

    {xy:yFq,xCi(e1,qs)}=q1e1gcd(qs1q1,e1)Ci(gcd(qs1)q1,e1),qs)

    其中q1e1gcd(qs1q1,e1)C(gcd(qs1q1,e1),qs)i表示C(gcd(qs1q1,e1),qs)i中的任意元素在集合xy:yFq,xC(e1,qs)i中出现的次数均为q1e1gcd(qs1q1,e1)

    N1是一个正整数且满足n=qs1N,N1=gcd(qs1q1,N)。令θ=wN,其中wFqs的本原元。Z(qs,β0)Z(qs,β1)分别为使等式Tr(αβ0θi)=0Tr(αβ1θj)=0成立的i,j的个数。对于C中任意的码字c,如果α0,则它的Hamming重量为WH(c)=2nZ(qs,β0)Z(qs,β1),其中

    Z(qs,β0)=|{0in1:Tr(αβ0θi)=0}|Z(qs,β1)=|{0jn1:Tr(αβ1θj)=0}|

    χFq上的加法特征,由引理4有

    Z(qs,β0)=nq+(q1)N1Nqηk(N1,qs)Z(qs,β1)=nq+(q1)N1Nqηl(N1,qs)

    其中,k,l=0,1,,N11

    下面给出当β0,β1Fqsθ,N1=2,3,4时,式(1)定义的指标为2的不可约拟循环码的重量分布。显然,C0(N1,qs)C0(N,qs),C1(N1,qs),,C(N1,qs)N11构成了Fqsθ的一个划分。定义SFqsθ中的一个分圆类。

    定理1 设Nqs1的正因数且使得N1=gcd(qs1q1,N)=2,则

    情形(1)如果β0,β1在相同的分圆类S中,则由式(1)定义的不可约拟循环码的重量分布如表1所示。

    表 1  情形(1):不可约拟循环码的重量分布
    重量(i)频数(Ai)
    01
    2Nq(q1)(qs12η(2,qs)0)qs12
    2Nq(q1)(qs12η(2,qs)1)qs12
    下载: 导出CSV 
    | 显示表格

    情形(2)如果β0,β1在不同的分圆类S中,则由式(1)定义的不可约拟循环码的重量分布如表2所示。

    表 2  情形(2):不可约拟循环码的重量分布
    重量(i)频数(Ai)
    01
    2N(q1)qs1qs1
    下载: 导出CSV 
    | 显示表格

    证明 对于情形(1)本文只给出当β0,β1C(2,qs)0C(N,qs)0时的证明过程。

    αC(2,qs)0时,因为β0,β1C(2,qs)0C(N,qs)0,所以αβ0,αβ1C(2,qs)0。此时,对于cC,有WH(c)=2n2(nq+2(q1)Nqη(2,qs)0),出现的频数为qs12。同理,当αC(2,qs)1时,对于cC,有WH(c)=2n2(nq+2(q1)Nqη(2,qs)1),出现的频数为qs12。因此,式(1)定义的指标为2的不可约拟循环码的重量分布如表1所示。

    对于情形(2),可采用相同的讨论方法。证毕

    推论1  如果β0,β1满足定理1中的情形(1),则式(1)定义的指标为2的不可约拟循环码是一个参数为[2(qs1)N,s,2(q1)(qsqs2)Nq]的2-重量线性码,其重量分布为1+qs12x2(q1)(qsqs2)Nq+qs12x2(q1)(qs+qs2)Nq

    证明 如果β0,β1满足定理1中的情形(1),则由引理1可知,当p1(mod4)时,η(2,qs)0=qs2+12,η(2,qs)1=qs212;当p3(mod4)时,η(2,qs)0=qs212,η(2,qs)1=qs2+12。将η(2,qs)0,η(2,qs)1的值分别代入表1中,可得到式(1)中定义的指标为2的不可约拟循环码是一个参数为[2(qs1)N,s,2(q1)(qsqs2)Nq]的2-重量线性码,其重量分布为1+qs12x2(q1)(qsqs2)Nq+qs12x2(q1)(qs+qs2)Nq。证毕

    类似于定理1,当N1=gcd(qs1q1,N)=3时,式(1)定义的指标为2的不可约拟循环码的重量分布由以下定理给出。

    定理2 设Nqs1的正因数且使得N1=gcd(qs1q1,N)=3,则

    情形(3) 如果β0,β1在相同的分圆类S中,则式(1)定义的不可约拟循环码的重量分布如表3所示。

    表 3  情形(3):不可约拟循环码的重量分布
    重量(i)频数(Ai)
    01
    2Nq(q1)(qs13η(3,qs)0)qs13
    2Nq(q1)(qs13η(3,qs)1)qs13
    2Nq(q1)(qs13η(3,qs)2)qs13
    下载: 导出CSV 
    | 显示表格

    情形(4) 如果β0,β1在不同的分圆类S中,则式(1)定义的不可约拟循环码的重量分布如表4所示。

    表 4  情形(4):不可约拟循环码的重量分布
    重量(i)频数(Ai)
    01
    1Nq(q1)[2(qs1)3(η(3,qs)0+η(3,qs)1)]qs13
    1Nq(q1)[2(qs1)3(η(3,qs)1+η(3,qs)2)]qs13
    1Nq(q1)[2(qs1)3(η(3,qs)2+η(3,qs)0)]qs13
    下载: 导出CSV 
    | 显示表格

    推论2 设Nqs1的正因子,N1=gcd(qs1q1,N)=3,p2(mod3)。若β0,β1满足定理2中情形(3)且sm0(mod4),则式(1)定义的不可约拟循环码是参数为[2(qs1)N,s,2(q1)(qsqs2)Nq]的2-重量线性码,其重量分布为

    1+2(qs1)3x2(q1)(qsqs2)Nq+qs13x2(q1)(qs+2qs2)Nq

    sm2(mod4),则式(1)定义的不可约拟循环码是参数为[2(qs1)N,s,2(q1)(qs2qs2)Nq]的2-重量线性码,其重量分布为

    1+qs13x2(q1)(qs2qs2)Nq+2(qs1)3x2(q1)(qs+qs2)Nq

    证明: 当β0,β1满足定理2中情形(3)时,如果sm0(mod4),此时sm2为偶数,由引理2可得η(3,qs)0=1+2qs23,η(3,qs)1=η(3,qs)1=12qs23,将其代入表3,可得到推论2中的第1种结果;如果 sm2(mod4), 此时sm2为奇数,则由引理2可得η(3,qs)0=12qs23,η(3,qs)1=η(3,qs)1=1+2qs23,将其代入表3,可得到推论2中的另一种结果。 证毕

    类似于推论2的证明方法,可得到以下3个结论。

    推论3 设Nqs1的正因子,N1=gcd(qs1q1,N)=3,p2(mod3)。若β0,β1满足定理2中的情形(4)且sm0(mod4),则式(1)定义的不可约拟循环码是参数为[2(qs1)N,s,2(q1)(qsqs2)Nq]的2-重量线性码,其重量分布为

    1+2(qs1)3x(q1)(2qs+qs2)Nq+qs13x2(q1)(qsqs2)Nq

    sm2(mod4),则式(1)定义的不可约拟循环码是参数为[2(qs1)N,s,(q1)(2qsqs2)Nq]的2-重量线性码,其重量分布为

    1+2(qs1)3x(q1)(2qsqs2)Nq+2(qs1)3x2(q1)(qs+qs2)Nq

    推论4 设Nqs1的正因子,N1=gcd(qs1q1,N)=3,p1(mod3)。若β0,β1满足定理2中的情形(3)且sm0(mod3),则式(1)定义的不可约拟循环码是一个参数为[2(qs1)N,s]的3-重量线性码,其重量分布为

    1+qs13x2(q1)(qsc1qs3)Nq+qs13x2(q1)(qs+12(c1+9d1)qs3)Nq+qs13x2(q1)(qs+12(c19d1)qs3)Nq

    推论5 设Nqs1的正因子,N1=gcd(qs1q1,N)=3,p1(mod3)。若β0,β1满足定理2中的情形(4)且sm0(mod3),则式(1)定义的不可约拟循环码是参数为[2(qs1)N,s]的3-重量线性码,其重量分布为

    1+qs13x(q1)(2qs+12(9d1c1)qs3)Nq+qs13x(q1)(2qs12(c1+9d1)qs3)Nq+qs13x(q1)(2qs+c1qs3)Nq

    类似定理1,当N1=gcd(qs1q1,N)=4时,式(1)定义的不可约拟循环码的重量分布由以下定理给出。

    定理3 设Nqs1的正因数使得N1=gcd(qs1q1,N)=4,则

    情形(5)如果β0,β1在相同的分圆类S中,则式(1)定义的不可约拟循环码的重量分布如表5所示。

    表 5  情形(5):不可约拟循环码的重量分布
    重量(i)频数(Ai)
    01
    2Nq(q1)(qs14η(4,qs)0)qs14
    2Nq(q1)(qs14η(4,qs)1)qs14
    2Nq(q1)(qs14η(4,qs)2)qs14
    2Nq(q1)(qs14η(4,qs)3)qs14
    下载: 导出CSV 
    | 显示表格

    情形(6)如果下列条件之一成立,则式(1)定义的不可约拟循环码的重量分布如表6所示。

    表 6  情形(6):不可约拟循环码的重量分布
    重量(i)频数(Ai)
    01
    2Nq(q1)[(qs1)2(η(4,qs)0+η(4,qs)1)]qs14
    2Nq(q1)[(qs1)2(η(4,qs)1+η(4,qs)2)]qs14
    2Nq(q1)[(qs1)2(η(4,qs)2+η(4,qs)3)]qs14
    2Nq(q1)[(qs1)2(η(4,qs)3+η(4,qs)0)]qs14
    下载: 导出CSV 
    | 显示表格

    (a) β0,β1{C(4,qs)0C(N,qs)0,C(4,qs)1}

    (b) β0,β1{C(4,qs)0C(N,qs)0,C(4,qs)3}

    (c) β0,β1{C(4,qs)1,C(4,qs)2}

    (d) β0,β1{C(4,qs)2,C(4,qs)3

    情形(7) 如果下列条件之一成立,则式(1)定义的不可约拟循环码的重量分布如表7所示。

    表 7  情形(7):不可约拟循环码的重量分布
    重量(i)频数(Ai)
    01
    2Nq(q1)[(qs1)2(η(4,qs)0+η(4,qs)2)]qs12
    2Nq(q1)[(qs1)2(η(4,qs)1+η(4,qs)3)]qs12
    下载: 导出CSV 
    | 显示表格

    (e) β0,β1{C(4,qs)0C(N,qs)0,C(4,qs)2}

    (f) β0,β1{C(4,qs)1,C(4,qs)3}

    由引理3、推论2和定理3,可得到以下推论。

    推论6 设Nqs1的正因子,N1=gcd(qs1q1,N)=4,p3(mod4)。若β0,β1满足定理3中的情形(5)且sm0(mod4),则式(1)定义的不可约拟循环码是参数为[2(qs1)N,s,2(q1)(qsqs2)Nq]的2-重量线性码,其重量分布为

    1+3(qs1)4x2(q1)(qsqs2)Nq+qs14x2(q1)(qs+3qs2)Nq

    sm2(mod4),则式(1)定义的不可约拟循环码是参数为[2(qs1)N,s,2(q1)(qs3qs2)Nq]的2-重量线性码,其重量分布为

    1+3(qs1)4x2(q1)(qs+qs2)Nq+qs14x2(q1)(qs3qs2)Nq

    推论7 设Nqs1的正因子,N1=gcd(qs1q1,N)=4p3(mod4)。如果β0,β1满足定理3中的情形(6)或情形(7),则式(1)定义的指标为2的不可约拟循环码是参数为[2(qs1)N,s,2(q1)(qsqs2)Nq]的2-重量线性码,其重量分布为1+qs12x2(q1)(qsqs2)Nq+qs12x2(q1)(qs+qs2)Nq

    推论8 设Nqs1的正因子,N1=gcd(qs1q1,N)=4,p1(mod4)。若β0,β1满足定理3中的情形(7)且sm0(mod4),则式(1)定义的指标为2的不可约拟循环码是参数为[2(qs1)N,s]的2-重量线性码,其重量分布为

    1+qs12x2(q1)(qs+(u1+2v1)qs4)Nq+qs12x2(q1)(qs(u1+2v1)qs4)Nq

    定理4 设Nqs1的正因子,N1=3,p2(mod3),sm2(mod4)。若β0,β1满足定理2中的情形(4),则有

    (1)如果p=2,s=2,N=3(2m1),则式(1)定义的指标为2的不可约拟循环码是一类2-重量MDS码,其参数为[23(2m+1),2,13(2m+11)]

    (2)如果s=2,2N=3(pm1),则式(1)定义的指标为2的不可约拟循环码是一类达到Griesmer界的最优码,其参数为[43(pm+1),2,23(2pm1)]

    证明 将情形(1)中所给的条件代入推论4, 可得式(1)定义的指标为2的不可约拟循环码是一类2-重量线性码且其参数满足n=k+d1。 由引理5可得到一类2-重MDS码,其参数为[23(2m+1),2,13(2m+11)]

    将情形(2)中所给的条件代入推论4,可得到式(1)定义的指标为2的不可约拟循环码是一类2-重量线性码且参数满足n=k1i=0dqi。由引理6可得到一类达到Griesmer界的最优码,其参数为[43(pm+1),2,23(2pm1)]。 证毕

    例1 令(N1,p,m,s,N)=(3,2,3,2,21)。由定理4,式(1)定义的指标为2的不可约拟循环码是F8上参数为[6,2,5]的2-重量MDS码,其重量分布为1+42x5+21x6

    例2 令(N1,p,m,s,N)=(3,11,1,2,15)。由定理4,式(1)定义的指标为2的不可约拟循环码是F11上参数为[16,2,14]的达到Griesmer界的2-重量最优码,其重量分布为1+80x14+40x16

    类似于定理4,可由以下定理得到另一类2-重量MDS码。

    定理5 设Nqs1的正因子,N1=4,p3(mod4)。若β0,β1满足定理3中的情形(6)或情形(7),s=2,N=4(pm1),则式(1)定义的指标为2的不可约拟循环码是一类2-重量MDS码,其参数为[12(pm+1),2,12(pm1)]

    例3 令(N1,p,m,s,N)=(4,3,3,2,104)。由定理5,式(1)定义的指标为2的不可约拟循环码是F27上参数为[14,2,13]的2-重量MDS码,其重量分布为1+364x13+364x14

    本文研究了有限域上指标为2的不可约拟循环码的重量分布,构造了几类2-重量线性码和3-重量线性码。特别地,本文得到了3类最优的2-重量线性码。如何利用ZpZp[v]-加性码[18,19]的重量分布构造有限域上的最优码是一个有意义的研究问题。

  • 图  1  不同可信卫星数目下时间欺骗导致的伪距拉偏距离与平均欺骗检测概率的关系曲线

    图  2  M=2时接收机钟差变化量对算法检测性能的影响曲线

    图  3  不同安全认证信号数目的检测性能

    图  4  观测时刻不同用户测距误差的标准差与不同可信卫星数目单点定位精度统计结果对比图

    图  5  不同卫星数目观测时刻检测统计量取值统计结果对比图

    图  6  伪距测量误差标准差与平均欺骗检测概率的关系

    图  7  欺骗导致的定位偏差r与平均欺骗检测概率之间的关系曲线

    图  8  Δd2与欺骗信号检测概率之间的曲线对比图(接收机位于北京,r=50 m, Pfa=0.01, σURE=5.9 m)

    图  9  算法检测盲区理论计算分布图

    表  1  只存在时间欺骗时,不同可信卫星数目下,算法可检测出的接收机钟差欺骗距离的理论值(s)

    可信卫星数目M 1 2 3
    算法可检测出的接收机钟差欺骗 7.7×108 5.9×108 5.0×108
    下载: 导出CSV

    表  2  仿真场景参数

    参数 描述
    仿真时段 2023年9月2日16:00-2023年9月3日16:00,观测时间间隔为10 min,观测历元共144个
    北斗星座 本文采用北斗三号卫星星座,包括:24颗MEO卫星,Walker24/3/1星座,卫星轨道高度为21 528 km,轨道倾角为55°;
    3颗GEO卫星,轨道高度为35 786 km,分别定位于东经80°, 110.5°, 140°;3颗IGSO卫星,轨道高度35 786 km,
    轨道倾角55°,相位间隔120°
    低轨星座 卫星播发认证信号,为近极地星座,为Walker类型,共有60颗卫星,6个轨道面,相位因子为1,轨道高度为1 175 km,
    轨道倾角为86.5°
    接收机位置 随机选择国内5个城市作为接收机所在位置(地心地固坐标),包括:
    哈尔滨(–2661.32,3576.94,4546.01) km、北京(–2176.85,4387.47,4071.96) km、兰州(–1231,5014.65,3736.14) km、
    成都(–1334.67,5326.66,3234.39) km、海口(–2083.66,5620.29,2172.46) km
    接收机天线 最低仰角为15°
    下载: 导出CSV
  • [1] CARROLL J V, VAN DYKE K, KRAEMER J H, et al. Vulnerability assessment of the U. S. transportation infrastructure that relies on GPS[C]. Proceedings of the 14th International Technical Meeting of the Satellite Division of the Institute of Navigation, Salt Lake City, USA, 2001: 975–981.
    [2] WU Zhijun, ZHANG Yun, YANG Yiming, et al. Spoofing and anti-spoofing technologies of global navigation satellite system: A survey[J]. IEEE Access, 2020, 8: 165444–165496. doi: 10.1109/ACCESS.2020.3022294.
    [3] YUAN Muzi, TANG Xiaomei, and OU Gang. Authenticating GNSS civilian signals: A survey[J]. Satellite Navigation, 2023, 4(1): 6. doi: 10.1186/s43020-023-00094-6.
    [4] NICOLA M, MOTELLA B, PINI M, et al. Galileo OSNMA public observation phase: Signal testing and validation[J]. IEEE Access, 2022, 10: 27960–27969. doi: 10.1109/ACCESS.2022.3157337.
    [5] FERNÁNDEZ I, RIJMEN V, ASHUR T, et al. Galileo navigation message authentication specification for signal-in-space testing - v1.0[R]. European Commission, 2016.
    [6] ANDERSON J M, CARROLL K L, DEVILBISS N P, et al. Chips-message robust authentication (chimera) for GPS civilian signals[C]. Proceedings of the 30th International Technical Meeting of the Satellite Division of the Institute of Navigation, Portland, USA, 2017: 2388–2416. doi: 10.33012/2017.15206.
    [7] CHAPMAN D C. Chips message robust authentication (chimera) enhancement for the L1C signal: Space segment/user segment interface[R]. IS-AGT-100, 2019.
    [8] MANANDHAR D and SHIBASAKI R. Signal authentication for anti-spoofing based on L1S[C]. Proceedings of the ION 2017 Pacific PNT Meeting, Honolulu, USA, 2017: 938–947. doi: 10.33012/2017.15029.
    [9] MANANDHAR D and SHIBASAKI R. Authenticating GALILEO open signal using QZSS signal[C]. Proceedings of the 31st International Technical Meeting of the Satellite Division of the Institute of Navigation, Miami, USA, 2018: 3995–4003. doi: 10.33012/2018.15872.
    [10] 王磊, 李德仁, 陈锐志, 等. 低轨卫星导航增强技术——机遇与挑战[J]. 中国工程科学, 2020, 22(2): 144–152. doi: 10.15302/J-SSCAE-2020.02.018.

    WANG Lei, LI Deren, CHEN Ruizhi, et al. Low earth orbiter (LEO) navigation augmentation: Opportunities and challenges[J]. Strategic Study of CAE, 2020, 22(2): 144–152. doi: 10.15302/J-SSCAE-2020.02.018.
    [11] SUN Tianyu, HU Min, and YUN Chaoming. Low-orbit large-scale communication satellite constellation configuration performance assessment[J]. International Journal of Aerospace Engineering, 2022, 2022: 4918912. doi: 10.1155/2022/4918912.
    [12] S·柳辛. 使用PVT解估算来检测和消除GNSS欺骗信号[P]. 中国, 110114695A, 2019.

    S·LIUXIN. Detection and elimination of GNSS spoofing signals with PVT solution estimation[P]. CN, 110114695A, 2019.
    [13] 张超, 吕志伟, 张伦东, 等. 欺骗干扰对GNSS/INS系统定位性能的影响[J]. 导航定位学报, 2022, 10(4): 20–28. doi: 10.3969/j.issn.2095-4999.2022.04.003.

    ZHANG Cao, LYU Zhiwei, ZHANG Lundong, et al. Influence analysis of spoofing interference on positioning performance of GNSS/INS system[J]. Journal of Navigation and Positioning, 2022, 10(4): 20–28. doi: 10.3969/j.issn.2095-4999.2022.04.003.
    [14] 谢钢. GPS原理与接收机设计[M]. 北京: 电子工业出版社, 2009.

    XIE G. Principles of GPS and Receiver Design[M]. Beijing: Publishing House of Electronics Industry, 2009.
    [15] PARKINSON B W and AXELRAD P. Autonomous GPS integrity monitoring using the pseudorange residual[J]. Navigation, 1988, 35(2): 255–274. doi: 10.1002/j.2161-4296.1988.tb00955.x.
    [16] FU Dong, PENG Jing, GONG Hang, et al. Impact analysis of meaconing attack on timing receiver[M]. YANG Changfeng and XIE Jun. China Satellite Navigation Conference (CSNC 2021) Proceedings. Singapore: Springer, 2021: 423–434. doi: 10.1007/978-981-16-3146-7_39.
  • 期刊类型引用(2)

    1. 李艳俊,王琦,项勇,谢惠琴. 动态密码组件对合MDS矩阵的设计与实现. 微电子学与计算机. 2024(07): 37-45 . 百度学术
    2. 李建佳,蔚保国,鲍亚川,杨梦焕,崔宋祚,田润泽,赵军. 一种适应复杂地下环境的导航增强型脉冲超宽带通导一体自组网体制设计研究. 电子与信息学报. 2023(11): 4179-4189 . 本站查看

    其他类型引用(1)

  • 加载中
图(9) / 表(2)
计量
  • 文章访问数:  207
  • HTML全文浏览量:  120
  • PDF下载量:  33
  • 被引次数: 3
出版历程
  • 收稿日期:  2024-01-29
  • 修回日期:  2024-09-05
  • 网络出版日期:  2024-09-09
  • 刊出日期:  2024-10-30

目录

/

返回文章
返回