CN110954923A - 一种gnss欺骗干扰信号的检测系统 - Google Patents

一种gnss欺骗干扰信号的检测系统 Download PDF

Info

Publication number
CN110954923A
CN110954923A CN201911261820.5A CN201911261820A CN110954923A CN 110954923 A CN110954923 A CN 110954923A CN 201911261820 A CN201911261820 A CN 201911261820A CN 110954923 A CN110954923 A CN 110954923A
Authority
CN
China
Prior art keywords
signal
output
module
filtering
receiving
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911261820.5A
Other languages
English (en)
Other versions
CN110954923B (zh
Inventor
申洋赫
闫金栋
皇甫松涛
张晓松
崔永军
余速
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Spacecraft System Engineering
Original Assignee
Beijing Institute of Spacecraft System Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Spacecraft System Engineering filed Critical Beijing Institute of Spacecraft System Engineering
Priority to CN201911261820.5A priority Critical patent/CN110954923B/zh
Publication of CN110954923A publication Critical patent/CN110954923A/zh
Application granted granted Critical
Publication of CN110954923B publication Critical patent/CN110954923B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • G01S19/215—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service issues related to spoofing

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Noise Elimination (AREA)

Abstract

本发明涉及一种GNSS欺骗干扰信号的检测系统及方法，它适用于GNSS欺骗干扰信号的快速检测技术领域。该方法为一种基于抗混叠滤波、抽取的快速检测方法，可降低欺骗干扰检测运算量。本发明的信号捕获模块与现有技术相比，在输入信号、本地载波发生器、本地伪码发生器后端增加了滤波抽取模块，可明显减少DFT变换的运算量，从而有效提高信号捕获速度。本发明的信号干扰检测模块与现有技术相比，在输入信号、本地载波发生器、本地伪码发生器后端增加了滤波抽取模块，可明显减少超前、即时、滞后多路相关运算量，从而有效提高干扰检测速度。

Description

一种GNSS欺骗干扰信号的检测系统
技术领域
本发明涉及一种GNSS欺骗干扰信号的检测系统，它适用于GNSS欺骗干扰信号的快速检测技术领域。
背景技术
当前导航卫星与地面用户常采用码分多址技术(CDMA)进行通信，尽管CDMA具备一定的抗干扰能力，但由于到达地面接收机的信号功率电平较低(约低于环境噪声电平20dB)，因此接收机抗干扰能力较弱。如今各国均在稳步推进卫星导航系统(GPS、GLONASS、Galileo、北斗、日本准天顶、印度区域卫星导航系统)的建设，使民用精密定位接收机经常面临无意干扰(多址信号)的影响；特别在战争区域，GNSS干扰机配合作战，导航接收机(如弹载接收机)容易受到各种有意(阻塞信号、欺骗信号)干扰的影响，无法正常工作。
阻塞干扰信号的检测难度较低，其干信比数值较大，简单观测接收机射频前端功率的变化即可对其有效检测；欺骗干扰信号较为复杂，它在多普勒频移、载波相位、码相位等特征上与实际信号相干，较为隐蔽，观测射频前端功率的方法不再奏效；需深入分析接收机跟踪环路的相干积分值，才可对其有效检测。然而观测相干积分值要经过复杂的捕获、跟踪环节；对于某些高精度接收机，其为提高接收信号的载噪比，还需额外增加相干积分时间，使跟踪环节更加耗费时间。可见，GNSS欺骗干扰信号的检测需要较长的时间，将严重影响导航接收机的实时性。
发明内容
本发明的技术解决问题是：针对传统GNSS欺骗干扰信号检测方法耗时较长的问题，提出一种GNSS欺骗干扰信号的检测系统，该系统为一种基于抗混叠滤波、抽取的快速检测系统，可降低欺骗干扰检测运算量。
本发明的技术方案：
一种GNSS欺骗干扰信号的检测系统，该系统包括射频前端模块、信号捕获模块和信号干扰检测模块；
所述的信号捕获模块包括本地载波发生器A、本地伪码发生器A、滤波抽取模块A、滤波抽取模块B、滤波抽取模块C、乘法器A、乘法器B、乘法器C、乘法器D、DFT模块A、DFT模块B、DFT模块C、IDFT模块A、IDFT模块B和捕获处理单元；
所述的信号干扰检测模块包括本地载波发生器B、本地伪码发生器B、滤波抽取模块D、滤波抽取模块E、滤波抽取模块F、混频器A、混频器B、相关器A、相关器B、相关器C、相关器D、相关器E、相关器F、中间积分器A、中间积分器B、滞后积分器A、滞后积分器B、超前积分器A、超前积分器B、鉴相器及滤波器A、鉴相器及滤波器B和干扰检测单元；
所述的射频前端模块用于将天线接收到的模拟射频信号经下变频混频、中频信号滤波放大、模数转换后，变为数字中频信号r(k)，并输出至信号捕获模块中的滤波抽取模块A；
滤波抽取模块A用于接收射频前端模块输入的数字中频信号r(k)，滤波抽取模块A对数字中频信号r(k)进行抗混叠滤波和M倍抽取，并将抽取结果
Figure BDA0002311803780000021
输出至乘法器A、乘法器B；
本地载波发生器A用于产生同相载波信号lI(k)和正交载波信号lQ(k)，并将产生的同相载波信号lI(k)和正交载波信号lQ(k)输出给滤波抽取模块B，本地载波发生器A还用于将产生的同相载波信号lI(k)和正交载波信号lQ(k)的载波频率发送给捕获处理单元；本地载波发生器A还用于接收捕获处理单元输出的相关结果；
滤波抽取模块B用于接收本地载波发生器A输出的同相载波信号lI(k)，对同相载波信号lI(k)进行抗混叠滤波和M倍抽取，输出抽取后的载波信号
Figure BDA0002311803780000031
至乘法器B，滤波抽取模块B同时接收本地载波发生器A输出的正交载波信号lQ(k)，对正交载波信号lQ(k)进行抗混叠滤波和M倍抽取，输出抽取后的载波信号
Figure BDA0002311803780000032
至乘法器A；
本地伪码发生器A用于产生伪随机码信号sL(k)，并将产生的伪随机码信号sL(k)输出给滤波抽取模块C，本地伪码发生器A还用于将产生的伪随机码信号sL(k)的伪码相位发送给捕获处理单元，本地伪码发生器A还用于接收捕获处理单元输出的相关结果；
滤波抽取模块C用于接收本地伪码发生器A产生的伪随机码信号sL(k)，并对伪随机码信号sL(k)进行抗混叠滤波和M倍抽取，输出抽取后的伪随机码信号
Figure BDA0002311803780000033
至DFT模块C；
乘法器A用于接收滤波抽取模块A输出的数字中频信号
Figure BDA0002311803780000034
和滤波抽取模块B输出的本地正交载波信号
Figure BDA0002311803780000035
将数字中频信号
Figure BDA0002311803780000036
和本地正交载波信号
Figure BDA0002311803780000037
相乘得到MulA，并将MulA输出至DFT模块A；
乘法器B用于接收滤波抽取模块A输出的数字中频信号
Figure BDA0002311803780000038
和滤波抽取模块B输出的本地同相载波信号
Figure BDA0002311803780000039
将数字中频信号
Figure BDA00023118037800000310
和本地同相载波信号
Figure BDA00023118037800000311
相乘得到MulB，并将MulB输出至DFT模块B；
DFT模块A用于接收乘法器A输出的信号MulA，对信号MulA进行频域DFT变换得到fA，输出至乘法器C；
DFT模块B用于接收乘法器B输出的信号MulB，对信号MulB进行频域DFT变换得到fB，输出至乘法器D；
DFT模块C用于接收滤波抽取模块C的输出信号
Figure BDA0002311803780000041
先对信号
Figure BDA0002311803780000042
进行频域DFT变换，随后进行共轭变换，得到fC，输出至乘法器C和乘法器D；
乘法器C用于分别接收DFT模块A输出的信号fA和DFT模块C输出的信号fC，将fA和fC相乘得到MulC，输出至IDFT模块A；
乘法器D用于接收DFT模块B输出的信号fB和DFT模块C输出的信号fC，将fB和fC相乘得到MulD，输出至IDFT模块B；
IDFT模块A用于接收乘法器C的输出信号MulC，对MulC进行频域IDFT变换，并将变换结果A输出至捕获处理单元；
IDFT模块B用于接收乘法器D的输出信号MulD，对MulD进行频域IDFT变换，并将变换结果B输出至捕获处理单元；
捕获处理单元用于接收IDFT模块A输出的变换结果A和IDFT模块B输出的变换结果B，还用于接收本地伪码发生器A发送的sL(k)的伪码相位和本地载波发生器A发送的载波频率，并对接收到的变换结果A和变换结果B的平方和取均方根后得到相关结果，将相关结果反馈给本地伪码发生器A和本地载波发生器A，在设定时间内，取相关结果的最大值对应的本地伪码发生器A的伪码相位
Figure BDA0002311803780000043
和本地载波发生器A的载波频率
Figure BDA0002311803780000044
并将本地伪码发生器A的伪码相位
Figure BDA0002311803780000045
输出至本地伪码发生器B，将本地载波发生器A的载波频率
Figure BDA0002311803780000046
输出至本地载波发生器B；
本地载波发生器B用于接收捕获处理单元输出的载波频率
Figure BDA0002311803780000047
作为本地载波信号频率的初始量，并根据鉴相器及滤波器A输出的载波相位差异作为每个迭代周期中载波信号频率和相位的调整参考量产生同相载波信号LI和正交载波信号LQ，并将产生的同相载波信号LI和正交载波信号LQ输出给滤波抽取模块E；
本地伪码发生器B用于接收捕获处理单元输出的伪码相位
Figure BDA0002311803780000048
作为本地伪码信号相位的初始量，根据鉴相器及滤波器B输出的伪码相位差异作为每个迭代周期中伪码信号相位的调整参考量，产生即时伪码信号xP、超前伪码信号xE和滞后伪码信号xL，每个迭代周期将xP、xE、xL输出至滤波抽取模块F；
滤波抽取模块D用于接收射频前端模块输入的数字中频信号r(k)，对数字中频信号r(k)进行抗混叠滤波和M倍抽取，并将抽取结果
Figure BDA0002311803780000051
输出至混频器A和混频器B；
滤波抽取模块E用于接收本地载波发生器B输出的同相载波信号LI，对同相载波信号LI进行抗混叠滤波和M倍抽取，输出抽取结果
Figure BDA0002311803780000052
至混频器A，同时接收本地载波发生器B输出的正交载波信号LQ，对正交载波信号LQ进行抗混叠滤波和M倍抽取，输出抽取结果
Figure BDA0002311803780000053
至混频器B；
滤波抽取模块F用于接收本地伪码发生器B输出的即时伪码信号xP、超前伪码信号xE、滞后伪码信号xL，对xP、xE、xL三路信号进行抗混叠滤波和M倍抽取，分别得到
Figure BDA0002311803780000054
将
Figure BDA0002311803780000055
输出至相关器B和相关器E，将
Figure BDA0002311803780000056
输出至相关器A和相关器D，将
Figure BDA0002311803780000057
输出至相关器C和相关器F；
混频器A用于接收滤波抽取模块D输出的数字中频信号
Figure BDA0002311803780000058
和滤波抽取模块E输出的
Figure BDA0002311803780000059
将
Figure BDA00023118037800000510
和
Figure BDA00023118037800000511
混频得到MixA，将MixA输出至相关器A、相关器B、相关器C；
混频器B用于接收滤波抽取模块D输出的数字中频信号
Figure BDA00023118037800000512
和滤波抽取模块E输出的
Figure BDA00023118037800000513
将
Figure BDA00023118037800000514
和
Figure BDA00023118037800000515
混频得到MixB，将MixB输出至相关器D、相关器E、相关器F；
相关器A用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure BDA00023118037800000516
将MixA和
Figure BDA00023118037800000517
相乘得到CA，将CA输出至滞后积分器A；
相关器B用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure BDA00023118037800000518
将MixA和
Figure BDA0002311803780000061
相乘得到CB，将CB输出至中间积分器A；
相关器C用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure BDA0002311803780000062
将MixA和
Figure BDA0002311803780000063
相乘得到CC，将CC输出至超前积分器A；
相关器D用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure BDA0002311803780000064
将MixB和
Figure BDA0002311803780000065
相乘得到CD，将CD输出至滞后积分器B；
相关器E用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure BDA0002311803780000066
将MixB和
Figure BDA0002311803780000067
相乘得到CE，将CE输出至中间积分器B；
相关器F用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure BDA0002311803780000068
将MixB和
Figure BDA0002311803780000069
相乘得到CF，将CF输出至超前积分器B；
滞后积分器A用于接收相关器A输出的CA，在相干积分时间内对其积分，得到IL，将IL输出至鉴相器及滤波器B，当IL收敛后，将IL输出至干扰检测单元；
中间积分器A用于接收相关器B输出的CB，在相干积分时间内对其积分，得到IP，将IP输出至鉴相器及滤波器A；当IP收敛后，将IP输出至干扰检测单元；
超前积分器A用于接收相关器C输出的CC，在相干积分时间内对其积分，得到IE，将IE输出至鉴相器及滤波器B；当IE收敛后，将IE输出至干扰检测单元；；
滞后积分器B用于接收相关器D输出的CD，在相干积分时间内对其积分，得到QL，将QL输出至鉴相器及滤波器B；
中间积分器B用于接收相关器E输出的CE，在相干积分时间内对其积分，得到QP，将QP输出至鉴相器及滤波器A；
超前积分器B用于接收相关器F输出的CF，在相干积分时间内对其积分，得到QE，将QE输出至鉴相器及滤波器B；
鉴相器及滤波器A用于接收超前积分器A输出IE、超前积分器B输出的QE、滞后积分器A输出的IL、滞后积分器B输出的QL，根据IE、QE、IL、QL得到载波相位差异，并滤波，将滤波结果反馈至本地载波发生器B；
鉴相器及滤波器B用于接收中间积分器A输出的IP、中间积分器B输出的QP，根据IP、QP得到伪码相位差异，并滤波，将滤波结果反馈至本地伪码发生器B；
干扰检测单元用于接收收敛后的IP、IE、IL，并根据收敛后的IP、IE、IL计算相干值比例RTC，若RTC超过检测门限，则接收机给出干扰告警。
使用该系统进行GNSS欺骗干扰信号的检测，包括如下步骤：
(1)射频前端模块将天线接收到的模拟射频信号经下变频混频、中频信号滤波放大、模数转换后，变为数字中频信号，并输出至信号捕获模块；
(2)信号捕获模块接收到数字中频信号后将数字中频信号的载波频率和伪码相位进行估计，并将本地伪码发生器A的伪码相位
Figure BDA0002311803780000071
和本地载波发生器A的载波频率
Figure BDA0002311803780000072
提供给信号干扰检测模块；
(3)干扰检测模块对收敛后的IP、IE、IL计算相干值比例RTC，若RTC超过检测门限，则接收机给出干扰告警，则存在欺骗干扰信号。
本发明与现有技术相比的有益效果是：
(1)本发明的信号捕获模块与现有技术相比，在输入信号、本地载波发生器、本地伪码发生器后端增加了滤波抽取模块，可明显减少DFT变换的运算量，从而有效提高信号捕获速度。
(2)本发明的信号干扰检测模块与现有技术相比，在输入信号、本地载波发生器、本地伪码发生器后端增加了滤波抽取模块，可明显减少超前、即时、滞后多路相关运算量，从而有效提高干扰检测速度。
附图说明
图1为本发明的系统组成示意图；
图2为信号捕获模块的组成示意图；
图3为信号干扰检测模块的组成示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
实施例
如图1-3所示，一种GNSS欺骗干扰信号的检测系统，该系统包括射频前端模块、信号捕获模块和信号干扰检测模块；
所述的射频前端模块用于将天线接收到的模拟射频信号经下变频混频、中频信号滤波放大、模数转换后，变为数字中频信号r(k)，并输出至信号捕获模块中的滤波抽取模块A；
滤波抽取模块A用于接收射频前端模块输入的数字中频信号r(k)，滤波抽取模块A对数字中频信号r(k)进行抗混叠滤波和M倍抽取，并将抽取结果
Figure BDA0002311803780000081
输出至乘法器A、乘法器B；
本地载波发生器A用于产生同相载波信号lI(k)和正交载波信号lQ(k)，并将产生的同相载波信号lI(k)和正交载波信号lQ(k)输出给滤波抽取模块B，本地载波发生器A还用于将产生的同相载波信号lI(k)和正交载波信号lQ(k)的载波频率发送给捕获处理单元；本地载波发生器A还用于接收捕获处理单元输出的相关结果；
滤波抽取模块B用于接收本地载波发生器A输出的同相载波信号lI(k)，对同相载波信号lI(k)进行抗混叠滤波和M倍抽取，输出抽取后的载波信号
Figure BDA0002311803780000082
至乘法器B，滤波抽取模块B同时接收本地载波发生器A输出的正交载波信号lQ(k)，对正交载波信号lQ(k)进行抗混叠滤波和M倍抽取，输出抽取后的载波信号
Figure BDA0002311803780000083
至乘法器A；
本地伪码发生器A用于产生伪随机码信号sL(k)，并将产生的伪随机码信号sL(k)输出给滤波抽取模块C，本地伪码发生器A还用于将产生的伪随机码信号sL(k)的伪码相位发送给捕获处理单元，本地伪码发生器A还用于接收捕获处理单元输出的相关结果；
滤波抽取模块C用于接收本地伪码发生器A产生的伪随机码信号sL(k)，并对伪随机码信号sL(k)进行抗混叠滤波和M倍抽取，输出抽取后的伪随机码信号
Figure BDA0002311803780000091
至DFT模块C；
乘法器A用于接收滤波抽取模块A输出的数字中频信号
Figure BDA0002311803780000092
和滤波抽取模块B输出的本地正交载波信号
Figure BDA0002311803780000093
将数字中频信号
Figure BDA0002311803780000094
和本地正交载波信号
Figure BDA0002311803780000095
相乘得到MulA，并将MulA输出至DFT模块A；
乘法器B用于接收滤波抽取模块A输出的数字中频信号
Figure BDA0002311803780000096
和滤波抽取模块B输出的本地同相载波信号
Figure BDA0002311803780000097
将数字中频信号
Figure BDA0002311803780000098
和本地同相载波信号
Figure BDA0002311803780000099
相乘得到MulB，并将MulB输出至DFT模块B；
DFT模块A用于接收乘法器A输出的信号MulA，对信号MulA进行频域DFT变换得到fA，输出至乘法器C；
DFT模块B用于接收乘法器B输出的信号MulB，对信号MulB进行频域DFT变换得到fB，输出至乘法器D；
DFT模块C用于接收滤波抽取模块C的输出信号
Figure BDA00023118037800000910
先对信号
Figure BDA00023118037800000911
进行频域DFT变换，随后进行共轭变换，得到fC，输出至乘法器C和乘法器D；
乘法器C用于分别接收DFT模块A输出的信号fA和DFT模块C输出的信号fC，将fA和fC相乘得到MulC，输出至IDFT模块A；
乘法器D用于接收DFT模块B输出的信号fB和DFT模块C输出的信号fC，将fB和fC相乘得到MulD，输出至IDFT模块B；
IDFT模块A用于接收乘法器C的输出信号MulC，对MulC进行频域IDFT变换，并将变换结果A输出至捕获处理单元；
IDFT模块B用于接收乘法器D的输出信号MulD，对MulD进行频域IDFT变换，并将变换结果B输出至捕获处理单元；
捕获处理单元用于接收IDFT模块A输出的变换结果A和IDFT模块B输出的变换结果B，还用于接收本地伪码发生器A发送的sL(k)的伪码相位和本地载波发生器A发送的载波频率，并对接收到的变换结果A和变换结果B的平方和取均方根后得到相关结果，将相关结果反馈给本地伪码发生器A和本地载波发生器A，在设定时间内，取相关结果的最大值对应的本地伪码发生器A的伪码相位
Figure BDA0002311803780000101
和本地载波发生器A的载波频率
Figure BDA0002311803780000102
并将本地伪码发生器A的伪码相位
Figure BDA0002311803780000103
输出至本地伪码发生器B，将本地载波发生器A的载波频率
Figure BDA0002311803780000104
输出至本地载波发生器B；
本地载波发生器B用于接收捕获处理单元输出的载波频率
Figure BDA0002311803780000105
作为本地载波信号频率的初始量，并根据鉴相器及滤波器A输出的载波相位差异作为每个迭代周期中载波信号频率和相位的调整参考量产生同相载波信号LI和正交载波信号LQ，并将产生的同相载波信号LI和正交载波信号LQ输出给滤波抽取模块E；
本地伪码发生器B用于接收捕获处理单元输出的伪码相位
Figure BDA0002311803780000106
作为本地伪码信号相位的初始量，根据鉴相器及滤波器B输出的伪码相位差异作为每个迭代周期中伪码信号相位的调整参考量，产生即时伪码信号xP、超前伪码信号xE和滞后伪码信号xL，每个迭代周期将xP、xE、xL输出至滤波抽取模块F；
滤波抽取模块D用于接收射频前端模块输入的数字中频信号r(k)，对数字中频信号r(k)进行抗混叠滤波和M倍抽取，并将抽取结果
Figure BDA0002311803780000107
输出至混频器A和混频器B；
滤波抽取模块E用于接收本地载波发生器B输出的同相载波信号LI，对同相载波信号LI进行抗混叠滤波和M倍抽取，输出抽取结果
Figure BDA0002311803780000111
至混频器A，同时接收本地载波发生器B输出的正交载波信号LQ，对正交载波信号LQ进行抗混叠滤波和M倍抽取，输出抽取结果
Figure BDA0002311803780000112
至混频器B；
滤波抽取模块F用于接收本地伪码发生器B输出的即时伪码信号xP、超前伪码信号xE、滞后伪码信号xL，对xP、xE、xL三路信号进行抗混叠滤波和M倍抽取，分别得到
Figure BDA0002311803780000113
将
Figure BDA0002311803780000114
输出至相关器B和相关器E，将
Figure BDA0002311803780000115
输出至相关器A和相关器D，将
Figure BDA0002311803780000116
输出至相关器C和相关器F；
混频器A用于接收滤波抽取模块D输出的数字中频信号
Figure BDA0002311803780000117
和滤波抽取模块E输出的
Figure BDA0002311803780000118
将
Figure BDA0002311803780000119
和
Figure BDA00023118037800001110
混频得到MixA，将MixA输出至相关器A、相关器B、相关器C；
混频器B用于接收滤波抽取模块D输出的数字中频信号
Figure BDA00023118037800001111
和滤波抽取模块E输出的
Figure BDA00023118037800001112
将
Figure BDA00023118037800001113
和
Figure BDA00023118037800001114
混频得到MixB，将MixB输出至相关器D、相关器E、相关器F；
相关器A用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure BDA00023118037800001115
将MixA和
Figure BDA00023118037800001116
相乘得到CA，将CA输出至滞后积分器A；
相关器B用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure BDA00023118037800001117
将MixA和
Figure BDA00023118037800001118
相乘得到CB，将CB输出至中间积分器A；
相关器C用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure BDA00023118037800001119
将MixA和
Figure BDA00023118037800001120
相乘得到CC，将CC输出至超前积分器A；
相关器D用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure BDA00023118037800001121
将MixB和
Figure BDA00023118037800001122
相乘得到CD，将CD输出至滞后积分器B；
相关器E用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure BDA00023118037800001123
将MixB和
Figure BDA00023118037800001124
相乘得到CE，将CE输出至中间积分器B；
相关器F用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure BDA0002311803780000121
将MixB和
Figure BDA0002311803780000122
相乘得到CF，将CF输出至超前积分器B；
滞后积分器A用于接收相关器A输出的CA，在相干积分时间内对其积分，得到IL，将IL输出至鉴相器及滤波器B，当IL收敛后，将IL输出至干扰检测单元；
中间积分器A用于接收相关器B输出的CB，在相干积分时间内对其积分，得到IP，将IP输出至鉴相器及滤波器A；当IP收敛后，将IP输出至干扰检测单元；
超前积分器A用于接收相关器C输出的CC，在相干积分时间内对其积分，得到IE，将IE输出至鉴相器及滤波器B；当IE收敛后，将IE输出至干扰检测单元；；
滞后积分器B用于接收相关器D输出的CD，在相干积分时间内对其积分，得到QL，将QL输出至鉴相器及滤波器B；
中间积分器B用于接收相关器E输出的CE，在相干积分时间内对其积分，得到QP，将QP输出至鉴相器及滤波器A；
超前积分器B用于接收相关器F输出的CF，在相干积分时间内对其积分，得到QE，将QE输出至鉴相器及滤波器B；
鉴相器及滤波器A用于接收超前积分器A输出IE、超前积分器B输出的QE、滞后积分器A输出的IL、滞后积分器B输出的QL，根据IE、QE、IL、QL得到载波相位差异，并滤波，将滤波结果反馈至本地载波发生器B；
鉴相器及滤波器B用于接收中间积分器A输出的IP、中间积分器B输出的QP，根据IP、QP得到伪码相位差异，并滤波，将滤波结果反馈至本地伪码发生器B；
干扰检测单元用于接收收敛后的IP、IE、IL，并根据收敛后的IP、IE、IL计算相干值比例RTC，若RTC超过检测门限，则接收机给出干扰告警。

Claims (10)

1.一种GNSS欺骗干扰信号的检测系统，其特征在于：该系统包括射频前端模块、信号捕获模块和信号干扰检测模块；
所述的信号捕获模块包括本地载波发生器A、本地伪码发生器A、滤波抽取模块A、滤波抽取模块B、滤波抽取模块C、乘法器A、乘法器B、乘法器C、乘法器D、DFT模块A、DFT模块B、DFT模块C、IDFT模块A、IDFT模块B和捕获处理单元；
所述的信号干扰检测模块包括本地载波发生器B、本地伪码发生器B、滤波抽取模块D、滤波抽取模块E、滤波抽取模块F、混频器A、混频器B、相关器A、相关器B、相关器C、相关器D、相关器E、相关器F、中间积分器A、中间积分器B、滞后积分器A、滞后积分器B、超前积分器A、超前积分器B、鉴相器及滤波器A、鉴相器及滤波器B和干扰检测单元。
2.根据权利要求1所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：所述的射频前端模块用于将天线接收到的模拟射频信号经下变频混频、中频信号滤波放大、模数转换后，变为数字中频信号r(k)，并输出至信号捕获模块中的滤波抽取模块A。
3.根据权利要求2所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：滤波抽取模块A用于接收射频前端模块输入的数字中频信号r(k)，滤波抽取模块A对数字中频信号r(k)进行抗混叠滤波和M倍抽取，并将抽取结果
Figure FDA0002311803770000011
输出至乘法器A、乘法器B；
本地载波发生器A用于产生同相载波信号lI(k)和正交载波信号lQ(k)，并将产生的同相载波信号lI(k)和正交载波信号lQ(k)输出给滤波抽取模块B，本地载波发生器A还用于将产生的同相载波信号lI(k)和正交载波信号lQ(k)的载波频率发送给捕获处理单元；本地载波发生器A还用于接收捕获处理单元输出的相关结果。
4.根据权利要求2所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：滤波抽取模块B用于接收本地载波发生器A输出的同相载波信号lI(k)，对同相载波信号lI(k)进行抗混叠滤波和M倍抽取，输出抽取后的载波信号
Figure FDA0002311803770000021
至乘法器B，滤波抽取模块B同时接收本地载波发生器A输出的正交载波信号lQ(k)，对正交载波信号lQ(k)进行抗混叠滤波和M倍抽取，输出抽取后的载波信号
Figure FDA0002311803770000022
至乘法器A。
5.根据权利要求2所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：本地伪码发生器A用于产生伪随机码信号sL(k)，并将产生的伪随机码信号sL(k)输出给滤波抽取模块C，本地伪码发生器A还用于将产生的伪随机码信号sL(k)的伪码相位发送给捕获处理单元，本地伪码发生器A还用于接收捕获处理单元输出的相关结果。
6.根据权利要求2所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：滤波抽取模块C用于接收本地伪码发生器A产生的伪随机码信号sL(k)，并对伪随机码信号sL(k)进行抗混叠滤波和M倍抽取，输出抽取后的伪随机码信号
Figure FDA0002311803770000023
至DFT模块C。
7.根据权利要求2所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：乘法器A用于接收滤波抽取模块A输出的数字中频信号
Figure FDA0002311803770000024
和滤波抽取模块B输出的本地正交载波信号
Figure FDA0002311803770000025
将数字中频信号
Figure FDA0002311803770000026
和本地正交载波信号
Figure FDA0002311803770000027
相乘得到MulA，并将MulA输出至DFT模块A；
乘法器B用于接收滤波抽取模块A输出的数字中频信号
Figure FDA0002311803770000028
和滤波抽取模块B输出的本地同相载波信号
Figure FDA0002311803770000029
将数字中频信号
Figure FDA00023118037700000210
和本地同相载波信号
Figure FDA0002311803770000031
相乘得到MulB，并将MulB输出至DFT模块B。
8.根据权利要求2所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：DFT模块A用于接收乘法器A输出的信号MulA，对信号MulA进行频域DFT变换得到fA，输出至乘法器C；
DFT模块B用于接收乘法器B输出的信号MulB，对信号MulB进行频域DFT变换得到fB，输出至乘法器D；
DFT模块C用于接收滤波抽取模块C的输出信号
Figure FDA0002311803770000032
先对信号
Figure FDA0002311803770000033
进行频域DFT变换，随后进行共轭变换，得到fC，输出至乘法器C和乘法器D。
9.根据权利要求2所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：乘法器C用于分别接收DFT模块A输出的信号fA和DFT模块C输出的信号fC，将fA和fC相乘得到MulC，输出至IDFT模块A；
乘法器D用于接收DFT模块B输出的信号fB和DFT模块C输出的信号fC，将fB和fC相乘得到MulD，输出至IDFT模块B；
IDFT模块A用于接收乘法器C的输出信号MulC，对MulC进行频域IDFT变换，并将变换结果A输出至捕获处理单元；
IDFT模块B用于接收乘法器D的输出信号MulD，对MulD进行频域IDFT变换，并将变换结果B输出至捕获处理单元；
捕获处理单元用于接收IDFT模块A输出的变换结果A和IDFT模块B输出的变换结果B，还用于接收本地伪码发生器A发送的sL(k)的伪码相位和本地载波发生器A发送的载波频率，并对接收到的变换结果A和变换结果B的平方和取均方根后得到相关结果，将相关结果反馈给本地伪码发生器A和本地载波发生器A，在设定时间内，取相关结果的最大值对应的本地伪码发生器A的伪码相位
Figure FDA0002311803770000034
和本地载波发生器A的载波频率
Figure FDA0002311803770000035
并将本地伪码发生器A的伪码相位
Figure FDA00023118037700000411
输出至本地伪码发生器B，将本地载波发生器A的载波频率
Figure FDA00023118037700000412
输出至本地载波发生器B。
10.根据权利要求2所述的一种GNSS欺骗干扰信号的检测系统，其特征在于：
本地载波发生器B用于接收捕获处理单元输出的载波频率
Figure FDA0002311803770000041
作为本地载波信号频率的初始量，并根据鉴相器及滤波器A输出的载波相位差异作为每个迭代周期中载波信号频率和相位的调整参考量产生同相载波信号LI和正交载波信号LQ，并将产生的同相载波信号LI和正交载波信号LQ输出给滤波抽取模块E；
本地伪码发生器B用于接收捕获处理单元输出的伪码相位
Figure FDA0002311803770000042
作为本地伪码信号相位的初始量，根据鉴相器及滤波器B输出的伪码相位差异作为每个迭代周期中伪码信号相位的调整参考量，产生即时伪码信号xP、超前伪码信号xE和滞后伪码信号xL，每个迭代周期将xP、xE、xL输出至滤波抽取模块F；
滤波抽取模块D用于接收射频前端模块输入的数字中频信号r(k)，对数字中频信号r(k)进行抗混叠滤波和M倍抽取，并将抽取结果
Figure FDA0002311803770000043
输出至混频器A和混频器B；
滤波抽取模块E用于接收本地载波发生器B输出的同相载波信号LI，对同相载波信号LI进行抗混叠滤波和M倍抽取，输出抽取结果
Figure FDA0002311803770000044
至混频器A，同时接收本地载波发生器B输出的正交载波信号LQ，对正交载波信号LQ进行抗混叠滤波和M倍抽取，输出抽取结果
Figure FDA0002311803770000045
至混频器B；
滤波抽取模块F用于接收本地伪码发生器B输出的即时伪码信号xP、超前伪码信号xE、滞后伪码信号xL，对xP、xE、xL三路信号进行抗混叠滤波和M倍抽取，分别得到
Figure FDA0002311803770000046
将
Figure FDA0002311803770000047
输出至相关器B和相关器E，将
Figure FDA0002311803770000048
输出至相关器A和相关器D，将
Figure FDA0002311803770000049
输出至相关器C和相关器F；
混频器A用于接收滤波抽取模块D输出的数字中频信号
Figure FDA00023118037700000410
和滤波抽取模块E输出的
Figure FDA0002311803770000051
将
Figure FDA0002311803770000052
和
Figure FDA0002311803770000053
混频得到MixA，将MixA输出至相关器A、相关器B、相关器C；
混频器B用于接收滤波抽取模块D输出的数字中频信号
Figure FDA0002311803770000054
和滤波抽取模块E输出的
Figure FDA0002311803770000055
将
Figure FDA0002311803770000056
和
Figure FDA0002311803770000057
混频得到MixB，将MixB输出至相关器D、相关器E、相关器F；
相关器A用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure FDA0002311803770000058
将MixA和
Figure FDA0002311803770000059
相乘得到CA，将CA输出至滞后积分器A；
相关器B用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure FDA00023118037700000510
将MixA和
Figure FDA00023118037700000511
相乘得到CB，将CB输出至中间积分器A；
相关器C用于接收混频器A的输出信号MixA和滤波抽取模块F输出的
Figure FDA00023118037700000512
将MixA和
Figure FDA00023118037700000513
相乘得到CC，将CC输出至超前积分器A；
相关器D用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure FDA00023118037700000514
将MixB和
Figure FDA00023118037700000515
相乘得到CD，将CD输出至滞后积分器B；
相关器E用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure FDA00023118037700000516
将MixB和
Figure FDA00023118037700000517
相乘得到CE，将CE输出至中间积分器B；
相关器F用于接收混频器B的输出信号MixB和滤波抽取模块F输出的
Figure FDA00023118037700000518
将MixB和
Figure FDA00023118037700000519
相乘得到CF，将CF输出至超前积分器B；
滞后积分器A用于接收相关器A输出的CA，在相干积分时间内对其积分，得到IL，将IL输出至鉴相器及滤波器B，当IL收敛后，将IL输出至干扰检测单元；
中间积分器A用于接收相关器B输出的CB，在相干积分时间内对其积分，得到IP，将IP输出至鉴相器及滤波器A；当IP收敛后，将IP输出至干扰检测单元；
超前积分器A用于接收相关器C输出的CC，在相干积分时间内对其积分，得到IE，将IE输出至鉴相器及滤波器B；当IE收敛后，将IE输出至干扰检测单元；；
滞后积分器B用于接收相关器D输出的CD，在相干积分时间内对其积分，得到QL，将QL输出至鉴相器及滤波器B；
中间积分器B用于接收相关器E输出的CE，在相干积分时间内对其积分，得到QP，将QP输出至鉴相器及滤波器A；
超前积分器B用于接收相关器F输出的CF，在相干积分时间内对其积分，得到QE，将QE输出至鉴相器及滤波器B；
鉴相器及滤波器A用于接收超前积分器A输出IE、超前积分器B输出的QE、滞后积分器A输出的IL、滞后积分器B输出的QL，根据IE、QE、IL、QL得到载波相位差异，并滤波，将滤波结果反馈至本地载波发生器B；
鉴相器及滤波器B用于接收中间积分器A输出的IP、中间积分器B输出的QP，根据IP、QP得到伪码相位差异，并滤波，将滤波结果反馈至本地伪码发生器B；
干扰检测单元用于接收收敛后的IP、IE、IL，并根据收敛后的IP、IE、IL计算相干值比例RTC，若RTC超过检测门限，则接收机给出干扰告警。
CN201911261820.5A 2019-12-10 2019-12-10 一种gnss欺骗干扰信号的检测系统 Active CN110954923B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911261820.5A CN110954923B (zh) 2019-12-10 2019-12-10 一种gnss欺骗干扰信号的检测系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911261820.5A CN110954923B (zh) 2019-12-10 2019-12-10 一种gnss欺骗干扰信号的检测系统

Publications (2)

Publication Number Publication Date
CN110954923A true CN110954923A (zh) 2020-04-03
CN110954923B CN110954923B (zh) 2020-11-10

Family

ID=69980718

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911261820.5A Active CN110954923B (zh) 2019-12-10 2019-12-10 一种gnss欺骗干扰信号的检测系统

Country Status (1)

Country Link
CN (1) CN110954923B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113204034A (zh) * 2021-06-15 2021-08-03 东南大学 一种基于Ratio的GNSS发转发式欺骗检测方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090128400A1 (en) * 2005-05-31 2009-05-21 Hans Hellsten Stepped frequency radar device
CN101489153A (zh) * 2009-02-27 2009-07-22 航天恒星科技有限公司 一种卫星转发器干扰信号的监测装置及监测方法
CN101833100A (zh) * 2010-03-29 2010-09-15 北京航空航天大学 一种全数字式gnss兼容导航接收机的构建方法
CN102928854A (zh) * 2012-11-23 2013-02-13 江苏东大集成电路系统工程技术有限公司 一种基于匹配滤波器的gps捕获单元设计方法
US20150264239A1 (en) * 2014-01-07 2015-09-17 Novatel Inc. Gnss jammer detection system with optical tracking and identification
CN105375937A (zh) * 2015-11-11 2016-03-02 中国电子科技集团公司第四十一研究所 一种数字中频可变带宽成形滤波器及滤波方法
CN108900209A (zh) * 2018-05-23 2018-11-27 中国电子科技集团公司第四十研究所 一种大动态范围的多通道中频接收机
CN109188469A (zh) * 2018-08-01 2019-01-11 南京航空航天大学 一种gnss信号接收机欺骗信号参数估计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090128400A1 (en) * 2005-05-31 2009-05-21 Hans Hellsten Stepped frequency radar device
CN101489153A (zh) * 2009-02-27 2009-07-22 航天恒星科技有限公司 一种卫星转发器干扰信号的监测装置及监测方法
CN101833100A (zh) * 2010-03-29 2010-09-15 北京航空航天大学 一种全数字式gnss兼容导航接收机的构建方法
CN102928854A (zh) * 2012-11-23 2013-02-13 江苏东大集成电路系统工程技术有限公司 一种基于匹配滤波器的gps捕获单元设计方法
US20150264239A1 (en) * 2014-01-07 2015-09-17 Novatel Inc. Gnss jammer detection system with optical tracking and identification
CN105375937A (zh) * 2015-11-11 2016-03-02 中国电子科技集团公司第四十一研究所 一种数字中频可变带宽成形滤波器及滤波方法
CN108900209A (zh) * 2018-05-23 2018-11-27 中国电子科技集团公司第四十研究所 一种大动态范围的多通道中频接收机
CN109188469A (zh) * 2018-08-01 2019-01-11 南京航空航天大学 一种gnss信号接收机欺骗信号参数估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SONGLIN LI ET AL.: "Studying on the power requirements of beidou deception jamming in two cases", 《2017 CHINESE AUTOMATION CONGRESS (CAC)》 *
张鑫: "卫星导航欺骗干扰信号检测技术综述", 《全球定位系统》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113204034A (zh) * 2021-06-15 2021-08-03 东南大学 一种基于Ratio的GNSS发转发式欺骗检测方法

Also Published As

Publication number Publication date
CN110954923B (zh) 2020-11-10

Similar Documents

Publication Publication Date Title
Ta et al. Partial differential postcorrelation processing for GPS L2C signal acquisition
CN104536016A (zh) 一种gnss新体制信号捕获装置及方法
CN102854516A (zh) 一种gnss接收机中的载噪比估计方法及其系统
WO2010065250A1 (en) Method and apparatus for multipath mitigation
CN114095070A (zh) 一种基于北斗卫星导航的箭体信息回传装置
US7852264B2 (en) Systems and methods for fast GNSS signals acquisition
US20070253471A1 (en) Gps Receiver Using Differential Correlation
US8477828B2 (en) Adaptive correlation for detection of a high-frequency signal
JP5202540B2 (ja) Ｂｏｃ変調無線航法信号を受信する方法および装置
CN110954923B (zh) 一种gnss欺骗干扰信号的检测系统
Savas et al. A comparative performance analysis of GPS L1 C/A, L5 acquisition and tracking stages under polar and equatorial scintillations
CN111399006A (zh) 一种高灵敏度gnss载波跟踪环路优化方法
Pany et al. On the state-of-the-art of real-time GNSS signal acquisition—A comparison of time and frequency domain methods
Shanmugam New enhanced sensitivity detection techniques for GPS L 1 C/A and modernized signal acquisition
CN105372678A (zh) 一种正弦boc调制信号的无模糊跟踪方法
US20070147483A1 (en) Bit synchronization in a communications device
CA2811830C (en) Apparatus and method
JP2007228424A (ja) 同期タイミング検出装置および受信機
CN114236578A (zh) 一种非全向单天线旋转条件下卫星导航信号跟踪方法
Hussain et al. The effects of using variable lengths for degraded signal acquisition in GPS receivers
Heinrichs et al. Galileo/GPS receiver architecture for high sensitivity acquisition
CN108732592B (zh) 基于GNU Radio的GPS软件接收机
Zhang et al. An Adaptive Method for BeiDou Dual-Frequency Joint Acquisition Based on SNR Estimation
Wang Research on BDS B1C Signal Synchronization Method
Kao et al. A DSP/FPGA design for the acquisition and tracking of GIOVE-A signals

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant