专利名称：基于重构信号和时间采样的空中机动目标检测方法
技术领域：
本发明属于雷达信号处理技术领域，特别是涉及一种基于重构信号和时间采样的空中机动目标检测方法。
背景技术：
由于载机平台的运动，机载相控阵雷达的地杂波谱在多普勒上表现出主杂波展宽和旁瓣杂波扩散，并且是空时二维耦合的，因此杂波抑制是机载预警雷达下视工作时的关键问题。空时二维自适应处理(Space-Time Adaptive Processing, STAP)作为一种有效的机载雷达地杂波抑制手段，已受到了越来越多的关注。它是在空间和时间上的二维联合滤波，同时具有空间和时间自由度，可以在二维谱平面上形成凹口，在对地杂波进行有效抑制的同时能够保证目标信号获得足够的增益。但是STAP是基于平稳信号模型的，即假定在相干处理时间内动目标回波多普勒频率恒定(即目标为匀速运动)；而当来袭目标具有很强的机动性时，其在一个相干处理时间(Coherent Processing Interval,CPI)内目标回波多普勒频率随时间发生变化，即发生多普勒走动，使得传统的STAP方法相参积累性能大大下降，从而导致目标检测能力下降。当机动目标做勻加速运动时目标回波信号为线性调频(Linear FrequencyModulation, LFM)信号。对LFM信号的参数检测与估计可以采用最大似然估计(MaximumLikelihood,ML)和时频分析等经典方法。其中ML方法是一种最有效的估计方法，其估计精度很高，对参数的估计方差接近于Cramer-Rao下界，但其缺点是运算量很大，并且该方法的参数估计模型是在高斯白噪声环境下进行的，而在实际情况中，噪声往往是色噪声，因此ML方法不利于在实际情况中进行实时处理和工程实现。近年来，基于各种时频分析工具的LFM信号的检测与参数估计方法不断出现，包括短时Fourier变换(STFT) ,Wigner-Ville分布(WVD)和分数阶 Fourier 变换(FRactional Fourier Transform, FRFT)等。其中 STFT是利用滑动窗函数对Fourier变换进行的简单扩展，即用窗函数对信号在时间上截取一段后再作Fourier变换，但STFT的时频分辨率会受到观察窗的影响，因而分辨率不高。WVD是一种时频分辨率最高的双线性时频分布。但由于WVD是非线性的，因此当信号包含多分量成分时，WVD会存在交叉项干扰，因此严重地干扰着人们对信号时频特性的解释。FRFT是一种一维的线性变换，在检测多个机动目标时它不会像WVD那样产生交叉项，而且其数值计算可借助快速傅立叶变换(FFT)快速实现，计算简单，易于实现，因此受到了广泛重视。但是，利用FRFT来估计机动目标参数时需要较多的采样点数，因此该类方法在地基雷达和合成孔径雷达(Synthetic Aperture Radar, SAR)中应用比较广泛(地基雷达和SAR发射脉冲数较多)。而当机载预警雷达的脉冲重复频率一定时，较多的采样点数意味着CPI加长，这会引起杂波和目标的距离走动，给后续处理带来更大困难，因此直接利用时频分析方法来检测空中机动目标会出现估计精度较差的问题
发明内容
为了解决上述问题，本发明的目的在于提供一种能够提高参数估计精度的机载雷达机动目标检测新方法。为了达到上述目的，本发明提供的基于重构信号和时间采样的空中机动目标检测方法包括按顺序进行的下列步骤I)对机载预警雷达接收到的总回波数据进行杂波抑制的SI阶段；2)估计目标回波幅度的S2阶段；3)确定目标参数二维搜索范围的S3阶段；4)对多阵元数据重构时间采样，构造代价函数，得到目标参数估计结果的S4阶段。在步骤I)中，所述的对机载预警雷达接收到的总回波数据进行杂波抑制的方法是利用子空间投影技术对机载预警雷达接收到的总回波数据进行杂波抑制，即将待检测单元数据投影到杂波子空间的正交子空间中，得到投影后的无杂波数据。在步骤2)中，所述的估计目标回波幅度的方法是利用特征值分解方法进行。在步骤3)中，所述的确定目标参数二维搜索范围的方法是利用快速解线调法粗略地估计出目标信号的初始速度和加速度，确定一个目标参数的二维搜索范围。在步骤4)中，所述的对多阵元数据重构时间采样，构造代价函数，得到目标参数估计结果的方法是利用重构时间采样的方法对多个阵元杂波抑制后的数据进行补偿相位，将补偿后的各阵元数据首尾拼接，即用空间采样来重构时间采样，等效于增加单个阵元脉冲点数的效果；同时，构造单个阵元直接增加脉冲点数的目标信号，利用非线性最小二乘(Nonlinear least squares,NLS)法将重构时间采样后的数据与单个阵元直接增加脉冲点数的目标信号相匹配，构造代价函数；以目标信号和拼接数据之间具有最小“平方和”时的参数作为目标参数的估计结果。本发明提供的基于重构信号和时间采样的空中机动目标检测方法是对多个阵元杂波抑制后的数据进行重构时间采样，等效于增加了单个阵元的脉冲点数，同时构造单个阵元直接增加脉冲点数时的目标信号，并利用非线性最小二乘法将上述两组数据进行匹配，通过计算最小的“平方和”距离来进行目标参数估计。本发明能够在脉冲点数有限的情况下得到精确的估计结果。
图I为本发明提供的重构信号和时间采样的空中机动目标检测方法流程图。
图2 (a)为N个阵元数据拼接前时频图(N == 4)。
图2(b)为N个阵元数据拼接后时频图(N == 4)。
图3 (a)为N个阵元数据拼接前相位图(N ==2)。
图3(b)为N个阵元数据拼接后相位图(N ==2)。
图4为杂波抑制前总回波的功率谱。
图5为常规方法杂波抑制后的功率谱。
图6为本发明方法处理后的功率谱。
图7为对加速度补偿前后功率谱。
图8(a)为初始速度均方根误差随信噪比变化曲线。
图8(b)为加速度均方根误差随信噪比变化曲线。
具体实施例方式下面结合附图和具体实例对本发明提供的基于重构信号和时间采样的空中机动目标检测方法进行详细说明。图I为本发明提供的基于重构信号和时间采样的空中机动目标检测方法流程图。如图I所示，本发明提供的重构信号和时间采样的空中机动目标检测方法包括按顺序进行的下列步骤I)对机载预警雷达接收到的总回波数据进行杂波抑制的SI阶段在此阶段中，利用子空间投影技术对机载预警雷达接收到的总回波数据进行杂波抑制，具体方法如下将待检测单元数据投影到杂波子空间的正交子空间中，得到投影后的无杂波数据TLpraj = V X(I)其中,Pf为杂波子空间的正交补空间的投影矩阵，X为机载预警雷达接收到的总回波数据，由于通常of <為=C1，<为估计得到的噪声方差，λ i为协方差矩阵的大特征值；而确定杂波子空间维数较复杂，所以我们用杂波加噪声的协方差的逆矩阵R-1代替1^抑制杂波，实际中杂波协方差矩阵R由参考距离单元数据估计出，因此，杂波抑制后的数据可写为Xproi=R-1X(2)2)估计目标回波幅度的S2阶段本阶段利用特征值分解方法估计出目标回波幅度 为估计得到的杂波加噪声的协方差矩阵，对其进行特征值分解R = E 為 W + Σ Ww
/=1十 I其中，X1S λ2彡…彡入5是协方差矩阵的P个大特征值，ΝΚ-ρ个特征值λρ+1> λρ+2彡…彡近似等于噪声方差为第i个特征向量。因此可由小特征值来估计噪声方差为O-I =---T Af(4)
"ΝΚ-Ρ + \& '通过矩阵求逆消杂波后可以得到信号加噪声的方差天+# 2,因此，用其减去估计得到的噪声方差式就可以得到信号方差€ ,此时目标回波幅度为bt = ^Ja2pmj-σ2η (5)3)确定目标参数二维搜索范围的S3阶段在此阶段中，为了降低本发明方法的计算量，首先对目标的初始速度和加速度进行粗估，来确定参数搜索范围，为下一步精估做准备。本发明方法采用快速解线调算法来确定参数搜索范围。4)对多阵元数据重构时间采样，构造代价函数，得到目标参数估计结果的S4阶段由背景技术分析可知，利用时频分析方法来估计机动目标参数时需要较多的脉冲点数，否则估计精度难以满足要求，而当机载预警雷达的脉冲重复频率一定时，较多的采样点数意味着CPI加长，这会引起杂波和目标的距离走动，给后续处理带来更大困难。当目标做匀加速运动时，对于机载预警雷达的每个阵元来说其回波信号均为一个LFM信号(杂波抑制后)，并且由目标的空时二维数据模型可知每个阵元的LFM信号只差空间相位，其时频图如图2(a)所示(图中N = 4)。因此可以利用干涉SAR中相位展开的思想，对空间中每个阵元的数据进行相位补偿后首尾拼接，使其等效于增加单个阵元时间采样点数的效果，拼接后时频图如图2(b)所示。图3为两个阵元数据拼接前后相位图。不失一般性，下面以两个阵元为例讨论将多阵元首尾拼接时每个阵元所需补偿的相位。由目标数据模型可知，当不考虑空间相位时，两个阵元接收到的目标数据为
rm/1 n; ''TT+^17 }2π·<Κ-\)^-+]π·(ΚΑ·ΓW_=[U 从 '"VA.人#]“(O)=[1 ej2;r'l'fd+j7r'uy "· eJ2/r<K-l>fci+J^<K-iy-ad ]τ其中，Xsl表示第I个阵元接收到的目标信号，Xs2表示第2个阵元接收到的目标信
O IO
号，K为脉冲点数，/d =IilV为初始频率，&为调频率。V和a分别为目标初始速度
Λ· /,,Ar Jr
和加速度。当一个阵元的时间采样点数由K增加到2Κ时，该阵元接收到的目标信号为
I = [I eJZn~'h~fl+^'U^ …eJ1πΛ K~l ■/(^-1 Γ ·( (+^·Λ~ Ud
^2π·[ K +1)-/, + JTTi K +1 Γ ...)■// + r ■ / j
=[1 ...
Ο [I I,.· I eJl7cKf'+-mKla'产 ρΚ' +ρπΚλ-α
、............. v·'..............
i 2A.V11
r Vft +7^--1^.JlwiK-Yyfd ^jWiK-Yf -a, .T
， ，CJ( 2AW1)
0jlKKf·,十jkK'a +j2^K-iK-l ha,
Ii Λ\1 I (7)
-￥jjzK~Qi + j2ttK-(AT—I )·α^
I
CA xl)
= Λι，χ 2。，2 Wi)其中，Θ为Hadamard积，T表示转置操作。由此可见，对第2个阵元的数据进行相位补偿后再跟第I个阵元的数据拼接就可以等效为一个阵元直接增加时间采样点数的效
果O由上式错误！未找到引用源。可以得出对第2个阵元应该补偿的相位为
= [./2冗人:/;/ +./-T人 ,2;τ人:/“. + ./>人 + ]2πΚ ■ I - ,(8).■ ■,./2πΚ/[, + ]πΚ% + βπΚ ·(/C - I).a,同理，第η个阵元应该补偿的相位为
权利要求
1.一种基于重构信号和时间采样的空中机动目标检测方法，其特征在于，所述的基于重构信号和时间采样的空中机动目标检测方法包括按顺序进行的下列步骤 1)对机载预警雷达接收到的总回波数据进行杂波抑制的Si阶段； 2)估计目标回波幅度的S2阶段； 3)确定目标参数二维搜索范围的S3阶段； 4)对多阵元数据重构时间采样，构造代价函数，得到目标参数估计结果的S4阶段。
2.根据权利要求I所述的基于重构信号和时间采样的空中机动目标检测方法，其特征在于在步骤I)中，所述的对机载预警雷达接收到的总回波数据进行杂波抑制的方法是利用子空间投影技术对机载预警雷达接收到的总回波数据进行杂波抑制，即将待检测单元数据投影到杂波子空间的正交子空间中，得到投影后的无杂波数据。
3.根据权利要求I所述的基于重构信号和时间采样的空中机动目标检测方法，其特征在于在步骤2)中，所述的估计目标回波幅度的方法是利用特征值分解方法进行。
4.根据权利要求I所述的基于重构信号和时间采样的空中机动目标检测方法，其特征在于在步骤3)中，所述的确定目标参数二维搜索范围的方法是利用快速解线调法粗略地估计出目标信号的初始速度和加速度，确定一个目标参数的二维搜索范围。
5.根据权利要求I所述的基于重构信号和时间采样的空中机动目标检测方法，其特征在于在步骤4)中，所述的对多阵元数据重构时间采样，构造代价函数，得到目标参数估计结果的方法是利用重构时间采样的方法对多个阵元杂波抑制后的数据进行补偿相位，将补偿后的各阵元数据首尾拼接，即用空间采样来重构时间采样，等效于增加单个阵元脉冲点数的效果；同时，构造单个阵元直接增加脉冲点数的目标信号，利用非线性最小二乘法将重构时间采样后的数据与单个阵元直接增加脉冲点数的目标信号相匹配，构造代价函数；以目标信号和拼接数据之间具有最小“平方和”时的参数作为目标参数的估计结果。
全文摘要
本发明公开了一种基于重构信号和时间采样的空中机动目标检测方法。包括1)对机载预警雷达接收到的总回波数据进行杂波抑制；2)估计目标回波幅度；3)确定目标参数二维搜索范围；4)对多阵元数据重构时间采样，构造代价函数，得到目标参数估计结果。本发明方法是对多个阵元杂波抑制后的数据进行重构时间采样，等效于增加了单个阵元的脉冲点数，同时构造单个阵元直接增加脉冲点数时的目标信号，并利用非线性最小二乘法将上述两组数据进行匹配，通过计算最小的“平方和”距离来进行目标参数估计。本发明能够在脉冲点数有限的情况下得到精确的估计结果。
文档编号G01S7/41GK102928826SQ20121040306
公开日2013年2月13日 申请日期2012年10月22日 优先权日2012年10月22日
发明者李海, 吴仁彪, 王小寒, 郑景忠 申请人:中国民航大学