CN112782647A - 信息联合的二次等式约束最小二乘辐射源定位方法 - Google Patents
信息联合的二次等式约束最小二乘辐射源定位方法 Download PDFInfo
- Publication number
- CN112782647A CN112782647A CN202011473715.0A CN202011473715A CN112782647A CN 112782647 A CN112782647 A CN 112782647A CN 202011473715 A CN202011473715 A CN 202011473715A CN 112782647 A CN112782647 A CN 112782647A
- Authority
- CN
- China
- Prior art keywords
- radiation source
- observation
- pseudo
- equation
- information
- 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.)
- Granted
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 127
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000013598 vector Substances 0.000 claims abstract description 60
- 239000005433 ionosphere Substances 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims description 6
- 230000004807 localization Effects 0.000 claims 3
- 239000007787 solid Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 230000014509 gene expression Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000004590 computer program Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001010 compromised effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Images
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/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/04—Position of source determined by a plurality of spaced direction-finders
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于辐射源定位领域,特别涉及一种信息联合的二次等式约束最小二乘辐射源定位方法,包含:利用若干测向站对地面短波辐射源进行二维测向,获取该辐射源方位角及仰角信息,同时构建两者的非线性观测方程;并获取电离层虚高观测模型;结合电离层虚高观测模型将方位角及仰角的非线性观测方程分别转换为伪线性观测方程并合并,形成二维角度伪线性观测方程;构建二次等式约束最小二乘优化模型,并将模型约束优化问题转化为关于拉格朗日乘子的非线性方程组问题;迭代求解并依据求解结果确定短波辐射源地心地固坐标系下的位置向量及其协方差矩阵,进而获取短波辐射源经纬度。本发明提升对短波辐射源的定位精度,获取更精确的辐射源目标定位信息。
Description
技术领域
本发明属于辐射源定位领域,特别涉及一种信息联合的二次等式约束最小二乘辐射源定位方法。
背景技术
无线信号定位技术广泛应用于通信、雷达、目标监测、导航遥测、地震勘测、射电天文、紧急救助、安全管理等领域,其在工业生产和国防安全中都发挥着重要作用。对目标进行定位(即位置参数估计)可以使用雷达、激光、声纳等有源设备来完成,该类技术称为有源定位技术,它具有全天候、高精度等优点。然而,有源定位系统通常需要依靠发射大功率电磁信号来完成,因此极易暴露自己位置,容易被对方发现,从而遭到对方电子干扰影响,导致定位性能急剧恶化,甚至会危及系统自身的安全性和可靠性。目标定位还可以利用目标(主动)辐射或者(被动)散射的无线电信号来实现,该类技术称为无源定位技术,它是指在观测站不主动发射电磁信号的情况下,通过接收目标辐射或者散射的无线电信号来估计目标的位置参数。与有源定位系统相比,无源定位系统具有生存能力强、侦察作用距离远等优点,从而受到国内外学者的广泛关注和深入研究。依据观测站数量可以将无源定位系统划分成单站无源定位系统和多站无源定位系统两大类,其中多站定位系统可以提供更多的观测量,从而提高目标定位精度。众所周知,短波测向交汇定位是一类十分重要的多站无源定位技术,其主要是针对远距离超视距辐射源进行定位,该定位技术的基本原理是利用每个测向站所提供的方位角信息以及地球椭圆模型确定短波辐射源的位置信息,该类定位方法的优势是能够快速确定短波辐射源的位置坐标,但不足之处在于仅能利用各个测向站所提供的方位角信息,难以直接利用各个测向站所提供的仰角信息进行交汇定位,因此其定位精度还存在进一步提升的空间。
发明内容
为此,针对短波辐射源定位问题,本发明提供一种信息联合的二次等式约束最小二乘辐射源定位方法,不仅可以利用方位角信息,还可以利用仰角信息和电离层虚高信息,进行联合定位,提高短波测向交汇定位精度。
按照本发明所提供的设计方案,一种信息联合的二次等式约束最小二乘辐射源定位方法,用于联合方位角、仰角和电离层虚高信息对短波辐射源进行测向交汇定位,包含如下内容:
利用若干测向站对地面短波辐射源进行二维测向,获取该辐射源方位角及仰角信息,同时构建两者的非线性观测方程;并通过测向站有源探测获取电离层虚高观测模型;
结合电离层虚高观测模型将方位角及仰角的非线性观测方程分别转换为伪线性观测方程并合并,形成二维角度伪线性观测方程;
依据二维角度伪线性观测方程并结合地球椭圆模型构建用于短波辐射源定位的二次等式约束最小二乘优化模型,并利用拉格朗日乘子法将模型约束优化问题转化为关于拉格朗日乘子的非线性方程组问题;
对非线性方程组进行迭代求解,依据拉格朗日乘子估计值确定短波辐射源地心地固坐标系下的位置向量及其协方差矩阵,进而获取短波辐射源经纬度。
作为本发明信息联合的二次等式约束最小二乘辐射源定位方法,进一步的,假设在地面放置M个测向站,利用该M个测向站对地面短波辐射源进行二维测向,依据地面短波辐射源经纬度及测向站经纬度来获取短波辐射源方位角及仰角。
作为本发明信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,依据测向站及短波辐射源在地心地固坐标系下位置向量来表示方位角非线性观测方程;根据测向站与辐射源之间地心角信息及短波辐射源信号到测向站传播路径所对应的电离层虚高来表示仰角非线性观测方程;根据测向站对应的电离层虚高真实值及电离层虚高观测值来表示电离层虚高观测模型。
作为本发明信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,利用三角函数性质将方位角非线性观测方程转化为方位角伪线性观测方程。
作为本发明信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,基于多项式求根原理和引入辅助变量将仰角非线性观测方程转化成仰角伪线性观测方程。
作为本发明信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,利用M个测向站的仰角观测值和电离层虚高观测值建立M个一元二次方程;利用该M个一元二次方程的正根并结合辐射源位置向量来建立仰角伪线性观测方程,其中,辐射源位置向量为通过引入辅助变量扩维后的辐射源位置向量。
作为本发明信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,二维角度伪线性观测方程表示为:其中, ξ1、分别表示方位角伪线性观测方程中方位角伪线性观测向量、观测误差向量、观测矩阵,ξ2分别表示仰角伪线性观测方程中仰角伪线性观测向量、观测矩阵、观测误差向量,OM×1表示M个测向站向量,表示辐射源位置向量。
作为本发明信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,二次等式约束最小二乘优化模型表示为:
其中,Q为二维角度伪线性观测误差向量ξ的协方差矩阵,Re为已知距离数值,Γ1、Γ2、η为预设单位向量矩阵。
其中,λ1和λ2表示两个未知的拉格朗日乘子,
作为本发明信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,利用Newton迭代算法求解非线性方程组;根据拉格朗日乘子估计值确定短波辐射源地心地固坐标系下的位置向量表示为:其中,和表示拉格朗日乘子估计值,I3、O3×1为预设单位向量。
本发明的有益效果:
本发明通过联合短波辐射源方位角、仰角和电离层虚高信息,将短波测向交汇定位问题转化成二次等式约束最小二乘优化问题,并利用拉格朗日乘子法进行求解,以获得短波辐射源的位置信息,相比已有短波测向交汇定位方法,能够进一步提升对短波辐射源的定位精度,获取更精确的辐射源目标定位信息,具有较好的应用前景。
附图说明:
图1为实施例中辐射源定位流程示意;
图2为实施例中短波辐射源定位均方根误差随着方位角估计误差标准差的变化曲线示意;
图3为实施例中短波辐射源定位均方根误差随着仰角估计误差标准差的变化曲线示意;
图4为实施例中短波辐射源定位均方根误差随着电离层虚高观测误差标准差的变化曲线示意。
具体实施方式:
为使本发明的目的、技术方案和优点更加清楚、明白,下面结合附图和技术方案对本发明作进一步详细的说明。
本发明实施例,提供一种信息联合的二次等式约束最小二乘辐射源定位方法,用于联合方位角、仰角和电离层虚高信息对短波辐射源进行测向交汇定位,包含如下内容:
利用若干测向站对地面短波辐射源进行二维测向,获取该辐射源方位角及仰角信息,同时构建两者的非线性观测方程;并通过测向站有源探测获取电离层虚高观测模型;
结合电离层虚高观测模型将方位角及仰角的非线性观测方程分别转换为伪线性观测方程并合并,形成二维角度伪线性观测方程;
依据二维角度伪线性观测方程并结合地球椭圆模型构建用于短波辐射源定位的二次等式约束最小二乘优化模型,并利用拉格朗日乘子法将模型约束优化问题转化为关于拉格朗日乘子的非线性方程组问题;
对非线性方程组进行迭代求解,依据拉格朗日乘子估计值确定短波辐射源地心地固坐标系下的位置向量及其协方差矩阵,进而获取短波辐射源经纬度。
过联合短波辐射源方位角、仰角和电离层虚高信息,将短波测向交汇定位问题转化成二次等式约束最小二乘优化问题,并利用拉格朗日乘子法进行求解,以获得短波辐射源的位置信息,相比已有短波测向交汇定位方法,能够进一步提升对短波辐射源的定位精度,提升定位方案在实际应用中的安全性和可靠性。
进一步,参见图1所示,首先利用多个测向站获得短波辐射源方位角、仰角以及信号传播路径所对应的电离层虚高信息。接着利用三角函数性质将方位角非线性观测方程转化成方位角伪线性观测方程,基于多项式求根原理和引入辅助变量将仰角非线性观测方程转化成仰角伪线性观测方程,并合并这两类伪线性观测方程,形成二维角度伪线性观测方程。然后结合地球椭圆模型和辅助变量的特点建立用于短波辐射源定位的二次等式约束最小二乘优化模型,利用拉格朗日乘子法将该约束优化问题转化成关于拉格朗日乘子的非线性方程组问题,并设计Newton迭代算法求解该方程组,以获得拉格朗日乘子估计值。最后利用拉格朗日乘子估计值确定短波辐射源地心地固坐标系下的位置向量及其协方差矩阵,并通过Gauss-Newton迭代方法进一步确定短波辐射源的经度和纬度。
作为本发明实施例中信息联合的二次等式约束最小二乘辐射源定位方法,进一步的,假设在地面放置M个测向站,利用该M个测向站对地面短波辐射源进行二维测向,依据地面短波辐射源经纬度及测向站经纬度来获取短波辐射源方位角及仰角。进一步地,依据测向站及短波辐射源在地心地固坐标系下位置向量来表示方位角非线性观测方程;根据测向站与辐射源之间地心角信息及短波辐射源信号到测向站传播路径所对应的电离层虚高来表示仰角非线性观测方程;根据测向站对应的电离层虚高真实值及电离层虚高观测值来表示电离层虚高观测模型。
式中us表示短波辐射源在地心地固坐标系下的位置向量(其为待求参量);um表示第m个测向站在地心地固坐标系下的位置向量(其为已知参量);εm1表示方位角估计误差,其服从零均值的独立高斯分布,并且方差为相关表达式为
式中Re=6378.160km和e=0.081819643716348。
式中Ro≈6370km;hm表示短波辐射源信号到达第m个测向站的传播路径所对应的电离层虚高;εm2表示仰角估计误差,其服从零均值的独立高斯分布,并且方差为θm表示第m个测向站与辐射源之间的地心角的一半,其表达式为
每个测向站利用有源探测(包括垂直探测、斜向探测以及返回散射等)手段获得电离层虚高观测值,其观测模型为
式中表示第m个测向站对应的电离层虚高观测值;hm(1≤m≤M)表示第m个测向站对应的电离层虚高真实值;δm(1≤m≤M)表示第m个测向站对应的电离层虚高观测误差,其服从零均值的独立高斯分布,并且方差为 表示由各个测向站对应的电离层虚高观测值所形成的列向量;h=[h1 h2 … hM]T表示由各个测向站对应的电离层虚高真实值所形成的列向量;δ=[δ1 δ2 … δM]T表示电离层虚高观测误差向量,其服从零均值的高斯分布,并且协方差矩阵为
作为本发明实施例中信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,利用三角函数性质将方位角非线性观测方程转化为方位角伪线性观测方程。
方位角伪线性观测方程如下式所示:
ξ1表示方位角伪线性观测误差向量,其近似服从零均值的高斯分布,并且协方差矩阵为
式中
作为本发明实施例中信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,基于多项式求根原理和引入辅助变量将仰角非线性观测方程转化成仰角伪线性观测方程。进一步地,利用M个测向站的仰角观测值和电离层虚高观测值建立M个一元二次方程;利用该M个一元二次方程的正根并结合辐射源位置向量来建立仰角伪线性观测方程,其中,辐射源位置向量为通过引入辅助变量扩维后的辐射源位置向量。
式中
ξ2表示仰角伪线性观测误差向量,其近似服从零均值的高斯分布,并且协方差矩阵为
式中
作为本发明实施例中信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,二维角度伪线性观测方程表示为:其中, ξ1、分别表示方位角伪线性观测方程中方位角伪线性观测向量、观测误差向量、观测矩阵, ξ2分别表示仰角伪线性观测方程中仰角伪线性观测向量、观测矩阵、观测误差向量,OM×1表示M个测向站向量,表示辐射源位置向量。ξ表示二维角度伪线性观测误差向量,其近似服从零均值的高斯分布,并且协方差矩阵为
Q=E[ξξT]=blkdiag{Q1,Q2}=diag[q11 q21 … qM1 q12 q22 … qM2]
作为本发明实施例中信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,二次等式约束最小二乘优化模型表示为:
其中,Q为二维角度伪线性观测误差向量ξ的协方差矩阵,Re为已知距离数值,Γ1、Γ2、η为预设单位向量矩阵。这些单位向量矩阵可表示为:
其中,λ1和λ2表示两个未知的拉格朗日乘子,f1(λ1,λ2)和f2(λ1,λ2)表示关于λ1和λ2的非线性函数,相应的表达式可表示为:
作为本发明实施例中信息联合的二次等式约束最小二乘辐射源定位方法,进一步地,利用Newton迭代算法求解非线性方程组;根据拉格朗日乘子估计值确定短波辐射源地心地固坐标系下的位置向量表示为:其中,和表示拉格朗日乘子估计值,I3、O3×1为预设单位向量。
μ(k)表示步长因子,其可以通过求解如下一维非线性优化问题获得
式中
基于以上定位估计结果,可利用Gauss-Newton迭代方法确定短波辐射源的经度和纬度,首先确定迭代初值,如下式所示
为验证本发明方案有效性,下面结合实验数据做进一步解释说明:
假设共有5个测向站利用方位角、仰角和电离层虚高信息对短波辐射源进行定位,测向站的经纬度和短波辐射源信号到达测向站的电离层虚高如表1所示,短波辐射源的经度为125.05度,纬度为27.47度。
表1测向站的经纬度和电离层虚高
首先将仰角估计误差标准差设为0.5度(即σm2=0.5度(1≤m≤M)),电离层虚高观测误差标准差设为5公里(即σm3=5公里(1≤m≤M)),图2给出了短波辐射源定位均方根误差随着方位角估计误差标准差σm1的变化曲线;然后将方位角估计误差标准差设为0.5度(即σm1=0.5度(1≤m≤M)),电离层虚高观测误差标准差设为5公里(即σm3=5公里(1≤m≤M)),图3给出了短波辐射源定位均方根误差随着仰角估计误差标准差σm2的变化曲线;最后将方位角估计误差标准差设为0.5度(即σm1=0.5度(1≤m≤M)),仰角估计误差标准差设为0.5度(即σm2=0.5度(1≤m≤M)),图4给出了短波辐射源定位均方根误差随着电离层虚高观测误差标准差σm3的变化曲线。
从图2至图4中可以看出,本案公开的定位方法确实可以提高短波测向交汇定位的精度,并且方位角估计误差越大,新方法的优势越明显;仰角估计误差和电离层虚高观测误差越小,新方法的优势越明显。
除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对步骤、数字表达式和数值并不限制本发明的范围。
基于上述的方法或系统,本发明实施例还提供一种网络设备,包括:一个或多个处理器;存储装置,用于存储一个或多个程序,当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现上述的系统或执行上述的方法。
基于上述的系统,本发明实施例还提供一种计算机可读介质,其上存储有计算机程序,其中,该程序被处理器执行时实现上述的系统。
本发明实施例所提供的装置,其实现原理及产生的技术效果和前述系统实施例相同,为简要描述,装置实施例部分未提及之处,可参考前述系统实施例中相应内容。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统和装置的具体工作过程,可以参考前述系统实施例中的对应过程,在此不再赘述。
在这里示出和描述的所有示例中,任何具体值应被解释为仅仅是示例性的,而不是作为限制,因此,示例性实施例的其他示例可以具有不同的值。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、系统和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述系统的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种信息联合的二次等式约束最小二乘辐射源定位方法,用于联合方位角、仰角和电离层虚高信息对短波辐射源进行测向交汇定位,其特征在于,包含如下内容:
利用若干测向站对地面短波辐射源进行二维测向,获取该辐射源方位角及仰角信息,同时构建两者的非线性观测方程;并通过测向站有源探测获取电离层虚高观测模型;
结合电离层虚高观测模型将方位角及仰角的非线性观测方程分别转换为伪线性观测方程并合并,形成二维角度伪线性观测方程;
依据二维角度伪线性观测方程并结合地球椭圆模型构建用于短波辐射源定位的二次等式约束最小二乘优化模型,并利用拉格朗日乘子法将模型约束优化问题转化为关于拉格朗日乘子的非线性方程组问题;
对非线性方程组进行迭代求解,依据拉格朗日乘子估计值确定短波辐射源地心地固坐标系下的位置向量及其协方差矩阵,进而获取短波辐射源经纬度。
2.根据权利要求1所述的信息联合的二次等式约束最小二乘辐射源定位方法,其特征在于,假设在地面放置M个测向站,利用该M个测向站对地面短波辐射源进行二维测向,依据地面短波辐射源经纬度及测向站经纬度来获取短波辐射源方位角及仰角。
3.根据权利要求1或2所述的信息联合的二次等式约束最小二乘辐射源定位方法,其特征在于,依据测向站及短波辐射源在地心地固坐标系下位置向量来表示方位角非线性观测方程;根据测向站与辐射源之间地心角信息及短波辐射源信号到测向站传播路径所对应的电离层虚高来表示仰角非线性观测方程;根据测向站对应的电离层虚高真实值及电离层虚高观测值来表示电离层虚高观测模型。
4.根据权利要求1所述的信息联合的二次等式约束最小二乘辐射源定位方法,其特征在于,利用三角函数性质将方位角非线性观测方程转化为方位角伪线性观测方程。
5.根据权利要求1所述的信息联合的二次等式约束最小二乘辐射源定位方法,其特征在于,基于多项式求根原理和引入辅助变量将仰角非线性观测方程转化成仰角伪线性观测方程。
6.根据权利要求5所述的信息联合的二次等式约束最小二乘辐射源定位方法,其特征在于,利用M个测向站的仰角观测值和电离层虚高观测值建立M个一元二次方程;利用该M个一元二次方程的正根并结合辐射源位置向量来建立仰角伪线性观测方程,其中,辐射源位置向量为通过引入辅助变量扩维后的辐射源位置向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011473715.0A CN112782647B (zh) | 2020-12-15 | 2020-12-15 | 信息联合的二次等式约束最小二乘辐射源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011473715.0A CN112782647B (zh) | 2020-12-15 | 2020-12-15 | 信息联合的二次等式约束最小二乘辐射源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112782647A true CN112782647A (zh) | 2021-05-11 |
CN112782647B CN112782647B (zh) | 2023-05-16 |
Family
ID=75750888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011473715.0A Active CN112782647B (zh) | 2020-12-15 | 2020-12-15 | 信息联合的二次等式约束最小二乘辐射源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112782647B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2764149C1 (ru) * | 2021-06-24 | 2022-01-13 | Акционерное общество "Концерн "Созвездие" | Угломерно-корреляционный способ определения местоположения наземных источников радиоизлучения |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7064702B1 (en) * | 2005-03-01 | 2006-06-20 | The Boeing Company | System, method and computer program product for reducing quadratic phase errors in synthetic aperture radar signals |
US20180149729A1 (en) * | 2015-06-17 | 2018-05-31 | Thales | Method and device for locating an electromagnetic emission source and system implementing such a method |
CN108761399A (zh) * | 2018-06-01 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 一种无源雷达目标定位方法与装置 |
CN109298388A (zh) * | 2018-08-21 | 2019-02-01 | 中国人民解放军战略支援部队信息工程大学 | 基于方位信息的超视距目标地理坐标直接估计方法 |
CN109975749A (zh) * | 2019-02-26 | 2019-07-05 | 中国人民解放军战略支援部队信息工程大学 | 一种校正源存在条件下的短波单站直接定位方法 |
CN109975755A (zh) * | 2019-02-26 | 2019-07-05 | 中国人民解放军战略支援部队信息工程大学 | 一种校正源存在条件下的短波多站直接定位方法 |
CN110568403A (zh) * | 2019-08-10 | 2019-12-13 | 中国人民解放军战略支援部队信息工程大学 | 一种协同短波与卫星系统的超视距目标无源定位方法 |
CN111199281A (zh) * | 2019-12-20 | 2020-05-26 | 中国人民解放军战略支援部队信息工程大学 | 基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法 |
CN111308530A (zh) * | 2020-02-17 | 2020-06-19 | 中国人民解放军战略支援部队信息工程大学 | 一种基于二维波达方向的短波多站和单星协同直接定位方法 |
CN111551895A (zh) * | 2020-04-25 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | 基于加权多维标度和拉格朗日乘子技术的运动辐射源tdoa和fdoa定位方法 |
-
2020
- 2020-12-15 CN CN202011473715.0A patent/CN112782647B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7064702B1 (en) * | 2005-03-01 | 2006-06-20 | The Boeing Company | System, method and computer program product for reducing quadratic phase errors in synthetic aperture radar signals |
US20180149729A1 (en) * | 2015-06-17 | 2018-05-31 | Thales | Method and device for locating an electromagnetic emission source and system implementing such a method |
CN108761399A (zh) * | 2018-06-01 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 一种无源雷达目标定位方法与装置 |
CN109298388A (zh) * | 2018-08-21 | 2019-02-01 | 中国人民解放军战略支援部队信息工程大学 | 基于方位信息的超视距目标地理坐标直接估计方法 |
CN109975749A (zh) * | 2019-02-26 | 2019-07-05 | 中国人民解放军战略支援部队信息工程大学 | 一种校正源存在条件下的短波单站直接定位方法 |
CN109975755A (zh) * | 2019-02-26 | 2019-07-05 | 中国人民解放军战略支援部队信息工程大学 | 一种校正源存在条件下的短波多站直接定位方法 |
CN110568403A (zh) * | 2019-08-10 | 2019-12-13 | 中国人民解放军战略支援部队信息工程大学 | 一种协同短波与卫星系统的超视距目标无源定位方法 |
CN111199281A (zh) * | 2019-12-20 | 2020-05-26 | 中国人民解放军战略支援部队信息工程大学 | 基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法 |
CN111308530A (zh) * | 2020-02-17 | 2020-06-19 | 中国人民解放军战略支援部队信息工程大学 | 一种基于二维波达方向的短波多站和单星协同直接定位方法 |
CN111551895A (zh) * | 2020-04-25 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | 基于加权多维标度和拉格朗日乘子技术的运动辐射源tdoa和fdoa定位方法 |
Non-Patent Citations (10)
Title |
---|
JIEXIN YIN 等: "Direct Localization of Multiple Stationary Narrowband Sources Based on Angle and Doppler", IEEE COMMUNICATIONS LETTERS * |
K.YANG 等: "Constrained total least-squares location algorithm using time-difference-of-arrival measurements", TRANSACTION ON VEHICULAR TECHNOLOGY * |
Y.H. HU等: "Optimal linear spectral unmixing", IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING * |
周恭谦 等: "改进的非完全约束加权最小二乘TDOA/FDOA无源定位方法", 系统工程与电子技术 * |
左燕;陈志猛;蔡立平;: "基于约束总体最小二乘的单站DOA/TDOA联合误差校正与定位算法", 电子与信息学报 * |
张杰;蒋建中;郭军利;: "利用约束条件的定位偏差改进算法", 信息工程大学学报 * |
曲长文;徐征;李炳荣;苏峰;: "一种新的基于角度和时差的稳健定位跟踪算法", 信号处理 * |
杨泽宇 等: "一种电离层虚高观测误差存在条件下的超视距直接定位方法", 信息工程大学学报 * |
王鼎 等: "基于外辐射源的约束总体最小二乘定位算法及其理论性能分析", 中国科学:信息科学 * |
王鼎;张莉;吴瑛;: "基于角度信息的约束总体最小二乘无源定位算法", 中国科学(E辑:信息科学) * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2764149C1 (ru) * | 2021-06-24 | 2022-01-13 | Акционерное общество "Концерн "Созвездие" | Угломерно-корреляционный способ определения местоположения наземных источников радиоизлучения |
Also Published As
Publication number | Publication date |
---|---|
CN112782647B (zh) | 2023-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106597364B (zh) | 一种单天线单站无源定位的目标辐射源初始位置估计方法 | |
CN108761399B (zh) | 一种无源雷达目标定位方法与装置 | |
Noroozi et al. | Efficient weighted least squares estimator for moving target localization in distributed MIMO radar with location uncertainties | |
Kazemi et al. | Efficient convex solution for 3-D localization in MIMO radars using delay and angle measurements | |
CN105137424A (zh) | 一种杂波背景下实波束扫描雷达角超分辨方法 | |
Kim et al. | Two‐step estimator for moving‐emitter geolocation using time difference of arrival/frequency‐difference of arrival measurements | |
CN111308530A (zh) | 一种基于二维波达方向的短波多站和单星协同直接定位方法 | |
He et al. | Bias compensation for AOA-geolocation of known altitude target using single satellite | |
Amiri et al. | Closed‐form positioning in MIMO radars with antenna location uncertainties | |
Xiong et al. | Efficient bias reduction approach of time‐of‐flight‐based wireless localisation networks in NLOS states | |
Yang et al. | Joint position and velocity estimation of a moving target in multistatic radar by bistatic range, TDOA, and Doppler shifts | |
CN112782647B (zh) | 信息联合的二次等式约束最小二乘辐射源定位方法 | |
Qu et al. | Recursive source localisation by time difference of arrival sensor networks with sensor position uncertainty | |
Kang et al. | An Efficient Hybrid Multi-Station TDOA and Single-Station AOA Localization Method | |
Li et al. | A closed-form solution for moving target localization with sphere constraint | |
Liu et al. | An efficient estimator for source localization using TD and AOA measurements in MIMO radar systems | |
CN108445445B (zh) | 一种分布式无源定位方法及装置 | |
Norouzi et al. | Angle of arrival‐based target localisation with low Earth orbit satellite observer | |
Jiexin et al. | Combination of land-based and satellite-based OTH geolocations using differentiable exact penalty method | |
Nuhoglu et al. | Iterative semidefinite relaxation for geolocation of uncooperative radars using Doppler frequency measurements | |
CN113835061B (zh) | 一种信号载波频率先验误差存在下单平台多普勒两阶段闭式定位方法 | |
Zhang et al. | Source localisation using TDOA and FDOA measurements under unknown noise power knowledge | |
Sun et al. | Efficient source positioning method based on two stations using GROA and AOA measurements | |
Congfeng et al. | A Joint TDOA/FDOA Localization Algorithm Using Bi‐iterative Method with Optimal Step Length | |
CN109884582B (zh) | 利用一维测向快速确定目标三维坐标的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |