CN110095766A - 基于非均匀重采样技术的机动目标相干积累检测方法 - Google Patents
基于非均匀重采样技术的机动目标相干积累检测方法 Download PDFInfo
- Publication number
- CN110095766A CN110095766A CN201910440085.8A CN201910440085A CN110095766A CN 110095766 A CN110095766 A CN 110095766A CN 201910440085 A CN201910440085 A CN 201910440085A CN 110095766 A CN110095766 A CN 110095766A
- Authority
- CN
- China
- Prior art keywords
- time
- signal
- radar
- interpolation
- signal sequence
- 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 62
- 238000009825 accumulation Methods 0.000 title claims abstract description 34
- 230000001427 coherent effect Effects 0.000 title claims abstract description 28
- 238000005516 engineering process Methods 0.000 title claims abstract description 19
- 238000005070 sampling Methods 0.000 title claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 49
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000005314 correlation function Methods 0.000 claims abstract description 20
- 238000007906 compression Methods 0.000 claims abstract description 11
- 238000012952 Resampling Methods 0.000 claims description 39
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 36
- 230000035485 pulse pressure Effects 0.000 claims description 11
- 230000006835 compression Effects 0.000 claims description 10
- 230000001934 delay Effects 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 238000009828 non-uniform distribution Methods 0.000 claims description 3
- 230000001502 supplementing effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 4
- 210000001367 artery Anatomy 0.000 abstract 1
- 210000003462 vein Anatomy 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 14
- 229910052704 radon Inorganic materials 0.000 description 11
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 11
- 230000005012 migration Effects 0.000 description 9
- 238000013508 migration Methods 0.000 description 9
- 230000033001 locomotion Effects 0.000 description 8
- 238000009826 distribution Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010420 art technique Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 235000020030 perry Nutrition 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 230000009466 transformation Effects 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/04—Systems determining presence of a target
-
- 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
- G01S7/411—Identification of targets based on measurements of radar reflectivity
-
- 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
- G01S7/418—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于非均匀重采样技术的机动目标相干积累检测方法，在快时间域对接收的雷达接收信号进行解调和脉冲压缩处理；对脉压缩信号做快速傅里叶变换；计算约束条件并生成空矩阵；判断所选时刻和所选时延是否同时满足约束条件，符合条件则对距离向频域信号进行非均匀重采样，计算对应元素值，不符合条件则对应元素值为0；依次选取所有的时刻和时延，将生成的元素值放入矩阵中；对相关函数矩阵数据执行快速傅里叶变换和快速傅里叶逆变换，完成相干积累；利用相干积累矩阵判定目标是否存在。本发明在检测过程中不受搜索方法的影响，计算复杂度低，检测性能优。
Description
技术领域
本发明属于雷达技术领域，更进一步涉及雷达信号处理技术领域中的一种基于非均匀重采样技术的机动目标相干积累检测方法。本发明可用于低信噪比环境下微弱机动目标检测。
背景技术
微弱目标检测由于其在实际应用中的需求而引起了越来越多的关注。这些微弱机动目标具有小RCS(雷达散射截面积)、回波能量小等典型特点，为了实现对微弱机动目标的有效检测，需要采用以时间换取能量的长时间积累技术对回波数据进行能量积累。传统的雷达检测方法面临着因过长的观察时间造成的高阶距离偏移和多普勒频率偏移等重大挑战。如何针对上述问题实现现有弱雷达机动目标的快速检测，是当前雷达探测领域亟需解决的问题。
中国人民解放军海军航空工程学院在其申请的专利文献“雷达动目标Radon-分数阶傅里叶变换长时间相参积累检测方法”(申请号201310053874.9，公开号CN103176178 A)中公开了一种基于Radon-分数阶傅里叶变换RFRFT的长时间相参积累检测方法。该方法首先将雷达回波距离向解调、脉压，完成脉内积累，然后采用RFRFT方法补偿距离和多普勒徙动，完成长时间脉间相参积累，同时遍历所有搜索参数，构建距离-RFRFT域检测单元图，从而提高目标检测性能，最后对检测单元图进行恒虚警检测，得到目标运动参数估计。但是，该方法仍然存在的不足之处是，在构建检测单元图时需要遍历所有搜索参数，因为该方法对雷达目标参数要进行多维搜索，使得运算复杂度高，导致检测速度慢。
李小龙在其发表的论文“高速机动目标长时间相参积累算法研究”(电子科技大学2017博士论文)中提出了一种基于相邻互相关函数(ACCF)的算法。该方法首先根据多个运动参数将目标的倾斜范围建模为多项式函数，然后采用迭代相邻互相关运算来消除范围偏移并减少多普勒频率偏移的顺序，最后通过傅立叶变换估计运动参数。该方法无需搜索过程即可估计参数，能在高信噪比(SNR)下获得紧密估计性能，计算成本低。但是，该方法仍然存在的不足之处是，该方法采用迭代相邻互相关运算来避免多维搜索，因为重复使用互相关函数会带来额外的噪声项，因此基于ACCF的算法只能在高信噪比环境中工作，不适用于低信噪比的工作环境。
发明内容
本发明的目的在于针对上述已有技术的不足，提出了一种基于非均匀重采样技术的机动目标相干积累检测方法。本发明利用非均匀重采样技术，消除了接收信号中的距离徙动和多普勒徙动，完成机动目标回波能量的相参积累，能更好的适应低信噪比的环境，具有更优的检测性能，本发明无需对运动参数进行搜索，有效地降低了计算复杂度，缩短了检测时间。
实现本发明目的的思路是，对雷达接收的回波信号在快时间域进行解调和脉冲压缩处理，对解调和脉冲压缩后的回波信号在快时间域做快速傅里叶变换，计算约束条件并生成空矩阵，从雷达回波信号慢时间中依次选取一个时刻和一个时延，判断所选时刻和所选时延是否同时满足约束条件，若是，则对快速傅里叶变换后的回波信号进行非均匀重采样，计算第所选时刻的序号行第所选时延的序号列的元素值，若否，将0作为空矩阵第所选时刻的序号行第所选时延的序号列的元素值，判断所有时刻和时延是否选取完毕，若是，将生成的元素值放入空矩阵中对应的位置得到相关函数矩阵，否则，继续选取，对相关函数矩阵数据执行快速傅里叶变换和快速傅里叶逆变换来对目标能量进行积分，完成相干积累，利用快时间域-多普勒域相干积累矩阵检测目标是否存在。
本发明的具体步骤包括如下：
(1)获取回波信号：
雷达接收其自身发射的线性调频信号的回波信号；
(2)获取解调信号：
利用解调公式，在快时间域内解调雷达接收的回波信号，得到解调信号；
(3)获取脉压信号：
利用脉冲压缩公式，在快时间域内对解调信号在距离维进行脉冲压缩，得到脉压信号；
(4)获取距离向频域信号：
利用快速傅里叶变换公式，在快时间域内对脉压信号做快速傅里叶变换，得到距离向频域信号；
(5)按照下式，计算约束条件：
其中，tu表示雷达回波信号慢时间的第u个时刻，u＝1,2,...,M，M表示雷达接收的回波信号的脉冲总数，y表示根据工程需要设定的雷达接收的回波信号的采样频率，z表示雷达发射的线性调频信号的载频，TM表示雷达波束在微弱机动目标上的驻留时间，τv表示雷达回波信号的第v个时延，v＝1,2,...,M；
(6)搭建一个的空矩阵，其行和列的取值均与雷达接收的回波信号的脉冲总数相等；
(7)从雷达回波信号慢时间中依次选取一个时刻和一个时延，判断所选时刻和所选时延是否同时满足约束条件，若是，则执行步骤(8)，否则，执行步骤(9)；
(8)利用非均匀重采样技术，计算元素值；
(8a)对步骤(4)得到的距离向频域信号做快速傅里叶变换，得到离散频域信号；
(8b)在离散频域信号的频域上大于采样频率y的位置之后补7N个零，N表示雷达接收信号的采样点数，得到补零后的离散频域信号；
(8c)将补零后的离散频域信号做快速逆傅里叶变换，得到重采样信号序列；
(8d)按照下式，计算第一次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置：
其中，xi,j表示第一次对重采样信号序列插值前，在重采样信号序列的时间轴上对应雷达回波信号慢时间的第i个时刻和雷达回波信号第j个时延非均匀分布的插值位置，ti表示第i个慢时间时刻，f表示距离频域内的距离频率，表示开平方根操作，τj表示雷达回波信号的第j个时延，i和j与所选时刻的序号和所选时延的序号的取值对应相等；
(8e)利用线性插值公式，计算第一次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置上待插入的值；
(8f)按照下式，计算第二次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置：
其中，yp,q表示第二次对重采样信号序列插值前，在重采样信号序列的时间轴上对应雷达回波信号慢时间的第p个时刻和雷达回波的第q个时延的非均匀分布的插值位置，其中p和q与所选时刻的序号和所选时延的序号的取值对应相等；
(8g)利用线性插值公式，计算第二次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置上待插入的值；
(8h)将两次待插入的值相乘，将乘积作为空矩阵第所选时刻的序号行第所选时延的序号列的元素值；
(9)将0作为空矩阵第所选时刻的序号行第所选时延的序号列的元素值；
(10)判断所有时刻和时延是否选取完毕，若是，将生成的元素值放入空矩阵中对应的位置得到相关函数矩阵，执行步骤(11)，否则，执行步骤(7)；
(11)生成快时间-多普勒域相干积累矩阵：
(11a)利用快速傅里叶变换公式，在快时间域-慢时间域，对相关函数矩阵做快速傅里叶变换，得到频域相关函数矩阵；
(11b)利用快速逆傅里叶变换公式，在距离频域，对频域相关函数矩阵做快速逆傅里叶变换，得到快时间-多普勒域相干积累矩阵；
(12)检测机动目标：
(12a)利用门限公式，计算检测门限值；
(12b)判断快时间-多普勒域相干积累矩阵的所有元素中绝对值的最大值是否大于或等于检测门限，若是，则将雷达接收的回波信号判定为含有微弱机动目标的目标回波信号，否则，将雷达接收的回波信号判定为无目标的高斯白噪声。
本发明与现有技术相比具有以下优点：
第一，由于本发明利用非均匀重采样技术，设计了一个用以消除接收信号中的距离徙动和多普勒徙动的相关函数矩阵，克服了现有技术因多次利用keystone变换来校正距离徙动和多普勒徙动，造成相参积累能量损失，机动目标的检测概率低的问题，使得本发明能够消除因机动目标的加速度带来的距离徙动和多普勒徙动，完成机动目标回波能量的相参积累，能更好的适应低信噪比的环境，具有更优的机动目标检测性能。
第二，由于本发明通过对相关函数矩阵数据执行快速傅里叶变换和快速傅里叶逆变换来对目标能量进行积分，积累了目标的能量，克服了现有技术需要搜索运动参数，造成计算复杂度高，机动目标的检测速度慢的问题，使得本发明避免了对运动参数的搜索，降低运算复杂度，提高了雷达对机动目标的检测速度，更加利于工程上的实时应用。
附图说明
图1为本发明的流程图；
图2为本发明的仿真图。
具体实施方式
下面结合附图，对本发明做进一步详细描述。
参照图1，对本发明的具体实施步骤做进一步详细描述。
步骤1，获取回波信号。
雷达接收其自身发射的线性调频信号的回波信号。
步骤2，获取解调信号。
利用解调公式，在快时间域内解调雷达接收的回波信号，得到解调信号。
所述的解调公式如下：
其中，h表示解调信号,r表示雷达接收的回波信号，exp表示以自然常数e为底的指数操作，j表示虚数单位符号，π表示圆周率，z表示雷达发射的线性调频信号的载频，tn表示雷达接收回波信号快时间域的采样时间点，tn的取值范围为Tr表示雷达发射的采样间隔为1/y的脉冲串信号重复周期，R0表示雷达和之间的初始径向距离，a1表示微弱机动目标的速度，a2表示微弱机动目标的加速度，tm表示雷达接收回波信号慢时间域的采样时间点，tm的取值范围为(0,TM)，TM表示雷达波束在微弱机动目标上的驻留时间，c表示光速。
步骤3，获取脉冲压缩信号。
利用脉冲压缩公式，在快时间域内对解调信号在距离维进行脉冲压缩，得到脉压信号。
所述的脉冲压缩公式如下：
其中，p表示脉压信号，表示卷积操作，rect(·)表示矩形函数，当时，该矩形函数的值为0，当时，该矩形函数的取值为1，|·|表示取绝对值操作，Tp表示雷达接收的回波信号的脉冲宽度，μ表示雷达接收的回波信号的调频斜率。
步骤4，获取距离向频域信号。
利用快速傅里叶变换公式，在快时间域内对脉压信号做快速傅里叶变换，得到距离向频域信号。
所述的快速傅里叶变换公式如下：
其中，s表示距离向频域信号，∫表示积分操作。
步骤5，按照下式，计算约束条件：
其中，tu表示雷达回波信号慢时间的第u个时刻，u＝1,2,...,M，M表示雷达接收的回波信号的脉冲总数，y表示根据工程需要设定的雷达接收的回波信号的采样频率，z表示雷达发射的线性调频信号的载频，TM表示雷达波束在微弱机动目标上的驻留时间，τv表示雷达回波信号的第v个时延，v＝1,2,...,M。
步骤6，搭建一个的空矩阵，其行和列的取值均与雷达接收的回波信号的脉冲总数相等。
步骤7，从雷达回波信号慢时间中依次选取一个时刻和一个时延，判断所选时刻和所选时延是否同时满足约束条件，若是，则执行步骤8，否则，执行步骤9。
步骤8，利用非均匀重采样技术，计算元素值。
对步骤4得到的距离向频域信号做快速傅里叶变换，得到离散频域信号。
在离散频域信号的频域上大于采样频率y的位置之后补7N个零，N表示雷达接收信号的采样点数，得到补零后的离散频域信号。
将补零后的离散频域信号做快速逆傅里叶变换，得到重采样信号序列。
按照下式，计算第一次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置：
其中，xi,j表示第一次对重采样信号序列插值前，在重采样信号序列的时间轴上对应雷达回波信号慢时间的第i个时刻和雷达回波信号第j个时延非均匀分布的插值位置，ti表示第i个慢时间时刻，f表示距离频域内的距离频率，表示开平方根操作，τj表示雷达回波信号的第j个时延，i和j与所选时刻的序号和所选时延的序号的取值对应相等。
利用线性插值公式，计算第一次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置上待插入的值。
所述的线性插值公式如下：
其中，c(xi,j)表示时间位置xi,j上待插入的值，floor(xi,j)表示重采样信号序列中位于插值位置xi,j后的第一个元素的位置，ceil(xi,j)表示重采样信号序列中位于插值位置xi,j前的第一个元素的位置，c(ceil(xi,j))表示重采样信号序列中位于插值位置xi,j前的第一个元素的值，c(floor(xi,j))表示重采样信号序列中位于插值位置xi,j后的第一个元素的值。
按照下式，计算对重采样信号序列进行第二次插值时，非均匀分布的时间位置：
其中，yp,q表示第二次对重采样信号序列插值前，在重采样信号序列的时间轴上对应雷达回波信号慢时间的第p个时刻和雷达回波的第q个时延的非均匀分布的插值位置，其中p和q与所选时刻的序号和所选时延的序号的取值对应相等。
利用线性插值公式，计算第二次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置上待插入的值。
所述的线性插值公式如下：
其中，c(yp,q)表示时间位置yp,q上待插入的值，floor(yp,q)表示重采样信号序列中位于插值位置yp,q后的第一个元素的位置，ceil(yp,q)表示重采样信号序列中位于插值位置yp,q前的第一个元素的位置，c(ceil(yp,q))表示重采样信号序列中位于插值位置yp,q前的第一个元素的值，c(floor(yp,q))表示重采样信号序列中位于插值位置yp,q后的第一个元素的值。
将两次待插入的值相乘，将乘积作为空矩阵第所选时刻的序号行第所选时延的序号列的元素值。
步骤9，将0作为空矩阵第所选时刻的序号行第所选时延的序号列的元素值。
步骤10，判断所有时刻和时延是否选取完毕，若是，将生成的元素值放入空矩阵中对应的位置得到相关函数矩阵，执行步骤11，否则，执行步骤7。
步骤11，生成快时间-多普勒域相干积累矩阵。
利用快速傅里叶变换公式，在快时间域-慢时间域，对相关函数矩阵做快速傅里叶变换，得到频域相关函数矩阵。
利用快速逆傅里叶变换公式，在距离频域，对频域相关函数矩阵做快速逆傅里叶变换，得到快时间-多普勒域相干积累矩阵。
步骤12，检测机动目标。
利用下述门限公式，计算检测门限值：
其中，η表示检测门限，ln表示以自然常数e为底的对数操作，Pfa表示系统给定的虚警概率，σ表示系统噪声方差。
判断快时间-多普勒域相干积累矩阵的所有元素中绝对值的最大值是否大于或等于检测门限，若是，则将雷达接收的回波信号判定为含有微弱机动目标的目标回波信号，否则，将雷达接收的回波信号判定为无目标的高斯白噪声。
下面结合仿真实验对本发明做进一步的描述。
1.仿真实验条件：
本发明的仿真实验的硬件平台为：处理器为Intel(R)Pentium(R)2CPU，主频为3.0GHz，内存8GB。
本发明的仿真实验的软件平台为：Window XP操作系统和MATLAB 7.5.0。
本发明仿真实验所使用的雷达系统参数如表1所示：
表1雷达系统参数
| 载频 | 带宽 | 采样频率 | 脉冲重复频率 | 脉冲持续时间 | 驻留时间 | 虚警概率 |
| 1.2GHz | 40MHz | 60MHz | 600Hz | 1μs | 0.8s | 10<sup>-4</sup> |
2.仿真内容及结果分析：
本发明的仿真实验是采用本发明的方法和四个现有技术(广义Radon傅里叶变换方法、Radon Lv分布方法、Radon傅里叶变换方法、传统动目标检测方法)分别对信噪比区间为[-20，20]dB的每一个信噪比进行500次Monte-Carlo仿真，将仿真后的结果绘制成机动目标检测概率随信噪比变化对比图如图2所示。
仿真实验中的四个现有技术是指：
现有技术广义Radon傅里叶变换方法是指，J.Xu等人在“Radar maneuveringtarget motion estimation based on generalized radon-Fourier transform”(IEEETrans.Signal Process.,vol.60,no.12,pp.6190–6201,Dec.2012)中提出的机动目标检测方法，简称广义Radon傅里叶变换方法。
现有技术Radon Lv分布方法是指，X.L.Li等人在“Coherent integration formaneuvering target detection based on Radon-Lv’s distribution”(IEEE SignalProcess.Lett.,vol.22,no.9,pp.1467–1471,Sep.2015)中提出的机动目标检测方法，简称Radon Lv分布方法。
现有技术Radon傅里叶变换方法是指，J.Xu等人在“Radon-Fourier transformfor radar target detection”(IEEE Trans.Aerosp.Electron.Syst.,vol.47,no.2,pp.1186–1202,Apr.2011)中提出的机动目标检测方法，简称Radon傅里叶变换方法。
现有技术传统动目标检测方法是指，R.P.Perry等人在“SAR imaging of movingtargets”(IEEE Trans.Aerosp.Electron.Syst.,vol.35,no.1,pp.188–200,Jan.1999)中提出的机动目标检测方法，简称传统动目标检测方法。
图2中的横轴为信噪比，单位dB，纵轴为检测概率，以圆圈标示的曲线表示使用广义Radon傅里叶变换方法进行仿真实验后得到的目标检测概率随信噪比变化的曲线。以方框标示的曲线表示使用Radon Lv分布方法进行仿真实验后得到的目标检测概率随信噪比变化的曲线。以菱形标示的曲线表示使用本发明方法进行仿真实验后得到的目标检测概率随信噪比变化的曲线。以正三角标示的曲线表示使用Radon傅里叶变换方法进行仿真实验后得到的目标检测概率随信噪比变化的曲线。以侧三角标示的曲线表示使用传统动目标检测方法进行仿真实验后得到的目标检测概率随信噪比变化的曲线。
从图2可以看出，在同一个信噪比值下，本发明的目标检测概率远高于现有技术的Radon傅里叶变换方法和传统动目标检测方法，略低于现有技术的广义Radon傅里叶变换方法和Radon Lv分布方法。然而由于本发明对运动参数无需搜索，而现有技术的广义Radon傅里叶变换方法和Radon Lv分布方法，均需要对运动参数进行二维搜索，因此本发明的计算复杂度低，对目标的检测速度快。
综上所述，本发明方法能在计算复杂度和检测概率之间实现了良好的平衡，是一种非常实用的雷达微弱机动目标检测方法。
Claims (7)
1.一种基于非均匀重采样技术的机动目标相干积累检测方法，其特征在于，对距离向频域信号进行非均匀重采样，利用在非均匀重采样后的距离向频域信号生成相关函数矩阵；该方法的具体步骤包括如下：
(1)获取回波信号：
雷达接收其自身发射的线性调频信号的回波信号；
(2)获取解调信号：
利用解调公式，在快时间域内解调雷达接收的回波信号，得到解调信号；
(3)获取脉压信号：
利用脉冲压缩公式，在快时间域内对解调信号在距离维进行脉冲压缩，得到脉压信号；
(4)获取距离向频域信号：
利用快速傅里叶变换公式，在快时间域内对脉压信号做快速傅里叶变换，得到距离向频域信号；
(5)按照下式，计算约束条件：
其中，tu表示雷达回波信号慢时间的第u个时刻，u＝1,2,...,M，M表示雷达接收的回波信号的脉冲总数，y表示根据工程需要设定的雷达接收的回波信号的采样频率，z表示雷达发射的线性调频信号的载频，TM表示雷达波束在微弱机动目标上的驻留时间，τv表示雷达回波信号的第v个时延，v＝1,2,...,M；
(6)搭建一个的空矩阵，其行和列的取值均与雷达接收的回波信号的脉冲总数相等；
(7)从雷达回波信号慢时间中依次选取一个时刻和一个时延，判断所选时刻和所选时延是否同时满足约束条件，若是，则执行步骤(8)，否则，执行步骤(9)；
(8)利用非均匀重采样技术，计算元素值；
(8a)对步骤(4)得到的距离向频域信号做快速傅里叶变换，得到离散频域信号；
(8b)在离散频域信号的频域上大于采样频率y的位置之后补7N个零，N表示雷达接收信号的采样点数，得到补零后的离散频域信号；
(8c)将补零后的离散频域信号做快速逆傅里叶变换，得到重采样信号序列；
(8d)按照下式，计算第一次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置：
其中，xi,j表示第一次对重采样信号序列插值前，在重采样信号序列的时间轴上对应雷达回波信号慢时间的第i个时刻和雷达回波信号第j个时延非均匀分布的插值位置，ti表示第i个慢时间时刻，f表示距离频域内的距离频率，表示开平方根操作，τj表示雷达回波信号的第j个时延，i和j与所选时刻的序号和所选时延的序号的取值对应相等；
(8e)利用线性插值公式，计算第一次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置上待插入的值；
(8f)按照下式，计算第二次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置：
其中，yp,q表示第二次对重采样信号序列插值前，在重采样信号序列的时间轴上对应雷达回波信号慢时间的第p个时刻和雷达回波的第q个时延的非均匀分布的插值位置，其中p和q与所选时刻的序号和所选时延的序号的取值对应相等；
(8g)利用线性插值公式，计算第二次对重采样信号序列插值前，在重采样信号序列的时间轴上非均匀分布的插值位置上待插入的值；
(8h)将两次待插入的值相乘，将乘积作为空矩阵第所选时刻的序号行第所选时延的序号列的元素值；
(9)将0作为空矩阵第所选时刻的序号行第所选时延的序号列的元素值；
(10)判断所有时刻和时延是否选取完毕，若是，将生成的元素值放入空矩阵中对应的位置得到相关函数矩阵，执行步骤(11)，否则，执行步骤(7)；
(11)生成快时间-多普勒域相干积累矩阵：
(11a)利用快速傅里叶变换公式，在快时间域-慢时间域，对相关函数矩阵做快速傅里叶变换，得到频域相关函数矩阵；
(11b)利用快速逆傅里叶变换公式，在距离频域，对频域相关函数矩阵做快速逆傅里叶变换，得到快时间-多普勒域相干积累矩阵；
(12)检测机动目标：
(12a)利用门限公式，计算检测门限值；
(12b)判断快时间-多普勒域相干积累矩阵的所有元素中绝对值的最大值是否大于或等于检测门限，若是，则将雷达接收的回波信号判定为含有微弱机动目标的目标回波信号，否则，将雷达接收的回波信号判定为无目标的高斯白噪声。
2.根据权利要求1所述的基于非均匀重采样技术的机动目标相干积累检测方法，其特征在于：步骤(2)中所述的解调公式如下：
其中，h表示解调信号,r表示雷达接收的回波信号，exp表示以自然常数e为底的指数操作，j表示虚数单位符号，π表示圆周率，tn表示雷达接收回波信号快时间域的采样时间点，tn的取值范围为Tr表示雷达发射的采样间隔为1/y的脉冲串信号重复周期，R0表示雷达和微弱机动目标之间的初始径向距离，a1表示微弱机动目标的速度，a2表示微弱机动目标的加速度，tm表示雷达接收回波信号慢时间域的采样时间点，tm的取值范围为(0,TM)，c表示光速。
3.根据权利要求1所述的基于非均匀重采样技术的机动目标相干积累检测方法，其特征在于：步骤(3)中所述的脉冲压缩公式如下：
其中，p表示脉压信号，表示卷积操作，rect(·)表示矩形函数，当时，该矩形函数的值为0，当时，该矩形函数的取值为1，|·|表示取绝对值操作，Tp表示雷达接收的回波信号的脉冲宽度，μ表示雷达接收的回波信号的调频斜率。
4.根据权利要求1所述的基于非均匀重采样技术的机动目标相干积累检测方法，其特征在于：步骤(4)中所述的快速傅里叶变换公式如下：
其中，s表示距离向频域信号，∫表示积分操作。
5.根据权利要求1所述的基于非均匀重采样技术的机动目标相干积累检测方法，其特征在于：步骤(8e)中所述的线性插值公式如下：
其中，c(xi,j)表示时间位置xi,j上待插入的值，floor(xi,j)表示重采样信号序列中位于插值位置xi,j后的第一个元素的位置，ceil(xi,j)表示重采样信号序列中位于插值位置xi,j前的第一个元素的位置，c(ceil(xi,j))表示重采样信号序列中位于插值位置xi,j前的第一个元素的值，c(floor(xi,j))表示重采样信号序列中位于插值位置xi,j后的第一个元素的值。
6.根据权利要求1所述的基于非均匀重采样技术的机动目标相干积累检测方法，其特征在于：步骤(8g)中所述的线性插值公式如下：
其中，c(yp,q)表示时间位置yp,q上待插入的值，floor(yp,q)表示重采样信号序列中位于插值位置yp,q后的第一个元素的位置，ceil(yp,q)表示重采样信号序列中位于插值位置yp,q前的第一个元素的位置，c(ceil(yp,q))表示重采样信号序列中位于插值位置yp,q前的第一个元素的值，c(floor(yp,q))表示重采样信号序列中位于插值位置yp,q后的第一个元素的值。
7.根据权利要求1所述的基于非均匀重采样技术的机动目标相干积累检测方法，其特征在于：步骤(12a)中所述的门限公式如下：
其中，η表示检测门限，ln表示以自然常数e为底的对数操作，Pfa表示系统给定的虚警概率，σ表示系统噪声方差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910440085.8A CN110095766B (zh) | 2019-05-24 | 2019-05-24 | 基于非均匀重采样技术的机动目标相干积累检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910440085.8A CN110095766B (zh) | 2019-05-24 | 2019-05-24 | 基于非均匀重采样技术的机动目标相干积累检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110095766A true CN110095766A (zh) | 2019-08-06 |
| CN110095766B CN110095766B (zh) | 2023-03-21 |
Family
ID=67449120
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910440085.8A Active CN110095766B (zh) | 2019-05-24 | 2019-05-24 | 基于非均匀重采样技术的机动目标相干积累检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110095766B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112114296A (zh) * | 2020-09-18 | 2020-12-22 | 王玉冰 | 用于无人机协同tdoa/fdoa复合定位的参数估计方法及系统 |
| CN113064131A (zh) * | 2021-03-18 | 2021-07-02 | 江苏蛮酷科技有限公司 | 雷达波处理方法、装置和计算机存储介质 |
| CN114690128A (zh) * | 2022-03-28 | 2022-07-01 | 西安电子科技大学 | 基于多普勒频率相关消除的等离子体鞘套包覆目标相参积累方法 |
| CN114817675A (zh) * | 2022-06-20 | 2022-07-29 | 中国电子科技集团公司第五十四研究所 | 利用广域随机稀布阵波束旁瓣特征的目标快速搜索方法 |
| CN116582155A (zh) * | 2023-05-26 | 2023-08-11 | 北京航天长征飞行器研究所 | 基于直扩的能量积累策略跳频信号捕获方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63284481A (ja) * | 1987-05-18 | 1988-11-21 | Oki Electric Ind Co Ltd | パツシブレンジングシステム |
| JP2011149805A (ja) * | 2010-01-21 | 2011-08-04 | Toshiba Corp | レーダ装置、飛翔体誘導装置及び目標検出方法 |
| JP2012098257A (ja) * | 2010-11-05 | 2012-05-24 | Mitsubishi Electric Corp | レーダ装置 |
| CN103323829A (zh) * | 2013-06-04 | 2013-09-25 | 中国人民解放军海军航空工程学院 | 基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法 |
| CN103399310A (zh) * | 2013-08-07 | 2013-11-20 | 中国人民解放军海军航空工程学院 | 基于相位差分Radon-Lv分布的雷达微弱动目标检测方法 |
| CN104375128A (zh) * | 2014-11-07 | 2015-02-25 | 电子科技大学 | 一种基于互相关函数的高机动目标快速积累检测方法 |
| CN104950295A (zh) * | 2015-05-05 | 2015-09-30 | 电子科技大学 | 一种基于相关函数与尺度变换的高速机动目标检测方法 |
| CN109521410A (zh) * | 2018-11-16 | 2019-03-26 | 西安电子科技大学 | 基于时间反转变换的高速机动目标相参积累检测方法 |
| CN109613507A (zh) * | 2018-12-21 | 2019-04-12 | 北京理工大学 | 一种针对高阶机动目标雷达回波的检测方法 |
-
2019
- 2019-05-24 CN CN201910440085.8A patent/CN110095766B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63284481A (ja) * | 1987-05-18 | 1988-11-21 | Oki Electric Ind Co Ltd | パツシブレンジングシステム |
| JP2011149805A (ja) * | 2010-01-21 | 2011-08-04 | Toshiba Corp | レーダ装置、飛翔体誘導装置及び目標検出方法 |
| JP2012098257A (ja) * | 2010-11-05 | 2012-05-24 | Mitsubishi Electric Corp | レーダ装置 |
| CN103323829A (zh) * | 2013-06-04 | 2013-09-25 | 中国人民解放军海军航空工程学院 | 基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法 |
| CN103399310A (zh) * | 2013-08-07 | 2013-11-20 | 中国人民解放军海军航空工程学院 | 基于相位差分Radon-Lv分布的雷达微弱动目标检测方法 |
| CN104375128A (zh) * | 2014-11-07 | 2015-02-25 | 电子科技大学 | 一种基于互相关函数的高机动目标快速积累检测方法 |
| CN104950295A (zh) * | 2015-05-05 | 2015-09-30 | 电子科技大学 | 一种基于相关函数与尺度变换的高速机动目标检测方法 |
| CN109521410A (zh) * | 2018-11-16 | 2019-03-26 | 西安电子科技大学 | 基于时间反转变换的高速机动目标相参积累检测方法 |
| CN109613507A (zh) * | 2018-12-21 | 2019-04-12 | 北京理工大学 | 一种针对高阶机动目标雷达回波的检测方法 |
Non-Patent Citations (6)
| Title |
|---|
| M. SIMEUNOVIC: "Non-uniform sampled cubic phase function", 《SIGNAL PROCESS》 * |
| PETER O"SHEA: "Improving Polynomial Phase Parameter Estimation by Using Nonuniformly Spaced Signal Sample Methods", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING 》 * |
| 关键等: "雷达高速高机动目标长时间相参积累检测方法", 《信号处理》 * |
| 周莹: "高速高机动雷达目标检测与参数估计算法研究", 《中国硕士学位论文全文数据库 信息科技辑》 * |
| 左衍琴等: "频率步进雷达高速目标运动补偿新方法", 《雷达科学与技术》 * |
| 田静: "雷达机动目标长时间积累信号处理算法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112114296A (zh) * | 2020-09-18 | 2020-12-22 | 王玉冰 | 用于无人机协同tdoa/fdoa复合定位的参数估计方法及系统 |
| CN112114296B (zh) * | 2020-09-18 | 2024-04-16 | 王玉冰 | 用于无人机协同tdoa/fdoa复合定位的参数估计方法及系统 |
| CN113064131A (zh) * | 2021-03-18 | 2021-07-02 | 江苏蛮酷科技有限公司 | 雷达波处理方法、装置和计算机存储介质 |
| CN113064131B (zh) * | 2021-03-18 | 2023-12-26 | 江苏蛮酷科技有限公司 | 雷达波处理方法、装置和计算机存储介质 |
| CN114690128A (zh) * | 2022-03-28 | 2022-07-01 | 西安电子科技大学 | 基于多普勒频率相关消除的等离子体鞘套包覆目标相参积累方法 |
| CN114817675A (zh) * | 2022-06-20 | 2022-07-29 | 中国电子科技集团公司第五十四研究所 | 利用广域随机稀布阵波束旁瓣特征的目标快速搜索方法 |
| CN114817675B (zh) * | 2022-06-20 | 2022-09-09 | 中国电子科技集团公司第五十四研究所 | 利用广域随机稀布阵波束旁瓣特征的目标快速搜索方法 |
| CN116582155A (zh) * | 2023-05-26 | 2023-08-11 | 北京航天长征飞行器研究所 | 基于直扩的能量积累策略跳频信号捕获方法 |
| CN116582155B (zh) * | 2023-05-26 | 2024-03-12 | 北京航天长征飞行器研究所 | 基于直扩的能量积累策略跳频信号捕获方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110095766B (zh) | 2023-03-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110095766B (zh) | 基于非均匀重采样技术的机动目标相干积累检测方法 | |
| CN107561508B (zh) | 一种用于匀加速运动目标的相参积累检测方法 | |
| CN109061589B (zh) | 随机跳频雷达的目标运动参数估计方法 | |
| CN107843892B (zh) | 一种基于最小二乘法的高速目标多普勒测速方法 | |
| CN107132534B (zh) | 一种高速雷达目标频域检测的优化方法 | |
| CN109613506B (zh) | 一种随机跳频重频捷变雷达目标回波信号检测方法 | |
| CN110161477B (zh) | 基于多变量重采样相关函数的机动目标检测方法 | |
| CN102353953A (zh) | 单通道sar地面运动目标成像方法 | |
| Huang et al. | Radar high speed small target detection based on keystone transform and linear canonical transform | |
| CN111505599B (zh) | 一种基于精确距离演化模型的相参积累检测方法 | |
| CN113504522B (zh) | 一种基于发射天线随机切换的时空解耦与超分辨测角方法 | |
| CN109001671B (zh) | 一种跳频信号的目标检测和参数估计方法及装置 | |
| CN111580063B (zh) | 基于广义解调频-楔形变换的雷达目标检测方法 | |
| CN110954885A (zh) | 基于sbl的捷变频雷达自适应目标重构方法 | |
| CN108919221A (zh) | 一种用于变加速运动目标的相参积累检测方法 | |
| CN109061625B (zh) | 一种距离速度匹配的运动目标检测方法 | |
| CN116449320A (zh) | 频率捷变雷达体制下的长时间积累与参数估计方法 | |
| CN110426701A (zh) | 一种基于速度维稀疏表示的捷变频雷达目标速度测量方法 | |
| CN117970255A (zh) | 车载毫米波雷达干扰抑制方法、装置及系统 | |
| CN114545351A (zh) | 基于距离频率轴反转变换和二阶wvd的机动目标相参检测方法及系统 | |
| CN109061626A (zh) | 一种步进频相参处理检测低信杂比动目标的方法 | |
| Guner et al. | Implementation aspects of Wigner-Hough Transform based detectors for LFMCW signals | |
| CN114384484B (zh) | 一种基于分段处理的匀加速运动目标快速相参积累方法 | |
| CN109541567A (zh) | 基于深度学习的高速机动目标检测方法 | |
| CN109901155A (zh) | 多普勒信息反馈调节的高超声速目标修正prpft检测方法 |
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 |