CN106597409A - 一种基于外辐射源雷达的时延‑多普勒频提取方法 - Google Patents
一种基于外辐射源雷达的时延‑多普勒频提取方法 Download PDFInfo
- Publication number
- CN106597409A CN106597409A CN201611180577.0A CN201611180577A CN106597409A CN 106597409 A CN106597409 A CN 106597409A CN 201611180577 A CN201611180577 A CN 201611180577A CN 106597409 A CN106597409 A CN 106597409A
- Authority
- CN
- China
- Prior art keywords
- delay
- frequency
- signal
- doppler frequency
- doppler
- 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
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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于信号处理领域,具体涉及一种基于外辐射源雷达的时延‑多普勒频提取方法。本发明的方法基于传统的时频二维相关方法,在对目标的时延和多普勒的提取过程中,主要利用参考信号和传统方法获得的初步时延和多普勒,通过对参考信号和目标信号的补偿及相关运算,再次获取更精确的时延和多普勒。本发明的有益效果为,相对于传统方法,本发明可以明显提高雷达探测中时延和多普勒频的检测精度同时能够在达到相同检测精度时,时间复杂度更低。
Description
技术领域
本发明属于信号处理领域,具体涉及一种基于外辐射源雷达的时延-多普勒频提取方法。
背景技术
外辐射源雷达利用地面广播电台、电视台、通信台等非合作的第三方辐射信号经目标散射后的回波来获取目标信息,具有生存能力强、抗干扰性能好等特点以及探测隐身目标的潜力,并具有探测低空目标的能力。在外辐射源雷达探测过程中,时延和多普勒频是最重要的提取参数。该参数的精度直接影响着后续对目标的定位。目前在对时延和多普勒频进行参数提取的主要方法是时频二维相关。
时频二维相关采用二维相干匹配滤波提取目标的时延和多普勒频移,其原理可用模糊函数表示:
其中,ei(n)为外辐射源雷达在经过DBF(数字波束形成)后的某一路主通道信号,Sref(n)为外辐射源的参考信号,M为波束的个数。(τ,f)为目标引起的时延和多普勒频率。当存在目标时,|ψi(τd,fd)|会产生一个峰值点,该峰值点对应的(τd,fd)即为目标的时延和多普勒。具体的原理框图如附图1所示。
当外辐射源信号(参考信号)为广播调频信号可表示为:
其中v(t)为加性高斯白噪声信号,fc为调频广播信号的载频,D1为外辐射源雷达接收广播信号的时间延迟。θ(t-D1)为音频信号调制项。
上式中,第一项为外辐射源雷达接收的广播信号探测到目标的回波信号,fd为由于目标运动导致多普勒频移,D2为雷达接收天线接收到目标信号的时间延迟。m(t)是高斯白噪声信号。因此,直达波信号的延时式为:
上式中延时长度Δt用采样周期,或者采样周期的整数倍。上式的共轭式与主通道信号相乘得:
由上式可以看出,当D2=D1+nΔt时,上式中音频调制项就出现相消的情况,此时,上式将只包含多普勒频移项和噪声项,即:
对上式做FFT运算,即可得到目标的时延和多普勒。
利用上述传统的时频二维相关的方法来获取高精度的目标参数时延和多普勒频需要很大的计算量。在时间维度要提高时延测量精度就需要做更精细的扫描,在频率维度上要提高精度则需要提高采样率,而这两者都会加大时间复杂度。正因如此,本发明提出了一种新的方法,在达到相同时延和多普勒频的测量精度的基础上可以大大减少计算量。
发明内容
本发明针对现有的时频二维相关方法去提取雷达探测过程中的时延-多普勒频参数往往误差较大的问题,提出了更为精确的时延-多普勒频参数提取方法。
本发明的技术方案为:
如图2所示,一种基于外辐射源雷达的时延-多普勒频提取方法,其特征在于,包括以下步骤:
a.采用时频二维相关方法和峰值搜索初步获取时延τd1和多普勒频fd1:设参考信号为Sref(t),目标探测信号为其中n(t)为高斯白噪声;时频二维相关方法的原理如下公式1:
公式1中,ei(n)为外辐射源雷达在经过数字波束形成后的某一路主通道信号,Sref(n)为外辐射源的参考信号,M为波束的个数,(τ,f)为目标引起的时延和多普勒频率;
根据峰值搜索获取|ψi(τd1,fd1)|产生的峰值点,该峰值点对应的(τd1,fd1)即为目标的时延和多普勒;
b.对目标探测信号进行多普勒频补偿:采用步骤a中获得的多普勒频移fd1对目标探测信号进行多普勒频补偿获得目标信号S1(t)如下公式2:
c.采用参考信号对步骤b中获得的目标信号S1(t)进行时延扫描:利用参考信号Sref(t),在步骤a中获得的时延τd1附近进行时延拟合获得信号S2(t)=Sref(t-τd2),将S2(t)与S1(t)做相关运算,相关值取得最大时信号S2(t)所对应的时延τd2,为相对于时延τd1更为准确的时延值;
d.采用步骤c中获得的时延τd2,对参考信号Sref(t)进行时延补偿,获得目标信号S3(t):S3(t)=Sref(t-τd2);
e.采用参考信号对步骤d中获得的目标信号S3(t)进行多普勒频扫描:采用目标信号S3(t),在步骤a中得到的多普勒频fd1附近进行多普勒频移拟合得到信号将S4(t)与S1(t)做相关运算,相关值取得最大时S4(t)所对应的多普勒频fd2,为相对于多普勒频fd1更为准确的多普勒频移值。
本发明的有益效果为,相对于传统方法,本发明可以明显提高雷达探测中时延和多普勒频的检测精度同时能够在达到相同检测精度时,时间复杂度更低。
附图说明
图1为时频二维相关原理示意图;
图2为本发明的时延-多普勒频参数提取方案示意图;
图3为巴克码相位编码信号时域图;
图4为时频二维相关图(SNR=10dB);
图5为时频二维相关图(SNR=15dB);
图6为时频二维相关图(SNR=20dB);
图7为时频二维相关图(SNR=25dB);
图8为时频二维相关图(SNR=30dB)。
具体实施方式
下面结合附图和实施例,详细描述本发明的技术方案:
实施例:
在本实施例中,采用巴克码线性调频信号作为参考信号,该信号时域图如附图3所示,该信号的仿真参数如表1所示。仿真平台操作系统是WIN7(64bit)CPU 3.20GHz,采用MATLAB R2012a(64bit)进行编译。仿真基本参数表如下表1所示:
表1仿真基本参数表
脉宽 | 信号持续时间 | 脉冲重复间隔 | 采样率 |
1(us) | 1(ms) | 0.1(ms) | 10MHz |
本次仿真的目标时延为2.11(us),多普勒频移5455Hz
本方法和传统方法时延和多普勒频检测结果如下表2所示:
表2不同信噪比检测结果对比表
图4~图7为信噪比10~30dB的时频二维相关所得结果图。
结论分析:由上表2可知当信噪比在10~30dB时,传统的时频二维相关方法和本方法都可以检测出目标。通过数值分析可知,利用时频二维相关方法得到时延误差为0.11us,而利用本方法所得的时延误差为0.01us,在传统的仅仅利用时频二维相关的方法的精度提高了一个数量级效果明显;利用时频二维相关方法得到多普勒频的偏差为而利用本方法所得的多普勒频偏差最大为通过改变目标的时延-多普勒频参数,利用传统时频二维相关得到的结果和利用本方法得到的结果如表3所示:
表3不同目标检测结果对比表
由表3分析可知,本方法侧得目标时延与多普勒更为精确。由表4可知:
表4不同采样率检测结果及计算量分析
为了提高传统时频二维相关的时延精度,需要提高采样率,但提高采样率并不能提高提取多普勒频的精度,同时当增大采样率时,时间复杂度大大提高。通过上述对比可知本方法相对于传统方法可以明显提高雷达探测中时延和多普勒频的检测精度同时能够在达到相同检测精度时,时间复杂度更低。
Claims (1)
1.一种基于外辐射源雷达的时延-多普勒频提取方法,其特征在于,包括以下步骤:
a.采用时频二维相关方法和峰值搜索初步获取时延τd1和多普勒频fd1:设参考信号为Sref(t),目标探测信号为其中n(t)为高斯白噪声;时频二维相关方法的原理如下公式1:
公式1中,ei(n)为外辐射源雷达在经过数字波束形成后的某一路主通道信号,Sref(n)为外辐射源的参考信号,M为波束的个数,(τ,f)为目标引起的时延和多普勒频率;
根据峰值搜索获取|ψi(τd1,fd1)|产生的峰值点,该峰值点对应的(τd1,fd1)即为目标的时延和多普勒;
b.对目标探测信号进行多普勒频补偿:采用步骤a中获得的多普勒频移fd1对目标探测信号进行多普勒频补偿获得目标信号S1(t)如下公式2:
c.采用参考信号对步骤b中获得的目标信号S1(t)进行时延扫描:利用参考信号Sref(t),在步骤a中获得的时延τd1附近进行时延拟合获得信号S2(t)=Sref(t-τd2),将S2(t)与S1(t)做相关运算,相关值取得最大时信号S2(t)所对应的时延τd2,为相对于时延τd1更为准确的时延值;
d.采用步骤c中获得的时延τd2,对参考信号Sref(t)进行时延补偿,获得目标信号S3(t):S3(t)=Sref(t-τd2);
e.采用参考信号对步骤d中获得的目标信号S3(t)进行多普勒频扫描:采用目标信号S3(t),在步骤a中得到的多普勒频fd1附近进行多普勒频移拟合得到信号将S4(t)与S1(t)做相关运算,相关值取得最大时S4(t)所对应的多普勒频fd2,为相对于多普勒频fd1更为准确的多普勒频移值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611180577.0A CN106597409B (zh) | 2016-12-19 | 2016-12-19 | 一种基于外辐射源雷达的时延-多普勒频移提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611180577.0A CN106597409B (zh) | 2016-12-19 | 2016-12-19 | 一种基于外辐射源雷达的时延-多普勒频移提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106597409A true CN106597409A (zh) | 2017-04-26 |
CN106597409B CN106597409B (zh) | 2019-02-12 |
Family
ID=58601890
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611180577.0A Active CN106597409B (zh) | 2016-12-19 | 2016-12-19 | 一种基于外辐射源雷达的时延-多普勒频移提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106597409B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108680917A (zh) * | 2018-05-18 | 2018-10-19 | 华瑞哲 | 一种急弯处汽车雷达不可视运动目标检测方法及装置 |
CN110609263A (zh) * | 2019-10-29 | 2019-12-24 | 电子科技大学 | 一种同时计算脉冲激光雷达目标回波时延和频偏的方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707272A (zh) * | 2012-06-13 | 2012-10-03 | 西安电子科技大学 | 基于gpu的外辐射源雷达信号实时处理系统及处理方法 |
CN103217670A (zh) * | 2013-03-29 | 2013-07-24 | 电子科技大学 | 一种基于pca的外辐射源微弱信号检测方法 |
JP2013167599A (ja) * | 2012-02-17 | 2013-08-29 | Mitsubishi Electric Corp | パッシブレーダ装置および目標位置検出方法 |
CN103293518A (zh) * | 2013-05-29 | 2013-09-11 | 电子科技大学 | 广播信号外辐射源定位检测方法 |
EP2677343A1 (en) * | 2011-02-18 | 2013-12-25 | Mitsubishi Electric Corporation | Passive radar device |
US20150145716A1 (en) * | 2013-11-22 | 2015-05-28 | Hobbit Wave | Radar using hermetic transforms |
CN105656825A (zh) * | 2016-02-26 | 2016-06-08 | 北京交通大学 | 一种用于轨道交通高速移动场景的多普勒频移校正方法 |
CN105929389A (zh) * | 2015-12-05 | 2016-09-07 | 中国人民解放军信息工程大学 | 一种基于外辐射源时延和多普勒频率的直接定位方法 |
US9489813B1 (en) * | 2006-09-22 | 2016-11-08 | Michael L. Beigel | System for location in environment and identification tag |
-
2016
- 2016-12-19 CN CN201611180577.0A patent/CN106597409B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9489813B1 (en) * | 2006-09-22 | 2016-11-08 | Michael L. Beigel | System for location in environment and identification tag |
EP2677343A1 (en) * | 2011-02-18 | 2013-12-25 | Mitsubishi Electric Corporation | Passive radar device |
JP2013167599A (ja) * | 2012-02-17 | 2013-08-29 | Mitsubishi Electric Corp | パッシブレーダ装置および目標位置検出方法 |
CN102707272A (zh) * | 2012-06-13 | 2012-10-03 | 西安电子科技大学 | 基于gpu的外辐射源雷达信号实时处理系统及处理方法 |
CN103217670A (zh) * | 2013-03-29 | 2013-07-24 | 电子科技大学 | 一种基于pca的外辐射源微弱信号检测方法 |
CN103293518A (zh) * | 2013-05-29 | 2013-09-11 | 电子科技大学 | 广播信号外辐射源定位检测方法 |
US20150145716A1 (en) * | 2013-11-22 | 2015-05-28 | Hobbit Wave | Radar using hermetic transforms |
CN105929389A (zh) * | 2015-12-05 | 2016-09-07 | 中国人民解放军信息工程大学 | 一种基于外辐射源时延和多普勒频率的直接定位方法 |
CN105656825A (zh) * | 2016-02-26 | 2016-06-08 | 北京交通大学 | 一种用于轨道交通高速移动场景的多普勒频移校正方法 |
Non-Patent Citations (2)
Title |
---|
李万春等: ""无参考通道的外辐射源雷达参数估计"", 《系统工程与电子技术》 * |
邢翠柳等: ""基于频率补偿的时延估计方法"", 《飞行器测控学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108680917A (zh) * | 2018-05-18 | 2018-10-19 | 华瑞哲 | 一种急弯处汽车雷达不可视运动目标检测方法及装置 |
CN110609263A (zh) * | 2019-10-29 | 2019-12-24 | 电子科技大学 | 一种同时计算脉冲激光雷达目标回波时延和频偏的方法 |
CN110609263B (zh) * | 2019-10-29 | 2022-11-04 | 电子科技大学 | 一种同时计算脉冲激光雷达目标回波时延和频偏的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106597409B (zh) | 2019-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Falcone et al. | Two‐dimensional location of moving targets within local areas using WiFi‐based multistatic passive radar | |
CN108490410B (zh) | 一种两坐标雷达对海目标联合检测跟踪方法 | |
Wang et al. | Direction finding in frequency-modulated-based passive bistatic radar with a four-element adcock antenna array | |
CN103217670B (zh) | 一种基于pca的外辐射源微弱信号检测方法 | |
CN103293518B (zh) | 广播信号外辐射源定位检测方法 | |
Radmard et al. | Data fusion in MIMO DVB-T-based passive coherent location | |
CN104007421B (zh) | 基于全变差和压缩感知的罗兰c被动雷达toa估计方法 | |
AU2002314766B2 (en) | System and method for detection and feature extraction in passive coherent location applications | |
CN110261819A (zh) | 基于时延补偿的多无人机协同定位方法 | |
CN105759241A (zh) | 基于时差频差的直接定位方法 | |
CN109839623B (zh) | 一种地外天体着陆测量雷达面目标回波信号测距处理方法 | |
CN108398676A (zh) | 一种外辐射源雷达微弱运动目标探测方法 | |
CN106597409B (zh) | 一种基于外辐射源雷达的时延-多普勒频移提取方法 | |
CN111948618B (zh) | 一种基于卫星外辐射源的前向散射目标探测方法及系统 | |
Chen et al. | Suppression method for main-lobe interrupted sampling repeater jamming in distributed radar | |
KR101908455B1 (ko) | 저전력 주파수 변조 연속파 레이더 시스템 및 그 제어 방법 | |
CN108415013A (zh) | 基于调频斜率精搜索的低信噪比信号参数提取方法 | |
CN107346022B (zh) | 基于微波干涉仪的高精度测船雷达及测速方法 | |
CN102778677B (zh) | 远距脉冲压缩与近距脉冲包络检波延时测量结合测距方法 | |
CN114280532B (zh) | 一种基于带内共轭点乘的雷达目标角度估计方法及系统 | |
CN109085568B (zh) | 一种基于二次混频的调频连续波多目标检测方法 | |
CN108919206B (zh) | 一种基于子载波处理的外辐射源雷达极化滤波方法 | |
CN108718223B (zh) | 一种非合作信号的盲频谱感知方法 | |
Canini et al. | Use of Digital-Television terrestrial (DTV) signals for passive radars | |
RU2589036C1 (ru) | Радиолокатор с непрерывным шумовым сигналом и способ расширения диапазона измеряемых дальностей в радиолокаторе с непрерывным сигналом |
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 |