CN104880698B - 基于距离频域多项式相位变换的空间机动目标检测方法 - Google Patents
基于距离频域多项式相位变换的空间机动目标检测方法 Download PDFInfo
- Publication number
- CN104880698B CN104880698B CN201510308318.0A CN201510308318A CN104880698B CN 104880698 B CN104880698 B CN 104880698B CN 201510308318 A CN201510308318 A CN 201510308318A CN 104880698 B CN104880698 B CN 104880698B
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- apart
- data
- slow time
- target
- 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.)
- Expired - Fee Related
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
一种基于距离频域多项式相位变换的空间机动目标检测方法。该方法沿接收到的机动目标回波数据的快时间向进行快速傅里叶变换以得到距离频域‑慢时间数据,然后在得到的距离频域‑慢时间数据基础上利用距离频域多项式相位变换方法对机动目标回波信号出现的多普勒走动现象进行校正,并在校正基础上对回波数据进行相参积累而用于目标检测。将仿真结果与动目标检测方法、Radon‑Fourier变换、分数阶傅里叶变换和Radon‑FRFT方法的处理结果进行比较;仿真结果表明本发明方法能在相干积累脉冲数有限且低信噪比的情况下有效地检测到目标,从而验证了所提方法的有效性。
Description
技术领域
本发明属于雷达信号处理和空间机动目标检测技术领域,特别是涉及一种基于距离频域多项式相位变换的空间机动目标检测方法。
背景技术
空间机动目标主要是指在轨运行的卫星、空间站和空间碎片等。近年来随着世界各国对地球空间资源利用的不断深入,空间机动目标的数量也呈快速上升的趋势。随着空间机动目标数量的日益增多,这些空间机动目标给航天器的在轨运行和国土防御都带来了较大威胁,所以天基雷达探测和识别空间机动目标技术将会变得越来越重要。空间机动目标通常具有较大加速度,如空间碎片的加速度可达100m/s2以上,此时天基雷达的回波信号会出现明显的多普勒走动现象,并且目标回波能量分散在不同的多普勒单元中;另外为了获得较长的预警时间,天基雷达对于目标的观测时间较短,这样就会对在积累脉冲数有限的情况下有效地检测到目标带来较大的困难。
目前针对空间机动目标检测的方法主要有分数阶Fourier变换(FRFT)、Radon-FRFT方法(RFRFT)、Radon-Fourier变换方法(RFT)和动目标检测(Moving TargetDetection,MTD)方法等。其中分数阶Fourier变换方法和Radon-FRFT方法在雷达积累脉冲数有限时校正多普勒走动的性能会受到较大影响,而Radon-Fourier变换方法和MTD方法不能对多普勒走动现象进行校正,从而无法有效地积累目标能量,因此检测性能较差。总之上述这些方法均不能直接用于空间机动目标检测。
发明内容
为了解决上述问题,本发明的目的在于提供一种基于距离频域多项式相位变换的空间机动目标检测方法。
为了达到上述目的,本发明提供的基于距离频域多项式相位变换的空间机动目标检测方法包括按顺序进行的下列步骤:
1)对天基雷达接收到的空间机动目标回波数据进行脉冲压缩处理;
2)沿上述经过脉冲压缩处理后数据的快时间向进行快速傅里叶变换以得到距离频域-慢时间数据并留一个备份数据,接着沿非备份数据的慢时间向进行时延处理;
3)将上述未经过时延的距离频域-慢时间数据与经过时延的距离频域-慢时间数据的共轭进行乘积处理,然后沿乘积结果的距离频域向进行逆快速傅里叶变换处理以得到经过多普勒走动校正后的距离单元-慢时间数据;
4)沿步骤3)中得到的经过多普勒走动校正后的距离单元-慢时间数据的慢时间向进行快速傅里叶变换处理以对目标回波进行相参积累,在相参积累结果基础上完成目标检测。
在步骤3)中,所述的将未经过时延的距离频域-慢时间数据与经过时延的距离频域-慢时间数据的共轭进行乘积处理,然后沿乘积结果的距离频域向进行逆快速傅里叶变换处理以得到经过多普勒走动校正后的距离单元-慢时间数据的方法是:
在步骤2)中得到的未经过时延的距离频域-慢时间数据表达式可写为:
其中A0为回波信号复幅度,f为距离频域,tm=mTr为慢时间,c为光速,为信号波长,M为相干积累脉冲数,m=1,...,M,Tr是脉冲重复周期,fc为载波频率,为目标与雷达之间的距离,r0为目标的初始距离,v0为目标运动速度,a0为目标加速度,且有:
留一份X(f,tm)的备份数据,并记为X'(f,tm);
则经过时延的距离频域-慢时间数据表达式可写为,:
τ为延时量;
将X'(f,tm)与式(4)的共轭进行相乘处理以实现距离频域多项式相位运算:
其中“*”表示共轭处理;
沿的距离频域向进行逆快速傅里叶变换以得到经过多普勒走动校正后的距离单元-慢时间数据表达式为:
其中,t'为快时间,sinc(x)=sin(πx)/(πx)。
在步骤4)中,所述的沿步骤3)中得到的经过多普勒走动校正后的距离单元-慢时间数据的慢时间向进行快速傅里叶变换处理以对目标回波进行相参积累,在相参积累结果基础上完成目标检测的方法是:沿距离单元-慢时间数据s'r(t',tm)的慢时间向进行快速傅里叶变换处理得到用于目标检测的距离单元-多普勒数据,即实现距离频域多项式相位变换,然后在距离单元-多普勒数据基础上利用单元平均恒虚警检测原理对空间机动目标进行检测。
本发明提供的基于距离频域多项式相位变换的空间机动目标检测方法首先是将天基雷达接收到的机动目标回波数据进行脉冲压缩处理,然后沿脉冲压缩处理后数据的快时间向进行快速傅里叶变换以得到距离频域-慢时间数据,沿距离频域-慢时间数据的慢时间向进行时延处理,接着将没有经过时延和经过时延的距离频域-慢时间数据的共轭相乘,沿乘积结果的距离频域向进行逆快速傅里叶变换,再沿慢时间向进行快速傅里叶变换以对目标回波进行相参积累,最后在能量相参积累结果基础上完成目标检测。将仿真结果与常规MTD、RFT、FRFT和Radon-FRFT方法进行对比可知本文所提方法能在低信噪比和积累脉冲数有限的情况下有效地检测出目标,从而验证了所提方法的有效性。
附图说明
图1为本发明提供的基于距离频域多项式相位变换的空间机动目标检测方法流程图。
图2为经过解调和脉冲压缩处理后的天基雷达接收数据分布图。
图3为接收到的数据进行时延后的结果图。
图4为经过距离频域多项式相位变换后的结果基础上得到的距离单元-多普勒数据的结果图。
图5为在虚警概率Pfa=10-4时本发明方法、RFRFT方法、FRFT方法、RFT方法和MTD方法的恒虚警检测性能对比曲线图。
图6为在虚警概率Pfa=10-6时本发明方法、RFRFT方法、FRFT方法、RFT方法和MTD方法的恒虚警检测性能对比曲线图。
具体实施方式
下面结合附图和具体实施例对本发明提供的基于距离频域多项式相位变换的空间机动目标检测方法进行详细说明。
如图1所示,本发明提供的基于距离频域多项式相位变换的空间机动目标检测方法包括按顺序进行的下列步骤:
1)对天基雷达接收到的空间机动目标回波数据进行脉冲压缩处理;
设天基雷达接收到的未经脉冲压缩处理的距离频域-慢时间域信号形式为:
其中A0为信号复幅度,TP为脉冲信号宽度,c为光速,为调频率,fc表示信号的载波频率,f表示距离频域,tm=mTr为慢时间,且m=1,...,M,为空间机动目标与天基雷达间的距离,r0为目标的初始距离,v0为目标运动速度,a0为目标加速度,且有:
对式(1)进行匹配滤波,即乘以P(f)的共轭,然后沿匹配滤波后的距离频域-慢时间回波数据的距离频域向进行逆快速傅里叶变换后得到经过脉冲压缩处理后的回波信号表达式为:
式(4)中t'为快时间,sinc(x)=sin(πx)/(πx)。
2)沿上述经过脉冲压缩处理后数据的快时间向进行快速傅里叶变换以得到距离频域-慢时间数据并留一个备份数据,接着沿非备份数据的慢时间向进行时延处理;
在步骤2)中,对步骤1)中式(4)的快时间向进行快速傅里叶变换后的得到距离频域-慢时间数据表达式为:
留一个与距离频域-慢时间数据X(f,tm)相同的备份数据,记为X'(f,tm);沿X(f,tm)的慢时间向进行时延(时延量τ为0.01秒)后得到:
3)将上述未经过时延的距离频域-慢时间数据与经过时延的距离频域-慢时间数据的共轭进行乘积处理,然后沿乘积结果的距离频域向进行逆快速傅里叶变换处理以得到经过多普勒走动校正后的距离单元-慢时间数据;
在步骤3)中,将未经时延和经过时延的距离频域-慢时间数据的共轭进行相乘处理以实现距离频域多项式相位运算处理:
其中“*”表示共轭处理。
沿的距离频域向进行逆快速傅里叶变换以得到经过多普勒走动校正后的距离单元-慢时间数据表达式为:
4)沿步骤3)中得到的经过多普勒走动校正后的距离单元-慢时间数据的慢时间向进行快速傅里叶变换处理以对目标回波进行相参积累,在相参积累结果基础上完成目标检测。
在步骤4)中,对步骤3)中的式(8)的慢时间向进行快速傅里叶变换处理以得到对应的距离单元-多普勒数据(即实现距离频域多项式相位变换):
其中fm为慢时间频域(多普勒频域),δ(·)表示狄拉克函数。
由s'r(t',tm)的表达式可以看出,经过距离频域多项式相位变换后得到的回波信号是一个单频信号,故沿s'r(t',tm)的慢时间进行快速傅里叶变换处理可以得到用于目标检测的距离单元-多普勒数据,然后在距离单元-多普勒数据基础上利用单元平均恒虚警检测原理对空间机动目标进行检测。单元平均恒虚警检测处理的具体方法是任选一个距离单元-多普勒数据中的一个数据作为检测统计量,然后在检测统计量周围选取一定数量的保护单元,在保护单元外选取20个参考单元并将计算出噪声平均功率,然后利用计算出来的噪声平均功率计算出门限,接着比较检测统计量是否高于门限,如果高于门限,则认为目标存在,如果低于门限,则认为目标不存在;重复上述操作,直到将距离单元-多普勒数据中所有数据检测完为止。
仿真结果及分析:
本发明提供的基于距离频域多项式相位变换的空间机动目标检测方法的效果可以通过以下仿真实验进一步说明。仿真参数设置:天基雷达发射信号波长为0.3m,脉冲重复频率为1500Hz,信号带宽为30MHz,载频为1GHz。目标参数设置为:空间机动目标初始在第200个距离单元中,速度v0=50m/s,加速度a0=180m/s2,信噪比为-11dB。
1、回波数据分布图
图2为经过脉冲压缩处理后的天基雷达接收数据图。从图中可以看出目标回波信号在雷达相干处理时间内没有出现距离徙动现象。
2、接收数据进行时延处理
图3给出了对接收到的数据进行时延后的结果图,其中时延量为0.01秒。
3、在经过距离频域多项式相位变换后的结果基础上得到距离单元-多普勒数据
图4给出了在经过距离频域多项式相位变换后的结果基础上得到距离单元-多普勒数据的结果图。从图中可以看出,经过距离频域多项式相位变换处理并得到距离单元-多普勒数据后,目标能量得到了有效地积累,能够在噪声中有效地检测出目标。
5、单元平均恒虚警检测性能
图5和图6分别在虚警概率为10-4、10-6情况下对本发明方法、RFRFT方法、FRFT方法、RFT方法和MTD方法的恒虚警检测性能对比曲线图。从图5可以看出在信噪比为-17dB时,本发明方法的检测概率Pd=0.9,而在信噪比分别为-5dB、-4dB、0dB和2B时,RFRFT方法、FRFT方法、RFT方法和MTD方法的检测概率Pd=0.9。故可以看出本发明方法可以在低信噪比情况下有效地检测出目标。从图6中也能够看出本发明方法能够在低信噪比情况下检测出目标,检测性能明显优于其他方法。
Claims (2)
1.一种基于距离频域多项式相位变换的空间机动目标检测方法,所述的方法包括按顺序进行下列的步骤:
1)对天基雷达接收到的空间机动目标回波数据进行脉冲压缩处理;
2)沿上述经过脉冲压缩处理后数据的快时间向进行快速傅里叶变换以得到距离频域-慢时间数据并留一个备份数据,接着沿非备份数据的慢时间向进行时延处理;
3)将上述未经过时延的距离频域-慢时间数据与经过时延的距离频域-慢时间数据的共轭进行乘积处理,然后沿乘积结果的距离频域向进行逆快速傅里叶变换处理以得到经过多普勒走动校正后的距离单元-慢时间数据;
4)沿步骤3)中得到的经过多普勒走动校正后的距离单元-慢时间数据的慢时间向进行快速傅里叶变换处理以对目标回波进行相参积累,在相参积累结果基础上完成目标检测;
其特征在于:在步骤3)中,所述的将未经过时延的距离频域-慢时间数据与经过时延的距离频域-慢时间数据的共轭进行乘积处理,然后沿乘积结果的距离频域向进行逆快速傅里叶变换处理以得到经过多普勒走动校正后的距离单元-慢时间数据的方法是:
在步骤2)中得到的未经过时延的距离频域-慢时间数据表达式可写为:
其中A0为回波信号复幅度,f为距离频域,tm=mTr为慢时间,c为光速,为信号波长,M为相干积累脉冲数,m=1,...,M,Tr是脉冲重复周期,fc为载波频率, 为目标与雷达之间的距离,r0为目标的初始距离,v0为目标运动速度,a0为目标加速度,且有:
留一份X(f,tm)的备份数据,并记为X'(f,tm);
则经过时延的距离频域-慢时间数据表达式可写为:
τ为延时量;
将X'(f,tm)与式(4)的共轭进行相乘处理以实现距离频域多项式相位运算:
其中“*”表示共轭处理;
沿的距离频域向进行逆快速傅里叶变换以得到经过多普勒走动校正后的距离单元-慢时间数据表达式为:
其中,t'为快时间,sinc(x)=sin(πx)/(πx)。
2.根据权利要求1所述的基于距离频域多项式相位变换的空间机动目标检测方法,其特征在于:在步骤4)中,所述的沿步骤3)中得到的经过多普勒走动校正后的距离单元-慢时间数据的慢时间向进行快速傅里叶变换处理以对目标回波进行相参积累,在相参积累结果基础上完成目标检测的方法是:沿距离单元-慢时间数据s'r(t',tm)的慢时间向进行快速傅里叶变换处理得到用于目标检测的距离单元-多普勒数据,即实现距离频域多项 式相位变换,然后在距离单元-多普勒数据基础上利用单元平均恒虚警检测原理对空间机动目标进行检测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510308318.0A CN104880698B (zh) | 2015-06-05 | 2015-06-05 | 基于距离频域多项式相位变换的空间机动目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510308318.0A CN104880698B (zh) | 2015-06-05 | 2015-06-05 | 基于距离频域多项式相位变换的空间机动目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104880698A CN104880698A (zh) | 2015-09-02 |
CN104880698B true CN104880698B (zh) | 2017-03-08 |
Family
ID=53948264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510308318.0A Expired - Fee Related CN104880698B (zh) | 2015-06-05 | 2015-06-05 | 基于距离频域多项式相位变换的空间机动目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104880698B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054153A (zh) * | 2016-05-23 | 2016-10-26 | 武汉大学 | 一种基于分数阶变换的海杂波区目标检测与自适应杂波抑制方法 |
CN108919221A (zh) * | 2018-07-17 | 2018-11-30 | 武汉大学 | 一种用于变加速运动目标的相参积累检测方法 |
CN109633595A (zh) * | 2019-01-23 | 2019-04-16 | 北京理工大学 | 星载脉冲多普勒雷达的匀加速运动目标参数快速估计方法 |
CN109946671B (zh) * | 2019-04-12 | 2020-11-03 | 哈尔滨工程大学 | 一种基于双门限判决的水下机动弱目标检测跟踪方法 |
CN110780293B (zh) * | 2019-12-11 | 2021-10-22 | 北京无线电测量研究所 | 搜索雷达的快速阵地优化方法及系统 |
CN114488023B (zh) * | 2021-12-31 | 2024-09-06 | 航天南湖电子信息技术股份有限公司 | 一种双重匹配滤波方法、装置和电子设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104502906A (zh) * | 2014-12-09 | 2015-04-08 | 中国民航大学 | 基于rmdcft的空间超高速机动目标检测方法 |
CN104849708A (zh) * | 2015-05-18 | 2015-08-19 | 中国民航大学 | 基于频域多项式相位变换的高速机动目标参数估计方法 |
-
2015
- 2015-06-05 CN CN201510308318.0A patent/CN104880698B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104502906A (zh) * | 2014-12-09 | 2015-04-08 | 中国民航大学 | 基于rmdcft的空间超高速机动目标检测方法 |
CN104849708A (zh) * | 2015-05-18 | 2015-08-19 | 中国民航大学 | 基于频域多项式相位变换的高速机动目标参数估计方法 |
Non-Patent Citations (3)
Title |
---|
A New Algorithm of ISAR Imaging for Maneuvering Target with Low SNR;YACHAO LI等;《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》;20130131;第49卷(第1期);第543-557页 * |
A Novel Dynamic Programming Algorithm for Track-Before-Detect in Radar Systems;Emanuele Grossi等;《IEEE TRANSACTIONS ON SIGNAL PROCESSING》;20130515;第61卷(第10期);第2608-2619页 * |
ISAR目标运动参数估计及成像技术研究;李亚超;《中国博士学位论文全文数据库 信息科技辑》;20110315;第26-35页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104880698A (zh) | 2015-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104880698B (zh) | 基于距离频域多项式相位变换的空间机动目标检测方法 | |
CN106970371B (zh) | 一种基于Keystone和匹配滤波的目标检测方法 | |
Suo et al. | Detection of high‐speed and accelerated target based on the linear frequency modulation radar | |
CN104076355B (zh) | 基于动态规划的强杂波环境中弱小目标检测前跟踪方法 | |
Shan et al. | Efficient architecture and hardware implementation of coherent integration processor for digital video broadcast‐based passive bistatic radar | |
CN103176168B (zh) | 一种机载非正侧视阵雷达近程杂波对消方法 | |
CN107450055B (zh) | 基于离散线性调频傅立叶变换的高速机动目标检测方法 | |
CN103499825B (zh) | 一种高轨航天器bds弱信号快速捕获方法 | |
CN104502906B (zh) | 基于rmdcft的空间超高速机动目标检测方法 | |
CN103399310A (zh) | 基于相位差分Radon-Lv分布的雷达微弱动目标检测方法 | |
CN103777190B (zh) | 一种双基地mimo雷达高速高机动目标的角度估计方法 | |
CN104020451A (zh) | 基于聚类的外辐射源雷达目标航迹处理方法 | |
CN103344949A (zh) | 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法 | |
CN106507958B (zh) | 外辐射源雷达信号实时相干积累的方法 | |
Huang et al. | Approach for space‐based radar manoeuvring target detection and high‐order motion parameter estimation | |
CN105306399A (zh) | 一种雷达通信一体化信号的优化方法 | |
CN103018727A (zh) | 一种基于样本训练的机载雷达非平稳杂波抑制方法 | |
CN102879766A (zh) | 一种高频地波雷达弱目标检测跟踪方法和装置 | |
CN104849708B (zh) | 基于频域多项式相位变换的高速机动目标参数估计方法 | |
CN103323829A (zh) | 基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法 | |
CN103954951A (zh) | 基于功率中值和归一化协方差估计的海面目标检测方法 | |
CN103529437A (zh) | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 | |
CN104007421A (zh) | 基于全变差和压缩感知的罗兰c被动雷达toa估计方法 | |
Oveis et al. | Coherent method for ground‐moving target indication and velocity estimation using Hough transform | |
CN114047502B (zh) | 一种针对高速机动目标的相参积累方法及相参积累系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170308 Termination date: 20170605 |