CN109521410A - 基于时间反转变换的高速机动目标相参积累检测方法 - Google Patents
基于时间反转变换的高速机动目标相参积累检测方法 Download PDFInfo
- Publication number
- CN109521410A CN109521410A CN201811372326.1A CN201811372326A CN109521410A CN 109521410 A CN109521410 A CN 109521410A CN 201811372326 A CN201811372326 A CN 201811372326A CN 109521410 A CN109521410 A CN 109521410A
- Authority
- CN
- China
- Prior art keywords
- maneuvering target
- speed maneuvering
- echo signal
- speed
- time domain
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 67
- 238000009825 accumulation Methods 0.000 title claims abstract description 62
- 230000009466 transformation Effects 0.000 title claims abstract description 46
- 230000001133 acceleration Effects 0.000 claims abstract description 50
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000007906 compression Methods 0.000 claims abstract description 13
- 230000001427 coherent effect Effects 0.000 claims description 54
- 239000011159 matrix material Substances 0.000 claims description 28
- 230000006835 compression Effects 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000012937 correction Methods 0.000 abstract description 14
- 230000008569 process Effects 0.000 abstract description 6
- 239000000523 sample Substances 0.000 abstract 1
- 238000013508 migration Methods 0.000 description 37
- 230000005012 migration Effects 0.000 description 37
- 238000004088 simulation Methods 0.000 description 18
- 238000012545 processing Methods 0.000 description 10
- 230000036461 convulsion Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000007306 turnover Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000035485 pulse pressure Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
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
技术领域
本发明涉及通信技术领域,更进一步涉及到雷达信号处理技术领域的一种基于时间反转变换的高速机动目标相参积累检测方法。本发明可用于现有雷达利用长时间相参积累技术实现对临近空间高超声速机动目标的远程探测。
背景技术
随着临近空间资源的不断开发利用,雷达对临近空间的高超声速机动目标的探测逐渐成为雷达领域关注的热点。这些高速机动目标具有高超音速、强机动性和小RCS(雷达散射截面积)等典型特点,为了实现对高速机动目标的有效检测,需要采用以时间换取能量的长时间积累技术对回波数据进行能量积累。传统的雷达检测方法面临着因目标高速度造成的距离徙动和多普勒徙动等重大挑战,尤其在强机动性下,多普勒徙动的复杂表现形式加剧了高速机动目标检测的难度。如何针对上述问题实现现有雷达对高速机动目标的快速检测,是当前雷达探测领域亟需解决的问题。
中国人民解放军海军航空工程学院在其申请的专利文献“基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法”(申请号201310240003.8,公开号CN103323829 A)中公开了一种基于Radon-分数阶模糊函数RFRAF的高速机动目标相参积累检测方法。该方法首先通过距离-速度-加速度-第二加速度-阶数域上的五维搜索抽取出目标回波信号,然后利用FRFT实现抽取信号能量的相参积累,同时完成目标的距离徙动校正和多普勒徙动补偿,从而提高目标检测性能。但是,该方法仍然存在的不足是,该方法需要同时搜索距离、速度、加速度、第二加速度以及分数阶傅立叶变换阶数五维的参数空间,运算复杂度高,导致检测速度慢。
李小龙在其发表的论文“高速机动目标长时间相参积累算法研究”(电子科技大学2017博士论文)中提出了一种基于广义KT与广义去调频处理GKTGDP的高速机动目标相参积累检测方法。该方法通过依次对第二加速度引起的三阶距离走动做三阶梯形变换KT,对加速度引起的二阶距离走动做六阶KT,对剩余一阶距离走动做二阶KT,校正和补偿距离徙动和多普勒徙动,实现目标回波信号能量的相参积累。但是,该方法仍然存在的不足是,该方法在补偿和校正过程中仍然需要加加速度、加速度及速度模糊数的参数搜索,运算量依旧很大,并且低信噪比情况下回波能量的积累性能会下降,导致目标的检测概率低。该论文中也提出了一种基于时间序列翻转TRT-SKT-LVD的变加速运动目标三阶距离徙动校正与相参积累检测方法。该方法首先利用时间反转变换TRT实现线性距离徙动与三阶距离徙动的校正,再利用二阶KTSKT实现二阶距离徙动的校正,最后利用吕氏分布LVD实现多普勒徙动的校正和目标能量的积累。但是,该方法仍然存在的不足是,运算流程过于复杂,运算复杂度过高,处理过程SKT会造成相参积累能量损失,导致目标的检测概率低,同时处理过程LVD会受到冗余信息和参数估计范围的限制。
综上所述,对于高速机动目标长时间相参积累检测方法在现有雷达远程探测领域的应用,目前已有的方法仍存在因运算复杂度过高而导致运算速度慢,因进行多次变换后相参积累能量损失而导致目标的检测概率低,因参数估计范围受限制而导致不利于工程应用等问题。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种基于时间反转变换的高速机动目标相参积累检测方法。本发明利用时间反转变换和构建的距离频域-慢时间域二次相位补偿函数同时校正距离徙动和多普勒徙动,只需对加速度参数进行搜索,有效地降低了计算复杂度,缩短了检测时间,仅进行一次变换,避免积累能量的损失,提高了对高速机动目标的检测概率。
实现本发明目的具体思路是:对接收的高速机动目标回波信号在快时间域进行解调和脉冲压缩处理,对解调和脉冲压缩后的高速机动目标回波信号在快时间域做快速傅立叶变换,对快速傅立叶变换后的高速机动目标回波信号做时间反转变换,利用时间反转变换对因高速机动目标的速度引起的线性距离徙动和因高速机动目标的加加速度引起的三阶距离徙动和多普勒徙动进行同时校正,构建距离频域-慢时间域二次相位补偿函数,利用该补偿函数补偿和校正因高速机动目标的加速度引起的距离徙动和多普勒徙动,对补偿后的信号完成相参积累,利用相参积累后的快时间域-多普勒域高速机动目标回波信号检测目标是否存在。
本发明的具体步骤包括如下:
(1)雷达接收一个调制形式为线性调频脉冲的高速机动目标回波信号;
(2)对高速机动目标回波信号在快时间域进行解调和脉冲压缩处理:
(2a)利用解调公式,在快时间域,对雷达接收的高速机动目标回波信号进行解调;
(2b)利用脉冲压缩公式,在快时间域,对解调后的高速机动目标回波信号进行脉冲压缩;
(3)对解调和脉冲压缩后的高速机动目标回波信号做快速傅立叶变换:
利用快速傅立叶变换公式,在快时间域,对解调和脉冲压缩后的高速机动目标回波信号做快速傅立叶变换;
(4)对快速傅立叶变换后的高速机动目标回波信号做时间反转变换:
(4a)将快速傅立叶变换后的高速机动目标回波信号中慢时间域的数据,存入到一个空二维矩阵的行中,将距离频域的数据存入到空二维矩阵的列中;
(4b)利用数据交换方式,对二维数据矩阵中的行数据进行慢时间域反转,得到反转后的高速机动目标回波二维数据矩阵;
(4c)将反转前的高速机动目标回波二维数据矩阵和反转后的高速机动目标回波二维数据矩阵做矩阵相乘,得到时间反转变换后的距离频域-慢时间域高速机动目标回波信号;
(5)确定高速机动目标加速度值的搜索范围和搜索间隔:
(5a)将高速机动目标加速度值的搜索范围设定为[-amax,amax],其中,amax表示高速机动目标最大加速度,amax=300m/s2;
(5b)按照下式,设定高速机动目标加速度值的搜索间隔:
其中,Δa表示高速机动目标加速度参数的搜索间隔,λ表示雷达接收的高速机动目标回波信号的波长,λ=c/fc,c表示光速,fc表示雷达接收的高速机动目标回波信号的载频,T表示雷达接收高速机动目标回波信号所需要的时间;
(6)构建距离频域-慢时间域二次相位补偿函数:
(6a)利用高速机动目标加速度值估计公式,在其搜索范围内,计算距离频域-慢时间域二次相位补偿函数中的高速机动目标加速度估计值;
(6b)按照下式,构建距离频域-慢时间域二次相位补偿函数:
其中,Cm表示距离频域-慢时间域二次相位补偿函数,exp表示以自然常数e为底的指数操作,j表示虚数单位符号,π表示圆周率,f表示距离频域内的距离频率,表示距离频域-慢时间域二次相位补偿函数中的高速机动目标加速度估计值,tm表示慢时间域的时间采样点,M表示雷达接收的高速机动目标回波信号的脉冲总数;
(7)补偿和校正时间反转变换后的高速机动目标回波二维数据矩阵:
将距离频域-慢时间域二次相位补偿函数与时间反转变换后的高速机动目标回波信号相乘,得到补偿和校正后的距离频域-慢时间域高速机动目标回波信号;
(8)完成高速机动目标回波信号能量的相参积累:
(8a)利用快速逆傅里叶变换公式,在距离频域,对补偿和校正后的距离频域-慢时间域高速机动目标回波信号做快速逆傅里叶变换,得到快时间域-慢时间域高速机动目标回波信号;
(8b)利用快速傅里叶变换公式,在慢时间域,对快时间域-慢时间域高速机动目标回波信号做快速傅里叶变换,得到相参积累后的快时间域-多普勒域高速机动目标回波信号;
(9)检测高速机动目标:
(9a)利用检测门限的计算公式,计算检测门限值;
(9b)对相参积累后的快时间域-多普勒域高速机动目标回波信号的所有数据取绝对值;
(9c)判断所有数据绝对值中的最大值是否大于等于检测门限,若是,判定回波信号中存在高速机动目标,否则,判定回波信号中不存在高速机动目标。
本发明与现有技术相比具有以下优点:
第一,由于本发明利用高速机动目标加速度估计值构建距离频域-慢时间域二次相位补偿函数,利用该补偿函数校正和补偿因高速机动目标的加速度引起的距离徙动和多普勒徙动,克服了现有技术因运算复杂度高,导致检测速度慢的问题,使得本发明只需对高速机动目标的加速度进行搜索,降低运算复杂度,提高了雷达对高速机动目标的检测速度,运算时间处于秒级,更加利于工程上的实时应用。
第二,由于本发明对距离频域-慢时间域的高速机动目标回波信号做时间反转变换,利用时间反转变换校正因高速机动目标的速度引起的线性距离徙动和因高速机动目标的加加速度引起的三阶距离徙动和多普勒徙动,克服了现有技术因利用多次梯形变换KT校正距离徙动和多普勒徙动,造成相参积累能量损失,高速机动目标的检测概率低的问题,使得本发明能够高效地补偿因高速机动目标的加加速度带来的三阶距离徙动和多普勒徙动,完成高速机动目标回波能量的相参积累,能更好的适应低信噪比的环境,具有更优的检测性能。
第三,本发明只需对高速机动目标的加速度进行搜索和估计,克服了现有技术在处理过程中参数估计范围受限制的问题,使得本发明能够不受参数估计范围的限制,具有更广泛的适用范围。
附图说明:
图1是本发明的流程图;
图2是本发明的仿真图。
具体实施方式
下面结合附图,对本发明做进一步详细描述。
参照图1,对本发明的具体实施步骤做进一步详细描述。
步骤1,雷达接收一个调制形式为线性调频脉冲的高速机动目标回波信号。
步骤2,对高速机动目标回波信号在快时间域进行解调和脉冲压缩处理。
利用解调公式,在快时间域,对雷达接收的高速机动目标回波信号进行解调。
所述的解调公式如下:
其中,sr'表示解调后的高速机动目标回波信号,sr表示雷达接收的高速机动目标回波信号,exp表示以自然常数e为底的指数操作,j表示虚数单位符号,π表示圆周率,fc表示雷达接收的高速机动目标回波信号的载频,表示快时间域的采样时间点,的范围是Tr表示雷达发射端发射的脉冲串的信号的脉冲重复周期,采样间隔为1/fs,fs表示采样频率。
利用脉冲压缩公式,在快时间域,对解调后的高速机动目标回波信号进行脉冲压缩。
所述的脉冲压缩公式如下:
其中,sp表示解调和脉冲压缩后的高速机动目标回波信号,表示卷积操作,rect(·)表示矩形窗因子,当当Tp表示雷达接收的高速机动目标回波信号的脉冲宽度,|·|表示取绝对值操作,γ表示雷达接收的高速机动目标回波信号的调频斜率。
步骤3,对解调和脉冲压缩后的高速机动目标回波信号做快速傅立叶变换。
利用快速傅立叶变换公式,在快时间域,对解调和脉冲压缩后的高速机动目标回波信号做快速傅立叶变换。
所述的快速傅里叶变换公式如下:
其中,sp'表示快速傅立叶变换后的距离频域-慢时间域高速机动目标回波信号,∫表示积分操作,f表示距离频域内的距离频率。
步骤4,对快速傅立叶变换后的高速机动目标回波信号做时间反转变换。
将快速傅立叶变换后的高速机动目标回波信号中慢时间域的数据存入到一个空二维矩阵的行中,将距离频域的数据存入到空二维矩阵的列中。
利用数据交换方式,对二维数据矩阵中的行数据进行慢时间域反转,得到反转后的高速机动目标回波二维数据矩阵。
所述的数据交换方式是指,将二维数据矩阵的每一行数据进行首尾依次交换,每一行第一个数据与每一行最后一个数据交换,每一行第二个数据与每一行倒数第二个数据交换,以此类推,得到反转后的高速机动目标回波二维数据矩阵。
将反转前的高速机动目标回波二维数据矩阵和反转后的高速机动目标回波二维数据矩阵做矩阵相乘操作,得到时间反转变换后的距离频域-慢时间域高速机动目标回波信号。
步骤5,确定高速机动目标加速度值的搜索范围和搜索间隔。
将高速机动目标加速度值的搜索范围设定为[-amax,amax],其中,amax表示高速机动目标最大加速度,amax=300m/s2。
按照下式,设定高速机动目标加速度值的搜索间隔:
其中,Δa表示高速机动目标加速度参数的搜索间隔,λ表示雷达接收的高速机动目标回波信号的波长,λ=c/fc,c表示光速,T表示雷达接收高速机动目标回波信号所需要的时间,实际搜索时搜索间隔可采用变步长法,即先采用大步长(10Δa)确定加速度的粗估计值,再采用小步长(0.1Δa)确定加速度的精估计值。
步骤6,构建距离频域-慢时间域二次相位补偿函数。
利用高速机动目标加速度值估计公式,在其搜索范围内,计算距离频域-慢时间域二次相位补偿函数中的高速机动目标加速度估计值。
所述的高速机动目标加速度值估计公式如下:
其中,表示距离频域-慢时间域二次相位补偿函数中的高速机动目标加速度估计值,表示当后接表达式有最大值时取对应的下标值的操作,tm表示慢时间域的时间采样点,M表示雷达接收的高速机动目标回波信号的脉冲总数,sn表示时间反转变换后的距离频域-慢时间域高速机动目标回波信号,ax表示待估计的高速机动目标加速度,ax的取值范围是[-amax,amax]。
在所述的高速机动目标加速度值估计公式中,称为能量积累函数,当搜索的高速机动目标的加速度估计值与高速机动目标的加速度真实值相匹配时,目标能量将集中于同一个距离单元中,此时,能量积累将达到最大值,将该估计值代入二次相位补偿函数模型中,即可完全补偿和校正由高速机动目标的加速度带来的距离徙动和多普勒徙动。
按照下式,构建距离频域-慢时间域二次相位补偿函数:
其中,Cm表示距离频域-慢时间域二次相位补偿函数。
步骤7,补偿和校正时间反转变换后的高速机动目标回波二维数据矩阵。
将距离频域-慢时间域二次相位补偿函数与时间反转变换后的高速机动目标回波信号相乘,得到补偿和校正后的距离频域-慢时间域高速机动目标回波信号。
步骤8,完成高速机动目标回波信号能量的相参积累。
利用快速逆傅里叶变换公式,在距离频域,对补偿和校正后的距离频域-慢时间域高速机动目标回波信号做快速逆傅里叶变换,得到快时间域-慢时间域高速机动目标回波信号。
所述的快速逆傅里叶变换公式如下:
其中,Sf'表示快速逆傅里叶变换后的快时间域-慢时间域高速机动目标回波信号,Sf表示补偿和校正后的距离频域-慢时间域高速机动目标回波信号。
利用快速傅里叶变换公式,在慢时间域,对快时间域-慢时间域高速机动目标回波信号做快速傅里叶变换,得到相参积累后的快时间域-多普勒域高速机动目标回波信号。
所述的快速傅里叶变换公式如下:
其中,SF表示相参积累后的快时间域-多普勒域高速机动目标回波信号,fm表示多普勒域的多普勒频率。
步骤9,检测高速机动目标。
利用检测门限的计算公式,计算检测门限值。
所述的检测门限的计算公式如下:
其中,η表示检测门限,表示开平方根操作,ln表示以自然常数e为底的对数操作,Pfa表示系统给定的虚警概率,σ表示系统噪声方差。
对相参积累后的快时间域-多普勒域高速机动目标回波信号的所有数据取绝对值。
判断所有数据绝对值中的最大值是否大于等于检测门限,若是,判定回波信号中存在高速机动目标,否则,判定回波信号中不存在高速机动目标。
下面结合仿真实验对本发明做进一步的描述。
1.仿真条件:
高速机动目标的运动参数设置为:初始斜距R0=250km,初始速度v0=1200m/s,加速度a0=300m/s2,加加速度g0=180m/s3。雷达的系统参数设置为:载频fc=10GHz,发射信号带宽B=20MHz,采样频率fs=40MHz,脉冲重复频率fr=200Hz,脉冲宽度Tp=10μs,发射的脉冲个数M=401。计算机主要配置为:CPU为Intel Core i7-4970,3.6GHz,RAM为8GB;操作系统为Windows 7;运行软件为Matlab R2012b。
2.仿真内容与结果分析:
本发明的仿真实验有三个。
仿真实验1:
本发明的仿真实验1是采用本发明和三个现有技术(基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法、基于广义KT与广义去调频处理的高速机动目标相参积累检测方法、基于时间序列翻转的变加速运动目标三阶距离徙动校正与相参积累检测方法),分别对上述仿真条件下的高速机动目标进行检测,信噪比为-5dB。
表1仿真结果对比表
本发明仿真实验1的仿真结果如表1所示,分别统计本发明和三个现有技术每种方法在计算机上运行时,从仿真实验开始到仿真实验结束所需的运算时间,分列表1中。从表1可以看出,基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法的运算时间巨大,基于广义KT与广义去调频处理的高速机动目标相参积累检测方法、基于时间序列翻转的变加速运动目标三阶距离徙动校正与相参积累检测方法消耗的时间也处于百秒级,与上述三种方法相比,本发明只需对一维运动参数搜索,运算时间处于秒级,运算时间大大缩短,更加利于工程上的实时应用。
仿真实验2:
本发明的仿真实验2是用本发明和一个现有技术(基于广义KT与广义去调频处理的高速机动目标相参积累检测方法),分别对上述实验条件下的高速机动目标进行检测,信噪比区间为[-10,20]dB,分别对两种方法在每一个信噪比下进行500次Monte-Carlo试验,得到的结果如图2(a)所示。
图2(a)中的横轴表示脉压后信噪比的值,单位是dB,纵轴表示检测概率值,以圆圈标示的曲线表示本发明的检测性能曲线,以星号标示的曲线表示基于广义KT与广义去调频处理的高速机动目标相参积累检测方法的检测性能曲线。从图2(a)可以看出,相比于基于广义KT与广义去调频处理的高速机动目标相参积累检测方法,本发明在低信噪比情况下,具有更优的检测性能。
仿真实验3:
本发明的仿真实验3是用本发明和一种现有技术(基于时间序列翻转的变加速运动目标三阶距离徙动校正与相参积累检测方法),分别对上述实验条件下的高速机动目标进行检测,信噪比为-5dB,得到结果如图2(b)和图2(c)所示。
图2(b)是采用本发明的方法对高速机动目标回波信号能量的相参积累图,x轴表示多普勒域的多普勒频率值,单位是Hz,y轴表示距离值,单位是Km,z轴表示幅度值。图2(c)是采用现有技术基于时间序列翻转的变加速运动目标三阶距离徙动校正与相参积累检测方法,对高速机动目标回波信号能量的相参积累图,x轴表示中心频率单元值,y轴表示调频斜率单元值,z轴表示幅度值。
从图2(b)可以看出,在本发明对高速机动目标的检测结果中,高速机动目标回波信号的能量积聚在同一个距离-多普勒单元中,形成比较清晰的峰值,从图2(c)可以看出,在基于时间序列翻转的变加速运动目标三阶距离徙动校正与相参积累检测方法对高速机动目标的检测结果中,高速机动目标回波信号的能量分散在不同的单元中。从而证明,在仿真条件设定的参数超出基于时间序列翻转的变加速运动目标三阶距离徙动校正与相参积累检测方法的参数估计范围时,基于时间序列翻转的变加速运动目标三阶距离徙动校正与相参积累检测方法无法实现回波信号能量的有效积累,检测性能较差;本发明可以不受上述方法的参数估计范围的限制,有效地实现相参积累,完成对高速机动目标的检测,具有更广泛的适用范围。
Claims (9)
1.一种基于时间反转变换的高速机动目标相参积累检测方法,其特征在于,对距离频域-慢时间域的高速机动目标回波信号做时间反转变换;利用高速机动目标的加速度估计值构建距离频域-慢时间域二次相位补偿函数,并补偿和校正时间反转变换后的距离频域-慢时间域高速机动目标回波信号;该方法的具体步骤包括如下:
(1)雷达接收一个调制形式为线性调频脉冲的高速机动目标回波信号;
(2)对高速机动目标回波信号在快时间域进行解调和脉冲压缩处理:
(2a)利用解调公式,在快时间域,对雷达接收的高速机动目标回波信号进行解调;
(2b)利用脉冲压缩公式,在快时间域,对解调后的高速机动目标回波信号进行脉冲压缩;
(3)对解调和脉冲压缩后的高速机动目标回波信号做快速傅立叶变换:
利用快速傅立叶变换公式,在快时间域,对解调和脉冲压缩后的高速机动目标回波信号做快速傅立叶变换;
(4)对快速傅立叶变换后的高速机动目标回波信号做时间反转变换:
(4a)将快速傅立叶变换后的高速机动目标回波信号中慢时间域的数据,存入到一个空二维矩阵的行中,将距离频域的数据存入到空二维矩阵的列中;
(4b)利用数据交换方式,对二维数据矩阵中的行数据进行慢时间域反转,得到反转后的高速机动目标回波二维数据矩阵;
(4c)将反转前的高速机动目标回波二维数据矩阵与反转后的高速机动目标回波二维数据矩阵做矩阵相乘,得到时间反转变换后的距离频域-慢时间域高速机动目标回波信号;
(5)确定高速机动目标加速度值的搜索范围和搜索间隔:
(5a)将高速机动目标加速度值的搜索范围设定为[-amax,amax],其中,amax表示高速机动目标最大加速度,amax=300m/s2;
(5b)按照下式,设定高速机动目标加速度值的搜索间隔:
其中,Δa表示高速机动目标加速度参数的搜索间隔,λ表示雷达接收的高速机动目标回波信号的波长,λ=c/fc,c表示光速,fc表示雷达接收的高速机动目标回波信号的载频,T表示雷达接收高速机动目标回波信号所需要的时间;
(6)构建距离频域-慢时间域二次相位补偿函数:
(6a)利用高速机动目标加速度值估计公式,在其搜索范围内,计算距离频域-慢时间域二次相位补偿函数中的高速机动目标加速度估计值;
(6b)按照下式,构建距离频域-慢时间域二次相位补偿函数:
其中,Cm表示距离频域-慢时间域二次相位补偿函数,exp表示以自然常数e为底的指数操作,j表示虚数单位符号,π表示圆周率,f表示距离频域内的距离频率,表示距离频域-慢时间域二次相位补偿函数中的高速机动目标加速度估计值,tm表示慢时间域的时间采样点,M表示雷达接收的高速机动目标回波信号的脉冲总数;
(7)补偿和校正时间反转变换后的高速机动目标回波二维数据矩阵:
将距离频域-慢时间域二次相位补偿函数与时间反转变换后的高速机动目标回波信号相乘,得到补偿和校正后的距离频域-慢时间域高速机动目标回波信号;
(8)完成高速机动目标回波信号能量的相参积累:
(8a)利用快速逆傅里叶变换公式,在距离频域,对补偿和校正后的距离频域-慢时间域高速机动目标回波信号做快速逆傅里叶变换,得到快时间域-慢时间域高速机动目标回波信号;
(8b)利用快速傅里叶变换公式,在慢时间域,对快时间域-慢时间域高速机动目标回波信号做快速傅里叶变换,得到相参积累后的快时间域-多普勒域高速机动目标回波信号;
(9)检测高速机动目标:
(9a)利用检测门限的计算公式,计算检测门限值;
(9b)对相参积累后的快时间域-多普勒域高速机动目标回波信号的所有数据取绝对值;
(9c)判断所有数据绝对值中的最大值是否大于等于检测门限,若是,判定回波信号中存在高速机动目标,否则,判定回波信号中不存在高速机动目标。
2.根据权利要求1所述的基于时间反转变换的高速机动目标相参积累检测方法,其特征在于:步骤(2a)中所述的解调公式如下:
其中,sr'表示解调后的高速机动目标回波信号,sr表示雷达接收的高速机动目标回波信号,表示快时间域的采样时间点,的范围是Tr表示雷达发射端发射的脉冲串的信号的脉冲重复周期,采样间隔为1/fs,fs表示采样频率。
3.根据权利要求1所述的基于时间反转变换的高速机动目标相参积累检测方法,其特征在于:步骤(2b)中所述的脉冲压缩公式如下:
其中,sp表示解调和脉冲压缩后的高速机动目标回波信号,表示卷积操作,rect(·)表示矩形窗因子,当时当时Tp表示雷达接收的高速机动目标回波信号的脉冲宽度,|·|表示取绝对值操作,γ表示雷达接收的高速机动目标回波信号的调频斜率。
4.根据权利要求1所述的基于时间反转变换的高速机动目标相参积累检测方法,其特征在于:步骤(3)中所述的快速傅里叶变换公式如下:
其中,sp'表示快速傅立叶变换后的距离频域-慢时间域高速机动目标回波信号,∫表示积分操作。
5.根据权利要求1所述的基于时间反转变换的高速机动目标相参积累检测方法,其特征在于:步骤(4b)中所述的数据交换方式是指,将二维数据矩阵的每一行数据进行首尾依次交换,每一行第一个数据与每一行最后一个数据交换,每一行第二个数据与每一行倒数第二个数据交换,以此类推,得到反转后的高速机动目标回波二维数据矩阵。
6.根据权利要求1所述的基于时间反转变换的高速机动目标相参积累检测方法,其特征在于:步骤(6a)中所述的高速机动目标加速度值估计公式如下:
其中,表示当后接表达式有最大值时取对应的下标值的操作,sn表示时间反转变换后的距离频域-慢时间域高速机动目标回波信号,ax表示待估计的高速机动目标加速度,ax的取值范围是[-amax,amax]。
7.根据权利要求1所述的基于时间反转变换的高速机动目标相参积累检测方法,其特征在于:步骤(8a)中所述的快速逆傅里叶变换公式如下:
其中,Sf'表示快速逆傅里叶变换后的快时间域-慢时间域高速机动目标回波信号,Sf表示补偿和校正后的距离频域-慢时间域高速机动目标回波信号。
8.根据权利要求1所述的基于时间反转变换的高速机动目标相参积累检测方法,其特征在于:步骤(8b)中所述的快速傅里叶变换公式如下:
其中,SF表示相参积累后的快时间域-多普勒域高速机动目标回波信号,fm表示多普勒域的多普勒频率。
9.根据权利要求1所述的基于时间反转变换的高速机动目标相参积累检测方法,其特征在于:步骤(9a)中所述的检测门限的计算公式如下:
其中,η表示检测门限,表示开平方根操作,ln表示以自然常数e为底的对数操作,Pfa表示系统给定的虚警概率,σ表示系统噪声方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811372326.1A CN109521410B (zh) | 2018-11-16 | 2018-11-16 | 基于时间反转变换的高速机动目标相参积累检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811372326.1A CN109521410B (zh) | 2018-11-16 | 2018-11-16 | 基于时间反转变换的高速机动目标相参积累检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109521410A true CN109521410A (zh) | 2019-03-26 |
CN109521410B CN109521410B (zh) | 2022-12-02 |
Family
ID=65778350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811372326.1A Active CN109521410B (zh) | 2018-11-16 | 2018-11-16 | 基于时间反转变换的高速机动目标相参积累检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109521410B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031809A (zh) * | 2019-04-15 | 2019-07-19 | 北京电子工程总体研究所 | 一种雷达多接收通道幅相一致性标校方法和可读存储介质 |
CN110095766A (zh) * | 2019-05-24 | 2019-08-06 | 西安电子科技大学 | 基于非均匀重采样技术的机动目标相干积累检测方法 |
CN110133661A (zh) * | 2019-05-15 | 2019-08-16 | 成都锦江电子系统工程有限公司 | 一种相位历程建模补偿相干积累信噪比损失方法 |
CN110161477A (zh) * | 2019-06-12 | 2019-08-23 | 西安电子科技大学 | 基于多变量重采样相关函数的机动目标检测方法 |
CN110398730A (zh) * | 2019-06-26 | 2019-11-01 | 中国人民解放军战略支援部队信息工程大学 | 基于坐标旋转和非均匀傅里叶变换机动目标相参检测方法 |
CN110865351A (zh) * | 2019-11-22 | 2020-03-06 | 北京理工大学 | 一种高速机动目标参数估计方法 |
CN113009446A (zh) * | 2021-03-02 | 2021-06-22 | 中国科学院空天信息创新研究院 | 一种基于最佳解调算子的悬停低慢小目标检测方法及装置 |
CN113030861A (zh) * | 2021-03-07 | 2021-06-25 | 西北工业大学 | 一种时间反转和分数阶小波变换的水声信号检测方法 |
CN113109779A (zh) * | 2020-02-28 | 2021-07-13 | 加特兰微电子科技(上海)有限公司 | 提升目标探测精度的方法、集成电路、无线电器件及电子设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120280851A1 (en) * | 2010-01-11 | 2012-11-08 | Bae Systems Plc | Pulse radar range profile motion compensation |
CN107450055A (zh) * | 2017-07-15 | 2017-12-08 | 西安电子科技大学 | 基于离散线性调频傅立叶变换的高速机动目标检测方法 |
CN108761404A (zh) * | 2018-03-23 | 2018-11-06 | 电子科技大学 | 一种基于二次相位函数参数估计及补偿的改进算法 |
-
2018
- 2018-11-16 CN CN201811372326.1A patent/CN109521410B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120280851A1 (en) * | 2010-01-11 | 2012-11-08 | Bae Systems Plc | Pulse radar range profile motion compensation |
CN107450055A (zh) * | 2017-07-15 | 2017-12-08 | 西安电子科技大学 | 基于离散线性调频傅立叶变换的高速机动目标检测方法 |
CN108761404A (zh) * | 2018-03-23 | 2018-11-06 | 电子科技大学 | 一种基于二次相位函数参数估计及补偿的改进算法 |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031809B (zh) * | 2019-04-15 | 2021-07-23 | 北京电子工程总体研究所 | 一种雷达多接收通道幅相一致性标校方法和可读存储介质 |
CN110031809A (zh) * | 2019-04-15 | 2019-07-19 | 北京电子工程总体研究所 | 一种雷达多接收通道幅相一致性标校方法和可读存储介质 |
CN110133661A (zh) * | 2019-05-15 | 2019-08-16 | 成都锦江电子系统工程有限公司 | 一种相位历程建模补偿相干积累信噪比损失方法 |
CN110095766A (zh) * | 2019-05-24 | 2019-08-06 | 西安电子科技大学 | 基于非均匀重采样技术的机动目标相干积累检测方法 |
CN110095766B (zh) * | 2019-05-24 | 2023-03-21 | 西安电子科技大学 | 基于非均匀重采样技术的机动目标相干积累检测方法 |
CN110161477A (zh) * | 2019-06-12 | 2019-08-23 | 西安电子科技大学 | 基于多变量重采样相关函数的机动目标检测方法 |
CN110161477B (zh) * | 2019-06-12 | 2022-12-06 | 西安电子科技大学 | 基于多变量重采样相关函数的机动目标检测方法 |
CN110398730A (zh) * | 2019-06-26 | 2019-11-01 | 中国人民解放军战略支援部队信息工程大学 | 基于坐标旋转和非均匀傅里叶变换机动目标相参检测方法 |
CN110398730B (zh) * | 2019-06-26 | 2021-07-06 | 中国人民解放军战略支援部队信息工程大学 | 基于坐标旋转和非均匀傅里叶变换机动目标相参检测方法 |
CN110865351A (zh) * | 2019-11-22 | 2020-03-06 | 北京理工大学 | 一种高速机动目标参数估计方法 |
CN113109779A (zh) * | 2020-02-28 | 2021-07-13 | 加特兰微电子科技(上海)有限公司 | 提升目标探测精度的方法、集成电路、无线电器件及电子设备 |
CN113109779B (zh) * | 2020-02-28 | 2024-03-22 | 加特兰微电子科技(上海)有限公司 | 提升目标探测精度的方法、集成电路、无线电器件及电子设备 |
CN113009446B (zh) * | 2021-03-02 | 2022-05-03 | 中国科学院空天信息创新研究院 | 一种基于最佳解调算子的悬停低慢小目标检测方法及装置 |
CN113009446A (zh) * | 2021-03-02 | 2021-06-22 | 中国科学院空天信息创新研究院 | 一种基于最佳解调算子的悬停低慢小目标检测方法及装置 |
CN113030861A (zh) * | 2021-03-07 | 2021-06-25 | 西北工业大学 | 一种时间反转和分数阶小波变换的水声信号检测方法 |
CN113030861B (zh) * | 2021-03-07 | 2023-06-23 | 西北工业大学 | 一种时间反转和分数阶小波变换的水声信号检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109521410B (zh) | 2022-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109521410B (zh) | 基于时间反转变换的高速机动目标相参积累检测方法 | |
CN106970371B (zh) | 一种基于Keystone和匹配滤波的目标检测方法 | |
CN104076351B (zh) | 一种用于高速高机动目标的相参积累检测方法 | |
CN107450055B (zh) | 基于离散线性调频傅立叶变换的高速机动目标检测方法 | |
CN107561508B (zh) | 一种用于匀加速运动目标的相参积累检测方法 | |
CN108089171B (zh) | 一种针对无人机目标的雷达快速检测方法 | |
CN106526568A (zh) | 基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法 | |
CN106896358A (zh) | 一种基于位置旋转变换的高速目标相参积累检测方法 | |
CN111736128A (zh) | 基于skt-siaf-mscft的相参积累方法 | |
CN102798855B (zh) | 一种基于数字电视信号的直升机目标识别方法 | |
CN110161477B (zh) | 基于多变量重采样相关函数的机动目标检测方法 | |
CN102353953A (zh) | 单通道sar地面运动目标成像方法 | |
CN102901956A (zh) | 一种雷达微弱目标检测方法 | |
CN111551922B (zh) | 三维空间双/多基地雷达高速目标检测方法 | |
CN111505599B (zh) | 一种基于精确距离演化模型的相参积累检测方法 | |
CN111045002A (zh) | 基于TRT和SNuFFT的机动目标相参积累方法 | |
CN108549066A (zh) | 一种基于尺度rft的宽带雷达高速目标积累检测方法 | |
CN106597425A (zh) | 一种基于机器学习的雷达目标定位方法 | |
CN112162246B (zh) | 一种基于贝叶斯网络雷达系统复杂电磁环境效应分析方法 | |
CN110398730A (zh) | 基于坐标旋转和非均匀傅里叶变换机动目标相参检测方法 | |
CN111580063B (zh) | 基于广义解调频-楔形变换的雷达目标检测方法 | |
CN104793194A (zh) | 基于改进的自适应多脉冲压缩的距离-多普勒估计方法 | |
CN110095766A (zh) | 基于非均匀重采样技术的机动目标相干积累检测方法 | |
CN114047502A (zh) | 一种针对高速机动目标的相参积累方法及相参积累系统 | |
CN110426701B (zh) | 一种基于速度维稀疏表示的捷变频雷达目标速度测量方法 |
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 |