联合无线电信号复包络和载波相位信息的超视距目标直接定位方法技术

本发明专利技术涉及一种联合无线电信号复包络和载波相位信息的超视距目标直接定位方法，首先基于信号超视距传播模型，建立到达信号复包络和载波相位关于目标位置参数的数学模型，获得多站阵列信号模型，接着利用基2‑FFT算法将多站阵列天线接收数据转化为频域数据，并基于电离层高度的先验观测以及最大似然估计准则建立联合估计目标位置参数和电离层高度的数学优化模型，最后根据矩阵特征值扰动理论设计出Newton型迭代算法进行超视距目标定位。本发明专利技术能够有效提高目标定位的精度，可以在电离层高度先验观测的基础上，进一步提高对电离层高度的估计精度，具有较快的收敛速度，无需高维搜索，性能可靠、运算高效。

【技术实现步骤摘要】

本专利技术涉及无线电信号定位领域，特别涉及一种联合无线电信号复包络和载波相位信息的超视距目标直接定位方法。
技术介绍
众所周知，无线电信号定位技术对于目标发现及其态势感知具有十分重要的意义，经过近几十年的发展，该技术在理论和工程应用中都取得了长足的进展。根据观测站的数目进行划分，可以将无线电信号定位体制划分为单站定位和多站定位两大类，前者具有系统简洁、灵活性高、无需信息同步和通信传输等优点，后者则能够提供更多的观测量，从而有助于取得更高的定位精度。本专利主要涉及到多站定位方式。在多站定位系统中，最重要的两类定位体制是多站测向交汇定位和多站测时差交汇定位。前者要求每个观测站安装天线阵列，各个观测站利用信号到达本站内不同天线的载波相位差估计信号方位，然后再在中心站进行交汇定位；第二种定位体制则要求估计信号复包络到达不同观测站的时延差，并利用时延差进行交汇定位。从所利用的信息来看，前者仅仅利用了到达信号的载波相位信息，而后者仅仅利用了到达信号的复包络信息。虽然这两类定位体制都有其自身优势，但是定位精度都存在较大的提升空间，为了大幅度提高多站定位精度，最好是能够同时利用到达信号的复包络和载波相位信息，并设计出合理的定位方法。对于远距离目标而言，目标信号往往是通过超视距传播的方式到达各个观测站，最常见的一种传播方式是信号经过电离层折射后入射至地面观测站，在这种定位场景下人们通常采用多站测向交汇法进行目标定位，单在实际应用中该方法的定位精度并不高。事实上，若想对超视距目标进行精确定位，要求对信号的传播路径建立有效的数学模型，并且需要一些先验知识(例如电离层高度的先验观测)，除此之外还应尽可能地同时利用信号复包络和载波相位信息。因此，如何在获得电离层高度先验观测的基础上，通过联合信号复包络和载波相位信息实现对超视距目标的精确定位是十分有意义的课题，这也是本专利重点要解决的问题。另一方面，现有的无线电信号定位过程都可以归纳为两步估计定位模式，即首先从信号数据中提取定位参数(例如方位、时延差、多普勒等)，然后再基于这些参数解算出目标位置。虽然这种两步定位模式在现代定位系统中已被广泛应用，但以色列学者A.J.Weiss和A.Amar却指出了其中所存在的诸多缺点，并提出了单步直接定位的思想，其基本理念是从信号采集数据域中直接估计目标的位置参数，而无需再估计其它中间定位参数。
技术实现思路
针对现有技术中的不足，本专利技术提供一种联合无线电信号复包络和载波相位信息的超视距目标直接定位方法，借鉴单步直接定位模式解决超视距目标的定位问题，定位精度要比传统的两步定位模式高，这也是提高超视距目标定位精度的一种重要手段，相比于传统的针对超视距目标的多站测向交汇定位法，能够有效提高目标定位精度，并且可以在电离层高度先验观测的基础上进一步提高对电离层高度的估计精度，具有较快的收敛速度，无需高维搜索，是一种性能可靠、运算高效的定位方法。按照本专利技术所提供的设计方案，一种联合无线电信号复包络和载波相位信息的超视距目标直接定位方法，具体包含如下步骤：步骤1.对N个观测站的M通道阵列天线接收系统进行时间同步，根据奈奎斯特采样定理采集目标辐射的无线电信号数据，获得阵列信号时域数据；步骤2.每个观测站将所采集到的Q个时域数据样本点做基2-FFT运算，得到阵列信号频域数据，其中，Q为2的整数次幂；步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站，中心站利用电离层高度的先验观测以及观测站的阵列信号频域数据建立最大似然参数估计准则；步骤4.在最大似然参数估计准则的基础上，通过数学推演建立联合估计目标位置参数和电离层高度的数学优化模型；步骤5.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法，利用数学优化模型中矩阵的最大特征值进行数值寻优，进行超视距目标的定位。上述的，步骤1中，第n个观测站的阵列天线所接收到的信号时域模型为xn(t)＝βnan(p)s(t-τn(p,hn)-t0)+εn(t)(1≤n≤N)，其中，p表示目标位置向量，hn表示目标信号经过超视距传播至第n个观测站所经历的电离层高度，t0表示目标发射信号时间，s(t)表示目标信号复包络，an(pd)表示目标信号相对于第n个天线阵列的阵列流形向量，βn表示目标信号传播至第n个观测站的损耗因子，εn(t)表示第n个观测站中天线阵列的阵元噪声向量，τn(p,hn)表示目标信号到达第n个观测站的传播时延，它同时是关于目标位置向量p和电离层高度hn的函数；步骤2中，第n个观测站的阵列天线所接收到的信号频域模型为其中，和分别表示s(t)和εn(t)的频域形式，ωq表示第q个数字频点，bn(p,hn,ωq)和的表达式分别为步骤3中，中心站所建立的最大似然估计准则为其中，向量h＝[h1h2…hN]T包含了每个观测站所对应的电离层高度，表示h的先验观测向量，其观测误差协方差矩阵记为P；步骤4中，所建立的联合估计目标位置向量p和电离层高度向量h的数学优化模型为其中，λmax[·]表示取矩阵的最大特征值，矩阵B(p,h)和的表达式分别为向量的表达式为上述的，步骤5中，所设计出的数值优化中的牛顿型迭代算法的实现步骤为：步骤5.1)利用多重信号分类估计算法和泰勒级数迭代定位算法获得目标位置向量的初始估计将电离层高度的先验观测向量作为h的初始估计形成初始迭代向量步骤5.2)进行牛顿型迭代的计算公式为式中k表示迭代次数，0＜μ＜1表示迭代步长因子，和分别表示目标函数的梯度向量和Hessian矩阵，相应的计算公式分别为其中， < g 1 ( η ^ ( k ) ) > i = u 0 H ( p ^ ( k ) , h ^ ( k ) ) Z · i ( a ) ( p ^ ( k ) , h ^ 本文档来自技高网...

【技术保护点】
一种联合无线电信号复包络和载波相位信息的超视距目标直接定位方法，具体包含如下步骤：步骤1.对N个观测站的M通道阵列天线接收系统进行时间同步，根据奈奎斯特采样定理采集目标辐射的无线电信号数据，获得阵列信号时域数据；步骤2.每个观测站将所采集到的Q个时域数据样本点做基2‑FFT运算，得到阵列信号频域数据，其中，Q为2的整数次幂；步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站，中心站利用电离层高度的先验观测以及观测站的阵列信号频域数据建立最大似然参数估计准则；步骤4.在最大似然参数估计准则的基础上，通过数学推演建立联合估计目标位置参数和电离层高度的数学优化模型；步骤5.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法，利用数学优化模型中矩阵的最大特征值进行数值寻优，进行超视距目标的定位。

【技术特征摘要】
2015.12.05 CN 20151088094321.一种联合无线电信号复包络和载波相位信息的超视距目标直接定位方法，具体包含如下步骤：步骤1.对N个观测站的M通道阵列天线接收系统进行时间同步，根据奈奎斯特采样定理采集目标辐射的无线电信号数据，获得阵列信号时域数据；步骤2.每个观测站将所采集到的Q个时域数据样本点做基2-FFT运算，得到阵列信号频域数据，其中，Q为2的整数次幂；步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站，中心站利用电离层高度的先验观测以及观测站的阵列信号频域数据建立最大似然参数估计准则；步骤4.在最大似然参数估计准则的基础上，通过数学推演建立联合估计目标位置参数和电离层高度的数学优化模型；步骤5.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法，利用数学优化模型中矩阵的最大特征值进行数值寻优，进行超视距目标的定位。2.根据权利要求1所述的联合无线电信号复包络和载波相位信息的超视距目标直接定位方法，其特征在于：步骤1中，第n个观测站的阵列天线所接收到的信号时域模型为xn(t)＝βnan(p)s(t-τn(p,hn)-t0)+εn(t)(1≤n≤N)，其中，p表示目标位置向量，hn表示目标信号经过超视距传播至第n个观测站所经历的电离层高度，t0表示目标发射信号时间，s(t)表示目标信号复包络，an(pd)表示目标信号相对于第n个天线阵列的阵列流形向量，βn表示目标信号传播至第n个观测站的损耗因子，εn(t)表示第n个观测站中天线阵列的阵元噪声向量，τn(p,hn)表示目标信号到达第n个观测站的传播时延，它同时是关于目标位置向量p和电离层高度hn的函数；步骤2中，第n个观测站的阵列天线所接收到的信号频域模型为其中，和分别表示s(t)和εn(t)的频域形式，ωq表示第q个数字频点，bn(p,hn,ωq)和的表达式分别为步骤3中，中心站所建立的最大似然估计准则为，其中，向量h＝[h1h2…hN]T包含了每个观测站所对应的电离层高度，表示h的先验观测向量，其观测误差协方差矩阵记为P；步骤4中，所建立的联合估计目标位置向量p和电离层高度向量h的数学优化模型为其中，λmax[·]表示取矩阵的最大特征值，矩阵B(p,h)和的表达式分别为，向量的表达式为(1≤n≤N)。3.根据权利要求2所述的联合无线电信号复包络和载波相位信息的超视距目标直接定位方法，其特征在于：步骤5中，所设计出的数值优化中的牛顿型迭代算法的实现步骤为：步骤5.1)利用多重信号分类估计算法和泰勒级数迭代定位算法获得目标位置向量的初始估计将电离层高度的先验观测向量作为h的初始估计形成初始迭代向量步骤5.2)进行Newton型迭代的计算公式为式中k表示迭代次数，0＜μ＜1表示迭代步长因子，和分别表示目标函数的梯度向量和Hessian矩阵，相应的计算公式分别为其中， < g 1 ( η ^ ( k ) ) > i = u 0 H ( p ^ ( k ) , h ^ ( k ) ) Z · i ( a ) ( p ^ ( k ) , h ^ ( k ) ) u 0 ( p ^ ( k ) , h ^ ( k ) ) ]]> < g 2 ( η ^ ( k ) ) > i = u 0 H ( p ^ ( k ) , h ^ ( k ) ) Z i ( b ) ( p ^ ( k ) , h ^ ( k ) ) u 0 ( p ^ ( k ) , h ^ ( k ) ) ]]> < G 11 ( η ^ ( k ) ) > i j = u 0 H ( p ^ ( k ) , h ^ ( k ) ) Z ·· i j ( a a ) ( p ^ ( k ) , h ^ ( k ) ) u 0 ( p ^ ( k ) , h ^ ( k ) ) + 2 u 0 H ( p ^ ( k ) , h ^ ( k ) ) Z · i ( a ) H ( p ^ ( k ) , h ^ ( k ) ) U 0 ( p ^ ( k ) , h ^ ( k ) ) Z · j ( a ) ( p ^ ( k ) , h ^ ( k ) ) u 0 ( p ^ ( k ) , h ^ ( k ) ) ]]> < G 12 ( η ^ ( k ) ) > i j = u 0 H ( p ^ ( k ) , h ^ ( k ) ) Z ·· i j ( a b ) ( p ^ ( k ) , h ^ ( k ) ) u 0 ( p ^ ( k ) , h ^ ( k ) ) + 2 u 0 H ( p ^ ( k ) , h ^ ( k ) ) Z · i ( a ) H ( p ^ ( k ) , h ^ ( k ) ) U 0 ( p ^ ( k ) , h ^ ( k ) ) Z · j ( b ) ( p ^ ( k ) , h ^ ( k ) ) u 0 ( p ^ ( k ) , h ^ ( k ) ) ]]>其中，表示厄米特矩阵最大特征值所对应的单位特征向量，矩阵和的计算公式分别为 Z · i ( a ) ( p ^ ( k ) , h ^ ( k ) ) = B · i ( a ) ( p ^ ( k ) , h ^ ( k ) ) X ~ X ~ H B H ( p ^ ( k ) , h ^ ( k ) ) + B ( p ^ ( k ) , h ^ ( k ) ) X ~ X ~ H B · i ( a ) H ( p ^ ( k ) , h ^ ( k ) ) ]]> Z · i ( b ) ( p ^ ( k ) , h ^ ( k ) ) = B · i ( b ) ( p ^ ( k ) , h ^ ( k ) ) X ~ X ~ H B H ( p ^ ( k ) , h ^ ( k ) ) + B ( p ^ ( k ) , h ^ ( k ) ) X ~ X ~ H B · i ( b ) H ( p ^ ( k ) , h ^ ( k ) ) ]]> Z · i j ( a a ) ( p ^ ( k ) , h ^ ( k ) ) = B ·· i j ( a a ) ( p ^ ( k ) , h ^ ( k ) ) X ~ X ~ H B H ( p ^ ( k ) , h ^ ( k ) ) + B · i ( a ) ( p ^ ( k ) , h ^ ( k ) ) X ~ X ~ H B · j ( a ) H ( p ^ ( k ) , h ^ ...

【专利技术属性】
技术研发人员：于宏毅，王鼎，吴瑛，杜剑平，杨宾，张莉，张刚，唐涛，吴江，
申请(专利权)人：中国人民解放军信息工程大学，
类型：发明
国别省市：河南;41