CN106842118A - 用于分布式多点定位监视***的时差定位方法 - Google Patents
用于分布式多点定位监视***的时差定位方法 Download PDFInfo
- Publication number
- CN106842118A CN106842118A CN201611204604.3A CN201611204604A CN106842118A CN 106842118 A CN106842118 A CN 106842118A CN 201611204604 A CN201611204604 A CN 201611204604A CN 106842118 A CN106842118 A CN 106842118A
- Authority
- CN
- China
- Prior art keywords
- target
- receiving station
- station
- positioning
- main website
- 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.)
- Pending
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/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/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- 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
- 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/16—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using electromagnetic waves other than radio waves
-
- 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/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
- G01S5/22—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Probability & Statistics with Applications (AREA)
- Electromagnetism (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于分布式多点定位监视***的实时定位技术领域,特别涉及一种用于分布式多点定位监视***的时差定位方法。本发明根据GDOP时差定位精度选择定位方法;依据定位方法确定接收站主站,建立观测模型;利用半正定松弛的时差定位算法解算出目标的位置信息,并将定位结果通过终端显示监控***呈现给用户,从而提高机场管制员对机场场面飞机和引导车监控管制能力。本发明的所有计算分析过程都是自发的,最大限度的减少了人为的因素对飞机和引导车定位结果的影响,而且所有的数据处理过程都是并行的,极大地提高了数据处理的效率,使用户执行操作时能最快的获得分析的结果。本发明定位准确,在测量误差较大时仍能对信号源的位置进行较好的估计。
Description
技术领域
本发明属于分布式多点定位监视***的实时定位技术领域,特别涉及一种用于分布式多点定位监视***的时差定位方法。
背景技术
随着我国经济的飞速发展,航班量也在飞速增长,对空管***管制员的指挥和设施设备的维护提出了更高的要求,机场场面监视雷达是一次雷达,主要用于监视机场场面的飞机及车辆,利用目标对电磁波的自身辐射或反射特性发现目标。设备维护人员在维护机场场面监视雷达的过程中包括如下三方面问题:飞机目标***、存在一定的盲区、容易出现假目标。这些问题的存在扰乱了管制员正常的指挥工作,因此,分布式多点定位***的引入必要而且必须,这项技术是机场场面监视雷达的升级和补充,并且能够通过增加接收站或者改变接收站的布局实现***的扩展。
传统的时差定位方法在进行定位精度解算时,考虑到各个目标接收站误差元素之间是相关的,降低了定位精度解算的效率,提高了定位精度解算复杂程度,而且通常使用泰勒级数算法、Chan算法、MDS算法来确定目标的位置信息,使用泰勒级数算法对迭代初始值的选取有要求,会出现局部收敛或者发散的现象;Chan算法和MDS算法在测量误差较大时导致对信号源的位置判断不准确,从而导致时差定位方法经常出现高度非线性的问题,使得定位不准确。
发明内容
本发明为了克服上述现有技术的不足,提供了用于分布式多点定位监视***的时差定位方法,本时差定位方法极大地降低了定位精度解算方法的复杂性,提高了定位精度解算的效率,而且对迭代初始值的选取没有要求,能够收敛到全局极小点,不会出现局部收敛或者发散的现象;在测量误差较大时仍能对信号源的位置进行较好的估计,定位准确。
为实现上述目的,本发明采用了以下技术措施:
用于分布式多点定位监视***的时差定位方法,包括以下步骤:
S1、分布式多点定位监视***中目标接收站主站和接收站辅站获取目标应答信号,对所述目标应答信号经过该分布式多点定位监视***处理得到解析后的测量参数数据;
S2、根据所述目标接收站主站和接收站辅站获取所述测量参数数据,对测量参数数据中相同目标应答码和大于50%的置信度进行配对处理,获取目标在同一时刻发出的测量参数数据;
S3、对配对处理后的测量参数数据进行主站位置选择,得到待定位区域的定位精度;
S4、根据目标的测量参数数据以及待定位区域的定位精度,选择目标的接收站主站,确定时差定位模型;
S5、利用接收站主站和接收站辅站位置信息测量参数数据,通过半正定松弛的时差定位算法解算出目标的位置信息;
S6、将所述目标的位置信息送至终端显示。
优选的,所述测量参数数据包括目标应答码、置信度、目标应答方式、目标应答时间、目标的幅度值。
优选的,步骤S3的具体步骤包括:
S31、确定目标的接收站主站的位置信息;
获取配对处理后的目标的接收站的测量数据,选取其中一个接收站作为目标的接收站主站位置,建立观测模型计算该观测模型条件下定位区域的GDOP值,实现对目标的接收站主站的位置信息确定;通过选取不同的接收站作为目标的接收站主站,分别计算不同观测模型条件下定位区域的GDOP;通过选取不同观测模型条件下GDOP最小值即为时差定位精度最小值,实现对分布式多点定位监视***时差定位模型的构建;
S32、对待定位区域的位置信息求微分;
通过目标的接收站主站的区域覆盖图、目标的接收站辅站的区域覆盖图,获取可定位区域的位置信息,利用所述可定位区域的位置信息以及到达目标的第i个接收站辅站与目标的接收站主站之间的时间差,通过时差定位方程,对待定位区域的位置信息求微分:
c·Δti=c·(ti-t0)=ri-r0(i=1,2)
其中,c为光速、Δti为到达目标的第i个接收站辅站与目标的接收站主站之间时间差、ti为到达目标的第i个接收站辅站的时间、t0为到达目标的接收站主站时间、ri为到达目标的第i个接收站辅站的距离、r0为到达目标的接收站主站距离、i为到达目标的接收站辅站的编号;
S33、求出所述待定位区域内目标与目标的接收站主站位置、目标与目标的接收站辅站位置的相关系数矩阵;
通过所述待定位区域内目标到目标的各接收站辅站与目标到目标的接收站主站的位置差,得出待定位区域内目标与站址位置的相关系数矩阵:
x,y均为目标的位置信息、x0,y0均为目标的接收站主站位置信息、xi、yi均为目标的第i个接收站辅站位置信息、r0为目标到目标的接收站主站距离、ri为目标到目标的第i个接收站辅站距离,
S34、利用伪逆法计算所述待定位区域的误差估计值;
定位误差方程式为:
其中,为定位误差在x方向上定位方差、为定位误差在y方向上定位方差、i为目标的第i个接收站辅站、j为目标的第j个接收站辅站,其中,ηij为Δti与Δtj之间的相关系数,为第i站的时间测量误差的标准差,为第j站的时间测量误差的标准差,σs为站址测量误差的标准差;
S35、求出待定位区域的定位精度;
根据时间测量误差的相关系数ηij,各站址测量误差的标准差σs,利用线性组合求出σij,然后求得x方向上定位方差y方向上定位方差通过公式,计算得到待定位区域的定位精度GDOP。
优选的,步骤S5的具体步骤包括:
S51、构建时差定位方程;
利用接收站主站位置信息测量参数数据建立相应的观测模型,确定的接收站主站的站址信息,利用目标到接收站主站的距离、目标到接收站辅站的距离差构建距离差定位方程:
其中,为到达接收站主站和接收站辅站时间差的测量值,c为光速,ri,k-r1,k为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差,ni1,k为测量误差,i为接收站辅站的索引,k为目标索引;
S52、对信号源位置即待定位的目标位置进行最大似然估计;
利用目标到接收站主站和第i个接收站辅站的距离差测量误差的协方差矩阵表示目标位置的最大似然估计:
其中,J(u)为信号源位置u的高度非线性、非凸代价函数,为所构成的距离差估计组成的向量,f(u)为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差组成的向量,Q为测量误差ni1,k服从均值为零的协方差矩阵,T为矩阵的转置,为使得代价函数J(u)最小时所需信号源位置u的值;
S53、引入辅助矢量,将距离差定位方程转化为约束最小二乘问题;
引入辅助矢量即为目标与接收站主站的位置和距离,转化为约束最小二乘问题,并在此基础上转换为凸的半正定规划问题,将原代价函数整理得到的矩阵如下:
其中,ηk=Bknk,Bk=diag{r2,k,...,rN,k},k=1,2,R=[s2-s1,...,sN-s1]T,O(N-1)×3为(N-1)×3的0矩阵,0N-1为N-1的0行向量,为目标1的到达接收站主站和接收站辅站距离差估计组成的向量,为目标2的到达接收站主站和接收站辅站距离差估计组成的向量;
S54、对所述辅助矢量进行加权最小二乘求解;
对原代价函数整理得到的矩阵进行加权最小二乘求解,得到的辅助矢量估计值如下:
其中,加权矩阵W=E[ηηT]-1=(BQBT)-1,B=diag{B1,B2}表示以B1和B2为对角块的矩阵,Q=E[nnT]为噪声协方差矩阵,n为到达接收站主站和接收站辅站距离差测量误差的矩阵形式, Gy-h=η;
S55、利用加权最小二乘解算的辅助矢量初始估计松弛等式约束,构造新的代价函数;
根据约束条件,构造新的代价函数:
s.t y(6+k)=||y(3(k-1)+1:3k)||2,k=1,2
将(Gy-h)TW(Gy-h)转化为
其中,Y=yyT,
由此可知,Y(6+k,6+k)=trace{Y(3k-2:3k,3k-2:3k)},其中,k=1,2;由Y=yyT,其中y=[(u1-s1)T,(u2-s1)T,r1,1,r1,2]T,可知Y(7,8)=r1,1r1,2;
S56、利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值,并通过特征值分解得到辅助矢量的值;
利用已构造的新的代价函数进行优化求解时,把相应的等式约束转化凸半正定规划优化问题,
待求解的公式为:其约束条件为:
s.t Y(6+k,6+k)=trace{Y(3k-2:3k,3k-2:3k)}
S57、根据求得的辅助矢量的值与目标位置之间的关系获取目标的位置信息;
对Y进行特征值分解可得:
其中,λi为Y的特征值,i=1,...,r,qi为相应的特征向量,设i=1,...,r已按从大到小排列,即λ1≥λ2≥...≥λr>0,根据Y的定义可得信号源位置估计为:
本发明的有益效果在于:
1)、本发明根据GDOP时差定位精度选择定位方法;依据定位方法确定接收站主站,建立观测模型;利用半正定松弛的时差定位算法解算出目标的位置信息,并将定位结果通过终端显示监控***呈现给用户,从而提高机场管制员对机场场面飞机和引导车监控管制能力。本发明的所有计算分析过程都是自发的,最大限度的减少了人为的因素对飞机和引导车定位结果的影响,而且所有的数据处理过程都是并行的,极大地提高了数据处理的效率,使用户执行操作时能最快的获得分析的结果。
2)、确定待定位区域的定位精度是利用多个目标的接收站的位置信息和内部噪声、环境噪声和时间测量的量化噪声信息对目标的定位估计值进行定位精度解算,通过不同的组合和不同的布站方式,得到测量目标区域覆盖的定位精度GDOP,从而与不同布站方式下的GDOP比对,选择最优布站方式,为后续时差定位解算模型的构建提供依据。
3)、本发明将高度非线性的时差定位问题转化为约束加权最小二乘估计问题,然后在此基础上通过半正定松弛转化为凸的半正定规划问题,进行优化求解,避免了传统迭代算法中出现的局部收敛和求解发散的问题,大大地提高了定位的精度。
附图说明
图1为本发明的流程图;
图2为本发明的步骤S3的工作流程图;
图3为本发明的步骤S5的工作流程图;
图4为本发明的Y形布站仿真图;
图5为本发明的倒三角形布站仿真图;
图6为本发明的菱形布站仿真图;
图7为利用本半正定松弛时差定位方法、Chan算法、泰勒级数法对信号源u1和u2位置估计随变化时的均方根误差的统计结果示意图;
图8为利用本半正定松弛时差定位方法、Chan算法、MDS算法、泰勒级数法对信号源u3位置估计随变化时的均方根误差的统计结果示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,用于分布式多点定位监视***的时差定位方法,包括以下步骤:
S1、分布式多点定位监视***中目标接收站主站和接收站辅站获取目标应答信号,对所述目标应答信号经过该分布式多点定位监视***处理得到解析后的测量参数数据;
S2、根据所述目标接收站主站和接收站辅站获取所述测量参数数据,对测量参数数据中相同目标应答码和大于50%的置信度进行配对处理,获取目标在同一时刻发出的测量参数数据;
S3、对配对处理后的测量参数数据进行主站位置选择,得到待定位区域的定位精度;
S4、根据目标的测量参数数据以及待定位区域的定位精度,选择目标的接收站主站,确定时差定位模型;
S5、利用接收站主站和接收站辅站位置信息测量参数数据,通过半正定松弛的时差定位算法解算出目标的位置信息;
S6、将所述目标的位置信息送至终端显示。
所述目标即为待定位目标。
所述测量参数数据包括目标应答码、置信度、目标应答方式、目标应答时间、目标的幅度值。
如图2所示,对配对处理后的测量参数数据进行主站位置选择,得到待定位区域的定位精度的具体操作步骤包括:
S31、确定目标的接收站主站的位置信息;
获取配对处理后的目标的接收站的测量数据,选取其中一个接收站作为目标的接收站主站位置,建立观测模型计算该观测模型条件下定位区域的GDOP值,实现对目标的接收站主站的位置信息确定;通过选取不同的接收站作为目标的接收站主站,分别计算不同观测模型条件下定位区域的GDOP;通过选取不同观测模型条件下GDOP最小值即为时差定位精度最小值,实现对分布式多点定位监视***时差定位模型的构建;
S32、对待定位区域的位置信息求微分;
通过目标的接收站主站的区域覆盖图、目标的接收站辅站的区域覆盖图,获取可定位区域的位置信息,利用所述可定位区域的位置信息以及到达目标的第i个接收站辅站与目标的接收站主站之间的时间差,通过时差定位方程,对待定位区域的位置信息求微分:
c·Δti=c·(ti-t0)=ri-r0(i=1,2)
其中,c为光速、Δti为到达目标的第i个接收站辅站与目标的接收站主站之间时间差、ti为到达目标的第i个接收站辅站的时间、t0为到达目标的接收站主站时间、ri为到达目标的第i个接收站辅站的距离、r0为到达目标的接收站主站距离、i为到达目标的接收站辅站的编号;
S33、求出所述待定位区域内目标与目标的接收站主站位置、目标与目标的接收站辅站位置的相关系数矩阵;
通过所述待定位区域内目标到目标的各接收站辅站与目标到目标的接收站主站的位置差,得出待定位区域内目标与站址位置的相关系数矩阵:
x,y均为目标的位置信息、x0,y0均为目标的接收站主站位置信息、xi、yi均为目标的第i个接收站辅站位置信息、r0为目标到目标的接收站主站距离、ri为目标到目标的第i个接收站辅站距离,
S34、利用伪逆法计算所述待定位区域的误差估计值;
定位误差方程式为:
其中,为定位误差在x方向上定位方差、为定位误差在y方向上定位方差、i为目标的第i个接收站辅站、j为目标的第j个接收站辅站,其中,ηij为Δti与Δtj之间的相关系数,为第i站的时间测量误差的标准差,为第j站的时间测量误差的标准差,σs为站址测量误差的标准差;
S35、求出待定位区域的定位精度;
根据时间测量误差的相关系数ηij,各站址测量误差的标准差σs,利用线性组合求出σij,然后求得x方向上定位方差y方向上定位方差通过公式,计算得到待定位区域的定位精度GDOP。
当主站的坐标取(0,0)时,三个辅站的坐标分别取不同的值时,如表1所示,形成的仿真图分别为Y形布站、倒三角布站、菱形布站。
表1:
从图4~图6可知,对于同一种布站形式来说,定位精度GDOP随着目标和接收站的分布区域之间的距离的增加而降低,对于不同的布站方式,全方位内定位误差有所不同,用于分布式多点定位监视***的定位精度方法,可以充分考虑Y形、倒三角形、菱形布站方法对不同区域定位精度的不同,来进行目标的接收站主站和辅站位置的选择。力求目标在不同的区域里,选择最优布站方案,实现目标的高精度定位。
如图3所示,利用接收站主站和接收站辅站位置信息测量参数数据,通过半正定松弛的时差定位算法解算出目标的位置信息的具体操作步骤为:
S51、构建时差定位方程;
利用接收站主站位置信息测量参数数据建立相应的观测模型,确定的接收站主站的站址信息,利用目标到接收站主站的距离、目标到接收站辅站的距离差构建距离差定位方程:
其中,为到达接收站主站和接收站辅站时间差的测量值,c为光速,ri,k-r1,k为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差,ni1,k为测量误差,i为接收站辅站的索引,k为目标索引。
S52、对信号源位置即待定位的目标位置进行最大似然估计;
利用目标到接收站主站和第i个接收站辅站的距离差测量误差的协方差矩阵表示目标位置的最大似然估计:
其中,J(u)为信号源位置u的高度非线性、非凸代价函数,为所构成的距离差估计组成的向量,f(u)为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差组成的向量,Q为测量误差ni1,k服从均值为零的协方差矩阵,T为矩阵的转置,为使得代价函数J(u)最小时所需信号源位置u的值;
S53、引入辅助矢量,将距离差定位方程转化为约束最小二乘问题;
引入辅助矢量即为目标与接收站主站的位置和距离,转化为约束最小二乘问题,并在此基础上转换为凸的半正定规划问题,将原代价函数整理得到的矩阵如下:
其中,ηk=Bknk,Bk=diag{r2,k,…,rN,k},k=1,2,R=[s2-s1,...,sN-s1]T,O(N-1)×3为(N-1)×3的0矩阵,0N-1为N-1的0行向量,为目标1的到达接收站主站和接收站辅站距离差估计组成的向量,为目标2的到达接收站主站和接收站辅站距离差估计组成的向量;
S54、对所述辅助矢量进行加权最小二乘求解;
对原代价函数整理得到的矩阵进行加权最小二乘求解,得到的辅助矢量估计值如下:
其中,加权矩阵W=E[ηηT]-1=(BQBT)-1,B=diag{B1,B2}表示以B1和B2为对角块的矩阵,Q=E[nnT]为噪声协方差矩阵,n为到达接收站主站和接收站辅站距离差测量误差的矩阵形式, Gy-h=η,y为辅助矢量;
S55、利用加权最小二乘解算的辅助矢量初始估计松弛等式约束,构造新的代价函数;
根据约束条件,构造新的代价函数:
s.t y(6+k)=||y(3(k-1)+1:3k)||2,k=1,2
将(Gy-h)TW(Gy-h)转化为
其中,Y=yyT,
由此可知,Y(6+k,6+k)=trace{Y(3k-2:3k,3k-2:3k)},其中,k=1,2;由Y=yyT,其中y=[(u1-s1)T,(u2-s1)T,r1,1,r1,2]T,可知Y(7,8)=r1,1r1,2;
S56、利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值,并通过特征值分解得到辅助矢量的值;
利用已构造的新的代价函数进行优化求解时,把相应的等式约束转化凸半正定规划优化问题,
待求解的公式为:其约束条件为:
s.t Y(6+k,6+k)=trace{Y(3k-2:3k,3k-2:3k)}
S57、根据求得的辅助矢量的值与目标位置之间的关系获取目标的位置信息;
对Y进行特征值分解可得:
其中,λi为Y的特征值,i=1,...,r,qi为相应的特征向量,设i=1,...,r已按从大到小排列,即λ1≥λ2≥...≥λr>0,根据Y的定义可得信号源位置估计为:
为了检验本半正定松弛时差定位方法对信号源位置估计的性能,我们将本半正定松弛时差定位方法同传统的Chan算法、泰勒级数法、MDS算法及克拉美罗界的仿真结果进行比较,5个接收站位置坐标分别为:s1=[300,100,150]T,s2=[400,150,100]T,s3=[300,500,200]T,s4=[350,200,100]T,s5=[-100,-100,-100]T,信号源坐标单位均为米。假设各个信号源的TDOA测量值服从均值为零、方差为的高斯分布,则Qk满足对角线元素为非对角线元素为为了对比的公平,利用对所述辅助矢量进行加权最小二乘求解的具体步骤中的所得初始定位结果计算利用加权最小二乘解算的辅助矢量初始估计松弛等式约束,构造新的代价函数的具体步骤中的F和Chan算法、MDS算法的权值以及作为泰勒级数法的迭代初始值。三个信号源位置坐标分别为u1=[314,483,209]T,u2=[600,650,550]T和u3=[285,325,275]T。
如图7所示,分别给出了本半正定松弛时差定位方法、Chan算法及泰勒级数法对信号源u1和u2位置估计随变化时的均方根误差的统计结果,较低的一组曲线对应信号源u1,较高的一组对应u2。可以看出,在噪声功率较小时,三种方法对两信号源位置的估计均方根误差均十分接近克拉美罗界。随着测量误差的增大,各算法的均方根误差均有所增加;其中泰勒级数法最先偏离克拉美罗界,这是因为在测量误差较大时,由步骤S4求得的初始值偏离真实值较远,泰勒级数法在进行迭代求解时,很容易陷入局部极小点或者发散;而本半正定松弛时差定位方法对两信号源位置的估计均方根误差均明显小于Chan算法,具有更为准确的定位性能。
如图8所示,分别给出了本半正定松弛时差定位方法、Chan算法、MDS算法及泰勒级数法对信号源u3位置估计随变化时的均方根误差的统计结果。可以看出,在噪声功率较小时,几种方法对信号源位置估计的均方根误差均十分接近克拉美罗界。随着噪声功率的增加,各算法的估计均方根误差都会有所增加。在测量误差较大时,泰勒级数法的估计均方根误差出现急剧增加的现象;MDS算法偏离克拉美罗界要明显晚于Chan算法,但是在测量误差较大时,出现MDS算法的均方根误差大于Chan算法的现象;而本半正定松弛时差定位方法的估计均方根误差一直保持最小。
本发明将定位结果通过终端显示监控***呈现给用户,从而提高机场管制员对机场场面飞机和引导车监控管制能力。本发明的所有计算分析过程都是自发的,最大限度的减少了人为的因素对飞机和引导车定位结果的影响,而且所有的数据处理过程都是并行的,极大地提高了数据处理的效率,使用户执行操作时能最快的获得分析的结果,本发明定位准确,在测量误差较大时仍能对信号源的位置进行较好的估计。
Claims (4)
1.用于分布式多点定位监视***的时差定位方法,其特征在于,包括以下步骤:
S1、分布式多点定位监视***中目标接收站主站和接收站辅站获取目标应答信号,对所述目标应答信号经过该分布式多点定位监视***处理得到解析后的测量参数数据;
S2、根据所述目标接收站主站和接收站辅站获取所述测量参数数据,对测量参数数据中相同目标应答码和大于50%的置信度进行配对处理,获取目标在同一时刻发出的测量参数数据;
S3、对配对处理后的测量参数数据进行主站位置选择,得到待定位区域的定位精度;
S4、根据目标的测量参数数据以及待定位区域的定位精度,选择目标的接收站主站,确定时差定位模型;
S5、利用接收站主站和接收站辅站位置信息测量参数数据,通过半正定松弛的时差定位算法解算出目标的位置信息;
S6、将所述目标的位置信息送至终端显示。
2.如权利要求1所述的用于分布式多点定位监视***的时差定位方法,其特征在于:所述测量参数数据包括目标应答码、置信度、目标应答方式、目标应答时间、目标的幅度值。
3.如权利要求1所述的用于分布式多点定位监视***的时差定位方法,其特征在于,步骤S3的具体步骤包括:
S31、确定目标的接收站主站的位置信息;
获取配对处理后的目标的接收站的测量数据,选取其中一个接收站作为目标的接收站主站位置,建立观测模型计算该观测模型条件下定位区域的GDOP值,实现对目标的接收站主站的位置信息确定;通过选取不同的接收站作为目标的接收站主站,分别计算不同观测模型条件下定位区域的GDOP;通过选取不同观测模型条件下GDOP最小值即为时差定位精度最小值,实现对分布式多点定位监视***时差定位模型的构建;
S32、对待定位区域的位置信息求微分;
通过目标的接收站主站的区域覆盖图、目标的接收站辅站的区域覆盖图,获取可定位区域的位置信息,利用所述可定位区域的位置信息以及到达目标的第i个接收站辅站与目标的接收站主站之间的时间差,通过时差定位方程,对待定位区域的位置信息求微分:
c·Δti=c·(ti-t0)=ri-r0(i=1,2)
其中,c为光速、Δti为到达目标的第i个接收站辅站与目标的接收站主站之间时间差、ti为到达目标的第i个接收站辅站的时间、t0为到达目标的接收站主站时间、ri为到达目标的第i个接收站辅站的距离、r0为到达目标的接收站主站距离、i为到达目标的接收站辅站的编号;
S33、求出所述待定位区域内目标与目标的接收站主站位置、目标与目标的接收站辅站位置的相关系数矩阵;
通过所述待定位区域内目标到目标的各接收站辅站与目标到目标的接收站主站的位置差,得出待定位区域内目标与站址位置的相关系数矩阵:
x,y均为目标的位置信息、x0,y0均为目标的接收站主站位置信息、xi、yi均为目标的第i个接收站辅站位置信息、r0为目标到目标的接收站主站距离、ri为目标到目标的第i个接收站辅站距离,
S34、利用伪逆法计算所述待定位区域的误差估计值;
定位误差方程式为:
其中,为定位误差在x方向上定位方差、为定位误差在y方向上定位方差、i为目标的第i个接收站辅站、j为目标的第j个接收站辅站,其中,ηij为Δti与Δtj之间的相关系数,为第i站的时间测量误差的标准差,为第j站的时间测量误差的标准差,σs为站址测量误差的标准差;
S35、求出待定位区域的定位精度;
根据时间测量误差的相关系数ηij,各站址测量误差的标准差σs,利用线性组合求出σij,然后求得x方向上定位方差y方向上定位方差通过公式,计算得到待定位区域的定位精度GDOP。
4.如权利要求3所述的用于分布式多点定位监视***的时差定位方法,其特征在于,步骤S5的具体步骤包括:
S51、构建时差定位方程;
利用接收站主站位置信息测量参数数据建立相应的观测模型,确定的接收站主站的站址信息,利用目标到接收站主站的距离、目标到接收站辅站的距离差构建距离差定位方程:
其中,为到达接收站主站和接收站辅站时间差的测量值,c为光速,ri,k-r1,k为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差,ni1,k为测量误差,i为接收站辅站的索引,k为目标索引;
S52、对信号源位置即待定位的目标位置进行最大似然估计;
利用目标到接收站主站和第i个接收站辅站的距离差测量误差的协方差矩阵表示目标位置的最大似然估计:
其中,J(u)为信号源位置u的高度非线性、非凸代价函数,为所构成的距离差估计组成的向量,f(u)为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差组成的向量,Q为测量误差ni1,k服从均值为零的协方差矩阵,T为矩阵的转置,为使得代价函数J(u)最小时所需信号源位置u的值;
S53、引入辅助矢量,将距离差定位方程转化为约束最小二乘问题;
引入辅助矢量即为目标与接收站主站的位置和距离,转化为约束最小二乘问题,并在此基础上转换为凸的半正定规划问题,将原代价函数整理得到的矩阵如下:
其中,ηk=Bknk,Bk=diag{r2,k,...,rN,k},k=1,2,R=[s2-s1,...,sN-s1]T,O(N-1)×3为(N-1)×3的0矩阵,0N-1为N-1的0行向量,为目标1的到达接收站主站和接收站辅站距离差估计组成的向量,为目标2的到达接收站主站和接收站辅站距离差估计组成的向量;
S54、对所述辅助矢量进行加权最小二乘求解;
对原代价函数整理得到的矩阵进行加权最小二乘求解,得到的辅助矢量估计值如下:
其中,加权矩阵W=E[ηηT]-1=(BQBT)-1,B=diag{B1,B2}表示以B1和B2为对角块的矩阵,Q=E[nnT]为噪声协方差矩阵,n为到达接收站主站和接收站辅站距离差测量误差的矩阵形式, Gy-h=η;
S55、利用加权最小二乘解算的辅助矢量初始估计松弛等式约束,构造新的代价函数;
根据约束条件,构造新的代价函数:
s.t y(6+k)=||y(3(k-1)+1:3k)||2,k=1,2
将(Gy-h)TW(Gy-h)转化为
其中,Y=yyT,
由此可知,Y(6+k,6+k)=trace{Y(3k-2:3k,3k-2:3k)},其中,k=1,2;由Y=yyT,其中y=[(u1-s1)T,(u2-s1)T,r1,1,r1,2]T,可知Y(7,8)=r1,1r1,2;
S56、利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值,并通过特征值分解得到辅助矢量的值;
利用已构造的新的代价函数进行优化求解时,把相应的等式约束转化凸半正定规划优化问题,
待求解的公式为:其约束条件为:
s.t Y(6+k,6+k)=trace{Y(3k-2:3k,3k-2:3k)}
S57、根据求得的辅助矢量的值与目标位置之间的关系获取目标的位置信息;
对Y进行特征值分解可得:
其中,λi为Y的特征值,i=1,...,r,qi为相应的特征向量,设i=1,...,r已按从大到小排列,即λ1≥λ2≥...≥λr>0,根据Y的定义可得信号源位置估计为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611204604.3A CN106842118A (zh) | 2016-12-23 | 2016-12-23 | 用于分布式多点定位监视***的时差定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611204604.3A CN106842118A (zh) | 2016-12-23 | 2016-12-23 | 用于分布式多点定位监视***的时差定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106842118A true CN106842118A (zh) | 2017-06-13 |
Family
ID=59135961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611204604.3A Pending CN106842118A (zh) | 2016-12-23 | 2016-12-23 | 用于分布式多点定位监视***的时差定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106842118A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107422340A (zh) * | 2017-09-05 | 2017-12-01 | 芜湖华创光电科技有限公司 | 一种多点定位***接收站位置误差校正方法 |
CN107504630A (zh) * | 2017-08-02 | 2017-12-22 | 广东美的制冷设备有限公司 | 空调器及其控制方法和装置 |
CN109633724A (zh) * | 2019-01-16 | 2019-04-16 | 电子科技大学 | 基于单星与多地面站联合测量的无源目标定位方法 |
CN109948791A (zh) * | 2017-12-21 | 2019-06-28 | 河北科技大学 | 利用遗传算法优化bp神经网络的方法及其在定位中的应用 |
CN110261820A (zh) * | 2019-07-18 | 2019-09-20 | 中电科仪器仪表有限公司 | 一种多测量站的时差定位方法及装置 |
CN111194083A (zh) * | 2020-04-09 | 2020-05-22 | 成都信息工程大学 | 一种无线电定位***及其定位方法 |
CN111273274A (zh) * | 2020-03-12 | 2020-06-12 | 四川九洲电器集团有限责任公司 | 多基协同定位方法、存储介质、雷达及雷达定位*** |
CN111505573A (zh) * | 2020-03-18 | 2020-08-07 | 中国民用航空总局第二研究所 | 一种分布式定位***的航迹生成方法及装置 |
CN112327251A (zh) * | 2020-11-05 | 2021-02-05 | 中国人民解放军32802部队 | 低速率通信条件下的无人值守式信号时差定位方法 |
CN112394372A (zh) * | 2020-11-06 | 2021-02-23 | 四川九洲空管科技有限责任公司 | 一种基于ads-b记录数据评估多点定位性能的方法及*** |
CN113050137A (zh) * | 2021-03-09 | 2021-06-29 | 江西师范大学 | 一种多点协测的空间信息获取方法 |
RU2790348C1 (ru) * | 2021-12-14 | 2023-02-16 | Общество с ограниченной ответственностью "Специальный Технологический Центр" | Разностно-дальномерный способ определения координат источника радиоизлучения |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102004235A (zh) * | 2010-10-22 | 2011-04-06 | 民航数据通信有限责任公司 | 多点定位***的接收站选取方法 |
CN201945685U (zh) * | 2010-11-04 | 2011-08-24 | 中国民用航空总局第二研究所 | 高精度分布式脉冲信号到达时间差检测*** |
WO2010138696A9 (en) * | 2009-05-27 | 2014-11-27 | Sensis Corporation | System and method for passive range-aided multilateration using time lag of arrival (tloa) measurements |
CN104640206A (zh) * | 2015-02-15 | 2015-05-20 | 中国民航大学 | 一种偏心式星形布站体制的多点定位方法 |
CN104833953A (zh) * | 2015-05-11 | 2015-08-12 | 中国民用航空总局第二研究所 | 一种机场非视距信道环境下的多点定位监视***和方法 |
EP2548041B1 (en) * | 2010-03-17 | 2016-07-27 | Honeywell International Inc. | Systems and methods for short baseline, low cost determination of airborne aircraft location |
-
2016
- 2016-12-23 CN CN201611204604.3A patent/CN106842118A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010138696A9 (en) * | 2009-05-27 | 2014-11-27 | Sensis Corporation | System and method for passive range-aided multilateration using time lag of arrival (tloa) measurements |
EP2548041B1 (en) * | 2010-03-17 | 2016-07-27 | Honeywell International Inc. | Systems and methods for short baseline, low cost determination of airborne aircraft location |
CN102004235A (zh) * | 2010-10-22 | 2011-04-06 | 民航数据通信有限责任公司 | 多点定位***的接收站选取方法 |
CN201945685U (zh) * | 2010-11-04 | 2011-08-24 | 中国民用航空总局第二研究所 | 高精度分布式脉冲信号到达时间差检测*** |
CN104640206A (zh) * | 2015-02-15 | 2015-05-20 | 中国民航大学 | 一种偏心式星形布站体制的多点定位方法 |
CN104833953A (zh) * | 2015-05-11 | 2015-08-12 | 中国民用航空总局第二研究所 | 一种机场非视距信道环境下的多点定位监视***和方法 |
Non-Patent Citations (3)
Title |
---|
刘刚: "分布式多站无源时差定位***研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
周建红: "基于MDS***的定位算法与同步技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
杨琳等: "一种利用半正定松弛的TDOA定位算法", 《万方》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107504630A (zh) * | 2017-08-02 | 2017-12-22 | 广东美的制冷设备有限公司 | 空调器及其控制方法和装置 |
CN107422340B (zh) * | 2017-09-05 | 2020-05-19 | 芜湖华创光电科技有限公司 | 一种多点定位***接收站位置误差校正方法 |
CN107422340A (zh) * | 2017-09-05 | 2017-12-01 | 芜湖华创光电科技有限公司 | 一种多点定位***接收站位置误差校正方法 |
CN109948791A (zh) * | 2017-12-21 | 2019-06-28 | 河北科技大学 | 利用遗传算法优化bp神经网络的方法及其在定位中的应用 |
CN109948791B (zh) * | 2017-12-21 | 2022-03-29 | 河北科技大学 | 利用遗传算法优化bp神经网络的方法及其在定位中的应用 |
CN109633724A (zh) * | 2019-01-16 | 2019-04-16 | 电子科技大学 | 基于单星与多地面站联合测量的无源目标定位方法 |
CN109633724B (zh) * | 2019-01-16 | 2023-03-03 | 电子科技大学 | 基于单星与多地面站联合测量的无源目标定位方法 |
CN110261820A (zh) * | 2019-07-18 | 2019-09-20 | 中电科仪器仪表有限公司 | 一种多测量站的时差定位方法及装置 |
CN111273274B (zh) * | 2020-03-12 | 2022-03-18 | 四川九洲电器集团有限责任公司 | 多基协同定位方法、存储介质、雷达及雷达定位*** |
CN111273274A (zh) * | 2020-03-12 | 2020-06-12 | 四川九洲电器集团有限责任公司 | 多基协同定位方法、存储介质、雷达及雷达定位*** |
CN111505573A (zh) * | 2020-03-18 | 2020-08-07 | 中国民用航空总局第二研究所 | 一种分布式定位***的航迹生成方法及装置 |
CN111194083B (zh) * | 2020-04-09 | 2020-07-07 | 成都信息工程大学 | 一种无线电定位***及其定位方法 |
CN111194083A (zh) * | 2020-04-09 | 2020-05-22 | 成都信息工程大学 | 一种无线电定位***及其定位方法 |
CN112327251A (zh) * | 2020-11-05 | 2021-02-05 | 中国人民解放军32802部队 | 低速率通信条件下的无人值守式信号时差定位方法 |
CN112394372A (zh) * | 2020-11-06 | 2021-02-23 | 四川九洲空管科技有限责任公司 | 一种基于ads-b记录数据评估多点定位性能的方法及*** |
CN112394372B (zh) * | 2020-11-06 | 2024-05-28 | 四川九洲空管科技有限责任公司 | 一种基于ads-b记录数据评估多点定位性能的方法及*** |
CN113050137A (zh) * | 2021-03-09 | 2021-06-29 | 江西师范大学 | 一种多点协测的空间信息获取方法 |
CN113050137B (zh) * | 2021-03-09 | 2022-04-26 | 江西师范大学 | 一种多点协测的空间信息获取方法 |
RU2790348C1 (ru) * | 2021-12-14 | 2023-02-16 | Общество с ограниченной ответственностью "Специальный Технологический Центр" | Разностно-дальномерный способ определения координат источника радиоизлучения |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842118A (zh) | 用于分布式多点定位监视***的时差定位方法 | |
CN106814357A (zh) | 用于分布式多点定位监视***的半正定松弛时差定位方法 | |
CN103945332B (zh) | 一种接收信号强度和多径信息联合神经网络室内定位方法 | |
CN105740203B (zh) | 多传感器无源协同测向定位方法 | |
CN106610483B (zh) | 基于张量子空间和谱峰搜索的mimo雷达角度估计算法 | |
CN103941233B (zh) | 一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法 | |
CN110174643A (zh) | 一种无需噪声功率信息的基于到达时间差的定位方法 | |
CN106872942A (zh) | 用于分布式多点定位监视***的定位精度解算方法 | |
CN104301999B (zh) | 一种基于rssi的无线传感器网络自适应迭代定位方法 | |
CN104502911B (zh) | 一种穿墙成像雷达的墙壁参数估计方法 | |
CN106353720B (zh) | 基于tdoa/groa的多站连续定位模型 | |
CN106612495B (zh) | 一种基于传播损耗学习的室内定位方法及*** | |
CN104619020A (zh) | 基于rssi和toa测距的wifi室内定位方法 | |
CN106706133B (zh) | 一种点斑状目标姿态估计方法及*** | |
CN112748397A (zh) | 一种非视距条件下基于自适应bp神经网络的uwb定位方法 | |
CN103596267A (zh) | 一种基于欧氏距离的指纹图匹配方法 | |
CN106941718A (zh) | 一种基于信号子空间指纹库的混合室内定位方法 | |
CN106443655B (zh) | 一种多输入多输出雷达近场定位算法 | |
CN103885030A (zh) | 一种无线传感器网络移动节点的定位方法 | |
CN104504728B (zh) | 多机动目标跟踪方法、***及其广义联合概率数据关联器 | |
CN109188355A (zh) | 一种多点定位***接收天线优化及最优布站方法 | |
TWI427313B (zh) | Rfid標籤定位演算方法 | |
CN104661232B (zh) | 基于费歇尔信息矩阵指纹定位精度基本限的ap布置方法 | |
CN110515037A (zh) | 非视距环境下能时频多域联合的被动定位方法 | |
CN108668358A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170613 |
|
RJ01 | Rejection of invention patent application after publication |