一种抗授时干扰方法与流程

技术特征：
1.一种抗授时干扰方法，其特征在于，包括以下步骤：利用授时干扰自适应检测子模型，对用时设备接收到的授时信号进行实时干扰检测；利用定时抗扰动模型，对用时设备接收到的授时干扰进行处理，其中当检测到无授时干扰时，用时设备采用无干扰定时子模型进行定时；当检测到有授时干扰时，用时设备采用有干扰定时子模型实现自主时间保持。2.根据权利要求1所述的抗授时干扰方法，其特征在于，所述定时抗扰动模型f(f
t
,f
s
)，包括：无干扰定时子模型和有干扰定时子模型，式(1)中，f
t
表示无干扰定时子模型，f
s
表示有干扰定时子模型，μ为授时干扰标识。3.根据权利要求1所述抗授时干扰方法，其特征在于，所述无干扰定时子模型f
t
：f
t
＝t
g
t
f
ꢀꢀ
(8)式(8)中，t
g
表示用时设备接收gnss授时信号进行定时解算的结果，t
f
表示用时设备接收gnss授时信号对自带的原子钟进行校频后的定时修正量。4.根据权利要求1所述抗授时干扰方法，其特征在于，所述有干扰定时子模型f
s
：f
s
＝t
c
t
b
ꢀꢀ
(9)式(9)中，t
c
表示用时设备利用自带的原子钟进行自主时间保持的定时结果，t
b
表示用时设备利用gnss以外的其他可用授时源对自带的原子钟进行校频后的定时修正量。5.根据权利要求1所述抗授时干扰方法，其特征在于，所述授时干扰自适应检测子模型f
r
(α,γ,ε)：f
r
(α,γ,ε)＝1-(1-α)(1-γ)(1-ε) (2)式(2)中，α表示gnss授时信号相关峰分布因子，γ表示gnss授时信号秒间隔突跳因子，ε表示gnss授时信号秒脉冲变速因子；当f
r
(α,γ,ε)＝0时，表示无授时干扰；当f
r
(α,γ,ε)＝1时，表示有授时干扰。6.根据权利要求5所述抗授时干扰方法，其特征在于，授时信号相关峰分布因子α：α＝1-(1-p(n))(1-p(τ)) (3)其中，p(n)表示授时信号相关峰数量变化权值，n为授时信号相关峰数量，p(τ)表示授时信号相关峰幅度变化权值，τ为授时信号相关峰幅度变化率。7.根据权利要求5所述抗授时干扰方法，其特征在于，所述构建的授时信号秒间隔突跳因子γ：γ＝1-isequal(δt,t) (6)式(6)中，isequal(δt,t)表示相等性判定，δt＝t时判定值为1，δt≠t时判定值为0；δt表示根据gnss授时信号导航电文得到的秒间隔，t表示秒间隔突跳检测周期。8.根据权利要求5所述抗授时干扰方法，其特征在于，所述授时信号秒脉冲变速因子ε：ε＝1-isless(dist(f
c
,f
p
),d) (7)式(7)中，isless(dist(f
c
,f
p
),d)表示小于性判定，dist(f
c
,f
p
)＜d时判定值为1，dist(f
c
,f
p
)≥d时判定值为0；dist(f
c
,f
p
)表示两条曲线f
c
和f
p
之间的距离，f
c
表示校准模式下用时设备输出秒脉冲与自带原子钟输出秒脉冲的偏差变化曲线，f
p
表示定时模式下用时设
备输出秒脉冲与自带原子钟输出秒脉冲的偏差变化曲线；d表示用时设备输出秒脉冲的偏移门限。9.根据权利要求5所述抗授时干扰方法，其特征在于，当且仅当gnss授时信号相关峰数量为1且gnss授时信号相关峰幅度变化率小于10％时，gnss授时信号相关峰分布因子α＝0，否则α＝1；当且仅当根据gnss授时信号导航电文得到的秒间隔δt等于秒间隔突跳检测周期t时，gnss授时信号秒间隔突跳因子γ＝0，否则γ＝1。10.根据权利要求5所述抗授时干扰方法，其特征在于，当且仅当校准模式下用时设备输出秒脉冲与自带原子钟输出秒脉冲的偏差变化曲线f
c
与定时模式下用时设备输出秒脉冲与自带原子钟输出秒脉冲的偏差变化曲线f
p
之间的距离小于用时设备输出秒脉冲的偏移门限d时，gnss授时信号秒脉冲变速因子ε＝0，否则ε＝1。

技术总结
本发明公开一种抗授时干扰方法，包括，利用授时干扰自适应检测子模型，对用时设备接收到的授时信号进行实时干扰检测；利用定时抗扰动模型，对用时设备接收到的授时干扰进行处理，其中当检测到无授时干扰时，用时设备采用无干扰定时子模型进行定时；当检测到有授时干扰时，用时设备采用有干扰定时子模型实现自主时间保持。本发明可实时检测GNSS授时干扰，并实现抗干扰定时效果。针对复杂的复合式GNSS授时干扰及极端条件下的GNSS授时中断，本发明可自适应地检测GNSS授时干扰，并实现不依赖GNSS授时的自主时间保持，支持用时设备在多场景及强博弈环境中实现可靠定时。强博弈环境中实现可靠定时。强博弈环境中实现可靠定时。


技术研发人员：郭慧杰 杨慧君 李丹丹 王学运 张升康
受保护的技术使用者：北京无线电计量测试研究所
技术研发日：2022.09.02
技术公布日：2022/12/16