CN106597402B - 一种基于统计规律的自适应动目标显示方法 - Google Patents
一种基于统计规律的自适应动目标显示方法 Download PDFInfo
- Publication number
- CN106597402B CN106597402B CN201611037810.XA CN201611037810A CN106597402B CN 106597402 B CN106597402 B CN 106597402B CN 201611037810 A CN201611037810 A CN 201611037810A CN 106597402 B CN106597402 B CN 106597402B
- Authority
- CN
- China
- Prior art keywords
- frequency
- result
- mti
- moment
- complex signal
- 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.)
- Active
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
- G01S7/414—Discriminating targets with respect to background clutter
-
- 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/415—Identification of targets based on measurements of movement associated with the target
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
技术领域
本发明属于脉冲压缩雷达技术领域,涉及自适应动目标显示,可用于雷达信号处理,涉及一种基于统计规律的自适应动目标显示方法。
背景技术
动目标显示技术是一种可以将目标信号显示杂波区抑制的技术。其主要是通过回波的多普勒频移来判断是静止杂波还是运动的真实目标。也就是说早期的动目标显示技术是只能对静止的杂波或者目标有抑制作用,对于具有一定运动速度的运动杂波没有效果,在频域上分析就是只消除零频附近的杂波,对于具有一定多普勒频移的运动杂波是没有办法的。
随着技术的进步,谱估计理论又被提出。人们设想首先设计出一组动目标显示(MTI)滤波器,使得该组滤波器可以均匀覆盖检测频率上的所有频率,然后对运动杂波谱完成估计,选择相应频率的滤波器,将杂波进行滤除。该技术整体处理过程中,滤波器系数可以提前设计好,最后影响系统性能的一般是谱估计的准确性。如果谱估计技术可以简单准确的估计运动杂波的频率,就可以使得整个雷达系统在性能良好的基础上,快速简单的进行实现。但运动杂波频点在一定距离范围内满足统计规律,但单个距离点上没有明确的规律,由此使得运动杂波中心频率的估计不易得到。权系数法就是以上描述的方法中具体的一种,但该方法并未使用统计规律对运动杂波进行估计,导致运动杂波谱中心估计误差较大。
发明内容
(一)发明目的
本发明的目的是:针对上述自适应动目标显示方法中存在的问题,提出一种基于统计规律的自适应动目标显示方法,以准确确定运动杂波频率,选择合适的滤波器完成运动杂波的滤除。
(二)技术方案
为了解决上述技术问题,本发明提供一种基于统计规律的自适应动目标显示方法,其包括以下步骤:
步骤1:根据雷达系统参数,设计一组零频MTI滤波器系数W(f0),f0=0Hz,W(f0)的阶数为n;
步骤2:通过旋转因子得到中心频率为fj的MTI滤波器组系数W(fj)=HW(f0),H=diag{exp[j2πfjT,j2πfj2T,j2πfj3T,...,j2πfjnT]},T为雷达的脉冲重复周期;
步骤3:对t1时刻及以前n-1个时刻进行脉冲压缩处理后的复信号进行零频MTI滤波处理得重复上述操作,得到t2,...,tm时刻零频MTI滤波处理结果 tm=t1+mT,m为雷达天线一个扇区所包含的脉冲回波个数;
其中,为脉冲压缩雷达信号处理过程中脉冲压缩处理后的复信号;
步骤4:对进行复信号求模处理,得到 对进行相同距离门的非相参积累求和处理得到i表示脉冲回波信号的距离门号,p为最大的距离门号;
步骤5:对进行单元平均恒虚警检测,得到检测结果目标距离门信息ND1,目标个数为
步骤6:将以步骤5中的目标距离门为分割点,分为段,并进行分段相邻两个脉冲回波信号的复相关运算;
步骤7:对复相关运算结果进行反正切运算;
步骤8:对第一段F11,F21,…,Fm1内的频率值按出现次数进行统计,记录出现次数最多的频率值fc1,同样可以得到各段内出现次数最多的频率值根据步骤2选择相应频率点上的滤波器系数
步骤9:对t1时刻及之后n-1个时刻进行零频MTI处理后的复信号进行分段二级MTI滤波处理;
步骤10:重复步骤9,得到分别对其进行复信号求模处理,并对求模结果进行相应距离门的非相参积累求和处理;
步骤11:对进行单元平均恒虚警检测,得到雷达会波信号中目标的距离们信息,即t1时刻自适应动目标显示方法的结果。
其中,所述步骤4中,复信号求模处理的结果为:
其中,所述步骤6中,所分得的段为:
其中,所述步骤6中,所分得的段中,进行分段相邻两个脉冲回波信号的复相关运算结果为:
其中,所述步骤7中,反正切运算的结果为:
其中,所述步骤9中,所述二级MTI滤波处理结果为:
表示t1时刻由步骤5得到的第一个目标的二级MTI处理后的复信号,表示t1时刻由步骤5得到的第一个目标的零频MTI处理后的复信号,为第段数据。
其中,所述步骤10中,复信号求模处理的结果为:
其中,所述步骤10中,对求模结果进行相同距离门的非相参积累求和处理得到 i表示脉冲回波信号的距离门号,p为最大的距离门号。
(三)有益效果
上述技术方案所提供的基于统计规律的自适应动目标显示方法,在经过零频滤波处理之后,零速目标被滤除;经过分段按统计规律进行频率估计并完成滤波后,两侧较高的杂波被滤除,且两个目标幅度未发生明显改变,该方法准确有效。
附图说明
图1是本发明的总流程图;
图2是本发明仿真时雷达回波信号脉冲压缩处理后的波形图;
图3是本发明仿真时雷达回波信号脉冲压缩处理后的三维图;
图4是本发明仿真时4阶零频滤波器的幅频响应图;
图5是本发明仿真时零频滤波处理后的波形图;
图6是本发明仿真时零频滤波处理后的三维图;
图7是本发明仿真时第一段按统计规律估计得到频点对应的4阶滤波器的幅频响应图;
图8是本发明仿真时第二段按统计规律估计得到频点对应的4阶滤波器的幅频响应图;
图9是本发明仿真时第三段按统计规律估计得到频点对应的4阶滤波器的幅频响应图;
图10是本发明仿真时雷达回波信号自适应动目标显示处理后的波形图;
图11是本发明仿真时雷达回波信号自适应动目标显示处理后的三维图。
具体实施方式
为使本发明的目的、内容和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
参照图1,本发明的实现步骤如下:
步骤1:根据雷达系统参数,设计一组零频MTI滤波器系数W(f0),f0=0Hz,W(f0)的阶数为n;
步骤2:通过旋转因子得到中心频率为fj的MTI滤波器组系数W(fj)=HW(f0),H=diag{exp[j2πfjT,j2πfj2T,j2πfj3T,...,j2πfjnT]},T为雷达的脉冲重复周期;
步骤3:对t1时刻及以前n-1个时刻进行脉冲压缩处理后的复信号进行零频MTI滤波处理得重复上述操作,得到t2,...,tm时刻零频MTI滤波处理结果 tm=t1+mT,m为雷达天线一个扇区所包含的脉冲回波个数;
其中,为脉冲压缩雷达信号处理过程中脉冲压缩处理后的复信号;
步骤4:对进行复信号求模处理:
对进行相同距离门的非相参积累求和处理得到i表示脉冲回波信号的距离门号,p为最大的距离门号;
步骤5:对进行单元平均恒虚警检测,得到检测结果目标距离门信息ND1,目标个数为
步骤6:将以步骤5中的目标距离门为分割点,分为段得到:
分段时注意把目标复信号剔除,进行分段相邻两个脉冲回波信号的复相关运算:
步骤7:对进行反正切运算得到:
步骤8:对第一段F11,F21,…,Fm1内的频率值按出现次数进行统计,记录出现次数最多的频率值fc1,同样可以得到各段内出现次数最多的频率值根据(2)选择相应频率点上的滤波器系数
步骤9:对t1时刻及之后n-1个时刻进行零频MTI处理后的复信号进行分段二级MTI滤波处理得到:
表示t1时刻由(5)得到的第一个目标的二级MTI处理后的复信号,表示t1时刻由(5)得到的第一个目标的零频MTI处理后的复信号,为第段数据;
步骤10:重复步骤(9),得到分别对其进行复信号求模处理:
对求模结果进行相同距离门的非相参积累求和处理得到i表示脉冲回波信号的距离门号,p为最大的距离门号;
步骤11:对进行单元平均恒虚警检测,得到雷达会波信号中目标的距离们信息,即t1时刻自适应动目标显示方法的结果。
本发明的效果通过以下仿真实验进一步说明:
(1)实验场景:雷达回波仿真信号使用线性调频信号,在回波中设定三个目标,第一个为800点处的静止目标,第二个目标为1000点处的运动目标,第三个目标为2750点处的运动目标。雷达的脉冲重复频率为400Hz,杂波区域限定在200-1800点距离范围内。
(2)仿真内容:
仿真1,根据以上要求,限定一个CPI内包含28个PRT,图2给出脉冲压缩后波形图,图3给出脉冲压缩结果的三维图;
仿真2,图4为4阶零频滤波器的幅频响应图;
仿真3,图5为零频滤波处理后的波形图,图6为零频滤波处理结果的三维图;
仿真4,图7为第一段用本发明的方法得到的估计频率f1对应的滤波器的频率响应图,图8为第二段用本发明的方法得到的估计频率f2对应的滤波器的频率响应图,图9为第三段用本发明的方法得到的估计频率f3对应的滤波器的频率响应图;
仿真5,图10为分段不同频率滤波处理后的波形图,图11为分段不同频率滤波处理后的三维图;
(3)仿真结果分析:
从图5可以看出,在经过零频滤波处理之后,零速目标被滤除,剩下两个运动目标,在1000点左右两侧有较高的杂波未被滤除;
从图10可以看出,经过分段按统计规律进行频率估计并完成滤波后,1000点左右两侧较高的杂波被滤除,且两个目标幅度未发生明显改变,由此证明了,本方法有效。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (1)
1.一种基于统计规律的自适应动目标显示方法,其特征在于,包括以下步骤:
步骤1:根据雷达系统参数,设计一组零频MTI滤波器系数W(f0),f0=0Hz,W(f0)的阶数为n;
步骤2:通过旋转因子得到中心频率为fj的MTI滤波器组系数W(fj)=HW(f0),H=diag{exp[j2πfjT,j2πfj2T,j2πfj3T,...,j2πfjnT]},T为雷达的脉冲重复周期;
步骤3:对t1时刻及以前n-1个时刻进行脉冲压缩处理后的复信号进行零频MTI滤波处理得重复上述操作,得到t2,...,tm时刻零频MTI滤波处理结果 tm=t1+mT,m为雷达天线一个扇区所包含的脉冲回波个数;
其中,为脉冲压缩雷达信号处理过程中脉冲压缩处理后的复信号;
步骤4:对进行复信号求模处理,得到 对进行相同距离门的非相参积累求和处理得到i表示脉冲回波信号的距离门号,p为最大的距离门号;
所述步骤4中,复信号求模处理的结果为:
…
步骤5:对进行单元平均恒虚警检测,得到检测结果目标距离门信息ND1,目标个数为
步骤6:将以步骤5中的目标距离门为分割点,分为段,并进行分段相邻两个脉冲回波信号的复相关运算;
所述步骤6中,所分得的段为:
…
所分得的段中,进行分段相邻两个脉冲回波信号的复相关运算结果为:
r11=Y(t1,0)Y(t2,0)T
r12=Y(t1,1)Y(t2,1)T
…
r21=Y(t2,0)Y(t3,0)T
…
…
步骤7:对复相关运算结果进行反正切运算;
所述步骤7中,反正切运算的结果为:
…
步骤8:对第一段F11,F21,…,Fm1内的频率值按出现次数进行统计,记录出现次数最多的频率值fc1,同样可以得到各段内出现次数最多的频率值fc2,…,根据步骤2选择相应频率点上的滤波器系数
步骤9:对t1时刻及之后n-1个时刻进行零频MTI处理后的复信号进行分段二级MTI滤波处理;
所述步骤9中,所述二级MTI滤波处理结果为:
…
表示t1时刻由步骤5得到的第一个目标的二级MTI处理后的复信号,表示t1时刻由步骤5得到的第一个目标的零频MTI处理后的复信号,为第段数据;
步骤10:重复步骤9,得到分别对其进行复信号求模处理,并对求模结果进行相应距离门的非相参积累求和处理;
所述步骤10中,复信号求模处理的结果为:
…
所述步骤10中,对求模结果进行相同距离门的非相参积累求和处理得到 i表示脉冲回波信号的距离门号,p为最大的距离门号;
步骤11:对进行单元平均恒虚警检测,得到雷达回波信号中目标的距离门信息,即t1时刻自适应动目标显示方法的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611037810.XA CN106597402B (zh) | 2016-11-23 | 2016-11-23 | 一种基于统计规律的自适应动目标显示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611037810.XA CN106597402B (zh) | 2016-11-23 | 2016-11-23 | 一种基于统计规律的自适应动目标显示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106597402A CN106597402A (zh) | 2017-04-26 |
CN106597402B true CN106597402B (zh) | 2019-03-19 |
Family
ID=58592874
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611037810.XA Active CN106597402B (zh) | 2016-11-23 | 2016-11-23 | 一种基于统计规律的自适应动目标显示方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106597402B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102387273A (zh) * | 2011-07-08 | 2012-03-21 | 歌尔声学股份有限公司 | 一种抑制残留回声的方法和装置 |
CN102447795A (zh) * | 2010-10-09 | 2012-05-09 | 杭州华三通信技术有限公司 | 一种回波消除的方法和装置 |
CN104502906A (zh) * | 2014-12-09 | 2015-04-08 | 中国民航大学 | 基于rmdcft的空间超高速机动目标检测方法 |
CN103616848B (zh) * | 2013-12-04 | 2016-08-17 | 固高科技(深圳)有限公司 | 改善速度前馈的自适应滤波器的滤波方法和系统 |
-
2016
- 2016-11-23 CN CN201611037810.XA patent/CN106597402B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102447795A (zh) * | 2010-10-09 | 2012-05-09 | 杭州华三通信技术有限公司 | 一种回波消除的方法和装置 |
CN102387273A (zh) * | 2011-07-08 | 2012-03-21 | 歌尔声学股份有限公司 | 一种抑制残留回声的方法和装置 |
CN102387273B (zh) * | 2011-07-08 | 2013-08-28 | 歌尔声学股份有限公司 | 一种抑制残留回声的方法和装置 |
CN103616848B (zh) * | 2013-12-04 | 2016-08-17 | 固高科技(深圳)有限公司 | 改善速度前馈的自适应滤波器的滤波方法和系统 |
CN104502906A (zh) * | 2014-12-09 | 2015-04-08 | 中国民航大学 | 基于rmdcft的空间超高速机动目标检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106597402A (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | Radar high-speed target detection based on the frequency-domain deramp-keystone transform | |
CN109814073B (zh) | 一种mtd雷达解模糊测速方法 | |
EP0932837B1 (en) | Procedure for the elimination of interference in a radar unit of the fmcw type | |
CN204495996U (zh) | 宽带声纳接收器 | |
CN106814353A (zh) | 一种雷达信号处理系统 | |
CN106597429A (zh) | 一种基于雷达信号处理系统的速度跟踪子系统 | |
CN106772302A (zh) | 一种复合高斯背景下的知识辅助stap检测方法 | |
CN111830480B (zh) | 一种雷达海杂波短时谱特征参数估计方法及系统 | |
CN108196241B (zh) | 一种基于Hough变换的高速动目标速度估计方法 | |
CN104215959A (zh) | 一种多机动目标径向初速度和径向加速度的估计方法 | |
RU2579353C1 (ru) | Способ сопровождения воздушной цели из класса "самолёт с турбореактивным двигателем" при воздействии уводящей по скорости помехи | |
CN107367715A (zh) | 基于稀疏表示的杂波抑制方法 | |
CN109100696A (zh) | 基于点状杂波图的慢速运动目标消除方法 | |
CN112782685B (zh) | 基于mimo雷达的多声源定位与声音重构方法及系统 | |
Varghese et al. | Characterization of tissue microstructure scatterer distribution with spectral correlation | |
CN104698446A (zh) | 基于动目标检测的雷达回波滤波方法 | |
CN107329115B (zh) | 基于gcrbf网络的lfm信号参数估计方法 | |
CN106199539A (zh) | 基于白化滤波器的地杂波抑制方法 | |
CN106597402B (zh) | 一种基于统计规律的自适应动目标显示方法 | |
KR100936173B1 (ko) | 표적 추적장치에서의 도플러 주파수를 추출하는 필터뱅크및 그 방법 | |
CN110596651B (zh) | 一种雷达检测的方法 | |
Tian et al. | Multiple target detection using split spectrum processing and group delay moving entropy | |
CN109061626A (zh) | 一种步进频相参处理检测低信杂比动目标的方法 | |
Iverson | Coherent processing of ultra-wideband radar signals | |
Altesmber et al. | Doppler-tolerant classification of distributed targets--a bionic sonar |
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 |