CN112051550A - 一种用于gps信号的外辐射源直达波抑制方法 - Google Patents
一种用于gps信号的外辐射源直达波抑制方法 Download PDFInfo
- Publication number
- CN112051550A CN112051550A CN202010694345.7A CN202010694345A CN112051550A CN 112051550 A CN112051550 A CN 112051550A CN 202010694345 A CN202010694345 A CN 202010694345A CN 112051550 A CN112051550 A CN 112051550A
- Authority
- CN
- China
- Prior art keywords
- signal
- signals
- direct wave
- gps
- code
- 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
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种用于GPS外辐射源***的直达波干扰抑制方法,属于信号处理领域。在获得直达波信号的捕获跟踪参数之后,首先根据GPS信号的特征构建初始GPS信号模型,在完成载波相位同步和码相位同步之后,选取时长为100ms的主通道信号,并构建100ms的初始GPS直达波信号,对两者进行预处理后,两者相除得到幅度相位比例差,最后重构出完整的主通道直达波干扰信号,完成对主通道内直达波的抑制消除。本发明抑制了主通道内的直达波信号,信号功率较小的目标回波信号信杂比得以提升,不会被淹没在其他干扰信号之中,为后续目标检测打下基础。
Description
技术领域
本发明属于信号处理领域,涉及一种外辐射源直达波抑制处理方法,具体涉及一种基于GPS信号的外辐射源直达波抑制处理方法。
背景技术
外辐射源雷达,也称作无源雷达,是一种不主动发射电磁波信号,而接收已有的第三方非合作辐射源的回波来工作的雷达***。这种雷达***利用广播、电视、卫星等民用机会照射源,根据对导弹、飞机等目标的反射信号进行微弱目标相干检测、传输等信号处理,估计目标反射信号的到达方向、多普勒频移以及反射信号与直接信号的时间延迟等参数,并以此来识别跟踪目标。相对于有源雷达收发同站的局限性,无源雷达***具有接收机隐蔽性好、抗隐身能力强、***成本低、可用信号源种类丰富等一系列优点。
以GPS信号为外辐射源的无源雷达在继承上述优点的基础上,又具有其独特的优势,主要表现为:1)GPS卫星***数量众多,可实现近地空间全天候全覆盖,有利于选择合适的单个或多个卫星作为信号源;2)信号传输路径大部分位于自由空间,受地面杂波等影响小;3)由于GPS信号本身具有一定的信息传递功能,便于解决双基地雷达***中的同步问题。
无源雷达***主要有两个通道:一个是主通道,接收目标回波信号;另一个则是辅助通道接收参考信号,即直达波通道。而GPS卫星作为连续波雷达,它的信号是以连续的方式而非脉冲方式输出的。那么在接收机接收目标回波的同时,还会接收到直达波信号。一般来说,直达波信号的功率较强,而目标回波的强度会远远弱于直达波信号,这样就会对目标的检测造成干扰,即直达波通道干扰(Direct Path Interference)。因此,如何抑制目标回波通道内的直达波,提高目标回波信号的信噪比,使得无源雷达在目标探测中发挥出更大的作用,也是无源雷达研究的热点之一。
发明内容
针对上述问题,本发明根据GPS信号的特征,提出一种用于外辐射源雷达的直达波抑制方法,主要侧重于以GPS作为信号源的无源雷达***,解决了外辐射源雷达中的直达波干扰问题,能够完成对主通道内直达波的抑制消除,提升待检测目标信号信杂比,为后续信号处理打下基础。
本发明一种用于GPS信号的外辐射源直达波抑制方法,具体步骤如下:
步骤一:读入辅助通道内GPS信号的捕获跟踪参数以及发射机与接收机参数。
步骤二:直达波初始信号构建。
a、构建GPS信号。
b、GPS信号调制,得到初始信号,并播发出去。
步骤三:C/A码相位同步。
步骤四:载波相位同步。
步骤五:直达波幅度估计。
1)、截取一段实际测量的主通道信号,进行中频解调、C/A码解调(距离压缩)、距离徙动矫正和数据码解调。
2)、构建相同时长的直达波信号,进行中频解调、C/A码解调(距离压缩)、距离徙动矫正和数据码解调。
3)、将构建的直达波信号与实际测量的主通道信号,取二者的中心列相除并计算平均值,即为两者的幅度相位比例差。
步骤六:根据前述步骤一~四方法,以及需要处理的主通道信号时长,构建同等时长的直达波信号,并乘以步骤五中得到的幅度相位比例差后,与主通道信号相减消除,完成直达波抑制消除。
本发明的优点在于:
1、效率高。本发明针对GPS信号结构进行重构,无需对整个辅助通道信号进行预处理,大大减少待处理信号的数据量,提高处理效率。
2、适用范围广。对于有多个卫星直达波信号干扰的场景,本发明可以一一重构每个卫星的直达波信号并消除主通道内的干扰信号,实现干扰抑制。本发明不仅适用于单个卫星***,也适用于多星***。
附图说明
图1是本发明基于信号重构的GPS直达波干扰抑制方法流程图。
图2是主通道信号的预处理过程。
图3是主通道信号抑制前的时域处理结果。
图4是主通道信号抑制后的时域处理结果。
图5是主通道信号抑制前的频域处理结果。
图6是主通道信号抑制后的频域处理结果。
图7是矩阵中心列(直达波位置)抑制前信号功率曲线。
图8是矩阵中心列(直达波位置)抑制后信号功率曲线。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明用于GPS信号的外辐射源直达波抑制方法,如图1所示,具体步骤如下:
步骤一:读入辅助通道内GPS信号的捕获跟踪参数以及发射机与接收机参数。
发射机与接收机参数包括卫星发射信号载波L1、频率f1、接收机本地振荡器频率fIF、接收机采样率fs。
步骤二:直达波初始信号构建。
a、构建GPS信号。
GPS信号从结构上分为三个层次:信号载波、C/A码和数据码。其中,信号载波和数据码通过捕获跟踪得到的参数可以读取,C/A码需要根据所选卫星的序号生成其独有的C/A码。如:选取PRN为30的卫星,通过两个10级反馈移位寄存器组合产生该卫星特有的C/A码C。
b、GPS信号调制。
首先,对数据码与C/A码进行扩频调制,得到调制码;
其次,将调制码与正弦载波信号进行BPSK调制,得到初始信号:
其中,C(ttrans)为ttrans时刻的C/A码;D(ttrans)为ttrans时刻数据码;为ttrans时刻的载波相位;ttrans为GPS卫星发射时间,卫星发射时间与接收机接收时间存在一定的时间延迟τ;j在MATLAB中等价于虚数单位。
最后,将初始信号播发出去。
步骤三:C/A码相位同步;
码相位差不仅包括接收信号首位码片位置差,还包括由于GPS卫星与接收机之间的相对运动,引起的时间延迟,导致每时刻码相位差。具体方法如下:
B、:由于GPS卫星与接收机之间存在相对运动,会产生距离徙动现象。通过计算距离徙动RCM可以得到所需的时间延迟τ。也就是说,考虑时间延迟τ之后,接收机接收到的信号为:
其中,t为接收机接收时间;t与GPS卫星发射时间ttrans的关系为:
t=ttrans-τ
计算由于卫星和接收机相对运动引起的多普勒中心频率fd和多普勒调频率fr:
fd=[fcarr(t)-fIF]|t=0
其中,fIF为接收机的本地振荡器频率。
根据公式计算距离徙动RCM(t):
其中,fd和fr为因卫星和接收机相对运动引起的多普勒中心频率和多普勒调频率,λ为GPS卫星发射信号的载波L1波长。
计算时间延迟τ(t):
其中,RCM(t)为距离徙动,c为光速。
步骤四:载波相位同步;
其中,t为接收机接收时间;fcen为为接收机接收信号所用时间的中心时刻载波频率。
步骤五:直达波幅度估计
完成步骤一至四后,得到一段重构的直达波信号,与主通道内干扰的直达波主要存在幅度和相位两方面的差别,即:
其中,s(t)为最终需要获得的直达波信号,幅度A和相位θ即为所求。
1)通道信号预处理;
为完成幅度与相位估计,需要先截取一小段实际测量的主通道信号,对其进行一系列的预处理,包括:中频解调、C/A码解调(距离压缩)、距离徙动矫正和数据码解调这五个步骤,如图3所示。
在本实施例中,截取100ms实际测量的主通道信号sr(t),进行以下预处理:
首先是中频解调,得到的解调信号sIF(t)为:
sIF(t)=sr(t)exp(-j2πfcentfs)
其中,tfs是以0时刻为中心、长度为100ms、采样率为fs的采样时间点。
其次是C/A码解调,也就是距离压缩。将信号sIF(t)以1ms为一行的标准变换为二维信号,并在其距离向(也就是横向)进行FFT傅里叶变换到频域。将步骤a中得到的C/A码同样进行FFT傅里叶变换到频域,并把两个频域信号相乘,然后把得到的信号通过IFFT傅里叶反变换到时域,得到的时域信号sc(t)为
sc(t)=IFFT{FFT{sIF(t)}×FFT{C(t)}}
接着是距离徙动矫正。将信号sc(t)变换到频域乘上距离徙动的补偿函数,得到经过距离徙动矫正之后的信号srcmc(t)为:
其中,frange为距离向频率系数。
最后是数据码解调,解调之后的信号sD(t)为:
sD(t)=srcmc(t)×D(t)
2)根据前述步骤一~四方法构建相同时长100ms的直达波信号,进行步骤501中的预处理。
3)由于直达波经过距离压缩之后聚焦于信号矩阵的中心列,那么重构的直达波信号与实际测量的主通道信号,取二者的中心列相除并计算平均值,即为两者的幅度相位比例差。
步骤六:根据前述步骤一~四方法,以及需要处理的主通道信号时长,构建同等时长的直达波信号,并乘以步骤五中得到的幅度相位比例差后,与主通道信号相减消除,完成直达波抑制消除。
本实施例的直达波抑制对比结果如图3~图6所示,图3和图4是对消前后的时域结果,图5和图6是对消前后的频域结果。从时域来看,图4中经过距离压缩后位于矩阵中心列的直达波信号,在图4中被消除;从频域来看,如图5和图6所示,位于频域中心的直达波干扰信号得到抑制。令抑制比D为:
其中,P1为未经过抑制的直达波信号功率强度,而P2为经过抑制后的直达波信号功率强度。在本实施例中,通过计算可得抑制比D为34.1dB,具体结果如图6、图7所示。
通过本发明方法抑制了主通道内的直达波信号,信号功率较小的目标回波信号信杂比得以提升,不会被淹没在其他干扰信号之中,为后续目标检测打下基础。
Claims (5)
1.一种用于GPS信号的外辐射源直达波抑制方法,其特征在于:具体步骤如下:
步骤一:读入辅助通道内GPS信号的捕获跟踪参数以及发射机与接收机参数;
步骤二:直达波初始信号构建;
a、构建GPS信号;
b、GPS信号调制,得到初始信号,并播发出去;
步骤三:C/A码相位同步;
步骤四:载波相位同步;
步骤五:直达波幅度估计;
1)、截取一段实际测量的主通道信号,进行中频解调、C/A码解调(距离压缩)、距离徙动矫正和数据码解调。
2)、构建相同时长的直达波信号,进行中频解调、C/A码解调(距离压缩)、距离徙动矫正和数据码解调;
3)、将构建的直达波信号与实际测量的主通道信号,取二者的中心列相除并计算平均值,即为两者的幅度相位比例差。
步骤六:根据前述步骤一~四方法,以及需要处理的主通道信号时长,构建同等时长的直达波信号,并乘以步骤五中得到的幅度相位比例差后,与主通道信号相减消除,完成直达波抑制消除。
2.如权利要求1所述一种用于GPS信号的外辐射源直达波抑制方法,其特征在于:步骤二中,GPS信号调制方法为:对数据码与C/A码进行扩频调制,得到调制码,再将调制码与正弦载波信号进行BPSK调制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010694345.7A CN112051550B (zh) | 2020-07-17 | 2020-07-17 | 一种用于gps信号的外辐射源直达波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010694345.7A CN112051550B (zh) | 2020-07-17 | 2020-07-17 | 一种用于gps信号的外辐射源直达波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112051550A true CN112051550A (zh) | 2020-12-08 |
CN112051550B CN112051550B (zh) | 2022-08-30 |
Family
ID=73601902
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010694345.7A Active CN112051550B (zh) | 2020-07-17 | 2020-07-17 | 一种用于gps信号的外辐射源直达波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112051550B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114152919A (zh) * | 2021-12-10 | 2022-03-08 | 哈尔滨工业大学 | 一种收发分置短波超视距雷达时域直达波抑制方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6693580B1 (en) * | 2002-09-04 | 2004-02-17 | Northrop Grumman Corporation | Multifunction millimeter-wave system for radar, communications, IFF and surveillance |
CN102033227A (zh) * | 2010-11-30 | 2011-04-27 | 哈尔滨工程大学 | 以gps导航卫星为外辐射源的无源雷达微弱目标检测方法 |
CN103399301A (zh) * | 2013-07-01 | 2013-11-20 | 北京航空航天大学 | 一种宽带sar信号的接收装置及接收方法 |
CN105866750A (zh) * | 2016-05-23 | 2016-08-17 | 西安电子科技大学 | 一种多个gps卫星微弱回波信号检测的方法 |
CN108896976A (zh) * | 2018-07-05 | 2018-11-27 | 电子科技大学 | 一种用于多频道外辐射源雷达的相参处理方法 |
CN110286354A (zh) * | 2019-04-23 | 2019-09-27 | 中国人民解放军63921部队 | 多目标检测和区分的方法、装置和计算机可读存储介质 |
CN110376563A (zh) * | 2018-07-27 | 2019-10-25 | 西安电子科技大学 | 基于gps辐射源的雷达运动目标检测方法 |
CN110632594A (zh) * | 2019-09-18 | 2019-12-31 | 北京航空航天大学 | 一种长波长星载sar成像方法 |
-
2020
- 2020-07-17 CN CN202010694345.7A patent/CN112051550B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6693580B1 (en) * | 2002-09-04 | 2004-02-17 | Northrop Grumman Corporation | Multifunction millimeter-wave system for radar, communications, IFF and surveillance |
CN102033227A (zh) * | 2010-11-30 | 2011-04-27 | 哈尔滨工程大学 | 以gps导航卫星为外辐射源的无源雷达微弱目标检测方法 |
CN103399301A (zh) * | 2013-07-01 | 2013-11-20 | 北京航空航天大学 | 一种宽带sar信号的接收装置及接收方法 |
CN105866750A (zh) * | 2016-05-23 | 2016-08-17 | 西安电子科技大学 | 一种多个gps卫星微弱回波信号检测的方法 |
CN108896976A (zh) * | 2018-07-05 | 2018-11-27 | 电子科技大学 | 一种用于多频道外辐射源雷达的相参处理方法 |
CN110376563A (zh) * | 2018-07-27 | 2019-10-25 | 西安电子科技大学 | 基于gps辐射源的雷达运动目标检测方法 |
CN110286354A (zh) * | 2019-04-23 | 2019-09-27 | 中国人民解放军63921部队 | 多目标检测和区分的方法、装置和计算机可读存储介质 |
CN110632594A (zh) * | 2019-09-18 | 2019-12-31 | 北京航空航天大学 | 一种长波长星载sar成像方法 |
Non-Patent Citations (2)
Title |
---|
HONG-YANZHAO 等: "Target detection based on F-test in passive multistatic radar", 《DIGITAL SIGNALPROCESSING》 * |
鲁振兴: "双基地噪声雷达中直达波干扰抑制技术研究", 《电磁场与微波》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114152919A (zh) * | 2021-12-10 | 2022-03-08 | 哈尔滨工业大学 | 一种收发分置短波超视距雷达时域直达波抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112051550B (zh) | 2022-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rao et al. | Axis rotation MTD algorithm for weak target detection | |
Wu et al. | Translational motion compensation in ISAR image processing | |
Li et al. | Robust autofocus algorithm for ISAR imaging of moving targets | |
US7151483B2 (en) | System and method for concurrent operation of multiple radar or active sonar systems on a common frequency | |
Li et al. | BeiDou-based passive multistatic radar maritime moving target detection technique via space–time hybrid integration processing | |
US20150145716A1 (en) | Radar using hermetic transforms | |
US6087981A (en) | Method for pulse compression with a stepped frequency waveform | |
He et al. | A hybrid integration method for moving target detection with GNSS-based passive radar | |
Li et al. | Passive radar array with low-power satellite illuminators based on fractional Fourier transform | |
Fang et al. | Improved passive SAR imaging with DVB-T transmissions | |
Sun et al. | Effective moving target deception jamming against multichannel SAR-GMTI based on multiple jammers | |
Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
CN110133646B (zh) | 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法 | |
CN112255596B (zh) | 一种基于空间频率插值的星载sar欺骗干扰信号生成方法 | |
CN109143235A (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN112051550B (zh) | 一种用于gps信号的外辐射源直达波抑制方法 | |
CN115616629B (zh) | 一种基于天基外辐射源信号的运动目标探测补偿方法 | |
CN115508867B (zh) | Gnss-r接收机双天线信号协同相关处理***及方法 | |
He et al. | Range resolution improvement of GNSS-based passive radar via incremental Wiener filter | |
Lu et al. | Range Migration Correction Method Based on Improved Keystone Transform for GPS-based Passive Radar | |
Tang et al. | A high-efficiency partitioned integration method for GNSS-based passive radar moving target detection | |
Wang et al. | Passive radar analysis using DTMB signal | |
CN115932921B (zh) | 一种多卫星多载波的非相参联合累积处理方法 | |
Kubica | Opportunistic radar imaging using a multichannel receiver | |
Sun et al. | A low complexity long-time coherent integration method for ubiquitous radar target detection |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Wang Pengbo Inventor after: Zhou Xinkai Inventor after: Chen Jie Inventor after: Zeng Hongcheng Inventor after: Pei Zengcan Inventor before: Wang Pengbo Inventor before: Pei Zengcan Inventor before: Chen Jie Inventor before: Zeng Hongcheng Inventor before: Zhou Xinkai |
|
GR01 | Patent grant | ||
GR01 | Patent grant |