联合无线电信号复包络和载波相位信息的超视距目标直接定位方法
【技术领域】
[0001]本发明涉及无线电信号定位领域，特别涉及一种联合无线电信号复包络和载波相 位信息的超视距目标直接定位方法。
【背景技术】
[0002] 众所周知，无线电信号定位技术对于目标发现及其态势感知具有十分重要的意 义，经过近几十年的发展，该技术在理论和工程应用中都取得了长足的进展。根据观测站的 数目进行划分，可以将无线电信号定位体制划分为单站定位和多站定位两大类，前者具有 系统简洁、灵活性高、无需信息同步和通信传输等优点，后者则能够提供更多的观测量，从 而有助于取得更高的定位精度。本专利主要涉及到多站定位方式。在多站定位系统中，最重 要的两类定位体制是多站测向交汇定位和多站测时差交汇定位。前者要求每个观测站安装 天线阵列，各个观测站利用信号到达本站内不同天线的载波相位差估计信号方位，然后再 在中心站进行交汇定位;第二种定位体制则要求估计信号复包络到达不同观测站的时延 差，并利用时延差进行交汇定位。从所利用的信息来看，前者仅仅利用了到达信号的载波相 位信息，而后者仅仅利用了到达信号的复包络信息。虽然这两类定位体制都有其自身优势， 但是定位精度都存在较大的提升空间，为了大幅度提尚多站定位精度，最好是能够同时利 用到达信号的复包络和载波相位信息，并设计出合理的定位方法。
[0003] 对于远距离目标而言，目标信号往往是通过超视距传播的方式到达各个观测站， 最常见的一种传播方式是信号经过电离层折射后入射至地面观测站，在这种定位场景下人 们通常采用多站测向交汇法进行目标定位，单在实际应用中该方法的定位精度并不高。事 实上，若想对超视距目标进行精确定位，要求对信号的传播路径建立有效的数学模型，并且 需要一些先验知识(例如电离层高度的先验观测），除此之外还应尽可能地同时利用信号复 包络和载波相位信息。因此，如何在获得电离层高度先验观测的基础上，通过联合信号复包 络和载波相位信息实现对超视距目标的精确定位是十分有意义的课题，这也是本专利重点 要解决的问题。另一方面，现有的无线电信号定位过程都可以归纳为两步估计定位模式，BP 首先从信号数据中提取定位参数(例如方位、时延差、多普勒等），然后再基于这些参数解算 出目标位置。虽然这种两步定位模式在现代定位系统中已被广泛应用，但以色列学者 A. J. WeiSS和A. Amar却指出了其中所存在的诸多缺点，并提出了单步直接定位的思想，其基 本理念是从信号采集数据域中直接估计目标的位置参数，而无需再估计其它中间定位参 数。
【发明内容】
[0004] 针对现有技术中的不足，本发明提供一种联合无线电信号复包络和载波相位信息 的超视距目标直接定位方法，借鉴单步直接定位模式解决超视距目标的定位问题，定位精 度要比传统的两步定位模式高，这也是提高超视距目标定位精度的一种重要手段，相比于 传统的针对超视距目标的多站测向交汇定位法，能够有效提高目标定位精度，并且可以在 电离层高度先验观测的基础上进一步提高对电离层高度的估计精度，具有较快的收敛速 度，无需高维搜索，是一种性能可靠、运算高效的定位方法。
[0005] 按照本发明所提供的设计方案，一种联合无线电信号复包络和载波相位信息的超 视距目标直接定位方法，具体包含如下步骤：
[0006] 步骤1.对N个观测站的Μ通道阵列天线接收系统进行时间同步，根据奈奎斯特采样 定理采集目标辐射的无线电信号数据，获得阵列信号时域数据；
[0007] 步骤2.每个观测站将所采集到的Q个时域数据样本点做基2-FFT运算，得到阵列信 号频域数据，其中，Q为2的整数次幂；
[0008] 步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站，中心站利用电离 层高度的先验观测以及观测站的阵列信号频域数据建立最大似然参数估计准则；
[0009] 步骤4.在最大似然参数估计准则的基础上，通过数学推演建立联合估计目标位置 参数和电离层高度的数学优化模型；
[0010]步骤5.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法，利用数学优 化模型中矩阵的最大特征值进行数值寻优，进行超视距目标的定位。
[0011]上述的，步骤1中，第η个观测站的阵列天线所接收到的信号时域模型为Xn(t)=i3 nan(p)s(t-Tn(p，hn)-t())+e n(t)(l < η < N)，其中，p表示目标位置向量，hn表示目标信号经过 超视距传播至第η个观测站所经历的电离层高度，to表示目标发射信号时间，s(t)表示目标 信号复包络，a n(pd)表示目标信号相对于第η个天线阵列的阵列流形向量，βη表示目标信号 传播至第η个观测站的损耗因子，e n(t)表示第η个观测站中天线阵列的阵元噪声向量，τη(ρ， hn)表示目标信号到达第η个观测站的传播时延，它同时是关于目标位置向量ρ和电离层高 度1的函数；步骤2中，第η个观测站的阵列天线所接收到的信号频域模型为
其中，？和分别表 示s(tWPen(t)的频域形式，coq表示第q个数字频点，bn(p，h n，coq)和F(q)的表达式分别为
步骤3中，中心站所建立的最大似然估计准则为
其中，向量 h=[hi h2…hN]T/&含了每个观测站所对应的电尚层尚度，石=[/7|/22.../7、] |表不11的先验 观测向量，其观测误差协方差矩阵记为P;步骤4中，所建立的联合估计目标位置向量P和电 离层高度向量h的数学优化模型，
其中，λΜΧ[ ·]表示取矩阵的最大特征值，Ρ = ，矩阵B(p，h)和文的表达式分别为
L m'」 L 向量（(pA.)的表达式为
[0012] 上述的，步骤5中，所设计出的数值优化中的牛顿型迭代算法的实现步骤为：
[0013] 步骤5.1)利用多重信号分类估计算法和泰勒级数迭代定位算法获得目标位置向 量的初始估计，将电离层高度的先验观测向量A作为h的初始估计左~，形成初始迭代向 i
c.
[0014] 步骤5.2)进行牛顿型迭代的计算公式为
式中k表示迭代次数，0<μ<1表示迭代步长因子，和分别表示目标函数的梯 度向量和H e s s i a η矩阵，相应的计算公式分别为
其中，《?以? 表示厄米特矩阵......最大特征值，^>
所对应的单位特征向量，矩阵
I 勺计算公式分别为 7
[0020]
[0021]
分别对应矩阵所除最大特征值以外的其余N-1个特征值及其对 应的单位特征向量。
[0029]本发明的有益效果：
[0030]本发明首先基于信号超视距传播模型，建立到达信号复包络和载波相位关于目标 位置参数的数学模型，获得多站阵列信号模型，接着利用基2-FFT算法将多站阵列天线接收 数据转化为频域数据，并基于电离层高度的先验观测以及最大似然估计准则建立联合估计 目标位置参数和电离层高度的数学优化模型，最后根据矩阵特征值扰动理论设计出Newton 型迭代算法，以实现对超视距目标的精确定位，相比于传统的针对超视距目标的多站测向 交汇定位法，能够有效提高目标定位的精度，并且可以在电离层高度先验观测的基础上，进 一步提高对电离层高度的估计精度;此外，本发明通过Newton型迭代来实现超视距定位，具 有较快的收敛速度，无需高维搜索，是一种性能可靠、运算高效的定位方法。
【附图说明】：
[0031 ]图1为本发明的多观测站超视距目标直接定位原理图；
[0032] 图2为本发明的流程示意图；
[0033] 图3为本发明的超视距目标定位实例场景示意图；
[0034] 图4为本发明的定位结果比对示意图。
【具体实施方式】：
[0035] 下面结合附图和技术方案对本发明作进一步详细的说明，并通过优选的实施例详 细说明本发明的实施方式，但本发明的实施方式并不限于此。
[0036]实施例一，参见图1~2所示，一种联合无线电信号复包络和载波相位信息的超视 距目标直接定位方法，具体包含如下步骤：
[0037]步骤1.对N个观测站的Μ通道阵列天线接收系统进行时间同步，根据奈奎斯特采样 定理采集目标辐射的无线电信号数据，获得阵列信号时域数据；
[0038] 步骤2.每个观测站将所采集到的Q个时域数据样本点做基2-FFT运算，得到阵列信 号频域数据，其中，Q为2的整数次幂；
[0039] 步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站，中心站利用电离 层高度的先验观测以及观测站的阵列信号频域数据建立最大似然参数估计准则；
[0040] 步骤4.在最大似然参数估计准则的基础上，通过数学推演建立联合估计目标位置 参数和电离层高度的数学优化模型；
[0041 ]步骤5.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法，利用数学优 化模型中矩阵的最大特征值进行数值寻优，进行超视距目标的定位。
[0042]实施例二，参见图1~2所示，一种联合无线电信号复包络和载波相位信息的超视 距目标直接定位方法，具体包含如下步骤：
[0043]步骤1.对Ν个观测站的Μ通道阵列天线接收系统进行时间同步，根据奈奎斯特采样 定理采集目标辐射的无线电信号数据，获得阵列信号时域数据，第η个观测站的阵列天线所 接收到的信号时域模型为xn(t) =0nan(p)s(t-Tn(p，hn)-tQ) + en(t)(l < η < Ν)，其中，ρ表示 目标位置向量，hn表示目标信号经过超视距传播至第η个观测站所经历的电离层高度，to表 示目标发射信号时间，s(t)表示目标信号复包络，a n(pd)表示目标信号相对于第η个天线阵 列的阵列流形向量，仏表示目标信号传播至第η个观测站的损耗因子，e n(t)表示第η个观测 站中天线阵列的阵元噪声向量，in(p，hn)表示目标信号到达第η个观测站的传播时延，它同 时是关于目标位置向量Ρ和电离层高度h n的函数；
[0044] 步骤2.每个观测站将所采集到的Q个时域数据样本点做基2-FFT运算，得到阵列信 号频域数据，其中，Q为2的整数次幂，第η个观测站的阵列天线所接收到的信号频域模型为
f 示s(tWPen(t)的频域形式，coq表示第q个数字频点，bn(p，h n，coq)和？(％)的表达式分别为 其中，5Υ?)和茗(6>)分别表
[0045] 步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站，中心站利用电离 层高度的先验观测以及观测站的阵列信号频域数据建立最大似然参数估计准则，中心站所
建立的最大似然估计准则戈 y 其中，向量1ι=[1η h 2··· hN]T包含了每个观测站所对应的电离层高度，1=?毛…4〇τ表 示h的先验观测向量，其观测误差协方差矩阵记为Ρ，在实际工程应用中，可以利用探测技术 对电离层高度进行预测，或利用公开渠道进行查询，作为电离层高度的先验知识加以利用；
[0046] 步骤4.在最大似然参数估计准则的基础上，通过数学推演建立联合估计目标位置 参数和电离层高度的数学优化模型，所建立的联合估计目标位置向量P和电离层高度向量h 的数学优化模型为
其中， Amax[ ·]表示取矩阵的最大特征值，，矩阵8(1)，1〇和#的表达式分别为
[0047] 步骤5.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法，利用数学优 化模型中矩阵的最大特征值进行数值寻优，所设计出的数值优化中的