-
目标检测及运动参数估计是雷达系统的一个基础任务,其精确性对战争的胜利起到重要的作用。对雷达信号进行长时间积累,以增加实际可利用的信号能量,是提高检测性能和参数估计精确性的有效手段,然而长时间信号积累用于高速、多目标的检测和运动参数估计会受到距离徙动、多普勒扩散以及多普勒模糊的影响[1]。改进Hough变换算法[2]利用非相参积累方式进行目标检测,其检测效果受高速目标距离弯曲程度和信噪比的影响严重。基于长时间的相参积累方式,文献[3-4]提出了Radon-Fourier变换算法,然而距离弯曲及多普勒扩散导致算法很难对高速目标的信号能量进行有效的积累。针对距离徙动及多普勒扩散,文献[5]提出了Keystone变换结合解线频调的算法,然而Keystone变换受制于多普勒模糊。文献[6]利用速度模糊数搜索的算法可以解决多普勒模糊,但未涉及到多普勒扩散问题。针对距离徙动、多普勒扩散及多普勒模糊,文献[7-8]提出了基于循环平稳理论的迭代算法,其计算复杂度低,相比传统循环平稳算法运动参数估计范围有所提高,然而工程应用中算法估计精度低和参数估计范围受限的缺陷限制了其应用。
针对以上问题,该文提出了一种新的高速多目标窄带雷达目标检测及参数估计算法。在多目标及存在距离徙动、多普勒扩散和多普勒模糊的情况下,该算法首先采用基于循环平稳的联合Keystone变换与模糊数搜索方法完成运动参数的粗估计,然后在粗估计基础上利用联合频域距离徙动补偿处理与分数阶傅里叶变换(fraction Fourier transform, FRFT)方法完成高速多目标的参数化检测。本文算法保持了循环平稳复杂度的优点,且克服了已有循环平稳算法工程应用中参数估计精度低和运动参数估计范围有限的缺陷。计算机仿真和实测数据验证了算法的有效性。
-
假设脉冲多普勒(PD)雷达发射的第$ m $个线性调频脉冲信号的基带形式如下:
$$ x(\hat t, {t_m}) = {\mathop{\rm rect}\nolimits} \left( {\frac{{\hat t}}{{{T_p}}}} \right)\exp ({\rm{j\pi }}\gamma {\hat t^2}) $$ (1) 式中,;$ \hat t $为发射脉冲的时间变量,为快时间;$ m $为发射脉冲序列索引;$ {t_m} $为慢时间;$ {T_p} $为发射脉冲宽度;$ \gamma $为调频斜率。
高速多目标的雷达回波信号基带形式为:
$$ \begin{array}{c} s(\hat t, {t_m}) = \sum\limits_i {{A_i}{\mathop{\rm rect}\nolimits} \left[{\frac{{\hat t-\frac{{2{R_i}({t_m})}}{c}}}{{{T_p}}}} \right]} \times \\ \exp \left[{{\rm{j\pi }}\gamma {{\left( {\hat t-\frac{{2{R_i}({t_m})}}{c}} \right)}^2}} \right] \times \\ \exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{R_i}({t_m})}}{c}} \right] + n(\hat t, {t_m}) \end{array} $$ (2) 式中,$ {A_i} $为第$ i $个目标的回波幅度;慢时间$ {t_m} $时刻目标与雷达的径向距离为$ {R_i}({t_m}) = {R_{0i}} + {v_{0i}}{t_m} + $ $ \frac{1}{2}{a_i}t_m^2 $,$ {R_{0i}} $,$ {v_{0i}} $以及$ {a_i} $分别为目标的径向初始距离,径向初始速度以及径向加速度;$ {f_c} $为雷达载频;$ n(\hat t, {t_m}) $为加性高斯白噪声。
式(2) 对应的频域表示形式为:
$$ \begin{array}{*{20}{c}} {S\left( {f, {t_m}} \right) = \sum\limits_i {{A_i}{\mathop{\rm rect}\nolimits} \left( {\frac{f}{{\gamma {T_p}}}} \right)\exp \left( { - {\rm{j\pi }}\frac{{{f^2}}}{\gamma }} \right) \times } }\\ {\exp \left( { - {\rm{j}}2{\rm{\pi }}f\frac{{2{R_{0i}}}}{c}} \right)\exp \left[{-{\rm{j}}2{\rm{\pi }}f\left( {\frac{{2{v_{0i}}{t_m}}}{c} + \frac{{{a_i}t_m^2}}{c}} \right)} \right] \times }\\ {\exp \left( { - {\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{R_{0i}}}}{c}} \right)\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{v_{0i}}{t_m}}}{c}} \right] \times }\\ {\exp \left( { - {\rm{j}}2{\rm{\pi }}{f_c}\frac{{{a_i}t_m^2}}{c}} \right) + n\left( {f, {t_m}} \right)} \end{array} $$ (3) 利用$ X(f, {t_m}) = {\mathop{\rm rect}\nolimits} \left( {\frac{f}{{\gamma {T_p}}}} \right)\exp \left( { - {\rm{j\pi }}\frac{{{f^2}}}{\gamma }} \right) $匹配滤波器对式(3) 信号进行脉冲压缩处理:
$$ \begin{array}{*{20}{c}} {{X_{RM}}\left( {f, {t_m}} \right) = \sum\limits_i {{A_i}{\mathop{\rm rect}\nolimits} \left( {\frac{f}{{\gamma {T_p}}}} \right)\exp \left( { - {\rm{j}}2{\rm{\pi }}f\frac{{2{R_{0i}}}}{c}} \right)} \times }\\ {\exp \left[{-{\rm{j}}2{\rm{\pi }}f\frac{{2{v_{0i}}{t_m} + {a_i}t_m^2}}{c}} \right]\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{R_{0i}}}}{c}} \right] \times }\\ {\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{v_{0i}}{t_m}}}{c}} \right]\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{{a_i}t_m^2}}{c}} \right] + n(f, {t_m})} \end{array} $$ (4) 式中,第一个指数项表示目标初始时刻位置;第二个指数项表示目标由于速度和加速度产生的距离徙动;第四个指数项表示目标的多普勒频移,第五个指数项表示目标的多普勒扩散。
-
循环平稳理论在信号处理领域的一个重要应用方向是目标检测以及信号参数和波形的估计,其信号与干扰分离效果好,计算复杂度低,然而容易受到多普勒模糊的影响[7-8]。对式(4) 信号进行慢时间维的自相关处理:
$$ \begin{array}{*{20}{c}} {U(u, f, {t_m}) = {X_{RM}}(f, {t_m}){X_{RM}}^ * (f, {t_m} + u) = }\\ {{U_{ii}}(u, f, {t_m}) + {U_{ij}}(u, f, {t_m}) + {U_n}(u, f, {t_m})} \end{array} $$ (5) 式中,
$$ \begin{array}{c} {U_{ii}}(u, f, {t_m}) = \sum\limits_i {{A_i}^2{\mathop{\rm rect}\nolimits} \left( {\frac{f}{{\gamma {T_p}}}} \right)} \times \\ \exp \left[{{\rm{j}}2{\rm{\pi }}(f + {f_c})\frac{{2{a_i}u{t_m}}}{c}} \right] \times \\ \exp \left[{{\rm{j}}2{\rm{\pi }}(f + {f_c})\frac{{2{v_{0i}}u + {a_i}{u^2}}}{c}} \right] \end{array} $$ (6) $$ \begin{array}{*{20}{c}} {{U_{ij}}(u, f, {t_m}) = \sum\limits_i {\sum\limits_{j, i \ne j} {{A_i}{A_j}{\mathop{\rm rect}\nolimits} \left( {\frac{f}{{\gamma {T_p}}}} \right)} } \times }\\ {\exp \left[{-{\rm{j\pi }}(f + {f_c})\frac{{4{R_{0i}}-4{R_{0j}}-2{v_{0j}}u - {a_j}{u^2}}}{c}} \right] \times }\\ {\exp \left[{{\rm{j}}2{\rm{\pi }}(f + {f_c})\frac{{2{v_{0j}}-2{v_{0i}} + 2{a_j}u}}{c}{t_m}} \right] \times }\\ {\exp \left[{{\rm{j}}2{\rm{\pi }}(f + {f_c})\frac{{{a_j}-{a_i}}}{c}{t_m}^2} \right]} \end{array} $$ (7) 式中,$ {U_{ii}}(u, f, {t_m}) $为目标信号自相关项,存在目标加速度所造成的距离徙动;$ {U_{ij}}(u, f, {t_m}) $为目标信号的互相关项,存在目标速度和加速度造成的距离徙动及目标加速度产生的多普勒扩散;$ u $代表循环延迟,根据文献[7-8]:延迟$ u = \frac{T}{{\rm{2}}} $时估计性能最优,其中$ T $为信号积累总时间;$ {U_n}(u, f, {t_m}) $为信号自相关后的噪声项,能量集中在零循环频率处[7-8]。对比$ {U_{ii}}(u, f, {t_m}) $与$ {U_{ij}}(u, f, {t_m}) $可以得出:即使$ {U_{ii}}(u, f, {t_m}) $中第$ i $个目标的距离徙动得到补偿,$ {U_{ij}}(u, f, {t_m}) $中每一个互相关项对应的距离徙动及多普勒扩散仍旧存在,且徙动与扩散严重,因此可以认为$ {U_{ij}}(u, f, {t_m}) $对于$ {U_{ii}}(u, f, {t_m}) $相当于干扰。基于上述分析,后续信号处理将针对信号自相关项$ {U_{ii}}(u, f, {t_m}) $距离徙动的补偿。
将$ {U_{ii}}(u, f, {t_m}) $中第$ i $个目标对应的多普勒频率$ {f_{di}} = \frac{{2{a_i}u}}{\lambda } = {f_{di0}} + {n_{ki}}{\rm{prf}} $代入式(6),可得:
$$ \begin{array}{c} {{U'}_{ii}}(u, f, {t_m}) = \sum\limits_i {{A_i}^2{\mathop{\rm rect}\nolimits} \left( {\frac{f}{{\gamma {T_p}}}} \right) \times } \\ \exp \left[{{\rm{j}}2{\rm{\pi }}\frac{{f + {f_c}}}{{{f_c}}}{f_{di0}}{t_m}} \right]\exp \left[{{\rm{j}}2{\rm{\pi }}\frac{f}{{{f_c}}}{n_{ki}}{\rm{prf}}{t_m}} \right] \times \\\exp \left[{{\rm{j}}2{\rm{\pi }}(f + {f_c})\frac{{2{v_{0i}}u + {a_i}{u^2}}}{c}} \right] \end{array} $$ (8) 式中,$ \lambda $为波长;$ {f_{di0}} $为$ {R_{ii}}(u, f, {t_m}) $对应的第$ i $个目标折叠以后的多普勒频率;$ {n_{ki}} $为$ {U_{ii}}(u, f, {t_m}) $对应的第$ i $个目标的多普勒模糊数;$ {\rm{prf}} $为脉冲重复频率。式(8) 信号存在距离徙动及多普勒模糊,为完成目标参数化检测,需要进行距离徙动补偿及解多普勒模糊。
Keystone变换不需要目标的先验信息,即可完成距离徙动补偿,在雷达信号处理中得到了广泛的应用,然而Keystone变换本身受限于多普勒模糊,如果目标发生模糊,则需要进行解模糊处理[6]。
因为$ {f_c} > > f $,则$ \frac{{{f_c}}}{{{f_c} + f}} \approx 1 $,那么进行Keystone变换消除$ {U'_{ii}}(u, f, {t_m}) $的距离徙动,令$ {t_m} = \frac{{{f_c}}}{{f + {f_c}}}{\tau _m} $,可得:
$$ {K_{ii}}(u, f, {t_m}) \approx \sum\limits_i {{A_i}^2{\mathop{\rm rect}\nolimits} \left( {\frac{f}{{\gamma {T_p}}}} \right)} \exp [{\rm{j}}2{\rm{\pi }}{f_{di0}}{\tau _m}] \times \\\exp \left( {{\rm{j}}2{\rm{\pi }}f\frac{{{n_{ki}}{\rm{prf}}}}{{{f_c}}}{\tau _m}} \right)\exp \left[{{\rm{j}}2{\rm{\pi }}(f + {f_c})\frac{{2{v_{0i}}u + {a_i}{u^2}}}{c}} \right] $$ (9) 根据式(9) 信号形式,定义解模糊操作函数:
$$ H(f, {\tau _m}, {n'_{ki}}) = \exp \left[{-{\rm{j}}2{\rm{\pi }}f\frac{{{{n'}_k}{\rm{prf}}}}{{{f_c}}}{\tau _m}} \right] $$ (10) 对式(9) 信号进行解模糊处理,并转换到时域:
$$ \begin{array}{*{20}{c}} {T({{n'}_k}, u, \hat t, {\tau _m}) = \sum\limits_i {{A_i}^2\gamma {T_p}\exp \left[{{\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{v_{0i}}u + {a_i}{u^2}}}{c}} \right] \times } }\\ {\sin c\left[{\gamma {T_p}\left( {\hat t + \frac{{2{v_{0i}}u + {a_i}{u^2}}}{c} + \frac{{({n_{ki}}-{{n'}_k}){\rm{prf}}{\tau _m}}}{{{f_c}}}} \right)} \right] \times } \end{array} $$ (11) 式中,$ {n'_k} $为模糊数范围中的取值,若$ {n_{ki}} - {n'_k} = 0 $,则目标不再存在距离徙动,反之目标存在距离徙动。
根据设定的模糊数范围,通过代价函数式(12) 可以计算出目标的运动参数:
$$ ({t'_i}, {f'_{di0}}, {n'_{ki}}) = \mathop {\arg \max }\limits_{{{\hat n}_k}} \left| {{\rm{FFT}}[T({{n'}_k}, u, \hat t, {\tau _m})]} \right| $$ (12) 式中,$ {t'_i} = \frac{{2{v_{0i}}u + {a_i}{u^2}}}{c} $为第$ i $个目标峰值的快时间时刻;$ {f'_{di0}} $和$ {n'_{ki}} $分别为第$ i $个目标的多普勒频移及模糊数;$ {\hat n_k} $为模糊数取值范围。$ {\rm{FFT}} $作用于信号慢时间维,则第$ i $个目标的加速度$ {a'_i} $及初始速度$ {v'_{0i}} $分别为:
$$ {a'_i} = \frac{{({{f'}_{di0}} + {{n'}_{ki}}{\rm{prf}})\lambda }}{{2u}},\;{v'_{0i}} = \frac{{{{t'}_i}c - {a_i}{u^2}}}{{2u}} $$ (13) -
根据文献[7-8]的参数所示,设雷达载频为10 $ {\rm{GHz}} $,带宽为10 $ {\rm{MHz}} $,重频为500 $ {\rm{Hz}} $,积累脉冲数为300,在不进行插值处理的情况下,文献[7-8]和式(13) 中的加速度估计精度为$ \frac{1}{6}{\rm{m/}}{{\rm{s}}^2} $,式(13) 中的速度估计精度为50$ {\rm{m/s}} $,文献[7-8]中的速度估计与加速度估计相关,则其速度估计精度更低。由PD雷达的设计原则可知:相参积累时间内目标距离徙动不能超过半个距离单元,多普勒扩散不能超过一个多普勒单元。根据所设定的雷达参数,积累时间内目标距离徙动不超过半个距离单元的最大目标速度$ v $和多普勒扩散不超过一个多普勒单元的最大目标加速度$ a $分别为:12.5 m/s和$ \frac{1}{{24}}{\rm{ m/}}{{\rm{s}}^2} $,可见文献[7-8]和式(13) 的参数估计精度不能满足相参积累的条件。针对此问题,同时考虑到实际工程需求,本文提出联合频域距离徙动补偿与FRFT方法。利用式(13) 得到的运动参数对式(4) 信号进行补偿:
$$ \begin{array}{*{20}{c}} {B(f, {t_m}) = \sum\limits_i {{A_i}{\mathop{\rm rect}\nolimits} \left( {\frac{f}{{\gamma {T_p}}}} \right)} \times }\\ {\exp \left[{-{\rm{j}}2{\rm{\pi }}f\left( {\frac{{2({v_{0i}}-{{v'}_{0i}}){t_m}}}{c} + \frac{{({a_i}-{{a'}_i})t_m^2}}{c}} \right)} \right] \times } \end{array} $$ $$ \begin{array}{*{20}{c}} {\exp \left[{-{\rm{j}}2{\rm{\pi }}(f + {f_c})\frac{{2{R_{0i}}}}{c}} \right]\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{v_{0i}}{t_m}}}{c}} \right] \times }\\ {\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{({a_i}-{{a'}_i})t_m^2}}{c}} \right] + n(f, {t_m})} \end{array} $$ (14) 根据式(14) 信号形式定义频域距离徙动补偿处理公式:
$$ M(f, \hat v, {t_m}) = \exp \left[{-{\rm{j}}2{\rm{\pi }}f\frac{{2\hat v}}{c}{t_m}} \right] $$ (15) 式中,$ \hat v $是为补偿距离徙动而假定的值。利用式(15) 对式(14) 信号进行距离徙动补偿,并转换到时域:
$$ \begin{array}{*{20}{c}} {{b_{rm}}(\hat t, \hat a, \hat v, {t_m}) = \sum\limits_i {{A_i}B\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{R_{0i}}}}{c}} \right] \times } }\\ {{\mathop{\rm sinc}\nolimits} \left[{B\left( {\hat t-\frac{{2{R_{\mathit{0i}}} + 2({v_{\mathit{0i}}}-{{v'}_{0i}}-\hat v){t_m} + ({a_i} - {{a'}_i})t_m^2}}{c}} \right)} \right] \times }\\ {\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{2{v_{0i}}{t_m}}}{c}} \right]\exp \left[{-{\rm{j}}2{\rm{\pi }}{f_c}\frac{{({a_i}-{{a'}_i})t_m^2}}{c}} \right] + }\\ {n(t, {t_m})} \end{array} $$ (16) 式中,$ {\mathop{\rm sinc}\nolimits} (a) = \sin ({\rm{\pi }}a)/{\rm{\pi }}a $。$ {a_i} - {a'_i} $造成的距离徙动可忽略,可见,只要$ {v_{\mathit{0i}}} - {v'_{0i}} - \hat v = 0 $目标距离徙动便可被补偿。
利用FRFT能够对线性调频信号进行能量聚集以及FRFT是线性的优势P[9-10],对距离徙动补偿后$ {b_{rm}}(\hat t, \hat a, \hat v, {t_m}) $的每个距离单元进行阶次变换的分数阶傅里叶变换并取模。本方法采用的FRFT的变换核[8]:
$$ \begin{array}{*{20}{c}} {{K_\alpha }({\tau _m}, w) = }\\ {\left\{ {\begin{array}{*{20}{c}} {\sqrt {\frac{{1 - {\rm{j}}\cot \alpha }}{{2{\rm{\pi }}}}} \exp \left( {{\rm{j}}\frac{{{\tau _m}^2 + {w^2}}}{2}\cot \alpha - {\rm{j}}\frac{{w{\tau _m}}}{{\sin \alpha }}} \right), {\rm{ }}\alpha \ne n\pi }\\ {\begin{array}{*{20}{c}} {\delta ({\tau _m} - w), {\rm{ }}\alpha = 2n\pi }\\ {\delta ({\tau _m} + w), {\rm{ }}\alpha = (2n + 1)\pi } \end{array}} \end{array}} \right.} \end{array} $$ (17) 式中,$ \alpha = \frac{{p\pi }}{2} $;$ p $为分数阶傅里叶变换的阶数。
根据$ {a'_i} $, $ v_{0i}^{'} $, $ a $及$ v $可以设定参数搜索范围及搜索范围中的参数取值,通过代价函数式(18) 可以计算出目标的运动参数:
$$ ({a''_i}, {v''_{0i}}) = \mathop {\arg \max }\limits_{(\hat a', \hat v')} \left| {{\rm{FRFT}}[{b_{rm}}(\hat t, \hat a, \hat v, {t_m}), p]} \right| $$ (18) 式中,$ (\hat a, \hat v) $为$ (\hat a', \hat v') $搜索范围中的参数取值;$ ({a''_i}, {v''_{0i}}) $是估计得到的目标运动参数。
-
为验证所提出方法的有效性,假设雷达仿真参数:Ka波段载频30 GHz雷达,带宽为10 $ {\rm{MHz}} $,脉宽为500 $ {\rm{\mu s}} $,采样频率为10 $ {\rm{MHz}} $,重复频率为200 $ {\rm{Hz}} $,积累脉冲数为200。
目标仿真参数的设定如表 1所示。
表 1 目标仿真参数
参数 目标A 目标B 目标C 目标距离/km 180 181.2 180.6 径向加速度/m·s2 110 121 101 径向速度/m·s-1 567 602 753 在雷达回波信号信噪比-12 dB条件下,噪声为加性高斯白噪声,距离维脉冲压缩的结果如图 1所示。可以看出目标存在距离徙动,并且由于目标的加速度较高,目标距离维脉冲压缩之后产生距离弯曲。对图 1脉冲压缩后信号的方位单元作傅里叶变换进行能量积累,结果如图 2所示。由图 2结果可以得出:由于距离徙动及多普勒扩散的存在,常规检测方法(目标回波做脉压后数据做纵向FFT进行信号能量积累)不适用于高速目标的检测及运动参数估计。
利用文献[2-8]以及本文算法分别对图 1距离维脉冲压缩后对应的目标进行检测。文献[2]的改进Hough变换算法的检测结果如图 3a所示,文献[3-4]的Radon-Fourier变换算法的检测结果如图 3b所示,文献[5]的Keystone变换结合解线频调算法的检测结果如图 3c所示,文献[6]的模糊数搜索算法的检测结果如图 3d所示,文献[7-8]的基于循环平稳理论的迭代算法的检测结果如图 3e所示,本文算法的检测结果如图 3f所示,其中,图 3c~图 3f的检测结果由3次并行处理结果直接相加而成。
文献[2]改进Hough变换算法为基于非相参积累方式的线积分算法,由检测结果图 3a可以看出在低信噪比和目标包络呈现弯曲的情况下,检测结果中出现虚假峰值。文献[3-4]中Radon-Fourier变换算法为基于相参积累方式的线积分算法,然而由于包络弯曲及多普勒扩散的原因,如图 3b所示,其信号能量不能够得到有效的积累。实际战场环境中,多普勒模糊数是未知的,由检测结果图 3c可以看出文献[5]中Keystone变换结合解线频调算法的信号积累效果不理想。文献[6]的算法未考虑多普勒扩散的影响,导致能量发生扩散。文献[7-8]中基于循环平稳理论的迭代算法参数估计精度低,而且在实际中对加速度的取值范围没有任何先验信息,为保证不出现相位模糊,第一步迭代中延时值选择最小的非零延时值,即$ u = {\rm{0}}{\rm{.002 s}} $,但是由于信噪比低,首次迭代的结果易出现模糊的结果,导致后续迭代不能收敛,参数估计错误,由检测结果图 3e可以看出信号积累能量发生扩散。比较图 3f与图 3a~图 3e可以看出,本文算法得到了比较理想的检测效果,原因在于低信噪比下目标距离徙动得到了校正,多普勒扩散得到了补偿。
利用本文提出的算法,在信噪比为-15 dB的情况下进行100次Monte Carlo实验估计目标的运动参数。目标A估计得到的平均初始速度和平均加速度分别为:562.85 $ {\rm{m/s}} $,81 $ {\rm{m/}}{{\rm{s}}^2} $;目标B估计得到的平均初始速度和平均加速度分别为:990.85 $ {\rm{m/s}} $,61 $ {\rm{m/}}{{\rm{s}}^2} $;目标C估计得到的平均初始速度和平均加速度分别为:759.85 $ {\rm{m/s}} $,101 $ {\rm{m/}}{{\rm{s}}^2} $。由估计得到的目标运动参数的精确性可以看出,本文算法在低信噪比下估计出满足相参积累的高速多目标的运动参数。
将本文提出的算法利用某型雷达的实测数据进行验证,雷达积累脉冲数为5 120。对目标进行常规检测,其检测结果如图 4a所示。可得到结论:通过常规检测方法,由于距离徙动及多普勒扩散的存在,目标能量发生明显扩散。通过本文检测算法,其检测结果如图 4b所示。由于本文所提出的算法中,目标的距离徙动及多普勒扩散得到一定的补偿,可以在低信噪比条件下估计出满足相参积累条件的目标的运动参数,从而保证目标能量得到有效的积累。
A Novel Parametric Detection Method for High-Speed Multi-Target
-
摘要: 该文提出一种新的高速、多目标窄带雷达目标检测和参数估计算法。首先采用基于循环平稳的联合Keystone变换与模糊数搜索方法完成目标运动参数粗估计,然后在粗估计基础上采用联合频域距离徙动补偿处理与分数阶傅里叶变换方法完成目标检测及参数估计。该算法适用于多目标及存在距离徙动、多普勒扩散和多普勒模糊的情况,其保持了循环平稳计算复杂度低的优点,且克服了已有循环平稳算法在工程应用中估计精度低和运动参数估计范围受限的缺陷。计算机仿真和实测数据验证了算法的有效性。
-
关键词:
- 循环平稳 /
- 多普勒模糊 /
- 分数阶傅立叶变换 /
- Keystone变换
Abstract: A novel algorithm for high-speed multi-target detection and parameters estimation with narrowband radar is proposed in this paper. Firstly, base on cyclostationarity, moving parameters are roughly estimated with the joint Keystone transform and ambiguity searching. And then, according to the coarse estimation, the joint compensating range migration in frequency domain and fractional Fourier transform are utilized to conduct the targets detection and parameters estimation. In the case of multi-target, range migration, Doppler spread and Doppler ambiguity, this algorithm is suitable, and retains the merit of low computational complexity of cyclostationarity. Compared with existing cyclostationarity based algorithms, the weaknesses of low estimation precision and limited range of moving parameters are conquered in engineering applications. The validity of the proposed algorithm is demonstrated by computer simulation and raw radar data results.-
Key words:
- cyclostationarity /
- Doppler ambiguity /
- fractional Fourier transform /
- Keystone transform
-
表 1 目标仿真参数
参数 目标A 目标B 目标C 目标距离/km 180 181.2 180.6 径向加速度/m·s2 110 121 101 径向速度/m·s-1 567 602 753 -
[1] 许述文. 窄带、宽带雷达机动目标检测技术研究[D]. 西安: 西安电子科技大学, 2011: 1-3. http://cdmd.cnki.com.cn/Article/CDMD-10701-1011200366.htm XU Shu-wen. Detection methods of manoeuvring targets in narrowband and wideband Radars[D]. Xi'an:Xidian University, 2011:1-3. http://cdmd.cnki.com.cn/Article/CDMD-10701-1011200366.htm [2] 曾建奎, 何子述, 刘红明.一种基于改进Hough变换的雷达检测方法[J].电波科学学报, 2008, 23(5):838-840. http://www.cnki.com.cn/Article/CJFDTOTAL-DBKX200805008.htm ZENG Jian-kui, HE Zi-shu, LIU Hong-ming. Improved detection algorithm for radar based on hough transform[J]. Chines Journal of Radio Science, 2008, 23(5):838-840. http://www.cnki.com.cn/Article/CJFDTOTAL-DBKX200805008.htm [3] JIA Xu, JI Yu, YING Ning-peng, et al. Radon-Fourier transform for radar target detection, (i):Generalized doppler filter bank[J]. IEEE Transaction on Aerospace and Electronic Systems, 2011, 47(4):2473-2489. doi: 10.1109/TAES.2011.6034645 [4] JIA Xu, JI Yu, YING Ning-peng, et al. Radon-Fourier transform for radar target detection, (ii):Blind speed sidelobe suppression[J]. IEEE Transaction on Aerospace and Electronic Systems, 2011, 47(2):1186-1202. doi: 10.1109/TAES.2011.5751251 [5] 吴兆平, 何学辉, 苏涛.带有距离走动和多普勒扩散的高速运动目标检测[J].哈尔滨工程大学学报, 2010, 31(4):476-480. http://mall.cnki.net/magazine/Article/HEBG201004014.htm WU Zhao-ping, HE Xue-hui, SU Tao. Detecting multiple high speed targets using range migration and doppler speed[J]. Journal of Harbin Engineering University, 2010, 31(4):476-480. http://mall.cnki.net/magazine/Article/HEBG201004014.htm [6] 曹宇飞, 屈晓光, 黄培康.基于keytone变换的宽带雷达精准测速[J].系统工程与电子技术, 2009, 31(1):1-4. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=xtyd200901002&dbname=CJFD&dbcode=CJFQ CAO Yu-fei, QU Xiao-guang, HUANG Pei-kang. Accurate velocity measurement method for wideband radar based on keytone transform[J]. Systems Engineering and Electronics, 2009, 31(1):1-4. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=xtyd200901002&dbname=CJFD&dbcode=CJFQ [7] JIN Yan, JI Hong-bing. A cyclostationarity based iterative algorithm for multi-component chirp signal parameter estimation[C]//Proceedings of 2007 International Symposium on Intelligent Signal Processing and Communication Systems.[S.l.]:IEEE Press, 2007. [8] JIN Ya, JI Hong-bing. Improved parameter estimation for chirp signals based on cyclostationarity[J]. Chinese Journal of Electronics, 2008, 17(2):309-313. http://www.cqvip.com/main/none.aspx?lngid=27546425 [9] 陶然, 邓兵, 王越.分数阶Fourier变换在信号处理领域的研究进展[J].中国科学E辑信息科学, 2006, 36(2):113-136. http://info.scichina.com:8084/sciF/CN/abstract/abstract405783.shtml TAO Ran, DENG Bing, WANG Yue. The research of the fractional Fourier transform in signal processing field[J]. Science in China Ser E Information Sciences, 2006, 36(2):113-136. http://info.scichina.com:8084/sciF/CN/abstract/abstract405783.shtml [10] ZHANG S, SUN S, ZHANG W, et al. High-Resolution bistatic ISAR image formation for high-speed and complex-motion targets[J]. IEEE Journal of Selected Topics in Applied Earth Observations & Remote Sensing, 2015, 8(7):3520-3531. http://ieeexplore.ieee.org/lpdocs/epic03/wrapper.htm?arnumber=7084575