一种多基地雷达多目标成像定位方法

一种多基地雷达多目标成像定位方法
【专利摘要】本发明涉及一种多基地雷达多目标成像定位方法，其包括：初始化系统参数，建立二维回波数据矩阵；脉冲压缩处理；构造回波信号测量矩阵；设定自适应OMP算法的初始参数；寻找回波信号测量矩阵与重构残余误差的最大相干项；更新索引集合；重构双程距离单元集合上的目标散射系数向量；计算重构残余误差；计算平均重构残余误差和平均散射系数向量；计算平均重构残余误差变化率和平均散射系数变化率；自适应迭代终止判定。本发明利用平均重构残余误差变化率和平均散射系数变化率作为自适应OMP算法迭代终止的判定条件，利用自适应OMP算法迭代终止次数估计探测区域目标个数，适用于实际情况下目标个数未知的多基地雷达多目标成像定位。
【专利说明】
一种多基地雷达多目标成像定位方法
技术领域
[0001 ]本发明属于多基地雷达技术领域，更具体地说是涉及一种基于自适应0ΜΡ的多基 地雷达多目标成像定位方法。
【背景技术】
[0002] 多基地雷达多目标定位中为实现对多个目标的可靠定位，需要首先完成数据关 联。然而，复杂的数据关联算法制约了多基地雷达多目标定位的实时性和准确性。为了解决 数据关联问题，多基地雷达多目标成像定位方法将多目标定位问题看作三维成像问题，通 过将多个接收机回波信号的幅度信息投影到三维图像空间，再从图像空间中逐一提取各目 标，实现对多个目标的定位。图像空间的三维分辨率分别由接收机阵列和线性调频信号决 定，其中，接收机阵列可看作二维稀疏阵列，获得二维分辨率;接收机对线性调频信号进行 脉冲压缩，获得距离向上的分辨率。
[0003] 多基地雷达多目标成像定位方法虽然解决了数据关联问题，但也存在以下两个不 足。第一个是，由于线性调频信号脉冲压缩后的输出波形近似为Sine函数(详见"皮亦鸣，杨 建宇，付毓生，杨晓波.合成孔径雷达成像原理.第一版.电子科技大学出版社.2007.3"），存 在主瓣展宽和旁瓣串扰问题，导致距离向分辨率下降，从而严重影响三维图像空间中的成 像质量，降低该定位方法的正确性和准确性。第二个不足是，该方法从图像空间提取目标 时，需已知探测区域目标个数，而这一条件在实际应用中通常难以实现。
[0004]压缩传感稀疏信号重构是一种近几年来提出的新兴的信号处理理论，其主要思想 是采用非自适应线性投影来保持信号的原始结构，通过数值最优化问题准确重构出原始信 号。Donoho指出只要信号存在稀疏性或可压缩性，就可用远低于Nyquist采样率的采样信号 恢复出原始信号（详见参考文献"D.Donoho.Compressed sensing. IEEE Trans · Inf ·Theory，vo 1 · 52，no ·4，pp: 1289-1306，April 2006"）。在多基地雷达多目标场景 中，由于三维空间中大量区域不包含散射点（目标），回波信号在三维空间中具有稀疏性，因 此，可以利用压缩传感稀疏重建的方法，估计和重构出空间单元上原始信号，抑制主瓣展宽 和旁瓣串扰问题，提高定位精度。正交匹配追踪(0ΜΡ)算法是稀疏信号重构理论中的经典算 法，详见参考文南犬"J.A.Tropp，A.C.Gilber.Signal Recovery from Random Measurements via Orthogonal Matching Pursuit. IEEE Transactions on Information Theory, vol·53，no· 12，pp:4655-4666，2007/' AMP算法基本思想是利用贪婪迭代追踪方法，在每一 次迭代中选定回波信号测量矩阵中与重构残余误差的最大相干项作为索引原子，并对所选 定原子矩阵进行正交化以保证迭代结果最优。0ΜΡ算法具有结构简单、计算复杂度低和运算 时间快等优点，在图像处理、医学成像和无线通信等领域中得到了成功应用。但是，0ΜΡ算法 的前提是需要已知原始稀疏信号中非零元素的个数(对应多基地雷达多目标成像定位中， 即需要已知探测区域中目标个数），如果原始稀疏信号中非零个数未知，该算法重构性能会 严重下降。也就是说，要采用0ΜΡ算法抑制线性调频信号脉冲压缩后的主瓣展宽和旁瓣串扰 问题，提高定位精度，前提条件也是假设探测区域目标个数已知。这一假设条件和多基地雷 达多目标成像定位方法从图像空间提取目标时假设目标个数已知的条件相同。然而，在多 基地雷达多目标定位的实际情况中，探测区域目标个数通常是未知的，因此，为了将多基地 雷达多目标成像定位方法应用于实际中的同时改善成像质量，提高定位精度，就必须克服 OMP算法及多基地雷达多目标成像定位方法对已知原始稀疏信号中非零元素个数的约束 (即：对已知探测区域中目标个数的约束）。

【发明内容】

[0005] 本发明是为了克服上述的多基地雷达多目标成像定位方法中探测区域目标个数 对正交匹配追踪算法及图像空间目标提取过程的约束，而提出一种基于自适应0ΜΡ的多基 地雷达多目标成像定位方法，该方法利用多基地雷达多目标场景中多个接收机的平均重构 残余误差变化率和平均散射系数变化率代替探测区域中目标个数作为自适应0ΜΡ算法迭代 终止的判定条件，利用自适应0ΜΡ算法迭代终止次数估计探测区域目标个数，适用于实际情 况中对目标个数未知的多基地雷达多目标成像定位。
[0006] 上述的多基地雷达多目标成像定位方法，其具体步骤包括：（1)初始化系统参数， 建立二维回波数据矩阵；（2)脉冲压缩处理；（3)构造回波信号测量矩阵；（4)设定自适应0ΜΡ 算法的初始参数；（5)寻找回波信号测量矩阵与重构残余误差的最大相干项；（6)更新索引 集合；（7)重构双程距离单元集合上的目标散射系数向量；（8)计算重构残余误差；（9)计算 平均重构残余误差和平均散射系数向量；（10)计算平均重构残余误差变化率和平均散射系 数变化率;（11)自适应迭代终止判定。
[0007] 所述多基地雷达多目标成像定位方法，其中，所述步骤（1)中初始化参数均为已 知，初始化参数具体如下：由一个发射机和S个接收机组成系统，其中S多3，S个接收机分散 布置在地面，以发射机位置为原点，在目标、发射机和接收机所在的地理空间建立X-Y-Z笛 卡尔直角坐标系，发射机坐标位置记为T= [0，0，0 ]τ，接收机的坐标位置分别为rj = [Xj，yj， Zj]T，j = l，2,…，S，[]T表示矩阵转置，系统在X方向、Y方向和Z方向的最远探测距离分别为 XXmax、YY_^PZZmax，最近探测距离分别为XXmin、YY min和ZZmin;&射机对探测区域宽波束泛光 发射线性调频信号，发射脉冲载频为，带宽B，脉冲宽度!^，距离分辨率P;发射机及各个接 收机在时间上保持同步，以发射机发射线性调频信号时刻为系统时间零点；接收机宽波束 泛光接收回波信号，接收机采样频率为f s，采样点数为Nrang，将S个接收机接收的回波信号 存储在二维数据矩阵echo = {echo (i，j)}中，其中，i = l，2，···，Nrang表示双程距离单元，j =1，2，…，S表示接收机，假设目标坐标为Pt，则该目标的回波信号将被存储在双程距离单 元1 =。6：[1(?)，」_ = 1，2，~，5中，其中，；1^=||1'-?*||2+||?1；-：1^||2为发射机到目标再到第 j个接收机的双程距离，ceil表示向上取整，I I I |2表示L2范数；
[0008] 所述步骤(2)脉冲压缩处理的具体步骤为:取出所述步骤（1)的数据矩阵echo,利 用脉冲压缩方法对echo的每一列数据echo」，j = 1，2，…，S进行脉冲压缩处理，得到脉冲压 缩处理后的列向量eeh% .，j = 1，2，…，S，将所有接收机脉冲压缩处理后的列向量_ech〇/, j = 1，2，."，S存储在二维数据矩阵ECHO = {eCHO(/,./)}中，其中i = l，2，."，Nrang;j = l，2，···， S；
[0009] 所述步骤(3)构造的回波信号测量矩阵为Θ = { Θ (m，c)}，其中m=l，2，-_，Nrang; c = l，2，."，Nrang，Θ (m，c)具体表达式为：
[0010]
[0011] 所述多基地雷达多目标成像定位方法，其中：所述步骤(4)的具体步骤为:设定自 适应0ΜΡ算法重构迭代处理的最大迭代次数，记作K;设定平均重构残余误差变化率门限，记 作％;设定平均散射系数变化率初始门限，记作;令Ω = {1，…，Nrang}表示双程距离单元 集合;设定第j，j = 1，…，S个接收机在双程距离单元集合Ω上的目标散射系数向量的初始 值，记为Xj，》_e CV_xl为Nrange行1列的复数向量，Xj的初始值选择为Xj = 0;设定第j，j = 1，…，S个接收机的重构残余误差初始迭代值记为为Nrange行1列的复 数向量，〃的值选择为，其中为所述步骤(2)获得的脉冲压缩处理后的 第j列向量;设定第j，j = 1，…，S个接收机索引集合的初始迭代值，记为Sf;;，的值选择为 =0，其中0表示空集;k表示自适应0ΜΡ算法中的第k次迭代，k的初始值设置为k = 1。
[0012] 所述多基地雷达多目标成像定位方法，其中，所述步骤(5)的具体步骤为:采用公 式
j计算得到自适应0ΜΡ算法第k次迭代过程中测量矩阵与第j个 接收机重构残余误差的最大相干项，其中ai"g (；)表示提取括号内向量中的最大值元素 对应的行序号m值，abs( ·)表示对括号内向量中的元素求模，Θ为所述步骤(3)中得到的回 波信号测量矩阵，右上角T为转置运算符号，^^广为算法第k-Ι次迭代过程中得到的重构残 余误差，k表示自适应0ΜΡ算法中的第k次迭代次数;若k= 1，的值为所述步骤(4)中得 到的初始迭代值否则〃 <广通过自适应0ΜΡ算法第k-Ι次迭代过程中所述步骤(8)的 重构残余误差计算提供。
[0013] 所述多基地雷达多目标成像定位方法，其中，所述步骤(6)的具体步骤为:采用公 式咩卜巧丨1) U/)计算得到自适应0ΜΡ算法第k次迭代过程中的第j个接收机索引集合，记为 ^丨，其中为自适应0ΜΡ算法第k-Ι次迭代过程中得到的第j个接收机索引集合，对为所 述步骤(5)计算得到的测量矩阵与第j个接收机重构残余误差的最大相干项，U表示并集运 算符合，k表示自适应0ΜΡ算法中的第k次迭代次数;若k = 1，Ξ $ U的值为所述步骤(4)中得到 的初始值Sg;k = 2时，Ξ^11的值为自适应0ΜΡ算法第1次迭代过程中得到的第j个接收机索 弓丨集合Ξ& ;k = 3时，Ξ&11的值为自适应0ΜΡ算法第2次迭代过程中得到的第j个接收机索引 集合Ξ丨;;依次类推，在第k次迭代中且k> 1时，Ξ&"的值为自适应0ΜΡ算法第k-Ι次迭代过程 中得到的第j个接收机索引集合。
[0014] 所述多基地雷达多目标成像定位方法，其中：所述步骤(7)的具体步骤为：重构第 j，j = 1，…，S个接收机在双程距离单元集合Ω上的目标散射系数向量，具体为:采用公式
计算得到自适应ΟΜΡ算法第k 次迭代过程中第j个接收机在双程距离单元集合Ω上的散射系数向量，记为Xj，其中卩丨为 第j个接收机在双程距离单元集合Ω上的散射系数向量中由索引集合对应的元素组 成的向量，为所述步骤(6)中第k次迭代过程中得到的第j个接收机索引集合，·丨为第j 个接收机在双程距离单元集合Ω上的散射系数向量中索引集合对应的元素组成的向 量，司丨;为索引集合在双程距离单元集合Ω上的补集，（.f表示对括号中矩阵求伪逆矩 阵，右上角?为矩阵伪逆运算符号"为测量矩阵Θ中由索引集合对应的列组成的矩 阵为所述步骤(2)获得的脉冲压缩处理后的第j列向量；
[0015] 所述步骤(8)的具体步骤为:采用公式^.7· - 计算得到自适应OMP算法 第k次迭代过程中第j个接收机的重构残余误差，记为〃sg，其中echo/为步骤(2)获得的脉 冲压缩处理后的第j列向量。
[0016] 所述多基地雷达多目标成像定位方法，其中，所述步骤(9)的具体步骤包括：
[0017] (9.1)采用公?!
:寸所述步骤(8)获得的S个接收机的重构 残余误差^，j = l，…，S求平均，计算得到平均重构残余误差向量，记为mean_res(k);
[0018] (9.2)采用公式801'1：_1」=801'1：(3&8(1」））对所述步骤(7)获得的第」个接收机在双 程距离单元集合Ω上的重构得到的散射系数向量中元素的模值按照升序从大到小排列，记 为S〇rt_Xj，其中sort( ·)表示将括号内向量中的元素根据其值的大小按照升序从大到小 排列，abs( ·)表示对括号内向量中的元素求模；
[0019] (9.3)采用公式
f所述步骤(9.2)获得的S个接收机排序后 Μ / 的散射系数向量sort_Xj，j = 1，…，S求平均，计算得到S个接收机在双程距离单元集合Ω上 的排序后散射系数的平均散射系数向量，记为mean_x。
[0020] 所述多基地雷达多目标成像定位方法，其中，所述步骤(10)的具体步骤为:采用公 式
-算得到第k次迭代过程中的平均重构残余误差变 - .2 化率，记为nk，其中mean_res(k)为所述步骤(9.1)计算得到的第k次迭代过程中S个接收机的 平均重构残余误差向量，I I · I U表示对向量求L2范数，| · |表示求绝对值，若k=l，mean_ 的值为
_中〃sg为步骤(4)中得到的初始迭代值;采用公 式 汁算得到第k次迭代过程中的平均散射系数变化率，记 为妒，其中mean_x (k)表示向量mean_x中的第k个元素，若k = 1，mean_x (k-Ι)的值为mean_x (0) =mean_x( 1) 〇
[0021 ]所述多基地雷达多目标成像定位方法，其中：所述步骤（11)的具体步骤为:如果nk 彡且k<K，贝ijk=k+i，返回所述步骤(5);否则自适应omp算法终止迭代，输出重构 的散射系数和目标个数，且此刻自适应0ΜΡ算法第k次迭代得到的散射系数向量&，j = 1，…，S即为S个接收机在双程距离单元集合Ω上最终重构获得的散射系数向量，第k次迭代 次数减去1即为自适应0ΜΡ算法对探测区域目标个数的估计值，即L = k-1，其中L表示探测区 域目标个数的估计值，n〇为所述步骤(4)中初始化得到的平均重构残余误差变化率门限，妒 为所述步骤(10)计算得到的第k次迭代过程中的平均散射系数变化率，β 1-1为所述步骤(10) 计算得到的第k-Ι次迭代过程中的平均散射系数变化率，若k=l，i3H的值为俨'βο,βο为所 述步骤(4)中初始化得到的平均散射系数变化率初始门限;将重构获得的散射系数向量 Xj， j = l，…，S存储在二维数据矩阵ECH0={ECH0(i，j)}中，其中i = l，2，".，Nrang;j = l，2，···， So
[0022] 有益效果：
[0023]本发明多基地雷达多目标成像定位方法构思合理，针对多基地雷达多目标场景在 三维空间上的稀疏特征，通过建立接收机回波信号与三维空间中目标散射系数的线性测量 矩阵，利用平均重构残余误差变化率和平均散射系数变化率作为正交匹配追踪(0ΜΡ)算法 迭代处理的迭代终止条件，提出了基于自适应0ΜΡ的多基地雷达多目标成像定位方法，克服 了传统正交匹配追踪算法在多基地雷达多目标成像定位中对探测区域目标个数的依赖，适 用于实际情况中未知目标个数时的多目标成像定位，同时提高了图像空间中的成像质量； 通过该方法的迭代终止次数能获得估计探测区域目标个数，该估计目标个数值可用于成像 定位中从图像空间中对目标的提取，克服了多基地雷达多目标成像定位方法目标提取过程 中对探测区域目标个数的依赖，更加适用于实际情况中对目标个数未知的多基地雷达多目 标成像定位。
[0024] 自适应0ΜΡ算法利用多基地雷达多目标场景中多个接收机的平均重构残余误差变 化率和平均散射系数变化率代替探测区域中目标个数作为0ΜΡ算法迭代终止的判定条件， 通过合理的设置平均重构残余误差变化率值门限，以及比较相邻两次迭代过程获得的平均 散射系数变化率值，自适应0ΜΡ在探测区域中目标个数未知时也能较精确地分离观测场景 中的主散射目标和弱散射背景，同时能较准确地估计出目标个数，因而可以应用于目标个 数未知时的多基地雷达多目标成像定位方法中。
【附图说明】
[0025] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可 以根据这些附图获得其他的附图。
[0026] 图1为本发明多基地雷达多目标成像定位方法的流程图；
[0027] 图2为本发明多基地雷达多目标成像定位方法的系统模型，其中，多个分散布置的 接收机可看作一个稀疏阵列；
[0028] 图3本发明多基地雷达多目标成像定位方法的投影成像示意图；（其中，A、B和F均 为回波信号经过脉冲压缩和稀疏重建后的幅度值，A和B为两个不同的目标，F为虚假目标，S 表示接收机个数，t表示时间轴，投影成像方法将同一目标在不同接收机中的幅度信息在图 像空间积累。）
[0029] 图4为本发明多基地雷达多目标成像定位方法的实施例1在监视区域中存在2个目 标情况下，第1个接收机脉冲压缩后获得的回波信号；
[0030] 图5为本发明多基地雷达多目标成像定位方法的实施例1在监视区域中存在2个目 标情况下，第1个接收机重构的目标散射系数；
[0031] 图6为本发明多基地雷达多目标成像定位方法的实施例1在监视区域中存在5个目 标情况下，第6个接收机脉冲压缩后获得的回波信号；
[0032] 图7为本发明多基地雷达多目标成像定位方法的实施例1在监视区域中存在5个目 标情况下，第6个接收机重构的目标散射系数；
【具体实施方式】
[0033] 如图1至7所示，本发明多基地雷达多目标成像定位方法，是基于自适应0ΜΡ的多基 地雷达多目标成像定位方法，其包括如下步骤：
[0034] S0010、初始化系统参数，建立二维回波数据矩阵
[0035] 初始化参数均为已知，初始化参数如下：系统由一个发射机和S个接收机组成，其 中S^：3，S个接收机分散布置在地面，以发射机位置为原点，在目标、发射机和接收机所在的 地理空间建立X-Y-Z笛卡尔直角坐标系，发射机坐标位置记为T = [ 0，0，0 ]τ，接收机的坐标 位置分别为0=[^^，^]^_ = 1，2，~，5，[]7表示矩阵转置，系统在父方向、￥方向和2方向 的最远探测距离分别为XXmax、YY_^PZZ max，最近探测距离分别为XXmin、YYmin和ZZmin;发射机 对探测区域宽波束泛光发射线性调频信号，发射脉冲载频为fo,带宽B，脉冲宽度T P，距离分 辨率ρ;发射机及各个接收机在时间上保持同步，以发射机发射线性调频信号时刻为系统时 间零点;接收机宽波束泛光接收回波信号，接收机采样频率为f s，采样点数为Nrang，将S个 接收机接收的回波信号存储在二维数据矩阵echo= {echo(i，j)}中，其中，i = 1，2,…， Nrang表示双程距离单元，j = 1，2，…，S表示接收机，假设目标坐标为Pt，则该目标的回波信 号将被存储在双程距离单元1 = 〇6：11(?)，」_ = 1，2，"，5中，其中，；1^=||1'-?1；||2+||?1；-1·」 2为发射机到目标再到第j个接收机的双程距离，ceil表示向上取整，| | | |2表示L2范数。
[0036] S0020、脉冲压缩处理
[0037]即取出步骤S0010的数据矩阵echo,利用脉冲压缩方法对echo的每一列数据 ech〇j，j = 1，2，…，S进行脉冲压缩处理，得到脉冲压缩处理后的列向量.ech〇_/，j = 1，2，…，S， 将所有接收机脉冲压缩处理后的列向量eck,,，j = l，2,…，S存储在二维数据矩阵 ECH() = -|ECHO(/，./)^A*i = l，2，."，Nrang;j = l，2，."，S。
[0038] S0030、构造回波信号测量矩阵，即 Θ = {Θ(ηι，(3)}，其中 m=l，2，···，Nrang;c = 1， 2，-"，Nrang，Θ (m，c)具体表达式为：
[0039]
[0040] S0040、设定自适应0ΜΡ算法（即正交匹配追踪算法)的初始参数
[0041 ]设定自适应0ΜΡ算法重构迭代处理的最大迭代次数，记作K;设定平均重构残余误 差变化率门限，记作no;设定平均散射系数变化率初始门限，记作;令Ω = {1，…，Nrang}表 示双程距离单元集合;设定第j，j = 1，…，S个接收机在双程距离单元集合Ω上的目标散射 系数向量的初始值，记为为Nrange行1列的复数向量，一般&的初始值选择为 Xj = 0;设定第j，j = 1，…，S个接收机的重构残余误差初始迭代值记为e 为Nrange行1列的复数向量，一般mg的值选择为f〇echo;·，其中为步骤S0020获得 的脉冲压缩处理后的第j列向量;设定第j，j = l，…，S个接收机索引集合的初始迭代值，记 为部;，見°;的值选择为'巧=0，其中0表示空集;k表示自适应0ΜΡ算法中的第k次迭代，k的 初始值设置为k=l。
[0042] S0050、寻找回波信号测量矩阵与重构残余误差的最大相干项
[0043] 即采用公另
I：得到自适应0ΜΡ算法第k次迭代过程 中测量矩阵与第j个接收机重构残余误差的最大相干项，其中表示提取括号内向 量中的最大值元素对应的行序号m值，abs( ·)表示对括号内向量中的元素求模，Θ为步骤 S0030中得到的回波信号测量矩阵，右上角T为转置运算符号，为算法第k-Ι次迭代过 程中得到的重构残余误差，k表示自适应0ΜΡ算法中的第k次迭代次数;若k= 1，的值为 步骤S0040中得到的初始迭代值〃#，否则通过自适应0ΜΡ算法第k-1次迭代过程中 步骤S0080的重构残余误差计算提供。
[0044] S0060、更新索引集合
[0045]即采用公式计算得到自适应0ΜΡ算法第k次迭代过程中的第j个接 收机索引集合，记为^;|，其中2?11为自适应0ΜΡ算法第k-Ι次迭代过程中得到的第j个接收 机索引集合，/；为步骤S0050计算得到的测量矩阵与第j个接收机重构残余误差的最大相干 项，U表示并集运算符合，k表示自适应0ΜΡ算法中的第k次迭代次数;若k = 1，11的值为步 骤S0040中得到的初始值S| ;k = 2时，Ξ&11的值为自适应0ΜΡ算法第1次迭代过程中得到的 第j个接收机索引集合424 = 3时，的值为自适应0ΜΡ算法第2次迭代过程中得到的第j 个接收机索引集合Ξ&依次类推，在第k次迭代中且k>l时，Ξ^1)的值为自适应0ΜΡ算法第 k-Ι次迭代过程中得到的第j个接收机索引集合。
[0046] S0070、重构双程距离单元集合上的目标散射系数向量
[0047]即重构第j，j = 1，…，S个接收机在双程距离单元集合Ω上的目标散射系数向量， 具体为:采用公式
丨丨计算得到自 适应0ΜΡ算法第k次迭代过程中第j个接收机在双程距离单元集合Ω上的散射系数向量，记 为&，其中'丨;丨为第j个接收机在双程距离单元集合Ω上的散射系数向量中由索引集合 〇寸应的元素组成的向量,?1为步骤S0060中第k次迭代过程中得到的第j个接收机索引 集合，为第j个接收机在双程距离单元集合Ω上的散射系数向量中索引集合对应 的元素组成的向量，为索引集合在双程距离单元集合Ω上的补集，表示对括号中 矩阵求伪逆矩阵，右上角f为矩阵伪逆运算符号，?sy为测量矩阵Θ中由索引集合对应 的列组成的矩阵，^,为步骤S0020获得的脉冲压缩处理后的第j列向量。
[0048] S0080、计算重构残余误差
[0049] 即采用公式〃$ -物^计算得到自适应0ΜΡ算法第k次迭代过程中第j个接 收机的重构残余误差，记为，其中为步骤S0020获得的脉冲压缩处理后的第j列向 量。
[0050] S0090、计算平均重构残余误差和平均散射系数向量
[0051] S0091、采用公?
开述步骤S0080获得的S个接收机的重 构残余误差，j = l,…，S求平均，计算得到平均重构残余误差向量，记为mean_res(k); [0052] 30092、采用公式8〇1'1:_1」=8〇1'1:(3&8(1」））对步骤30070获得的第」个接收机在双程 距离单元集合Ω上的重构得到的散射系数向量中元素的模值按照升序从大到小排列，记为 sort_xj，其中sort (·)表示将括号内向量中的元素根据其值的大小按照升序从大到小排 列，abs( ·)表示对括号内向量中的元素求模；
[0053] S0093、采用公￥
?所述步骤S0092获得的S个接收机排序 后的散射系数向量sort_Xj，j = 1，…，S求平均，计算得到S个接收机在双程距离单元集合Ω 上的排序后散射系数的平均散射系数向量，记为mean_x。
[0054] S0100、计算平均重构残余误差变化率和平均散射系数变化率
[0055] 即采用公另
-计算得到第k次迭代过程中的平 一 2. 均重构残余误差变化率，记为nk，其中mean_res(k)为步骤S0091计算得到的第k次迭代过程 中S个接收机的平均重构残余误差向量，U · | |2表示对向量求L2范数，| · |表示求绝对值， 若k=l，mean_res(k-1〉的值为
，其中为步骤S0040中得到的初 始迭代值;采用公式 计算得到第k次迭代过程中的平均散 射系数变化率，记为妒，其中mean_x (k)表示向量mean_x中的第k个元素，若k = 1，mean_x (k-1)的值为mean_x(0) =mean_x(l) 〇
[0056] S0110、自适应迭代终止判定
[0057] 如果nk彡1且k<K，贝ljk = k+l，返回步骤S0050;否则自适应0ΜΡ算法终止 迭代，输出重构的散射系数和目标个数，且此刻自适应0ΜΡ算法第k次迭代得到的散射系数 向量&，j = 1，…，S即为S个接收机在双程距离单元集合Ω上最终重构获得的散射系数向 量，第k次迭代次数减去1即为自适应0ΜΡ算法对探测区域目标个数的估计值，即L = k-1，其 中L表示探测区域目标个数的估计值，％为步骤S0040中初始化得到的平均重构残余误差变 化率门限，为步骤S0100计算得到的第k次迭代过程中的平均散射系数变化率，β 1-1为步骤 S0100计算得到的第k-Ι次迭代过程中的平均散射系数变化率，若k=l，i^的值为βΗ = β〇， β〇为步骤S0040中初始化得到的平均散射系数变化率初始门限；将重构获得的散射系数向 量xj，j = 1，…，S存储在二维数据矩阵ECH0= {ECH0( i，j)}中，其中 i = l，2，···，Nrang; j = 1， 2，.-.，S〇
[0058] 下面结合具体实施例，对本发明做进一步详述：
[0059] 实施例1
[0000 ] S1010、初始化系统参数，建立二维回波数据矩阵
[0061 ]准备发射机和接收机，初始化系统参数:系统由一个发射机和S个接收机组成，其 中S = 20;接收机分散布置在地面，以发射机位置为原点，在目标、发射机和接收机所在的地 理空间建立X-Y-Z笛卡尔直角坐标系，发射机坐标位置记为Τ=[0,0,0] τ，接收机的坐标位 置分别记为1'」=|^」，7」，2」]1'，」=1，2，~，20，其中，接收机1、2、3分别位于1'1=[-25,0,0] Tkm，Γ2 = [ 25，0，0 ] Tkm和Γ3 = [ 0，43，0 ] Tkm处，接收机4~20随机散布于由接收机1、2、和3决定 的三角形区域内；系统在X方向、Y方向和Z方向的最远探测距离分别为XX max = 50.5km、YYmax =50 · 5km和ZZmax= 10 · 5km，最近探测距离分别为XXmin = 49 · 5km、YYmin = 49 · 5km和ZZmin = 9.5km;发射机对探测区域宽波束泛光发射线性调频信号，发射脉冲载频为fQ= 1GHz，带宽B =15MHz，脉冲宽度!^ = lys，距离分辨率P = 10m;发射机及各个接收机在时间上保持同步， 以发射机发射线性调频信号时刻为系统时间零点；接收机宽波束泛光接收回波信号，接收 机采样频率为f s = 30MHz，采样点数为Nrang=10000,将S = 20个接收机接收的回波信号存 储在二维数据矩阵echo = {echo(i，j)}中，其中，i = l，2,…，10000表示双程距离单元，j = 1，2，…，20表示接收机，假设目标坐标为Pt，则该目标的回波信号将被存储在双程距离单元 1 =。6：[1(；1^/10)，」_ = 1，2，"_，20中，其中，；1^=||1'-?1；||2+||?1；-：1^||2为发射机到目标再到第 j个接收机的双程距离，ceil表示向上取整，I I I |2表示L2范数。
[0062] S1020、回波数据进行脉冲压缩处理
[0063]对每个接收机的回波数据进行脉冲压缩处理:取出步骤S1010数据矩阵echo,利用 脉冲压缩方法对echo的每一列数据echo」，j = 1，2，…，20进行脉冲压缩处理，得到脉冲压缩 处理后的列向量echOj，j = 1，2，…，20，将所有接收机脉冲压缩处理后的列向量echo./ :，j = 1， 2，···，20存储在二维数据矩阵 ECHO = {ECHO(/../)}中，其中i = l，2,…，10000;j = l，2,…， 20〇
[0064] S1030、构造回波信号测量矩阵
[0065] 回波信号测量矩阵0={0(111，(；）}，其中111=1，2，~，1〇〇〇〇;。= 1，2，."，1〇〇〇〇，0 〇11，(3)具体表达式为：
[0066]
[0067] S1040、设定自适应正交匹配追踪(简称0ΜΡ)算法的初始参数
[0068]初始化自适应0ΜΡ算法的参数包括：自适应0ΜΡ算法重构迭代处理的最大迭代次 数，记作κ = 10;平均重构残余误差变化率门限，记作％ = 〇. 03;平均散射系数变化率初始门 限，记作= 5;令Ω={1，···，10000}表示双程距离单元集合;第j，j = 1，…，20个接收机在 双程距离单元集合Ω上的目标散射系数向量的初始值，记为，*；； e Ciw(wxl为10000行1列的 复数向量，一般xj的初始值选择为xj = 〇;第j，j = 1，…，20个接收机的重构残余误差初始迭 代值记为e 为 1 0 0 0 0行1列的复数向量，一般的值选择为 ，:其中ecii&为步骤S1020获得的脉冲压缩处理后的第j列向量;第j，j = 1，…，20 个接收机索引集合的初始迭代值，记为Ξ丨;;，Ξ;:;的值选择为Ξ;;;=0，其中0表示空集;k表示 自适应0ΜΡ算法中的第k次迭代，k的初始值设置为k= 1。
[0069] S1050、寻找回波信号测量矩阵与重构残余误差的最大相干项
[0070] 采用公式
|计算得到自适应0ΜΡ算法第k次迭代过程中 测量矩阵与第j个接收机重构残余误差的最大相干项，其中at_gmaX(_)表示提取括号内向量 t'n 中的最大值元素对应的行序号m值，abs( ·)表示对括号内向量中的元素求模，Θ为步骤 S1030中得到的回波信号测量矩阵，右上角T为转置运算符号，为算法第k-Ι次迭代过 程中得到的重构残余误差，k表示自适应0ΜΡ算法中的第k次迭代次数;若k= 1，的值为 步骤S1040中得到的初始迭代值〃s((j，否则通过自适应0ΜΡ算法第k-1次迭代过程中 步骤S1080的重构残余误差计算提供。
[0071] S1060、更新索引集合
[0072] 采用公式与;卜<)"|」/;计算得到自适应〇1?算法第1^欠迭代过程中的第」个接收 机索引集合，记为其中为自适应0ΜΡ算法第k-Ι次迭代过程中得到的第j个接收机 索引集合，#为步骤S1050计算得到的测量矩阵与第j个接收机重构残余误差的最大相干 项，U表示并集运算符合，k表示自适应0ΜΡ算法中的第k次迭代次数;若k= 1，的值为步 骤S1040中得到的初始值：^;k = 2时，的值为自适应0ΜΡ算法第1次迭代过程中得到的 第j个接收机索引集合Sg;k = 3时，Ξ&11的值为自适应OMP算法第2次迭代过程中得到的第j 个接收机索引集合Ξ;;;;依次类推，在第k次迭代中且k>l时，Ξ;;-丨11的值为自适应0MP算法第 k-Ι次迭代过程中得到的第j个接收机索引集合。
[0073] S1070、重构第j，j = 1，…，S个接收机在双程距离单元集合Ω上的目标散射系数向 量，具体为：
[0074]采用公式
：十算得到自 适应0ΜΡ算法第k次迭代过程中第j个接收机在双程距离单元集合Ω上的散射系数向量，记 为^，其中为第j个接收机在双程距离单元集合Ω上的散射系数向量中由索引集合 对应的元素组成的向量，为步骤S1060中第k次迭代过程中得到的第j个接收机索引 集合，％丨|为第j个接收机在双程距离单元集合Ω上的散射系数向量中索引集合对应 的元素组成的向量，为索引集合在双程距离单元集合Ω上的补集，表示对括号中 矩阵求伪逆矩阵，右上角?为矩阵伪逆运算符号，为测量矩阵Θ中由索引集合Ξ|对应 的列组成的矩阵，eciio,·为步骤S1020获得的脉冲压缩处理后的第j列向量。
[0075] S1080、计算重构残余误差
[0076] 采用公式〃计算得到自适应0ΜΡ算法第k次迭代过程中第j个接收 机的重构残余误差，记为〃<|，其中^,为步骤S1020获得的脉冲压缩处理后的第j列向 量。
[0077] S1090、计算平均重构残余误差和平均散射系数向量，具体包括：
[0078] S1091、采用公￥
#斤述步骤S1080获得的S个接收机的重 构残余误差^4^.，j = l，…，S求平均，计算得到平均重构残余误差向量，记为mean_res(k); [0079] 31092、采用公式8〇1'1:_1」=8〇1'1:(3&8(1」））对步骤31070获得的第」个接收机在双程 距离单元集合Ω上的重构得到的散射系数向量中元素的模值按照升序从大到小排列，记为 sort_xj，其中sort (·)表示将括号内向量中的元素根据其值的大小按照升序从大到小排 列，abs( ·)表示对括号内向量中的元素求模；
[0080] S1093、采用公另
对所述步骤S1092获得的S个接收机排序 后的散射系数向量sort_Xj，j = 1，…，S求平均，计算得到S个接收机在双程距离单元集合Ω 上的排序后散射系数的平均散射系数向量，记为mean_x。
[0081 ] S1100、计算平均重构残余误差变化率和平均散射系数变化率
[0082] 即采用公；E
卜算得到第k次迭代过程中的平 均重构残余误差变化率，记为nk，其中mean_resw为步骤S1091计算得到的第k次迭代过程 中S个接收机的平均重构残余误差向量，U · | |2表示对向量求L2范数，| · |表示求绝对值， 若k=l mearurestn的值〉
车中叹^为步骤S1040中得到的初 始迭代值;采用公?
卜算得到第k次迭代过程中的平均散 射系数变化率，记为妒，其中mean_x (k)表示向量mean_x中的第k个元素，若k = 1，mean_x (k-1)的值为mean_x(0) =mean_x(l) 〇
[0083] S1110、自适应迭代终止判定
[0084] 如果nk彡1且k<K，贝ljk = k+l，返回步骤S1050;否则自适应0ΜΡ算法终止 迭代，输出重构的散射系数和目标个数，且此刻自适应0ΜΡ算法第k次迭代得到的散射系数 向量&，j = 1，…，S即为S个接收机在双程距离单元集合Ω上最终重构获得的散射系数向 量，第k次迭代次数减去1即为自适应0ΜΡ算法对探测区域目标个数的估计值，即L = k-1，其 中L表示自适应0ΜΡ算法对探测区域目标个数的估计值，％为步骤S1040中初始化得到的平 均重构残余误差变化率门限，栌为步骤S1110计算得到的第k次迭代过程中的平均散射系数 变化率，β 1-1为步骤S1110计算得到的第k-Ι次迭代过程中的平均散射系数变化率，若k= 1， 俨1的值为俨1 = β〇，β〇为步骤S1040中初始化得到的平均散射系数变化率初始门限;将重构 获得的散射系数向量xj，j = 1，…，S存储在二维数据矩阵ECH0= {ECH0( i，j)}中，其中i = 1， 2，…，Nrang; j = 1，2，…，S 〇
[0085] 表1为本发明多基地雷达多目标成像定位方法估计的探测区域目标个数和真实目 标个数对照表；（其中，估计的探测区域目标个数是对50次蒙特卡洛仿真获得的估计目标个 数值求平均获得。）
[0086]
[0087]本发明实施例1主要采用仿真实验的方法进行验证，所有步骤、结论都在 MATLABR2012a上验证正确。
[0088] 本发明利用多基地雷达多目标场景中多个接收机的平均重构残余误差变化率和 平均散射系数变化率代替探测区域中目标个数作为自适应0ΜΡ算法迭代终止的判定条件， 利用自适应0ΜΡ算法迭代终止次数估计探测区域目标个数，适用于实际情况中对目标个数 未知的多基地雷达多目标成像定位。
[0089] 以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精 神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种多基地雷达多目标成像定位方法，其特征在于，具体步骤包括： (1) 初始化系统参数，建立二维回波数据矩阵； (2) 脉冲压缩处理； (3) 构造回波信号测量矩阵； (4) 设定自适应OMP算法的初始参数； (5) 寻找回波信号测量矩阵与重构残余误差的最大相干项； (6) 更新索引集合； (7) 重构双程距离单元集合上的目标散射系数向量； (8) 计算重构残余误差； (9) 计算平均重构残余误差和平均散射系数向量； (10) 计算平均重构残余误差变化率和平均散射系数变化率； (11) 自适应迭代终止判定。2. 如权利要求1所述的多基地雷达多目标成像定位方法，其特征在于，所述步骤（1)中 初始化参数均为已知，初始化参数具体如下：由一个发射机和S个接收机组成系统，其中 3,S个接收机分散布置在地面，W发射机位置为原点，在目标、发射机和接收机所在的地理 空间建立X-Y-Z笛卡尔直角坐标系，发射机坐标位置记为Τ=[0,0,0]τ，接收机的坐标位置 分别为rj = [xj，yj，zj]τ，j = l，2，…，S，[]τ表示矩阵转置，系统在X方向、Y方向和Z方向的最 远探测距离分别为XXmax、YYmax和ZZmax，最近探测距离分别为XXmin、YYmin和ZZmin;发射机对探 测区域宽波束泛光发射线性调频信号，发射脉冲载频为fo,带宽B，脉冲宽度Τρ，距离分辨率 Ρ;发射机及各个接收机在时间上保持同步，W发射机发射线性调频信号时刻为系统时间零 点；接收机宽波束泛光接收回波信号，接收机采样频率为fs，采样点数为Nrang，将S个接收 机接收的回波信号存储在二维数据矩阵echo = {echoQ , j)}中，其中，i = l ,2,…，化ang表 示双程距离单元，j = 1，2，…，S表示接收机，假设目标坐标为Pt，则该目标的回波信号将被 存储在双程距离单元i = ceil (r j/p)，j = 1，2，…，S中，其中，η = II T-Pt II2+II P广η I 12为发 射机到目标再到第j个接收机的双程距离，ceil表示向上取整，I I I h表示L2范数； 所述步骤(2)脉冲压缩处理的具体步骤为:取出所述步骤（1)的数据矩阵echo,利用脉 冲压缩方法对echo的每一列数据echoj，j = 1，2，…，S进行脉冲压缩处理，得到脉冲压缩处 理后的列向量ediOjj = 1，2，···，S，将所有接收机脉冲压缩处理后的列向量echoy .，j = 1， 2，...，5存储在二维数据矩阵玩帘3 = {瓦丽(/。4中，其中1 = 1，2，...，化曰雌^ = 1，2，...，5; 所述步骤(3)构造的回波信号测量矩阵为Θ = { Θ (m，c)}，其中m=l，2,…，Nrang;c = 1，2，...，化日雌，0(111，(3)具体表达式为：3. 如权利要求2所述的多基地雷达多目标成像定位方法，其特征在于:所述步骤(4)的 具体步骤为:设定自适应0MP算法重构迭代处理的最大迭代次数，记作K;设定平均重构残余 误差变化率口限，记作η〇;设定平均散射系数变化率初始口限，记作阶;令Ω = {1，…，化ang} 表示双程距离单元集合;设定第j，j = 1，…，S个接收机在双程距离单元集合Ω上的目标散 射系数向量的初始值，记为为化ange行1列的复数向量，xj的初始值选择为xj =0;设定第j，j = 1，…，S个接收机的重构残余误差初始迭代值记为r份,(;,V份说erwnmgxi为 化ange行1列的复数向量，的值选择为1*巧<^0|=品品,:，其中ec品/为所述步骤(2)获得的脉 冲压缩处理后的第j列向量；设定第j，j = l，…，S个接收机索引集合的初始迭代值，记为 邱！，邱；的值选择为己為=0，其中0表示空集;k表示自适应OMP算法中的第k次迭代，k的初 始值设置为k=l。4. 如权利要求3所述的多基地雷达多目标成像定位方法，其特征在于，所述步骤(5)的 具体步骤为:采用公式计算得到自适应OMP算法第k次迭代过 程中测量矩阵与第j个接收机重构残余误差的最大相干项，其中表示提取括号内 向量中的最大值元素对应的行序号m值，abs( ·)表示对括号内向量中的元素求模，Θ为所 述步骤(3)中得到的回波信号测量矩阵，右上角T为转置运算符号，算法第k-1次迭 代过程中得到的重构残余误差，k表示自适应0MP算法中的第k次迭代次数;若k = 1， 的值为所述步骤(4)中得到的初始迭代值re哉j，否则喃过自适应OMP算法第k-1次迭 代过程中所述步骤(8)的重构残余误差计算提供。5. 如权利要求4所述的多基地雷达多目标成像定位方法，其特征在于，所述步骤(6)的 具体步骤为:采用公式己瓦"U/;计算得到自适应0MP算法第k次迭代过程中的第j个接 收机索引集合，记为瑞-|，其中为自适应0MP算法第k-1次迭代过程中得到的第j个接收 机索引集合，/?为所述步骤(5)计算得到的测量矩阵与第j个接收机重构残余误差的最大相 干项，U表示并集运算符合，k表示自适应0MP算法中的第k次迭代次数;若k= 1，SjS户的值为 所述步骤(4)中得到的初始值S;^;;k = 2时，Ξ打)的值为自适应0MP算法第1次迭代过程中得 到的第j个接收机索引集合Ξ打;k = 3时，Ξ祈U的值为自适应0MP算法第2次迭代过程中得到 的第j个接收机索引集合百拐;依次类推，在第k次迭代中且k> 1时，旦汾，的值为自适应0MP算 法第k-1次迭代过程中得到的第j个接收机索引集合。6. 如权利要求5所述的多基地雷达多目标成像定位方法，其特征在于:所述步骤(7)的 具体步骤为:重构第j，j = 1，…，S个接收机在双程距离单元集合Ω上的目标散射系数向量， 具体为:采用公??开算得到自 适应0ΜΡ算法第k次迭代过程中第j个接收机在双程距离单元集合Ω上的散射系数向量，记 为xj，其中X啤,1为第j个接收机在双程距离单元集合Ω上的散射系数向量xj中由索引集合 畔;財应的元素组成的向量，Ξ;;;为所述步骤(6)中第k次迭代过程中得到的第j个接收机索 弓煤合，为第j个接收机在双程距离单元集合Ω上的散射系数向量X冲索引集合畐;对 应的元素组成的向量，宜货为索引集合邱洽双程距离单元集合Ω上的补集，(丫表示对括号 中矩阵求伪逆矩阵，右上角?为矩阵伪逆运算符号戸瑕1为测量矩阵Θ中由索引集合邱;对 应的列组成的矩阵，品品，为所述步骤(2)获得的脉冲压缩处理后的第j列向量； 所述步骤(8)的具体步骤为:采用公式^心',';;=6〇^1^;-(^,.计算得到自适应01?算法第1^ 次迭代过程中第j个接收机的重构残余误差，记为其中品而,·为步骤(2)获得的脉冲压 缩处理后的第j列向量。7. 如权利要求6所述的多基地雷达多目标成像定位方法，其特征在于，所述步骤(9)的 具体步骤包括： (9.1) 采用公式对所述步骤(8)获得的S个接收机的重构残余 误差，j = 1，…，S求平均，计算得到平均重构残余误差向量，记为mean_res W . (9.2) 采用公式soパ_Xj = soパ(abs(Xj))对所述步骤(7)获得的第j个接收机在双程距 离单元集合Ω上的重构得到的散射系数向量中元素的模值按照升序从大到小排列，记为 S〇d_Xj，其中sod( ·)表示将括号内向量中的元素根据其值的大小按照升序从大到小排 列，abs( ·)表示对括号内向量中的元素求模； (9.3) 采用公；1^对所述步骤(9.2)获得的S个接收机排序后的散 射系数向量sod_Xj，j = 1，…，S求平均，计算得到S个接收机在双程距离单元集合Ω上的排 序后散射系数的平均散射系数向量，记为mean_x。8. 如权利要求7所述的多基地雷达多目标成像定位方法，其特征在于，所述步骤(10)的 具体步骤为： 采用公式1^算得到第4次迭代过程中的平均重构 残余误差变化率，记为扣其中mean_resW为所述步骤(9.1)计算得到的第k次迭代过程中S 个接收机的平均重构残余误差向量，I I · I U表示对向量求L2范数，I · I表示求绝对值，若k = l，mean_res(W的值夫其中reSfl;;为步骤(4)中得到的初始迭代 值;采用公式计算得到第k次迭代过程中的平均散射系数 变化率，记为护，其中mean_x化)表示向量mean_x中的第k个元素，若k = 1，mean_x化-1)的值 为mean_x(0)=mean_x(l)D9.如权利要求8所述的多基地雷达多目标成像定位方法，其特征在于:所述步骤(11)的 具体步骤为:如果，护且k<K，则k=k+i，返回所述步骤(5);否则自适应OMP算法 终止迭代，输出重构的散射系数和目标个数，且此刻自适应0MP算法第k次迭代得到的散射 系数向量xj，j = 1，…，S即为S个接收机在双程距离单元集合Ω上最终重构获得的散射系数 向量，第k次迭代次数减去1即为自适应0MP算法对探测区域目标个数的估计值，目化= k-l， 其中L表示探测区域目标个数的估计值，no为所述步骤(4)中初始化得到的平均重构残余误 差变化率口限，0k为所述步骤（10)计算得到的第k次迭代过程中的平均散射系数变化率， βΚ为所述步骤(10)计算得到的第k-1次迭代过程中的平均散射系数变化率，若k=l，i]W的 值为护-1=阶，抗为所述步骤(4)中初始化得到的平均散射系数变化率初始口限;将重构获得 的散射系数向量xj，j = 1，…，S存储在二维数据矩阵ECH0=巧畑0( i，j)}中，其中i = 1， 2，...，Nrang; j = l，2，...，S。
【文档编号】G01S13/89GK105974386SQ201610291839
【公开日】2016年9月28日
【申请日】2016年5月5日
【发明人】樊玲
【申请人】乐山师范学院