发明名称 一种基于追踪滤波器的雷达数据动目标检测处理方法
摘要 本发明公开了一种基于追踪滤波器的雷达数据动目标检测处理方法,包括以下步骤:(1)对脉压后信号经过信号预处理;(2)对步骤1得到的信号进行滤波,分为检测和追踪两种状态:检测状态采用M路MTD滤波器输出M路多普勒频率通道的目标信号,追踪状态采用追踪滤波器对特定多普勒频率信号进行滤波;(3)对步骤(2)的输出数据,进行杂波图处理;(4)对步骤(3)的结果,进行恒虚警处理;(5)根据杂波分布,对杂波特性进行分析计算,对杂波区分类,对应各种杂波区采用对应的检测门限,进行目标检测;(6)将各个通道相同距离元的检测信号进行合成,获取距离单元的目标检测结果。
申请公布号 CN103885053B 申请公布日期 2016.07.06
申请号 CN201410120474.X 申请日期 2014.03.27
申请人 中国电子科技集团公司第二十八研究所 发明人 龙超;林学;黄绍斌;赵春光
分类号 G01S7/41(2006.01)I;G01S7/36(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 江苏圣典律师事务所 32237 代理人 胡建华
主权项 一种基于追踪滤波器的雷达数据动目标检测处理方法,其特征在于,包括以下步骤:(1)对脉压信号预处理,去除信号中的干扰和杂波;(2)对步骤1得到的信号进行滤波,分为检测和追踪两种状态:检测状态采用M路MTD滤波器输出M路多普勒频率通道的目标信号,追踪状态采用追踪滤波器对特定多普勒频率信号进行滤波;(3)对步骤(2)的输出数据,进行杂波图处理;(4)对步骤(3)的结果,进行恒虚警处理;(5)根据杂波分布,对杂波特性进行分析计算,对杂波区分类,对应各种杂波区采用对应的检测门限,进行目标检测;目标检测数据再经过滑窗检测,根据区域类型,分别设置滑窗检测门限;(6)将各个通道相同距离元的检测信号进行合成,获取距离单元的目标检测结果;步骤(2)中,追踪滤波器对特定多普勒频率信号进行滤波包括,将追踪滤波器的最大零点对准目标信号,并估计目标信号的多普勒频率中心;追踪滤波器估计目标的多普勒频率中心时:如果追踪滤波器处于追踪状态,利用目标检测的速度变化估计多普勒频率中心;如果追踪滤波器处于检测状态,利用经典谱估计算法周期图法,或质心算法估计多普勒频率中心;步骤(2)中,追踪状态下采用以下公式估计多普勒频率中心:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>f</mi><mi>c</mi></msub><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>V</mi><mi>r</mi></msub></mrow><mi>&lambda;</mi></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000938136370000011.GIF" wi="214" he="126" /></maths>其中,f<sub>c</sub>是多普勒频率,V<sub>r</sub>是目标径向速度,λ是雷达的工作波长;步骤(2)中,检测状态下质心算法估计多普勒频率中心包括以下步骤:(a)对于步骤(1)得到的脉压信号p<sub>m</sub>(n):p<sub>m</sub>(n),n=0,1,…,N‑1,m=m<sub>1</sub>,m<sub>1</sub>+1,…,m<sub>2</sub>,m是参与信号功率谱估计的距离单元,m<sub>1</sub>是距离起始单元,m<sub>2</sub>是距离终止单元,N为脉压信号的长度,n表示对应脉压信号;(b)计算功率谱S<sub>m</sub>(n):S<sub>m</sub>(n)=|FFT[s<sub>m</sub>(n)]|<sup>2</sup>;(c)计算所有距离单元的功率谱平均值S(n):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>S</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>m</mi><mn>2</mn></msub><mo>-</mo><msub><mi>m</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><msub><mi>m</mi><mn>1</mn></msub><msub><mi>m</mi><mn>2</mn></msub></munderover><msub><mi>S</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000938136370000021.GIF" wi="661" he="158" /></maths>(d)计算质心矢量对应的相位φ<sub>0</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&phi;</mi><mn>0</mn></msub><mo>=</mo><mi>p</mi><mi>h</mi><mi>a</mi><mi>s</mi><mi>e</mi><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><msub><mi>m</mi><mn>1</mn></msub><msub><mi>m</mi><mn>2</mn></msub></munderover><mrow><mo>(</mo><msub><mi>S</mi><mi>m</mi></msub><mo>(</mo><mi>n</mi><mo>)</mo><mi>exp</mi><mo>(</mo><mrow><mi>j</mi><mo>&CenterDot;</mo><mn>2</mn><mi>&pi;</mi><mi>n</mi><mo>/</mo><mi>N</mi></mrow><mo>)</mo><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo></mrow>]]></math><img file="FDA0000938136370000022.GIF" wi="1078" he="207" /></maths>φ<sub>0</sub>是质心矢量对应的相位,N是脉压数据长度,Phase是对复信号求相位,j是虚部单位,exp表示求指数运算;(e)估计频率中心f<sub>c</sub>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>f</mi><mi>c</mi></msub><mo>=</mo><msub><mi>&phi;</mi><mn>0</mn></msub><mo>&CenterDot;</mo><mfrac><msub><mi>f</mi><mi>p</mi></msub><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000938136370000023.GIF" wi="350" he="173" /></maths>其中,f<sub>p</sub>是雷达脉冲重复频率,是采集信号的雷达系统参数;对信号先进行加权自相关再进行FFT处理。
地址 210007 江苏省南京市苜蓿园东街1号1406信箱07分箱