专利名称：一种空时自适应处理中的自适应权值迭代计算方法
技术领域：
本发明属于机载相控阵雷达杂波抑制技术领域，涉及一种空时自适应处理中的自适应权值迭代计算方法。
背景技术：
机载相控阵雷达能实现对地面运动目标的有效检测，不过处于下视工作状态的机载相控阵雷达将面临比地基雷达更为严重的地/海杂波问题。地/海杂波不仅分布范围广、强度大，而且呈现出很强的空时耦合特性。空时自适应处理(STAP)技术能充分利用空域和时域信息，对目标信号进行相干积累的同时，将自适应空域处理和自适应多普勒处理两者优势结合起来，在空时域联合自适应滤除杂波，能获得更好的主瓣杂波抑制性能，改善对慢速目标的检测；同时可以对受旁瓣杂波干扰的小目标进行有效的检测。
常规的全维STAP算法需要大量满足独立同分布(I. I. D)条件的训练样本对协方差矩阵进行估计，在非均匀杂波环境中此条件尤其得不到满足；而且当系统自由度很高时，全维协方差矩阵直接求逆运算(SMI)在现有计算水平下几乎无法实现。虽然降维和非均匀STAP算法能够降低STAP自适应权值计算的运算量和改善非均匀杂波环境中STAP算法杂波抑制性能，但是常规降维STAP算法在求解自适应权值时依然面临着对协方差矩阵直接求逆的运算，这将耗费系统很大的运算量和设备量，使得STAP技术难以满足系统实时性的要求。另外，虽然基于协方差矩阵逆更新的SMI算法无需估计采样协方差矩阵，其迭代计算的次数等于训练样本的数目，但是因为零矩阵不存在逆矩阵，该算法将难以设置初始逆矩阵来等价求解SMI算法的自适应权矢量，只能获得近似的解。
发明内容
本发明针对常规STAP算法自适应权值计算需对空时协方差矩阵直接求逆，耗费系统很大运算量和设备量，使得STAP技术难以满足实时性要求的问题，根据协方差矩阵为Hermitian矩阵的特性，利用脉冲数据阶数分块递推，提出一种空时自适应处理中的自适应权值迭代计算方法。步骤一、建立空时接收数据模型；假设雷达天线阵元数目为N，发射脉冲数目为M，阵元间距为d，载机飞行速度为V，高度为h，脉冲重复频率PRF为f；，Tr = l/fr为脉冲重复时间；若将斜距为R。处的杂波距离环在方位角度Q上分成N。个间隔为A 0 = 2 /N。的杂波散射单兀，0和—是杂波散射单元的方位角和俯仰角，/s =c/cos0cos识/A和ft = ^ fs分别为归一化空间频率和多普勒频率，^ = 2vTr/d表示杂波谱的斜率；那么第i个杂波散射单元的NX I维空间导向矢量c(fs；i)和MX I维时间导向矢量c (ft，i)即表示为
ciL.,) = [lexp(./2^r,/；v.),exp(；2^-(N ^ l)fs i)]'4T(5)
C(Zlj) = [l，exp(./2;r.4)”..，exp(,/2;r(M -1),/；,,)]
斜距R。处的杂波回波为N。个在空间上相互独立的杂波散射源响应之和
权利要求
1.一种空时自适应处理中的自适应权值迭代计算方法，其特征在于，包括以下步骤 步骤一、建立空时接收数据模型； 假设雷达天线阵元数目为N，发射脉冲数目为M，阵元间距为d，载机飞行速度为V，高度为h，脉冲重复频率PRF为f；，Tr = l/fr为脉冲重复时间；若将斜距为R。处的杂波距离环在方位角度9上分成N。个间隔为A 0 = 2 /N。的杂波散射单兀，0和是杂波散射单兀的方位角和俯仰角，X =rfcos0cosp//^P ft = ^ fs分别为归一化空间频率和多普勒频率，3=2vTyd表示杂波谱的斜率；那么第i个杂波散射单元的NX I维空间导向矢量c (U和MX I维时间导向矢量c(ft，i)即表示为
2.如权利要求I所述的一种空时自适应处理中的自适应权值迭代计算方法，其特征在于，步骤三中所述的按照脉冲阶数进行递归计算空时协方差矩阵的逆采用以下方法 在接收到第I个脉冲数据时，利用距离训练样本计算协方差矩阵R(I)，因为该矩阵是Hermitian矩阵,利用Hermitian矩阵的顺序主子矩阵迭代完成第I个脉冲协方差矩阵的求逆；接着在接收第2个直至第M个脉冲数据时，便利用以上迭代方法计算前2个直至前M个脉冲接收数据协方差矩阵的逆；递推过程中，中间变量矩阵P(2)直至矩阵P(M)的逆，同样根据Hermitian矩阵的性质,利用Hermitian矩阵的顺序主子矩阵迭代完成上述NXN维 的变量矩阵逆的计算。
全文摘要
本发明针对常规STAP算法自适应权值计算需对空时协方差矩阵直接求逆，耗费系统很大运算量和设备量，使得STAP技术难以满足实时性要求的问题，提出一种空时自适应处理中的自适应权值迭代计算方法。该方法首先根据Hermitian矩阵性质，递推得到第一个脉冲协方差矩阵的逆矩阵，然后按照脉冲阶数逐级嵌套递推得到最终的空时协方差矩阵的逆，这样极大地降低了计算STAP自适应权值的运算量。本发明能够获得和协方差矩阵直接求逆STAP算法同样的杂波抑制性能，但是由于避免了协方差矩阵直接求逆的运算，求解自适应权值的计算量只有协方差矩阵直接求逆的50%左右，因此更利于工程实现。
文档编号G01S7/36GK102866388SQ201210222659
公开日2013年1月9日 申请日期2012年6月28日 优先权日2012年6月28日
发明者杨小鹏, 刘永旭, 龙腾, 曾涛 申请人:北京理工大学