正交频分复用雷达超分辨距离与角度值联合估计方法
【技术领域】
[0001] 本发明属于通信技术领域，更进一步涉及雷达阵列信号处理技术领域中的一种正 交频分复用（orthogonal-frequency-divisionmultiplexing,OFDM)雷达超分辨距离与角 度值联合估计方法。采用本发明估计得到的距离值和角度值，可用于正交频分复用雷达对 空间动目标进行定位与跟踪。
【背景技术】
[0002] 波达角（directionofarrival，D0A)估计是利用处于空间中不同位置的多个传 感器阵列估计信号的空域参数。基于傅立叶波束形成的方法是由一组阵元接收得到信号样 本，对样本进行傅立叶变换来估计波达角。然而，由于傅里叶变换方法受到瑞利限的限制， 导致其波束有高旁瓣并且分辨率也较低。基于数据二阶统计特性的超分辨方法能够突破阵 列孔径的限制，其中，多重信号分类方法在理论上能够得到空间目标渐进无偏的角度估计， 在一定的信噪比门限上，该算法的性能接近于最大似然（maximumlikelihood，最大似然方 法）的方法。它利用信号子空间与噪声子空间的正交性，构造空间的谱函数，通过搜索与噪 声子空间正交的阵列流行矢量估计D0A。
[0003] 西安电子科技大学在其申请的专利"基于实数域广义多重信号分类算法的目标低 仰角估计方法"（专利申请号201310296530. 0,公开号CN103364772A)中公开了一种基于实 数域广义多重信号分类算法的目标低仰角估计方法。该方法的具体步骤是：（1)对雷达接 收回波进行采样；（2)利用采样数据计算采样协方差矩阵；（3)对采样协方差矩阵进行空间 平滑和酉变换，得到实数域协方差矩阵；(4)对实数域协方差矩阵进行特征值分解，获得噪 声投影矩阵；（5)构造实数域导向矢量流行；(6)利用噪声投影矩阵和实数域导向矢量流行 构造空间谱，对空间谱二维角度搜索，获得初次角度估计值；（7)利用初次角度估计值估计 多径衰减系数，构造二次空间谱，通过二维角度搜索得到角度估计值；(8)比较估计值中两 个角度大小，将最大角度作为目标仰角值。该方法降低了运算复杂度，提高了雷达在低信噪 比下的角度估计性能。但是，该方法仍然存在的不足之处是，在步骤（1)和步骤（2)中是利 用多次快拍的采样数据估计采样协方差矩阵，不利于目标的实时跟踪。
[0004] 西安电子科技大学大学在其申请的专利"基于空间平滑协方差矩阵稀疏表示的波 达方向角估计方法"（专利申请号201410280744. 3,公开号CN104020439A)中公开了一种基 于空间平滑协方差矩阵稀疏表示的波达方向角估计方法。该方法的具体步骤是：（1)采用 天线接收机形成均匀线阵；（2)采用空间平滑技术计算阵列输出的空间平滑协方差矩阵； (3)将空间平滑协方差矩阵矢量化，得到稀疏模型向量；(4)将空域网格划分，构造超完备 基；（5)根据稀疏模型向量和超完备基的稀疏表示关系，建立约束优化方程；(6)采用凸优 化方法求解约束优化方程得到最优估计；（7)根据最优估计值绘制幅度谱图，获得波达方 向角度值。该方法提高了无源测向的运算速度及低信噪比下对相干信号源的估计性能。但 是，该方法仍然存在的不足之处是，步骤（1)中接收阵列接收到的信号载频均相同，故不能 利用信号的频率分集，导致信号的距离信息被忽略。
【发明内容】
[0005] 本发明的目的在于克服上述已有技术的不足，提出一种正交频分复用雷达超分辨 距离与角度值联合估计方法。本发明是基于在窄带信号条件下，正交频分复用雷达的回波 信号中目标的距离和角度值具有无耦合特性，基于直接数据域方法在距离角度二维域进行 二维平滑，实现了单快拍条件下的目标超分辨距离角度值估计。
[0006] 实现本发明目的的基本思路是：首先，正交频分复用雷达发射正交频分复用波形 信号；其次，正交频分复用雷达接收经空间动目标反射后的回波信号，对回波信号进行向 量化表示；然后，利用二维空间平滑方法取样；最后，联合估计空间动目标的距离值和角度 值。
[0007] 为实现上述目的，本发明具体实现步骤如下：
[0008] (1)发射正交频分复用波形信号：
[0009] 正交频分复用雷达的一维等距线阵中的阵元，同时发射正交频分复用波形信号， 每个正交频分复用波形信号包含相同的子载波个数，每个相邻的子载波频率间隔相同；
[0010] (2)向量化表示：
[0011] (2a)正交频分复用雷达的一维等距线阵，接收由正交频分复用波形信号经空间动 目标反射后的回波信号；
[0012] (2b)正交频分复用雷达对回波信号进行一次快拍，得到回波数据；
[0013] (2c)使用滤波器，滤除回波数据中的载波；
[0014] (2d)将回波数据中的时间取样离散化后以向量形式表示；
[0015] ⑶相位补偿：
[0016] 使用乘法器，对回波数据中的每一列进行相位补偿，得到含有空间动目标的距离 信息和角度信息的回波矩阵；
[0017] (4)二维空间平滑取样：
[0018] (4a)提取回波矩阵的第一列数据重新排列成一个矩阵，重新排列后的矩阵的行数 为正交频分复用波形信号子载波的个数，列数为正交频分复用雷达的一维等距线阵阵元的 个数；
[0019] (4b)采用二维循环滑窗取样方法，对重新排列后的矩阵进行取样，得到样本子 阵；
[0020] (4c)将样本子阵重新排列为一个矩阵，将该矩阵作为米样样本；
[0021] (5)按照下式，计算空间动目标的二维搜索矢量：
[0022] a", (OM) =aWv ((R)
[0023] 其中，aw(0，R)表示空间动目标的二维搜索矢量，W表示采样样本的协方差矩阵特 征值的个数，且W=NsXMs，Ns表不滑窗所用窗的行数，Ms表不滑窗所用窗的列数，0表 示空间动目标的波达角度，R表示空间动目标与正交频分复用雷达的距离，aMs(0)表示空 间动目标的方向导向矢量，?表示做克罗内克相乘操作，aNs (R)表示空间动目标的距离导 向矢量；
[0024] (6)计算包含空间动目标的距离信息和角度信息的谱强度：