The Spoofing Detection Method of Navigation Terminal Using Partial Authenticated Signals
-
摘要: 导航信号认证服务处于初步部署阶段,认证信号对地覆盖重数无法满足独立定位授时需求,现有研究对这一阶段利用部分通过认证的信号,即可信信号,实现欺骗检测的方法关注度较低。针对这一现状,该文根据欺骗攻击原理,提出以可信信号为基准,基于可信信号伪距残差的欺骗检测方法,建立该场景下的欺骗检测模型,并分析影响所提方法检测性能的因素。经过仿真,在可信卫星数目为3颗、用户定位精度约10 m条件下,当欺骗导致的定位偏差为100 m时,该方法的平均欺骗检测概率可达0.96。此外,该文对算法欺骗检测盲区进行了分析,证明所提算法对于绝大部分欺骗导致的定位结果均有效。
-
关键词:
- 抗欺骗 /
- 导航认证信号 /
- 接收机自主完好性监测 /
- 伪距残差
Abstract: The navigation signal authentication service is in the initial stage. The coverage multiple numbers of the authentication signal to ground can not meet the requirement of independent positioning and timing. The existing research has paid little attention to the deception detection method based on partially trusted signals at this stage. Aiming at the status quo, according to the principle of spoofing attack, a spoofing detection method is proposed based on the pseudo-distance residual of the authentication signal, and the spoofing detection model is established in this scenario, and the factors that affect the detection performance of the proposed method are analyzed. After simulation, the average deception detection probability of the algorithm can reach 0.96 when the positioning deviation is 100 m, the positioning accuracy is about 10 m, and the number of trusted satellites is 3. In addition, the blind area of the algorithm is analyzed, and it is proved that the algorithm is effective for most of the deception positions. -
1. 引言
线性码的重量分布是编码理论中的一个重要研究问题。近几年,线性码的重量分布,尤其是循环码的重量分布,被国内外的编码学者广泛关注与研究。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-重量线性码。
2. 基本知识
令
p 是一个素数,q=pm ,r=qt ,其中m 和t 是正整数。令Fq 表示q 元有限域。Fql 的线性子空间称为有限域Fq 上码长为l 的线性码。Ai 表示线性码C 中Hamming重量为i 的码字个数。 定义码长为l 的线性码C 的重量分布多项式为1+A1x+⋯+Alxl 。令
l=m0+m1 。如果C 中的任意码字c=(c1,0,c1,1,⋯,c1,m0−1|c2,0,c2,1,⋯,c2,m1−1) 经循环移位T 作用后有T(c)=(c1,m0−1,c1,0,⋯,c1,m0−2|c2,m1−1,c2,0,⋯ ,c2,m1−2)∈C 则称线性码C 是分块长度为(m0,m1) 的双循环码。若m0=m1 ,则称线性码C 是指标为2的拟循环码。令
Fr=Fqt 表示qt 元有限域且ζ 为其本原元。设θ∈Fr∗=Fr∖{0} ,ord(θ)=n ,h(x) 是θ−1 在Fq 上的极小多项式,deg(h(x))=s ,则Fq[x]/Fq[x]⟨h(x)⟩⟨h(x)⟩=Fqs 是一个qs 元有限域且θ∈Fqs∗=Fqs∖{0} 。因为deg(h(x))=s ,所以s|t ,(qs−1)|(qt−1) 。令Fqs∗=⟨w⟩ ,其中w=ζqt−1qs−1 ,则w 为Fqs 的本原元。假设
θ=wN ,其中nN=qs−1 。定义指标为2的不可约拟循环码为C={(Tr(αβ0),Tr(αβ0θ),⋯,Tr(αβ0θn−1)|Tr(αβ1),Tr(αβ1θ),⋯,Tr(αβ1θn−1))|α∈Fqs} (1) 其中,
β0,β1∈Fqs∗∖⟨θ⟩ ,Tr(⋅) 表示从Fqs 到Fq 的迹映射。此外,C 在Fq 上的校验多项式为θ−1 的极小多项式h(x) 且维数为dim(C)=s 。设
Trq/qpp 表示从Fq 到Fp 的迹映射,χ 是从Fq 到模长为1的复数组成的乘法群的映射,对任意的x,y∈Fq 有χ(x+y)=χ(x)χ(y) 。设b∈Fq ,定义χb(c)=e2π√−1Trq/p(bc)/p,∀c∈Fq (2) 为
Fq 上的加法特征。如果b=0 ,则χ0(c)=1 ,χ0 称为Fq 上的平凡加法特征;如果b=1 ,则χ1 称为Fq 上的标准加法特征。定义
Ci(N,qs)=wi⟨wN⟩,i=0,1,⋯,N ,其中⟨wN⟩ 表示Fqs∗ 的一个子群,陪集Ci(N,qs) 称为Fqs 的阶为N 的分圆类。根据加法特征和分圆类,定义Fqs 上的高斯周期为ηi(N,qs)=∑x∈Ci(N,qs)χ(x),i=0,1,⋯,N 其中,
χ 是Fqs 上的标准加法特征。一般情况下,高斯周期的值很难计算,但是可以通过高斯多项式ψ(N,qs)(X) [17]得到一些特殊的值,其中高斯多项式定义为ψ(N,qs)(X)=N−1∏i=0(X−ηi(N,qs)) 引理1[1] 当
N=2 时,高斯周期的值为η0(2,qs)={−1+(−1)s⋅m−1qs22,p≡1(mod4)−1+(−1)s⋅m−1(√−1)smqs22,p≡3(mod4) 和
η1(2,qs)=−1−η0(2,qs) 。引理2[1] 当
N=3 时,高斯多项式ψ(3,qs)(X) 的分解如下:(1)如果
p≡2(mod3) ,则s⋅m 为偶数并且ψ(3,qs)(X)={3−3(3X+1+2qs2)(3X+1−qs2)2,sm2为偶数3−3(3X+1−2qs2)(3X+1+qs2)2,sm2为奇数 (2)如果
p≡1(mod3),s⋅m≡0(mod3) ,则ψ(3,qs)(X)=127(3X+1−c1qs3)⋅(3X+1+12(c1+9d1)qs3)⋅(3X+1+12(c1−9d1)qs3) 其中,
c1 和d1 满足4qs3=c21+27d21,c1≡1(mod3),gcd(c1,p)=1 。引理3[1] 当
N=4 时,高斯多项式ψ(4,qs)(X) 的分解如下:(1)如果
p≡3(mod4) ,则s⋅m 为偶数并且ψ(3,qs)(X)={4−4(4X+1+3qs2)(4X+1−qs2)3,sm2为偶数4−4(4X+1−3qs2)(4X+1+qs2)3,sm2为奇数 (2)如果
p≡1(mod4),s⋅m≡0(mod4) ,则ψ(3,qs)(X)=4−4((4X+1)+qs2+2u1qs4)⋅((4X+1)+qs2−2u1qs4)3×4−4((4X+1)−qs2+4v1qs4)⋅((4X+1)−qs2−4v1qs4) 其中,
u1,v1 满足qs2=u21+4v21,u1≡1(mod4),gcd(u1,p)=1 。为了确定指标为2的不可约拟循环码的重量分布,本文还需要以下引理。
引理 4[1] 设
e1 是qs−1 的正因数并且0≤i≤e1 ,则{xy:y∈Fq∗,x∈Ci(e1,qs)}=q−1e1gcd(qs−1q−1,e1)⋅Ci(gcd(qs−1)q−1,e1),qs), 其中
q−1e1gcd(qs−1q−1,e1)⋅C(gcd(qs−1q−1,e1),qs)i 表示C(gcd(qs−1q−1,e1),qs)i 中的任意元素在集合{ xy:y∈F∗q,x∈C(e1,qs)i} 中出现的次数均为q−1e1gcd(qs−1q−1,e1) 。3. 几类指标为2的不可约拟循环码的重量分布
设
N>1 是一个正整数且满足n=qs−1N,N1=gcd(qs−1q−1,N) 。令θ=wN ,其中w 为Fqs 的本原元。Z(qs,β0) 和Z(qs,β1) 分别为使等式Tr(αβ0θi)=0 和Tr(αβ1θj)=0 成立的i,j 的个数。对于C 中任意的码字c ,如果α≠0 ,则它的Hamming重量为WH(c)=2n−Z(qs,β0)−Z(qs,β1) ,其中Z(qs,β0)=|{0≤i≤n−1:Tr(αβ0θi)=0}|,Z(qs,β1)=|{0≤j≤n−1:Tr(αβ1θj)=0}| 令
χ′ 为Fq 上的加法特征,由引理4有Z(qs,β0)=nq+(q−1)N1Nq⋅ηk(N1,qs),Z(qs,β1)=nq+(q−1)N1Nq⋅ηl(N1,qs) 其中,
k,l=0,1,⋯,N1−1 。下面给出当
β0,β1∈F∗qs∖⟨θ⟩,N1=2,3,4 时,式(1)定义的指标为2的不可约拟循环码的重量分布。显然,C0(N1,qs)∖C0(N,qs),C1(N1,qs),⋯,C(N1,qs)N1−1 构成了F∗qs∖⟨θ⟩ 的一个划分。定义S 为F∗qs∖⟨θ⟩ 中的一个分圆类。定理1 设
N 是qs−1 的正因数且使得N1=gcd(qs−1q−1,N)=2 ,则情形(1)如果
β0,β1 在相同的分圆类S 中,则由式(1)定义的不可约拟循环码的重量分布如表1所示。表 1 情形(1):不可约拟循环码的重量分布重量(i) 频数(Ai) 0 1 2Nq(q−1)(qs−1−2η(2,qs)0) qs−12 2Nq(q−1)(qs−1−2η(2,qs)1) qs−12 情形(2)如果
β0,β1 在不同的分圆类S 中,则由式(1)定义的不可约拟循环码的重量分布如表2所示。表 2 情形(2):不可约拟循环码的重量分布重量(i) 频数(Ai) 0 1 2N(q−1)qs−1 qs−1 证明 对于情形(1)本文只给出当
β0,β1∈C(2,qs)0∖C(N,qs)0 时的证明过程。当
α∈C(2,qs)0 时,因为β0,β1∈C(2,qs)0∖C(N,qs)0 ,所以αβ0,αβ1∈C(2,qs)0 。此时,对于c∈C ,有WH(c)=2n−2(nq+2(q−1)Nq⋅η(2,qs)0) ,出现的频数为qs−12 。同理,当α∈C(2,qs)1 时,对于c∈C ,有WH(c)=2n−2(nq+2(q−1)Nq⋅η(2,qs)1) ,出现的频数为qs−12 。因此,式(1)定义的指标为2的不可约拟循环码的重量分布如表1所示。对于情形(2),可采用相同的讨论方法。证毕
推论1 如果
β0,β1 满足定理1中的情形(1),则式(1)定义的指标为2的不可约拟循环码是一个参数为[2(qs−1)N,s,2(q−1)(qs−qs2)Nq] 的2-重量线性码,其重量分布为1+qs−12x2(q−1)(qs−qs2)Nq+qs−12x2(q−1)(qs+qs2)Nq 。证明 如果
β0,β1 满足定理1中的情形(1),则由引理1可知,当p≡1(mod4) 时,η(2,qs)0=−qs2+12,η(2,qs)1=qs2−12 ;当p≡3(mod4) 时,η(2,qs)0=qs2−12,η(2,qs)1=−qs2+12 。将η(2,qs)0,η(2,qs)1 的值分别代入表1中,可得到式(1)中定义的指标为2的不可约拟循环码是一个参数为[2(qs−1)N,s,2(q−1)(qs−qs2)Nq] 的2-重量线性码,其重量分布为1+qs−12x2(q−1)(qs−qs2)Nq+qs−12x2(q−1)(qs+qs2)Nq 。证毕类似于定理1,当
N1=gcd(qs−1q−1,N)=3 时,式(1)定义的指标为2的不可约拟循环码的重量分布由以下定理给出。定理2 设
N 是qs−1 的正因数且使得N1=gcd(qs−1q−1,N)=3 ,则情形(3) 如果
β0,β1 在相同的分圆类S 中,则式(1)定义的不可约拟循环码的重量分布如表3所示。表 3 情形(3):不可约拟循环码的重量分布重量(i) 频数(Ai) 0 1 2Nq(q−1)(qs−1−3η(3,qs)0) qs−13 2Nq(q−1)(qs−1−3η(3,qs)1) qs−13 2Nq(q−1)(qs−1−3η(3,qs)2) qs−13 情形(4) 如果
β0,β1 在不同的分圆类S 中,则式(1)定义的不可约拟循环码的重量分布如表4所示。表 4 情形(4):不可约拟循环码的重量分布重量(i) 频数(Ai) 0 1 1Nq(q−1)[2(qs−1)−3(η(3,qs)0+η(3,qs)1)] qs−13 1Nq(q−1)[2(qs−1)−3(η(3,qs)1+η(3,qs)2)] qs−13 1Nq(q−1)[2(qs−1)−3(η(3,qs)2+η(3,qs)0)] qs−13 推论2 设
N 是qs−1 的正因子,N1=gcd(qs−1q−1,N)=3,p≡2(mod3) 。若β0,β1 满足定理2中情形(3)且s⋅m≡0(mod4) ,则式(1)定义的不可约拟循环码是参数为[2(qs−1)N,s,2(q−1)(qs−qs2)Nq] 的2-重量线性码,其重量分布为1+2(qs−1)3x2(q−1)(qs−qs2)Nq+qs−13x2(q−1)(qs+2qs2)Nq。 若
s⋅m≡2(mod4) ,则式(1)定义的不可约拟循环码是参数为[2(qs−1)N,s,2(q−1)(qs−2qs2)Nq] 的2-重量线性码,其重量分布为1+qs−13x2(q−1)(qs−2qs2)Nq+2(qs−1)3x2(q−1)(qs+qs2)Nq。 证明: 当
β0,β1 满足定理2中情形(3)时,如果s⋅m≡0(mod4) ,此时s⋅m2 为偶数,由引理2可得η(3,qs)0=−1+2qs23,η(3,qs)1=η(3,qs)1=−1−2qs23 ,将其代入表3,可得到推论2中的第1种结果;如果s⋅m≡2(mod4) , 此时s⋅m2 为奇数,则由引理2可得η(3,qs)0=−1−2qs23,η(3,qs)1=η(3,qs)1=−1+2qs23 ,将其代入表3,可得到推论2中的另一种结果。 证毕类似于推论2的证明方法,可得到以下3个结论。
推论3 设
N 是qs−1 的正因子,N1=gcd(qs−1q−1,N)=3,p≡2(mod3) 。若β0,β1 满足定理2中的情形(4)且s⋅m≡0(mod4) ,则式(1)定义的不可约拟循环码是参数为[2(qs−1)N,s,2(q−1)(qs−qs2)Nq] 的2-重量线性码,其重量分布为1+2(qs−1)3x(q−1)(2qs+qs2)Nq+qs−13x2(q−1)(qs−qs2)Nq。 若
s⋅m≡2(mod4) ,则式(1)定义的不可约拟循环码是参数为[2(qs−1)N,s,(q−1)(2qs−qs2)Nq] 的2-重量线性码,其重量分布为1+2(qs−1)3x(q−1)(2qs−qs2)Nq+2(qs−1)3x2(q−1)(qs+qs2)Nq。 推论4 设
N 是qs−1 的正因子,N1=gcd(qs−1q−1,N)=3,p≡1(mod3) 。若β0,β1 满足定理2中的情形(3)且s⋅m≡0(mod3) ,则式(1)定义的不可约拟循环码是一个参数为[2(qs−1)N,s] 的3-重量线性码,其重量分布为1+qs−13x2(q−1)(qs−c1qs3)Nq+qs−13x2(q−1)(qs+12(c1+9d1)qs3)Nq+qs−13x2(q−1)(qs+12(c1−9d1)qs3)Nq。 推论5 设
N 是qs−1 的正因子,N1=gcd(qs−1q−1,N)=3,p≡1(mod3) 。若β0,β1 满足定理2中的情形(4)且s⋅m≡0(mod3) ,则式(1)定义的不可约拟循环码是参数为[2(qs−1)N,s] 的3-重量线性码,其重量分布为1+qs−13x(q−1)(2qs+12(9d1−c1)qs3)Nq+qs−13x(q−1)(2qs−12(c1+9d1)qs3)Nq+qs−13x(q−1)(2qs+c1qs3)Nq。 类似定理1,当
N1=gcd(qs−1q−1,N)=4 时,式(1)定义的不可约拟循环码的重量分布由以下定理给出。定理3 设
N 是qs−1 的正因数使得N1=gcd(qs−1q−1,N)=4 ,则情形(5)如果
β0,β1 在相同的分圆类S 中,则式(1)定义的不可约拟循环码的重量分布如表5所示。表 5 情形(5):不可约拟循环码的重量分布重量(i) 频数(Ai) 0 1 2Nq(q−1)(qs−1−4η(4,qs)0) qs−14 2Nq(q−1)(qs−1−4η(4,qs)1) qs−14 2Nq(q−1)(qs−1−4η(4,qs)2) qs−14 2Nq(q−1)(qs−1−4η(4,qs)3) qs−14 情形(6)如果下列条件之一成立,则式(1)定义的不可约拟循环码的重量分布如表6所示。
表 6 情形(6):不可约拟循环码的重量分布重量(i) 频数(Ai) 0 1 2Nq(q−1)[(qs−1)−2(η(4,qs)0+η(4,qs)1)] qs−14 2Nq(q−1)[(qs−1)−2(η(4,qs)1+η(4,qs)2)] qs−14 2Nq(q−1)[(qs−1)−2(η(4,qs)2+η(4,qs)3)] qs−14 2Nq(q−1)[(qs−1)−2(η(4,qs)3+η(4,qs)0)] qs−14 (a)
β0,β1∈{C(4,qs)0∖C(N,qs)0,C(4,qs)1} ;(b)
β0,β1∈{C(4,qs)0∖C(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) 0 1 2Nq(q−1)[(qs−1)−2(η(4,qs)0+η(4,qs)2)] qs−12 2Nq(q−1)[(qs−1)−2(η(4,qs)1+η(4,qs)3)] qs−12 (e)
β0,β1∈{C(4,qs)0∖C(N,qs)0,C(4,qs)2} ;(f)
β0,β1∈{C(4,qs)1,C(4,qs)3} 。由引理3、推论2和定理3,可得到以下推论。
推论6 设
N 是qs−1 的正因子,N1=gcd(qs−1q−1,N)=4,p≡3(mod4) 。若β0,β1 满足定理3中的情形(5)且s⋅m≡0(mod4) ,则式(1)定义的不可约拟循环码是参数为[2(qs−1)N,s,2(q−1)(qs−qs2)Nq] 的2-重量线性码,其重量分布为1+3(qs−1)4x2(q−1)(qs−qs2)Nq+qs−14x2(q−1)(qs+3qs2)Nq。 若
s⋅m≡2(mod4) ,则式(1)定义的不可约拟循环码是参数为[2(qs−1)N,s,2(q−1)(qs−3qs2)Nq] 的2-重量线性码,其重量分布为1+3(qs−1)4x2(q−1)(qs+qs2)Nq+qs−14x2(q−1)(qs−3qs2)Nq。 推论7 设
N 是qs−1 的正因子,N1=gcd(qs−1q−1,N)=4,p≡3(mod4) 。如果β0,β1 满足定理3中的情形(6)或情形(7),则式(1)定义的指标为2的不可约拟循环码是参数为[2(qs−1)N,s,2(q−1)(qs−qs2)Nq] 的2-重量线性码,其重量分布为1+qs−12x2(q−1)(qs−qs2)Nq+qs−12x2(q−1)(qs+qs2)Nq 。推论8 设
N 是qs−1 的正因子,N1=gcd(qs−1q−1,N)=4,p≡1(mod4) 。若β0,β1 满足定理3中的情形(7)且s⋅m≡0(mod4) ,则式(1)定义的指标为2的不可约拟循环码是参数为[2(qs−1)N,s] 的2-重量线性码,其重量分布为1+qs−12x2(q−1)(qs+(u1+2v1)qs4)Nq+qs−12x2(q−1)(qs−(u1+2v1)qs4)Nq。 4. 最优码的构造
定理4 设
N 是qs−1 的正因子,N1=3,p≡2(mod3),s⋅m≡2(mod4) 。若β0,β1 满足定理2中的情形(4),则有(1)如果
p=2,s=2,N=3(2m−1) ,则式(1)定义的指标为2的不可约拟循环码是一类2-重量MDS码,其参数为[23(2m+1),2,13(2m+1−1)] 。(2)如果
s=2,2N=3(pm−1) ,则式(1)定义的指标为2的不可约拟循环码是一类达到Griesmer界的最优码,其参数为[43(pm+1),2,23(2pm−1)] 。证明 将情形(1)中所给的条件代入推论4, 可得式(1)定义的指标为2的不可约拟循环码是一类2-重量线性码且其参数满足
n=k+d−1 。 由引理5可得到一类2-重MDS码,其参数为[23(2m+1),2,13(2m+1−1)] 。将情形(2)中所给的条件代入推论4,可得到式(1)定义的指标为2的不可约拟循环码是一类2-重量线性码且参数满足
n=∑k−1i=0⌈dqi⌉ 。由引理6可得到一类达到Griesmer界的最优码,其参数为[43(pm+1),2,23(2pm−1)] 。 证毕例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 设
N 是qs−1 的正因子,N1=4,p≡3(mod4) 。若β0,β1 满足定理3中的情形(6)或情形(7),s=2,N=4(pm−1) ,则式(1)定义的指标为2的不可约拟循环码是一类2-重量MDS码,其参数为[12(pm+1),2,12(pm−1)] 。例3 令
(N1,p,m,s,N)=(4,3,3,2,104) 。由定理5,式(1)定义的指标为2的不可约拟循环码是F27 上参数为[14,2,13] 的2-重量MDS码,其重量分布为1+364x13+364x14 。5. 结论
本文研究了有限域上指标为2的不可约拟循环码的重量分布,构造了几类2-重量线性码和3-重量线性码。特别地,本文得到了3类最优的2-重量线性码。如何利用
ZpZp[v] -加性码[18,19]的重量分布构造有限域上的最优码是一个有意义的研究问题。 -
表 1 只存在时间欺骗时,不同可信卫星数目下,算法可检测出的接收机钟差欺骗距离的理论值(s)
可信卫星数目M 1 2 3 算法可检测出的接收机钟差欺骗 ≥7.7×10−8 ≥5.9×10−8 ≥5.0×10−8 表 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° -
[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)
-