CN105954712B - 联合无线电信号复包络和载波相位信息的多目标直接定位方法 - Google Patents
联合无线电信号复包络和载波相位信息的多目标直接定位方法 Download PDFInfo
- Publication number
- CN105954712B CN105954712B CN201610255410.XA CN201610255410A CN105954712B CN 105954712 B CN105954712 B CN 105954712B CN 201610255410 A CN201610255410 A CN 201610255410A CN 105954712 B CN105954712 B CN 105954712B
- Authority
- CN
- China
- Prior art keywords
- signal
- target
- matrix
- array
- frequency domain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/0009—Transmission of position information to remote stations
- G01S5/0018—Transmission from mobile station to base station
- G01S5/0027—Transmission from mobile station to base station of actual mobile position, i.e. position determined on mobile
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/46—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being of a radio-wave signal type
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0278—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves involving statistical or probabilistic considerations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Probability & Statistics with Applications (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Mobile Radio Communication Systems (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明涉及一种联合无线电信号复包络和载波相位信息的多目标直接定位方法,首先建立到达信号的复包络和载波相位关于目标位置参数的解析模型,通过基2‑FFT算法将多站阵列信号时域数据转化成频域数据,在频域利用子空间正交原理构造直接估计目标位置参数的数学优化模型,基于矩阵特征值扰动理论设计出Newton型迭代算法,进行多目标的精确定位。本发明相比传统的多站测向交汇定位和多站测时差交汇定位方法,能够利用到达信号的复包络和载波相位信息,通过在底层数据域直接进行目标位置估计,能够取得更高的定位精度,避免两步定位模式中存在的门限效应和“目标‑量测”匹配问题,具有较快的收敛速度,无需高维搜索,性能可靠、运算高效。
Description
技术领域
本发明涉及无线电信号的定位领域,特别涉及一种联合无线电信号复包络和载波相位信息的多目标直接定位方法。
背景技术
众所周知,无线电信号定位对于目标发现及其态势感知具有重要意义,其在通信信号侦察、电子信息对抗、无线电监测、遥测与导航等诸多领域具有十分重要的应用,近几十年来一直受到国内外相关学者和工程技术人员的广泛关注和研究。根据观测站的数目进行划分可将无线电信号定位体制划分为单站定位和多站定位两大类,这两类体制各有其自身优势。具体来说,单站定位系统具有灵活性高、机动性好、系统简洁、无需信息同步和信息传输等优点,而多站定位系统则能够提供更多的观测信息量,有助于获得更高的定位精度。本专利主要涉及多站定位体制。在多站定位系统中,最具有代表性的两种定位方式分别为多站测向交汇定位和多站测时差交汇定位。第一种定位系统要求每个观测站安装天线阵列,每个观测站利用信号到达本站内不同天线的载波相位差信息估计出信号方位,然后在中心站进行交汇定位;第二种定位体制则要求估计信号到达不同观测站的复包络时延差,并利用时延差信息进行交汇定位。从所利用的信息来看,前者仅仅利用了信号的载波相位信息,而后者仅仅利用了信号的复包络信息。虽然这两种定位系统都有其自身优势,但定位精度都存在较大提升空间,为了大幅度提高多站定位精度,最好能同时利用到达信号的复包络和载波相位信息,并设计出新型定位方法。
需要指出的是,现有的无线电信号定位方法都可以归纳成两步估计定位模式,即先从信号数据中估计定位参数(例如方位,时延差等),然后再基于这些参数估计目标位置。虽然这种两步估计定位模式应用广泛,但存在诸多缺点。首先,根据信息处理的理论可知,两步估计定位模式难以获得统计最优的定位精度,这是因为从原始信号数据到最终估计结果之间每增加一步信息处理环节,就会引入一些不确定性,从而损失掉部分信息,并且影响最终的定位精度(尤其在低信噪比和小样本数条件下该问题尤为突出)。其次,两步估计定位模式存在门限效应,以测向交汇定位体制为例,当两个目标相对于某个观测站的方位十分邻近时(小于其角度分辨门限时),该观测站会将其误判为同一信号,若该测向站将这一错误信息传送至中心站则可能导致中心站误判。最后,当有多个目标同时存在时,两步估计定位模式存在“目标-量测”匹配问题,即如何将信号测量参数与目标进行正确关联,从而完成后续的多目标定位。为了克服两步定位模式的缺点,一种较好的措施是采用单步定位(亦称直接定位)模式,即从信号采集数据域中直接提取目标的位置参数,而无需估计其它的中间定位参数,这一定位思想最早是由以色列学者A.J.Weiss和A.Amar所提出,然而他们所提出的定位算法都需要高维搜索运算,可能会导致庞大的运算量,需要设计出可以避免多维搜索的高效算法。
发明内容
针对现有技术中多站测向交汇定位和多站测时差交汇定位的不足,本发明提供一种联合无线电信号复包络和载波相位信息的多目标直接定位方法,通过在底层数据域直接进行目标位置估计,取得更高的定位精度,同时避免两步定位模式中存在的门限效应和“目标-量测”匹配问题,具有较快的收敛速度、无需高维搜索,性能稳定、可靠。
按照本发明所提供的设计方案,一种联合无线电信号复包络和载波相位信息的多目标直接定位方法,具体包含如下步骤:
步骤1.对N个观测站的M通道阵列天线接收系统做时间同步,根据奈奎斯特采样定理采集目标辐射的无线电信号数据,获得阵列信号时域数据;
步骤2.对每个观测站的阵列信号时域数据在时域上划分成K个子段,每个子段内均包含有Q个采集数据点,对每个子段内的Q个数据样本做基2-FFT运算,得到阵列信号频域数据;
步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站,中心站对每个观测站传输的阵列信号频域数据按照观测站的顺序堆栈排列,构造高维阵列信号频域数据;
步骤4.针对每个采集数据点,利用K个子段数据构造高维阵列协方差矩阵的最大似然估计值;
步骤5.对每个采集数据点所构造的高维阵列协方差矩阵进行特征值分解,利用其MN-D个最小特征值对应的单位特征向量构造特征向量矩阵,并计算特征向量矩阵列空间的正交投影矩阵,D为目标个数;
步骤6.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法进行多目标定位。
上述的,步骤1中第n个观测站的阵列天线所接收到的信号时域模型表示为:其中,pd表示第d个目标的位置向量,表示第d个目标信号的发射信号时间,sd(t)表示第d个目标信号的复包络,an(pd)表示第d个目标信号相对于第n个观测站的天线阵列流形向量,τn(pd)表示第d个目标信号到达第n个观测站的传播时延,βnd表示第d个目标信号传播至第n个观测站的损耗因子,εn(t)表示第n个观测站中天线阵列的阵元噪声向量;第n个观测站的阵列天线接收信号在第k个子段内的频域模型为:其中,和分别表示和的频域形式,ωq表示第q个采集数据点;步骤3中高维阵列信号频域模型表示为:其中,步骤4中第q个采集数据点的高维阵列协方差矩阵的最大似然估计值表示为:步骤5中,特征值分解表示为:其中,中的列向量是矩阵的MN-D个最小特征值对应的单位特征向量,该矩阵列空间的正交投影矩阵为
上述的,步骤6中的数值优化中的牛顿型迭代算法进行多目标定位具体包含如下内容:
步骤6.1)利用多重信号分类估计算法和Taylor级数迭代定位算法获得目标位置向量的初始估计
步骤6.2)依次对D个目标信号的位置向量进行牛顿迭代,其迭代公式为:其中,m表示迭代次数,0<μ<1表示迭代步长因子,和分别表示目标函数的梯度向量和Hessian矩阵,相应的计算公式分别表示为:
表示矩阵最小特征值所对应的单位特征向量,矩阵和的计算公式分别表示为:
,
和分别对应矩阵除最小特征值以外的其余N-1个特征值及其对应的单位特征向量。
本发明的有益效果:
本发明首先建立到达信号的复包络和载波相位关于目标位置参数的解析模型,通过基2-FFT算法将多站阵列信号时域数据转化成频域数据,并在频域利用子空间正交原理构造直接估计目标位置参数的数学优化模型,然后基于矩阵特征值扰动理论设计出Newton型迭代算法,以实现对多目标的精确定位;相比于传统的多站测向交汇定位和多站测时差交汇定位方法,本发明能够同时利用到达信号的复包络和载波相位信息,通过在底层数据域直接进行目标位置估计,能够取得更高的定位精度,同时还可以避免两步定位模式中存在的门限效应和“目标-量测”匹配问题;此外,本发明公开的定位方法是通过Newton型迭代来实现,具有较快的收敛速度,无需高维搜索,性能可靠、运算高效。
附图说明:
图1为本发明的多观测站多目标直接定位原理图;
图2为本发明的流程示意图;
图3为本发明的定位实例场景示意图;
图4为本发明的定位结果对比示意图。
具体实施方式:
下面结合附图和技术方案对本发明作进一步详细的说明,并通过优选的实施例详细说明本发明的实施方式,但本发明的实施方式并不限于此。
实施例一,参见图1~2所示,一种联合无线电信号复包络和载波相位信息的多目标直接定位方法,具体包含如下步骤:
步骤1.对N个观测站的M通道阵列天线接收系统做时间同步,根据奈奎斯特采样定理采集目标辐射的无线电信号数据,获得阵列信号时域数据;
步骤2.对每个观测站的阵列信号时域数据在时域上划分成K个子段,每个子段内均包含有Q个采集数据点,对每个子段内的Q个数据样本做基2-FFT运算,得到阵列信号频域数据;
步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站,中心站对每个观测站传输的阵列信号频域数据按照观测站的顺序堆栈排列,构造高维阵列信号频域数据;
步骤4.针对每个采集数据点,利用K个子段数据构造高维阵列协方差矩阵的最大似然估计值;
步骤5.对每个采集数据点所构造的高维阵列协方差矩阵进行特征值分解,利用其MN-D个最小特征值对应的单位特征向量构造特征向量矩阵,并计算该矩阵列空间的正交投影矩阵,D为目标信号个数;
步骤6.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法进行多目标定位。
实施例二,参见图1~2,一种联合无线电信号复包络和载波相位信息的多目标直接定位方法,具体包含如下步骤:
步骤1.对N个观测站的M通道阵列天线接收系统做时间同步,根据奈奎斯特采样定理采集目标辐射的无线电信号数据,获得阵列信号时域数据,第n个观测站的阵列天线所接收到的信号时域模型表示为:其中,pd表示第d个目标的位置向量,表示第d个目标信号的发射信号时间,sd(t)表示第d个目标信号的复包络,an(pd)表示第d个目标信号相对于第n个观测站的天线阵列流形向量,τn(pd)表示第d个目标信号到达第n个观测站的传播时延,βnd表示第d个目标信号传播至第n个观测站的损耗因子,εn(t)表示第n个观测站中天线阵列的阵元噪声向量;
步骤2.对每个观测站的阵列信号时域数据在时域上划分成K个子段,每个子段内均包含有Q个采集数据点,对每个子段内的Q个数据样本做基2-FFT运算,得到阵列信号频域数据,第n个观测站的阵列天线接收信号在第k个子段内的频域模型为:其中,和分别表示和的频域形式,ωq表示第q个采集数据点;
步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站,中心站对每个观测站传输的阵列信号频域数据按照观测站的顺序堆栈排列,构造高维阵列信号频域数据,高维阵列信号频域模型表示为:其中,
步骤4.针对每个采集数据点,利用K个子段数据构造高维阵列协方差矩阵的最大似然估计值,第q个采集数据点的高维阵列协方差矩阵的最大似然估计值表示为:
步骤5.对每个采集数据点所构造的高维阵列协方差矩阵进行特征值分解,利用其MN-D个最小特征值对应的单位特征向量构造特征向量矩阵,并计算该矩阵列空间的正交投影矩阵,D为目标信号个数,特征值分解表示为:其中,中的列向量是矩阵的MN-D个最小特征值对应的单位特征向量,该矩阵列空间的正交投影矩阵为
步骤6.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法进行多目标定位,具体包含如下内容:
步骤6.1)利用多重信号分类估计算法和Taylor级数迭代定位算法获得目标位置向量的初始估计
步骤6.2)依次对D个目标信号的位置向量进行牛顿迭代,其迭代公式为:其中,m表示迭代次数,0<μ<1表示迭代步长因子,和分别表示目标函数的梯度向量和Hessian矩阵,相应的计算公式分别表示为:
表示矩阵最小特征值所对应的单位特征向量,矩阵和的计算公式分别表示为:
和分别对应矩阵除最小特征值以外的其余N-1个特征值及其对应的单位特征向量。
参见图3~4所示,结合具体的试验数据对本发明做进一步解释说明:
如图3所示,假设两个待定位目标,两目标的位置坐标分别为(2km,2km)和(4km,4km),现有三个测向站对其进行定位,其位置坐标分别为(6km,6km)、(6km,-6km)和(-6km,6km),每个测向站均安装5元均匀线阵,信号复包络服从零均值圆高斯分布,信号到达三个测向站的复传播系数分别为0.9848+0.1736j、0.6428+0.7660j和0.1736+0.9848j(由不同的传播路径所产生),对信号采用基2-FFT算法的点数为512点,每个频点累积的样本点数为20。
下面将本专利公开的直接定位方法与传统先测向再交汇定位方法的性能进行比较,这里测向采用多重信号分类估计算法,交汇定位采用Taylor级数迭代定位算法。
首先,将信噪比固定为-5dB,图4-1为两种定位方法进行500次蒙特卡洛实验的定位结果,图4-2为目标1定位均方根误差随着信噪比的变化曲线,图4-3为目标2定位均方根误差随着信噪比的变化曲线,图4-1给出了两种方法进行了500次蒙特卡洛实验的定位结果;然后,改变信噪比,图4-2和图4-3分别给出了两个目标的定位均方根误差随着信噪比的变化曲线。
由图4-1至图4-3可以明显看出,相比于传统的先测向再交汇定位方法,本专利公开的直接定位方法可以显著提升定位精度,并且信噪比越低,其优势愈加明显。
本发明并不局限于上述具体实施方式,本领域技术人员还可据此做出多种变化,但任何与本发明等同或者类似的变化都应涵盖在本发明权利要求的范围内。
Claims (1)
1.一种联合无线电信号复包络和载波相位信息的多目标直接定位方法,具体包含如下步骤:
步骤1.对N个观测站的M通道阵列天线接收系统做时间同步,根据奈奎斯特采样定理采集目标辐射的无线电信号数据,获得阵列信号时域数据;
步骤2.对每个观测站的阵列信号时域数据在时域上划分成K个子段,每个子段内均包含有Q个采集数据点,对每个子段内的Q个数据样本做基2-FFT运算,得到阵列信号频域数据;
步骤3.每个观测站将所获得的阵列信号频域数据传输至中心站,中心站对每个观测站传输的阵列信号频域数据按照观测站的顺序堆栈排列,构造高维阵列信号频域数据;
步骤4.针对每个采集数据点,利用K个子段数据构造高维阵列协方差矩阵的最大似然估计值;
步骤5.对每个采集数据点所构造的高维阵列协方差矩阵进行特征值分解,利用其MN-D个最小特征值对应的单位特征向量构造特征向量矩阵,并计算特征向量矩阵列空间的正交投影矩阵,D为目标个数;
步骤6.基于矩阵特征值扰动公式设计数值优化中的牛顿型迭代算法进行多目标定位;
其中,步骤1中第n个观测站的阵列天线所接收到的信号时域模型表示为:其中,1≤n≤N,pd表示第d个目标的位置向量,表示第d个目标信号的发射信号时间,sd(t)表示第d个目标信号的复包络,an(pd)表示第d个目标信号相对于第n个观测站的天线阵列流形向量,tn(pd)表示第d个目标信号到达第n个观测站的传播时延,βnd表示第d个目标信号传播至第n个观测站的损耗因子,εn(t)表示第n个观测站中天线阵列的阵元噪声向量;步骤2中第n个观测站的阵列天线接收信号在第k个子段内的频域模型为:其中,和分别表示和的频域形式,ωq表示第q个采集数据点;步骤3中高维阵列信号频域模型表示为:其中,步骤4中第q个采集数据点的高维阵列协方差矩阵的最大似然估计值表示为:步骤5中,特征值分解表示为:其中,中的列向量是矩阵的MN-D个最小特征值对应的单位特征向量,该矩阵列空间的正交投影矩阵为步骤6中的数值优化中的牛顿型迭代算法进行多目标定位具体包含如下内容:
步骤6.1)利用多重信号分类估计算法和Taylor级数迭代定位算法获得目标位置向量的初始估计
步骤6.2)依次对D个目标信号的位置向量进行牛顿迭代,其迭代公式为:其中,1≤d≤D,m表示迭代次数,0<μ<1表示迭代步长因子,和分别表示目标函数的梯度向量和Hessian矩阵,相应的计算公式分别表示为:
表示矩阵最小特征值所对应的单位特征向量,矩阵和的计算公式分别表示为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<mo>(</mo>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<msup>
<mi>A</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<mover>
<mo>&Pi;</mo>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>A</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mover>
<mi>A</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>l</mi>
<mi>H</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<mover>
<mo>&Pi;</mo>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>X</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<mo>(</mo>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mover>
<mi>A</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>H</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<mover>
<mo>&Pi;</mo>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>A</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<mover>
<mo>&Pi;</mo>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>A</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msubsup>
<mover>
<mi>A</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<mi>H</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<mover>
<mo>&Pi;</mo>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>A</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mover>
<mi>A</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<mi>H</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<mover>
<mo>&Pi;</mo>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>A</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>q</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mo>)</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
<mo>(</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mi>n</mi>
</msub>
<mo>(</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>u</mi>
<mi>n</mi>
</msub>
<mo>(</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>u</mi>
<mi>n</mi>
<mi>H</mi>
</msubsup>
<mo>(</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
和分别对应矩阵除最
小特征值以外的其余N-1个特征值及其对应的单位特征向量。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510887968.5A CN105510948A (zh) | 2015-12-05 | 2015-12-05 | 联合无线电信号复包络和载波相位信息的多目标直接定位方法 |
CN2015108879685 | 2015-12-05 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105954712A CN105954712A (zh) | 2016-09-21 |
CN105954712B true CN105954712B (zh) | 2018-02-02 |
Family
ID=55719056
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510887968.5A Withdrawn CN105510948A (zh) | 2015-12-05 | 2015-12-05 | 联合无线电信号复包络和载波相位信息的多目标直接定位方法 |
CN201610255799.8A Active CN105738865B (zh) | 2015-12-05 | 2016-04-22 | 信号波形已知条件下的多目标直接定位方法 |
CN201610255410.XA Active CN105954712B (zh) | 2015-12-05 | 2016-04-22 | 联合无线电信号复包络和载波相位信息的多目标直接定位方法 |
Family Applications Before (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510887968.5A Withdrawn CN105510948A (zh) | 2015-12-05 | 2015-12-05 | 联合无线电信号复包络和载波相位信息的多目标直接定位方法 |
CN201610255799.8A Active CN105738865B (zh) | 2015-12-05 | 2016-04-22 | 信号波形已知条件下的多目标直接定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (3) | CN105510948A (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106483496A (zh) * | 2016-10-19 | 2017-03-08 | 河南城建学院 | 基于chan算法与改进牛顿迭代的联合时差定位方法 |
CN107132505B (zh) * | 2017-05-19 | 2019-11-08 | 中国人民解放军信息工程大学 | 直达与非直达混合场景中的多目标直接定位方法 |
CN108182474B (zh) * | 2017-12-27 | 2020-03-13 | 中国人民解放军战略支援部队信息工程大学 | 基于未校正阵列和神经网络的多目标直接定位方法 |
CN108872932B (zh) * | 2018-05-03 | 2021-04-02 | 中国人民解放军战略支援部队信息工程大学 | 基于神经网络的超视距目标直接定位结果纠偏方法 |
CN108770059A (zh) * | 2018-05-31 | 2018-11-06 | 清华大学 | 基于波束空间的低通信量高精度定位方法及系统 |
CN109031196A (zh) * | 2018-07-19 | 2018-12-18 | 中国人民解放军战略支援部队信息工程大学 | 基于运动观测站对多信号源的最大似然直接定位方法 |
CN109298388B (zh) * | 2018-08-21 | 2020-09-18 | 中国人民解放军战略支援部队信息工程大学 | 基于方位信息的超视距目标地理坐标直接估计方法 |
CN109581291B (zh) * | 2018-12-11 | 2023-08-01 | 哈尔滨工程大学 | 一种基于人工蜂群的直接定位方法 |
CN110632555B (zh) * | 2019-08-26 | 2022-12-13 | 中国人民解放军战略支援部队信息工程大学 | 一种基于矩阵特征值扰动的tdoa直接定位方法 |
CN111199280B (zh) * | 2019-12-20 | 2022-10-14 | 中国人民解放军战略支援部队信息工程大学 | 短波信道模型误差存在下联合信号复包络和载波相位信息的多站目标源地理坐标估计方法 |
CN112578335B (zh) * | 2020-11-29 | 2023-05-02 | 西北工业大学 | 一种强弱信号的直接定位与信源数估计方法 |
CN113075649B (zh) * | 2021-03-30 | 2022-07-22 | 电子科技大学 | 一种适用于分布式网络化雷达的信号级直接定位方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7450064B2 (en) * | 2005-03-22 | 2008-11-11 | Qualcomm, Incorporated | Methods and systems for deriving seed position of a subscriber station in support of unassisted GPS-type position determination in a wireless communication system |
CN101750600B (zh) * | 2008-12-18 | 2012-02-29 | 同济大学 | 一种用于测绘与导航的实时高精度定位方法 |
KR101767924B1 (ko) * | 2012-12-05 | 2017-08-14 | 한화테크윈 주식회사 | 다중 목표물 위치 추정 시스템 및 방법 |
CN103096465B (zh) * | 2013-01-11 | 2015-06-24 | 南京师范大学 | 一种环境自适应的多目标直接定位方法 |
CN103338516B (zh) * | 2013-07-19 | 2016-06-01 | 大连理工大学 | 一种基于总体最小二乘的无线传感器网络两步定位方法 |
CN103713277B (zh) * | 2013-12-19 | 2016-02-10 | 中国航天科工集团八五一一研究所 | 一种基于位置信息场的辐射源定位方法 |
CN103826298B (zh) * | 2014-03-05 | 2017-02-08 | 浙江农林大学 | 一种协作式迭代优化的无线传感器网络定位计算方法 |
CN104537257B (zh) * | 2015-01-12 | 2017-09-26 | 电子科技大学 | 一种基于时差的分布式自适应直接定位方法 |
CN104656104B (zh) * | 2015-02-27 | 2017-06-13 | 清华大学 | 基于最大似然估计的卫星导航欺骗信号识别方法及系统 |
-
2015
- 2015-12-05 CN CN201510887968.5A patent/CN105510948A/zh not_active Withdrawn
-
2016
- 2016-04-22 CN CN201610255799.8A patent/CN105738865B/zh active Active
- 2016-04-22 CN CN201610255410.XA patent/CN105954712B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105954712A (zh) | 2016-09-21 |
CN105738865A (zh) | 2016-07-06 |
CN105738865B (zh) | 2018-03-27 |
CN105510948A (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105954712B (zh) | 联合无线电信号复包络和载波相位信息的多目标直接定位方法 | |
CN109298388B (zh) | 基于方位信息的超视距目标地理坐标直接估计方法 | |
CN105911521B (zh) | 联合无线电信号复包络和载波相位信息的超视距目标直接定位方法 | |
CN107132505B (zh) | 直达与非直达混合场景中的多目标直接定位方法 | |
CN103713277B (zh) | 一种基于位置信息场的辐射源定位方法 | |
CN110290491B (zh) | 一种基于多径辅助的室内目标定位方法 | |
CN108872932B (zh) | 基于神经网络的超视距目标直接定位结果纠偏方法 | |
CN104519571B (zh) | 一种基于rss的室内定位方法 | |
CN104619020A (zh) | 基于rssi和toa测距的wifi室内定位方法 | |
CN106597363A (zh) | 一种室内wlan环境下的行人定位方法 | |
CN109061554A (zh) | 一种基于空间离散网格动态更新的目标到达角度估计方法 | |
CN104469937A (zh) | 压缩感知定位技术中的传感器高效部署方法 | |
CN108872971B (zh) | 一种基于运动单阵列的目标定位方法与装置 | |
CN109975754A (zh) | 联合角度、时延与多普勒信息的运动多站直接定位方法 | |
CN111199280B (zh) | 短波信道模型误差存在下联合信号复包络和载波相位信息的多站目标源地理坐标估计方法 | |
CN108182474A (zh) | 基于未校正阵列和神经网络的多目标直接定位方法 | |
CN110657806B (zh) | 一种基于CKF、chan解算和Savitzky-Golay平滑滤波的位置解算方法 | |
CN113281701B (zh) | 协同短波多站角度与三星时差的超视距目标直接定位方法 | |
CN102833016B (zh) | 一种在多径环境下针对同源信号时差的测定方法 | |
CN109001687A (zh) | 基于广义旁瓣相消结构的机载雷达空时自适应滤波方法 | |
CN105652235A (zh) | 基于线性回归算法的wlan室内定位多用户rss融合方法 | |
CN103338024B (zh) | 天线组阵中时延的互补卡尔曼滤波装置与方法 | |
CN110413939B (zh) | 一种基于原子范数的到达角估计方法 | |
CN103096465A (zh) | 一种环境自适应的多目标直接定位方法 | |
CN110531309A (zh) | 存在幅相误差时基于原子范数的相关信号到达角估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |