本发明属于雷达技术领域,特别涉及机载雷达杂波抑制方面,具体是一种基于知识辅助稀疏迭代协方差估计的杂波抑制方法,适用于非均匀环境下机载预警雷达抑制强地杂波以及检测地面慢速运动目标。
背景技术:
对于机载预警雷达而言,波束指向一般是水平或稍下方向,此时雷达不可避免的要观测到许多强地杂波。复杂强地杂波严重影响了机载预警雷达对地面运动目标的检测,因此在复杂强地杂波背景下检测运动目标必须首先抑制地杂波。空时自适应处理技术(stap)能够从空域和时域两维上同时区分运动目标和杂波,能够有效改善机载预警雷达的杂波抑制和运动目标检测性能,尤其是对慢速目标的检测更为有效,因而得到了广泛的关注。
由于待检测单元杂波信号的统计特性通常事先是不知道的,传统stap需要选取待检测单元附近的回波数据作为训练样本对其进行估计。为获得较优的stap性能(相对最优处理的性能损失不超过3db),所需要的训练样本数不应少于系统自由度的两倍。此处假设所选用的训练样本与待检测单元数据满足独立同分布条件,即这些训练样本是均匀的且和待检测单元数据的统计特性相同。然而,由于实际场景中存在复杂的地表类型、目标污染等情况,机载雷达通常都工作在非均匀杂波环境中。另外,当雷达天线旋转时,这也会导致机载雷达回波数据非均匀,使得上述独立同分布假设在实际环境中很难满足,进而影响待检测单元杂波协方差矩阵估计精度,导致stap的杂波抑制性能显著下降。
在非均匀环境下,林肯实验室ward报告所提出的传统stap方法不能有效的抑制地杂波。常见的非均匀杂波抑制方法有降秩方法和结构化方法,它们均能提高非均匀环境下杂波抑制的性能。降秩方法(rr)充分利用回波中杂波的分布特性,选取完备的杂波空间来形成自适应权值,对消杂波分量,改善杂波抑制性能,这主要包括主分量(pc)法,互谱法(csm)和多级维纳滤波(mwf),但是无论是pc,csm还是mwf,都需要预先估计出杂波子空间维数,而求取降秩变换阵和确定杂波子空间维数比较复杂,也因此增加了处理的运算量。结构化方法需要通过特征分解进行协方差矩阵的最大似然估计,利用已知的杂波协方差矩阵的结构特性对估计得到的杂波协方差矩阵进行修正,提高了杂波协方差矩阵估计的精度,但由于要进行最大似然估计同样也存在运算量较大的问题。
技术实现要素:
针对上述现有技术存在的问题,本发明的目的在于提出一种精度更高,运算量更小的基于知识辅助稀疏迭代协方差估计的杂波抑制方法。
本发明是一种基于知识辅助稀疏迭代协方差估计的机载雷达杂波抑制方法,其特征是:包括有以下步骤:
步骤1,获取待检测单元数据,计算空时导向矢量矩阵v:机载雷达的检测范围内存在杂波散射体s,分别获取nmax个距离门雷达回波数据和nmax个待检测单元数据,nmax为距离门的个数;分别计算杂波散射体s的归一化多普勒频率和杂波散射体的归一化空域频率,确定杂波脊,进而计算空时导向矢量矩阵v。
步骤2,确定初始杂波功率矩阵与空域协方差矩阵,构造中间变量:确定第l个待检测单元数据xl在杂波脊上的初始杂波功率矩阵
步骤3,计算杂波功率矩阵
(3a)计算第二中间变量ρ(i-1)并由第一中间变量d和第二中间变量ρ(i-1)得到第i-1次迭代后第l个待检测单元在杂波脊上的杂波功率矩阵
(3b)由第i次迭代后第l个待检测单元在杂波脊上的杂波功率矩阵
步骤4,对杂波功率矩阵
步骤5,确定第l个待检测单元数据xl重构的空时二维协方差矩阵
步骤6,确定nmax个待检测单元数据重构的空时二维协方差矩阵
本发明的有益效果:
第一,本发明方法利用基于知识辅助稀疏迭代协方差估计方法的优势,通过将待检测单元数据分解成空域快拍数据,再利用空域快拍数据重构待检测单元的空域协方差矩阵,最后利用杂波脊先验信息和估计的杂波块功率矩阵得到待检测单元的空时协方差矩阵,因此有效解决了基于采样协方差估计待检测单元杂波协方差矩阵存在的样本非均匀问题,相比传统的杂波谱重构方法谱估计精度更高,从而使得本发明方法能够有效抑制复杂强地杂波,显著改善了对运动目标的检测性能。
第二,本发明方法通过利用知识辅助稀疏迭代协方差估计谱重构技术重构待检测单元的杂波协方差矩阵,这样的重构过程使得杂波协方差矩阵的估算运算量更小,并且这种方法具有全局收敛性,因此该方法的实时性更好,实用价值更高。
附图说明
图1为本发明的基于知识辅助稀疏迭代协方差估计的杂波抑制方法流程图;
图2为机载雷达几何构型示意图;
图3(a)为最优处理方法得到的正侧阵时第100号距离门空时二维谱示意图;
图3(b)为本发明方法得到的正侧阵时第100号距离门空时二维谱示意图;
图4(a)为最优处理方法得到的30度安装角时第100号距离门空时二维谱示意图;
图4(b)为本发明方法得到的30度安装角时第100号距离门空时二维谱示意图;
图5(a)为正侧阵时第100号距离门本发明方法与最优处理方法改善因子对比示意图;
图5(b)为30度安装角时第100号距离门本发明方法与最优处理方法改善因子对比示意图。
具体实施方式
下面结合附图和具体实例对本发明作详细描述。
实施例1
传统stap需要选取待检测单元附近的回波数据作为训练样本对其进行估计,所选用的训练样本与待检测单元数据必须满足独立同分布条件,然而事实上机载雷达通常都工作在非均匀杂波环境中,独立同分布假设在实际环境中很难满足,导致stap的杂波抑制性能显著下降。而针对上述非均匀问题,出现了各种非均匀杂波的抑制方法,对传统的stap进行了改进,常见的有降秩方法和结构化方法,这两种方法均克服了非均匀问题,提高了杂波抑制的精度,但也加大了运算量和复杂度。因此对于传统stap的改进,能有效抑制地杂波,明显改善对微弱慢速运动的目标检测性能的同时,又能保证运算量不大,使工程上易于实现成为了主要的研究目的。
为达到上述技术目的,本发明采用如下技术方案予以实现:
本发明是一种基于知识辅助稀疏迭代协方差估计的杂波抑制方法,包括有以下步骤,参见图1:
步骤1,为得到杂波脊的先验信息,必须先获取待检测单元数据,进而计算空时导向矢量矩阵v:机载雷达的检测范围内存在杂波散射体s,分别获取nmax个距离门雷达回波数据和nmax个待检测单元数据;分别计算杂波散射体的归一化多普勒频率和杂波散射体的归一化空域频率,确定杂波脊,进而计算空时导向矢量矩阵v。具体有几下步骤实现:
1a确定机载雷达几何构型参数,机载雷达的检测范围内存在杂波散射体s,参见图2,图2为机载雷达几何构型示意图,其中,x-y平面表示地平面,机载雷达载机以高度h、速度v沿着x轴方向飞行,s为相对于机载雷达载机平台斜距为rl的杂波散射体,机载雷达发射信号并接收雷达回波数据,雷达回波数据包含nmax个距离门雷达回波数据,nmax为距离门的个数,将每个距离门的雷达回波数据分别记为一个待检测单元数据,进而得到nmax个待检测单元数据。
1b利用设定的雷达几何构型参数分别计算杂波散射体的归一化多普勒频率和杂波散射体的归一化空域频率,确定杂波脊,进而计算空时导向矢量矩阵v。
步骤2,确定初始杂波功率矩阵与空域协方差矩阵,构造中间变量:传统方法求得的空时协方差矩阵
2a利用第l个待检测单元数据xl确定xl在杂波脊上的初始杂波功率矩阵为
2b为求出第i次迭代的杂波功率矩阵
步骤3,计算杂波功率矩阵
3a计算第二中间变量ρ(i-1)并由第一中间变量d和第二中间变量ρ(i-1)得到第i-1次迭代后第l个待检测单元在杂波脊上的杂波功率矩阵
3b由3a求得的杂波功率矩阵
步骤4,对杂波功率矩阵
若i<γ,则令i的值加1,返回步骤3,继续计算杂波功率矩阵
若i≥γ,则迭代停止,并将迭代停止时对应的第i次迭代后第l个待检测单元数据xl在杂波脊上的杂波功率矩阵
步骤5,计算一个待检测单元的重构协方差矩阵和与其对应的权:以第l个待检测单元为例,确定第l个待检测单元数据xl重构的空时二维协方差矩阵
步骤6,确定nmax个待检测单元数据
为解决非均匀杂波抑制的技术问题,本发明采用基于知识辅助稀疏迭代协方差估计的杂波抑制的技术。基于知识辅助稀疏迭代协方差估计的杂波抑制方法实际上就是利用待检测单元的数据重构待检测单元的杂波协方差矩阵,该方法可以有效避免训练样本的非均匀性,并且不需要人为设置参数,大大减小了运算量和运算时间,在工程上比较容易实现。
本发明通过计算多个待检测单元数据在杂波脊上的初始杂波功率矩阵,构造中间变量,利用杂波脊先验信息与估计的空域协方差矩阵,通过稀疏迭代方法求得最终杂波功率矩阵,进而得到重构的空时协方差矩阵,重构的空时协方差矩阵克服了基于采样待检测单元杂波空时协方差矩阵存在的样本非均匀问题,提高待检测单元杂波协方差矩阵的估计精度,工程上易于实现,同时又能取得良好的目标检测性能。
实施例2
基于知识辅助稀疏迭代协方差估计的杂波抑制方法同实施例1,步骤2中所述的构造第一中间变量d和第二中间变量ρ(i),具体为:
为对初始估计的杂波功率矩阵
其中,a(fs,m)表示归一化空域频率为fs,m对应的第m个杂波块的空域导向矢量,上标h表示复共轭运算,m∈{1,2,…,nc},nc表示等距离环上独立杂波块数,tr表示求矩阵的迹运算,
计算第i-1次迭代后第二中间变量ρ(nc×1)的值ρ(i-1),第一中间变量d对于每一个距离门来说为一固定值,第二中间变量ρ(i)与迭代次数i有关,其表达式为:
其中,||·||表示矩阵的2范数,a(fs,m)表示归一化空域频率为fs,m对应的第m个杂波块的空域导向矢量,
本发明用杂波脊的先验信息与空域协方差矩阵
实施例3
基于知识辅助稀疏迭代协方差估计的杂波抑制方法同实施例1-2,步骤3中所述的第i次迭代后第l个待检测单元数据xl在杂波脊上的杂波功率矩阵
经过第i次迭代后第l个待检测单元数据xl在杂波脊上第m个杂波块的功率
其中,a(fs,m)表示归一化空域频率为fs,m对应的第m个杂波块的空域导向矢量,
实施例4
基于知识辅助稀疏迭代协方差估计的杂波抑制方法同实施例1-3,步骤5中所述的确定第l个待检测单元数据xl重构的空时二维协方差矩阵
已得到最终修正的杂波功率矩阵
其中,v表示空时二维导向矢量矩阵,
传统stap方法利用公式
下面给出一个更加详细和完整的例子,对本发明进一步说明:
实施例5
基于知识辅助稀疏迭代协方差估计的杂波抑制方法同实施例1-4,参照图1,图1为本发明的一种基于知识辅助稀疏迭代协方差估计的杂波抑制方法流程图,本发明的实现包括以下步骤:
步骤1,确定机载雷达,机载雷达的检测范围内存在杂波散射体s,所述杂波散射体s为机载雷达波束照射到地面上时能够反射机载雷达不需要的回波的物体。
机载雷达发射信号并接收雷达回波数据,雷达回波数据包含nmax个距离门雷达回波数据,依次记为第1个距离门的雷达回波数据、第2个距离门的雷达回波数据、……、第nmax个距离门的雷达回波数据。
参照图2,图2为机载雷达几何构型图;其中,x-y平面表示地平面,所述地平面未考虑地表的曲率问题,并且垂直于x轴方向的平面在水平面的交线为y轴,垂直于地平面方向定义为z轴方向;机载雷达载机以高度h、速度v沿着x轴方向飞行;s为相对于机载雷达载机平台斜距为rl的杂波散射体;将机载雷达载机速度v的方向与机载雷达天线轴向的夹角,记为安装角α;将杂波散射体s和机载雷达天线中心之间的连线与水平方向的夹角,记为俯仰角
其中,将机载雷达平台到杂波散射体的直线距离记为斜距,第l个距离门的雷达回波数据对应斜距为rl,rl=h+(l-1)δr,δr表示机载雷达的距离分辨率,
l为待检测单元的序号,当l=1时,rl=h,l=nmax,对应第nmax个距离门的雷达回波数据;当l=nmax时,rl=h+(nmax-1)δr,对应第nmax个距离门的雷达回波数据,且取值与机载雷达的最大不模糊距离ru取值相等,
将每个距离门的雷达回波数据分别记为一个待检测单元数据,进而得到nmax个待检测单元数据。
杂波散射体s的空域频率表示为
杂波散射体s的归一化空域频率表示为fs:
其中,d为机载雷达的阵元间距,λ为机载雷达发射电磁波的波长,
杂波散射体s的多普勒频率表示为
其中,v表示机载雷达载机飞行速度。
杂波散射体s的归一化多普勒频率表示为fd:
其中,v为机载雷达载机飞行速度,ψ为机载雷达载机速度方向与杂波散射体s之间的夹角,λ为机载雷达发射电磁波的波长,fr表示机载雷达的脉冲重复周期,
杂波脊的定义如下:由载机平台运动引起的反映多普勒频率和空域频率之间耦合关系的轨迹。杂波脊由雷达系统参数和杂波块相对于载机平台的角度唯一确定,其与雷达天线布置形式、地表覆盖内型、目标污染等因素无关,因此可以将杂波脊当做一种先验信息来运用。
定义杂波脊上归一化多普勒频率为fd,m的杂波散射体的时域导向矢量为b(fd,m),其表达式为:
其中,m∈{1,2,…,nc},fd,m表示杂波脊上第m个杂波散射体的归一化多普勒频率,m表示一个相干处理间隔内机载雷达发射的脉冲个数,[·]t表示矩阵转置运算,e表示自然对数的底数;j表示虚数单位,
定义杂波脊上归一化空域频率为fs,m的杂波散射体的空域导向矢量为a(fs,m),其表达式为:
其中,n表示机载雷达的天线阵列包含的阵元总数,[·]t表示转置运算,m∈{1,2,…,nc},fs,m表示杂波脊上第m个杂波散射体的归一化空域频率。
定义杂波脊上归一化多普勒频率为fd,m,归一化空域频率为fs,m的杂波散射体的空时二维导向矢量为v(fd,m,fs,m),其表达式为:
其中,fd,m表示杂波脊上第m个杂波散射体的归一化多普勒频率,fs,m表示杂波脊上第m个杂波散射体的归一化空域频率,m∈{1,2,…,nc},b(fd,m)表示杂波脊上归一化多普勒频率为fd,m的杂波散射体的时域导向矢量,a(fs,m)表示杂波脊上归一化空域频率为fs,m的杂波散射体的的空域导向矢量,
空域导向矢量矩阵a的表达式为:
空时导向矢量矩阵v的表达式为:
其中,a(fs,m)表示杂波脊上归一化空域频率为fs,m杂波散射体对应的空域导向矢量,v(fd,m,fs,m)表示杂波脊上归一化多普勒频率为fd,m、归一化空域频率为fs,m杂波散射体对应的空时二维导向矢量,m∈{1,2,…,nc},nc表示距离环上划分的杂波块数。
步骤2,确定第l个待检测单元数据xl在杂波脊上的初始杂波功率矩阵为
第l个距离门数据xl(mn×1)可以表示如下:
其中,yk表示第k个脉冲空域采样数据,其中t表示矩阵转置运算。
根据待检测单元的空域采样数据可以计算待检测单元的空域协方差矩阵
其中,yk表示第k个脉冲空域采样数据,上标h表示矩阵复共轭转置运算。
确定第l个待检测单元数据xl在杂波脊上第m个杂波块的初始杂波功率为
初始化:令i表示第i次迭代,i的初始值为1,l∈{1,2,…,nmax},l的初始值为1;nmax表示雷达回波数据包含的nmax个距离门雷达回波数据总数。
进而得到第l个待检测单元数据xl在杂波脊上的初始杂波功率矩阵为
其中,rnc×nc表示nc×nc维实矩阵,m∈{1,2,…,nc},nc表示等距离环上独立杂波块数。
计算中间变量d(nc×1)的值,其表达式为:
其中,a(fs,m)表示归一化空域频率为fs,m对应的第m个杂波块的空域导向矢量,上标h表示复共轭运算,tr表示求矩阵的迹运算。m∈{1,2,…,nc},nc表示等距离环上独立杂波块数。
计算第i-1次迭代后中间变量ρ(nc×1)的值ρ(i-1),其表达式为:
其中,||·||表示矩阵的2范数,a(fs,m)表示归一化空域频率为fs,m对应的第m个杂波块的空域导向矢量,
步骤3,分别计算第i次迭代后第l个待检测单元数据xl在杂波脊上的空域协方差矩阵杂波功率矩阵
经过第i次迭代后第l个待检测单元数据xl在杂波脊上第m个杂波块的功率
其中,a(fs,m)表示归一化空域频率为fs,m对应的第m个杂波块的空域导向矢量,m∈{1,2,…,nc},
经过第i次迭代后第l个待检测单元数据xl在杂波脊上第m个杂波块的功率矩阵
其中,
进而计算得到第i次迭代后第l个待检测单元数据xl的空域重构协方差矩阵
其中,a表示n×nc维空域导向矢量矩阵,
步骤4,若i<γ,则令i的值加1,返回步骤3;其中,γ为预先设定的最大迭代次数,本例中取10。
若i≥γ,则迭代停止,并将迭代停止时对应的第i次迭代后第l个待检测单元数据xl在杂波脊上的杂波功率矩阵
步骤5,计算第l个待检测单元数据xl重构的空时二维协方差矩阵
其中,v表示空时二维导向矢量矩阵,
步骤6,令l的值加1,返回步骤3,直到得到第nmax个待检测单元数据
至此,本发明的基于知识辅助稀疏迭代协方差估计的杂波抑制的方法结束。
以下通过仿真实验对本发明效果作进一步验证说明。
实施例6
基于知识辅助稀疏迭代协方差估计的机载雷达杂波抑制方法同实施例1-5。
(一)仿真实验数据说明
为了验证本发明方法的准确性,通过仿真实验予以证明;仿真实验参数如表1所示:仿真实验参数
(二)仿真结果及分析
本发明的仿真结果分别如图3(a)、图3(b)、图4(a)、图4(b)、图5(a)和图5(b)所示,图3(a)为正侧阵时最优处理方法得到的第100号距离门空时二维谱示意图,图3(b)为正侧阵时本发明方法得到的第100号距离门空时二维谱示意图,图4(a)为30度安装角时最优处理方法得到的第100号距离门空时二维谱示意图,图4(b)为30度安装角时本发明方法得到的空时二维谱示意图,其中横坐标都为归一化多普勒频率,纵坐标都为归一化空域频率;图5(a)为正侧阵时本发明方法与最优处理方法改善因子对比示意图,图5(b)为30安装角时本发明方法与最优处理方法改善因子对比示意图;其中横坐标都为归一化多普勒频率,纵坐标为都改善因子(db)。
从图3(a)与图3(b)的对比和图4(a)与图4(b)的对比可以看出,无论机载雷达天线以正侧阵方式还是30度安装角方式安置时,运用本发明方法估计出来的待检测单元杂波谱都与最优杂波谱很相似,图3(a)与图3(b)中杂波谱从功率分布上与长度上基本保持一致,图4(a)与图4(b)中杂波谱的拐角坐标与杂波脊的曲线弧度也基本相近,说明本发明方法能够有效估计待检测单元的杂波协方差矩阵,杂波谱的估计精度很高,从而保证了机载雷达对杂波的抑制能力,证明了本发明在杂波协方差矩阵的估计过程中,运算量的减小并没有影响到杂波谱的估计精度,反而使得本方法与最优处理方法两者的杂波谱十分接近,保证了运算量的同时也保证了机载雷达对杂波的抑制能力。
实施例7
基于知识辅助稀疏迭代协方差估计的机载雷达杂波抑制方法同实施例1-6,仿真条件和内容同实施例6,为了更直观明显的看到杂波抑制的效果,并显示出对慢速运动目标的检测能力,仿真对本发明方法的改善因子和最优处理方法的改善因子作出了对比,从图5(a)和图5(b)可以看出,无论机载雷达天线以正侧阵方式还是30度安装角方式安置时,本发明方法的改善因子和最优处理方法的改善因子基本重合,这说明采用本发明方法的机载雷达能获得良好的杂波抑制性能,对待检测单元的杂波协方差矩阵的估计精度很高,并且本发明方法的改善因子的凹口很窄,也因此证明本发明显著提高了对慢速运动目标的检测能力。
以上仿真实验验证了本发明在减小运算量与运算时间的同时,也提高了机载雷达良好的杂波抑制性能,保证了本发明的正确性,有效性和可靠性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围;这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
综上所述,本发明公开的基于知识辅助稀疏迭代协方差估计的杂波抑制方法,解决了传统空时自适应处理技术(stap)由于机载雷达所工作的杂波环境的非均匀性而使杂波抑制性能显著下降的问题,实现步骤为:计算空时导向矢量矩阵;确定初始杂波功率矩阵并构造中间变量;利用中间变量计算迭代中的杂波功率矩阵;进行迭代并得到最终杂波功率矩阵;确定一个待检测单元数据重构的空时二维协方差矩阵和与其对应的权;重复上述过程,直到得到所有待检测单元数据的权并计算空时自适应处理结果。本发明利用待检测单元的数据重构待检测单元的杂波协方差矩阵,可以避免训练样本的非均匀性,有效抑制复杂强地杂波,显著改善了对运动目标的检测性能,同时保证了运算量小,实时性更好,在工程上比较容易实现,适用于非均匀环境下机载预警雷达抑制强地杂波以及检测地面慢速运动目标。