，在目标的先验信息未知的情况下，最大后验概率（MAP)准则就 是寻找最合适的f满足下式：
[0131]
(16)
[0132] 其中，/为目标信息在最大后验概率准则下得到的解；在先验信息最少的情况下， 使用均匀分布函数描述目标在场景中的分布情形；在扫描雷达体系中，方位向回波数据中 与目标直接相关的目标数较少，可以作为统计学中的大观察数据下的小样本事件，因而需 要使用泊松分布函数来描述这种现象的统计特性，得出下式：
[0133]
(1Π
[0134] 对公式（17)两边进行取对数运算，得到：
[0135]
(IS)
[0136] 为了得到公式（18)的最大值，就需要对公式（18)进行梯度运算，并令结果为0,表 达式为：
[0137]
(If)
[0138] 有公式（19)可以推出如下表达式：
[0139]
C 20)
[0140] 这里，k表示迭代次数。
[0141] 步骤62、设置门限threshold，控制天线方向图的截断位置；
[0142] 步骤63、改变天线方向图插值的点数Neha，控制波束的宽度，得到不同的天线方向 图；
[0143] 步骤64、将不同的天线方向图代入步骤61推导的公式中进行卷积反演，得到不同 的反演结果图；
[0144] 步骤65、计算图像熵，基于最小图像熵准则确定合适的天线方向图。所述步骤65 中的图像熵定义为：
[0145]
(20
[0146] 其中，^为像素的灰度，P(X1)为各灰度级出现的概率。
[0147] 步骤66、计算得出最终结果。
[0148] 在原始天线方向图的基础上，我们首先设置门限范围为threshold = 0. 1+(κ-1)*0. 01，控制天线方向图的截断位置；然后改变天线方向图插值的点数N^= Nbeta-50+(J_l)*5,控制波束的宽度，得到不同的天线方向图；将不同的天线方向图应用 反卷积方法进行卷积反演，得到不同的反演结果图；再利用图像熵的性质，得到在K = 55, threshold = 0. 64和J = 18, Ncha= 435时，反卷积结果的熵最小，为H = 2. 066,如图 7所示。得到最终的天线方向图，如附图8所示。最后，将得到的最优的天线方向图应用到 反卷积方法中进行卷积反演，实现扫描雷达角超分辨成像。最后得到的最终结果图9所示
[0149] 本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发 明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的 普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各 种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。
【主权项】
1. 一种最优天线方向图选取的雷达角超分辨方法，其特征在于，包括如下步骤： 步骤1、前视扫描雷达方位回波建模； 步骤2、回波数据距离向脉冲压缩； 步骤3、距离走动判断； 步骤4、距离走动校正； 步骤5、扫描雷达角方位向回波卷积模型； 步骤6、基于最优天线方向图选取进行的反卷积。2. 如权利要求1所述的最优天线方向图选取的雷达角超分辨方法，其特征在于，所述 步骤6包括如下步骤： 步骤61、公式推导； 步骤62、设置门限threshold，控制天线方向图的截断位置； 步骤63、改变天线方向图插值的点数^ha，控制波束的宽度，得到不同的天线方向图； 步骤64、将不同的天线方向图代入步骤61推导的公式中进行卷积反演，得到不同的反 演结果图； 步骤65、计算图像熵，基于最小图像熵准则确定合适的天线方向图； 步骤66、计算得出最终结果。3. 如权利要求2所述的最优天线方向图选取的雷达角超分辨方法，其特征在于：所述 步骤1中前视扫描雷达方位回波建模过程为，由于基于机载雷达，所以载机平台速度为V ; 雷达天线的扫描速度为波束俯仰角为0 ;目标方位角为承；场景中目标到雷达天线的 初始斜距记为R。。经过时间t，载机平台与场景中位于（x，y)点处目标的距离，记为R(t); 此时目标到雷达的斜距表示为：记发射信号的载频为f。，脉冲重复时间为PRI ;假设雷达发射信号为：其中，rect( ?)表示矩形信号，其定义为Tp为发射脉冲持续时 间，k为调频斜率，T是快时间； 对于成像区域Q，回波可以表示为发射信号与目标的卷积加上噪声的结果，为了保证 理论与实际验证情况相符，就需要对距离向和方位向进行了离散处理，则回波的解析表达 式可以写成：其中，E为求和运算，（x，y)为场景中目标的位置；f(x，y)为点（x，y)处目标的散射函 数；为慢时间域的窗函数，表示天线方向图函数在方位向的调制；n a。为天线方位角初 始时刻；Tfi是目标在3dB天线波束宽度的驻留时间；c为电磁波传播速度；e Jt，n)表示 回波中的噪声。4. 如权利要求2所述的最优天线方向图选取的雷达角超分辨方法，其特征在于：所述 步骤2回波数据距离向脉冲压缩过程为，按照脉冲压缩原理，构造距离向脉冲参考信号的 公式为：这里，Traf表示距离向参考时间，k表示发射信号调频斜率。将Sg与回波数据 gl( T，n)进行最大自相关运算，实现回波信号在距离向的脉冲压缩；将公式（4)代入公式 (3)实现回波信号在距离向的脉冲压缩，脉冲压缩后的信号可以表示为：其中，B为发射信号带宽，e2( T，q)为yi( T，q)在进行脉冲压缩操作后引入系统的 噪声。5. 如权利要求2所述的最优天线方向图选取的雷达角超分辨方法，其特征在于：所述 步骤3距离走动判断的过程为，根据步骤1中的公式（1)，可得距离走动量公式为： Ag =Vl\cos(9COS^ (6)为波束扫描驻留时间，9 teta为天线波束宽度。判断其是否跨越距离单元其中fs为距离向采样率； 若满足A R〈 A r，直接进行步骤5 ;若A R> A r，则进入步骤4。6. 如权利要求2所述的最优天线方向图选取的雷达角超分辨方法，其特征在于：所述 步骤4距离走动校正的过程为，对数据&(1，n)进行尺度变换，得到数据平面内，消除距 离走动后的信号表达式为：其中，e3( T，q)为y2( T，q)进行距离走动校正操作后引入系统的总噪声。7. 如权利要求2所述的最优天线方向图选取的雷达角超分辨方法，其特征在于：所述 步骤5扫描雷达角方位向回波卷积模型的过程为，将扫描雷达方位向回波建立为天线方向 图和目标散射系数的卷积模型；为了数学计算上的方便，将公式（5)和公式（7)转化成矩阵 与向量的运算形式，即得： y = Af+e ； (8) 其中： y = [y3(Ti, rI1)-y3(T!, nM)；y3(T2, n2) -y3(t2, nM)； ;y3 (T n,rI i) ??? g4(T n, nM)]T； (9) f = !^(X1J1X(X1Jk) (X2J2) ...f(x2，yK) 〇(XnJ1) ...f(xN，yK)]T; (10) e - [ejTj，riM) ;e3(T2，ri2)...e3(T2，ri M) ; ； e3 ( x N, rI j) *** e3(xN, nM)]T； (11) 其中，M是方位向回波的采样点数，N是距离向回波的采样点数，K是目标的数目，有K =M+L-l，L为天线方向图的采样点数目。上标T表示转置运算；另有以下已知公式：其中，E{ ?}表示期望值运算；因此，（12)中的卷积矩阵A结构如下：其中，矩阵A是一个导向矩阵，/4 (A1,/?:,…為)；，.h# 0, 1 =1，…，L为天线向量； 因此，实波束扫描雷达角超分辨成像可转化为：给定公式（8)中y和A，求解f ;至此，扫 描雷达角超分辨成像问题转换化卷积反演问题；解卷积方法可以表述为寻找线性因子K = A\使得下式成立：这里，F(取),F(ST) J(W),#(取）分别对应⑷中f，y，A，e的傅里叶变换。8.如权利要求2所述的最优天线方向图选取的雷达角超分辨方法，其特征在于：所述 步骤61公式推导过程为，对于公式（8)，回波数据的后验概率可以表示为：其中，P (f I y)，P (y I f)和P (f)分别代表回波数据的后验概率，似然函数和先验概率；由 于回波数据是已知的，在目标的先验信息未知的情况下，最大后验概率（MAP)准则就是寻 找最合适的f满足下式：其中，j为目标信息在最大后验概率准则下得到的解；在先验信息最少的情况下，使 用均匀分布函数描述目标在场景中的分布情形；在扫描雷达体系中，方位向回波数据中与 目标直接相关的目标数较少，可以作为统计学中的大观察数据下的小样本事件，因而需要 使用泊松分布函数来描述这种现象的统计特性，得出下式：为了得到公式（18)的最大值，就需要对公式（18)进行梯度运算，并令结果为0,表达式 为：这里，k表示迭代次数。9.如权利要求2所述的最优天线方向图选取的雷达角超分辨方法，其特征在于，所述 步骤65中的图像熵定义为：其中，^为像素的灰度，P (X J为各灰度级出现的概率。
【专利摘要】本发明公开了一种最优天线方向图选取的雷达角超分辨方法，包括前视扫描雷达方位回波建模；回波数据距离向脉冲压缩；距离走动判断；距离走动校正；扫描雷达角方位向回波卷积模型；基于最优天线方向图选取进行的反卷积。最后一步包括公式推导；设置threshold，控制截断位置；改变点数Ncha，控制宽度，得到天线方向图；将天线方向图代入推导的公式中卷积反演，得到反演结果图；计算图像熵，确定天线方向图；计算得出最终结果。本发明的有益效果：首先是从代数学的角度分析了天线方向图误差对卷积反演结果的影响，然后基于图像熵的方法确定了最优的天线方向图，将最优天线方向图带入发明提出的反卷积方法中，实现卷积反演问题的求解，实现了扫描雷达角超分辨成像。
【IPC分类】G01S13/89, G01S7/41, G01S13/28
【公开号】CN105137408
【申请号】CN201510435711
【发明人】黄钰林, 王月, 周小军, 任建宇, 张永超, 杨建宇
【申请人】电子科技大学
【公开日】2015年12月9日
【申请日】2015年7月23日