CN111308418B - 一种对高度未知的目标进行二维定位的稳健方法 - Google Patents
一种对高度未知的目标进行二维定位的稳健方法 Download PDFInfo
- Publication number
- CN111308418B CN111308418B CN202010162061.3A CN202010162061A CN111308418B CN 111308418 B CN111308418 B CN 111308418B CN 202010162061 A CN202010162061 A CN 202010162061A CN 111308418 B CN111308418 B CN 111308418B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- determining
- target
- matrix
- height
- 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/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/12—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 by co-ordinating position lines of different shape, e.g. hyperbolic, circular, elliptical or radial
-
- 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
-
- 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/28—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 by co-ordinating position lines of different shape, e.g. hyperbolic, circular, elliptical or radial
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于电子信息技术领域,具体涉及一种对高度未知的目标进行二维定位的稳健方法。本发明的目的在于针对多站时差定位体制下观测站不足导致原定位系统失效,以及目标的实际高度与假设高度之间的误差导致定位精度下降的问题,使用与目标具有一定高度差的、不少于三个的观测站的时差测量,通过对目标位置坐标进行三维网格搜索并构建三阶张量,利用张量分解的方法进行降维,从而稳健的确定高度未知的目标的二维坐标。
Description
技术领域
本发明属于电子信息技术领域,具体涉及一种利用多个观测站的时差测量对高度未知目标进行二维定位的稳健方法。
背景技术
时差定位技术利用目标信号到达多个观测站的时间差来实现对目标的定位。在二维平面,两个观测站之间的时间差可以确定一组双曲线,故至少需要三个观测站才可能对二维平面上的目标进行二维定位;在三维空间,两个观测站之间的时间差可以确定一组双曲面,故至少需要四个观测站测才可能对三维空间里的目标进行三维定位。在实际的时差定位场景中,由于观测站往往同目标不在一个平面,利用四个及以上的观测站对目标进行三维定位是必要的。但是,当一部分观测站发生故障、被遮挡或测量误差过大导致实际可用的观测站数目不满足要求时,原先的定位系统便会失效。因此,如何在较少的时差测量下仍能准确地确定目标的坐标参数具有重要意义。
当目标高度已知时,可以利用三个观测站的时差测量对三维空间中的目标进行二维定位,例如,基于WGS_84椭球模型的三星时差定位通常假设目标位于地球表面。当目标不是位于地球表面时,只要目标的实际高度相对于卫星高度足够小而且与假设高度之间的误差足够小,则零高度假设带来的定位误差就可以忽略不计。但是,当目标的实际高度与假设高度之间的误差越来越大时,对目标进行二维定位的误差也会越来越差。由于目标的实际高度与假设高度之间的误差是未知的,对目标进行二维定位的性能好坏也是未知的,这给实际应用带来了不确定性。
为此,有必要发展对高度未知的目标进行二维定位的稳健方法。
发明内容
本发明的目的在于针对多站时差定位体制下观测站不足导致原定位系统失效,以及目标的实际高度与假设高度之间的误差导致定位精度下降的问题,使用与目标具有一定高度差的、不少于三个的观测站的时差测量,通过对目标位置坐标进行三维网格搜索并构建三阶张量,利用张量分解的方法进行降维,从而稳健的确定高度未知的目标的二维坐标。
为实现上述技术目的,本发明采用如下技术方案进行实现:
步骤1:初始化,设置观测站数目及对目标水平面二维及高度一维的搜索个数,将观测站的位置坐标、目标信号到达多个观测站的时差测量、对目标水平面二维及高度一维的搜索个数、三维网格点数据写入主机内存;
步骤2:由观测站的位置坐标和时差测量确定每个三维网格点的三维代价函数值,从而得到三阶张量;
步骤3:对三阶张量进行降维处理,并由降维的三阶张量确定与二维位置坐标相关的矩阵;
步骤4:确定与二维位置坐标相关的矩阵的最小元素对应的二维网格点为目标的二维坐标估计。
进一步地,所述步骤1的具体执行过程如下:
初始化,设置观测站数目为N,目标水平面二维及高度一维的搜索个数分别为Nx、Ny、Nz,将观测站的位置坐标(αi,βi,γi),i=1,…,N、目标信号到达多个观测站的时差测量ti1,i=2,…,N、三维网格点nx=1,2,...,Nx、ny=1,2,...,Ny、nz=1,2,...,Nz,写入主机内存;
进一步地,所述步骤2的具体执行过程如下:
其中t=[t21 … tN1]T为时差测量值,T表示转置操作,nx=1,2,...,Nx、ny=1,2,...,Ny、nz=1,2,...,Nz,
进一步地,所述步骤3的具体执行过程如下:
确定门限
从而确定降维的三阶张量为
由降维的三阶张量F′确定与二维位置坐标相关的矩阵为
进一步地,所述步骤4的具体执行过程如下:
与现有技术相比,本发明方案的有益效果是:
本发明利用不少于三个观测站的时差测量和观测站的位置坐标,即可在目标高度未知的情况下确定目标的平面二维坐标,一方面可实现在三维空间中利用更少的时差测量对高度未知的目标进行平面二维定位,另一方面可降低目标假设高度与目标实际高度之间的误差带来的定位精度的恶化。
具体实施方式
下面结合实施例,对本发明作进一步的详细描述:
实施例
本例中的观测站数目为N=3,其位置坐标分别为(300,0,1402)、(-300,0,1509)、(0,300,1829),其单位均为米,目标位置坐标为(1224.7,866.6,208.3)米。目标信号到达多个观测站的时差真实值为[1.4750,1.2265]T微秒。
本发明的具体实施方式流程如下:
步骤1:初始化设置观测站数目为N=3,观测站均距离地面有一定高度,将观测站的位置坐标(α1,β1,γ1)=(300,0,1402)、(α2,β2,γ2)=(-300,0,1509)、(α3,β3,γ3)=(0,300,1829)(单位:米)以及目标信号到达多个观测站的时差测量t21=1.4632、t31=1.2166(单位:微妙)等数据写入主机内存;
步骤2:对目标位置坐标范围x∈[0,2000]、y∈[0,2000]、z∈[0,1000](单位:米)进行三维网格点划分,得到三维网格点其中,nx=1,2,...,Nx、ny=1,2,...,Ny、nz=1,2,...,Nz,Nx=41、Ny=41、Nz=11分别表示对目标水平面二维及高度一维的搜索个数。对于网格点计算三维代价函数值为
其中t=[t21 … tN1]T为时差测量值,T表示转置操作,
c表示光速,nx=1,2,...,Nx、ny=1,2,...,Ny、nz=1,2,...,Nz。
步骤3:确定门限
λ1=1.2228×10-8、λ2=2.2281×10-10、λ3=1.7062×10-10、λ4=7.6722×10-12
从而确定
令F和U(1)做模-1乘积得到
从而确定
从而确定
其中
从而确定降维的三阶张量为
由降维的三阶张量F′确定与二维位置坐标相关的矩阵为
定义定位误差为目标的二维坐标位置估计与目标的实际二维坐标之间的距离。在本实例中,在时差测量误差分别为11.7692纳秒、9.9053纳秒时,本发明确定的目标二维坐标为(1250,850)米,可见本发明的定位误差为30.3047米。然而,在目标的假设高度分别为0米、200米、400米、600米时,确定的目标二维位置分别为(1388.8,992.2)、(1216.7,880.9)、(1045.3,769.1)、(875.2,656.5)(单位:米),定位误差分别为206.5857米、16.3402米、204.1483米、407.7642米,由于目标的实际高度为208.3米,可见,随着目标高度假设误差增大,定位误差也随之增加。因此,本发明不要求已知目标的实际高度,定位误差明显小于目标的假设高度为0米、400米、600米时的定位误差,可见本发明具有定位的稳健性。
Claims (2)
1.一种对高度未知的目标进行二维定位的稳健方法,其特征在于,包括以下步骤:
S1、设置观测站数目为N,目标水平面二维及高度一维的搜索个数分别为Nx、Ny、Nz,定义观测站的位置坐标为(αi,βi,γi),i=1,…,N、目标信号到达多个观测站的时差测量为ti1,i=2,…,N、三维网格点为
S2、由观测站的位置坐标和时差测量确定每个三维网格点的三维代价函数值,从而得到三阶张量;具体方法为:
其中t=[t21…tN1]T为时差测量值,T表示转置操作,nx=1,2,...,Nx、ny=1,2,...,Ny、nz=1,2,...,Nz,
c表示光速;
S3、对三阶张量进行降维处理,并由降维的三阶张量确定与二维位置坐标相关的矩阵;具体方法为:
设置门限
从而确定降维的三阶张量为
由降维的三阶张量F′确定与二维位置坐标相关的矩阵为F′::i
S4、确定与二维位置坐标相关的矩阵的最小元素对应的二维网格点为目标的二维坐标估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010162061.3A CN111308418B (zh) | 2020-03-10 | 2020-03-10 | 一种对高度未知的目标进行二维定位的稳健方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010162061.3A CN111308418B (zh) | 2020-03-10 | 2020-03-10 | 一种对高度未知的目标进行二维定位的稳健方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111308418A CN111308418A (zh) | 2020-06-19 |
CN111308418B true CN111308418B (zh) | 2021-11-23 |
Family
ID=71157103
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010162061.3A Active CN111308418B (zh) | 2020-03-10 | 2020-03-10 | 一种对高度未知的目标进行二维定位的稳健方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111308418B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112083406B (zh) * | 2020-09-17 | 2022-05-17 | 电子科技大学 | 一种对目标高度稳健的外辐射源二维定位方法 |
CN113391266B (zh) * | 2021-05-28 | 2023-04-18 | 南京航空航天大学 | 基于非圆多嵌套阵降维子空间数据融合的直接定位方法 |
CN113359086B (zh) * | 2021-06-25 | 2023-05-12 | 南京航空航天大学 | 基于增广互质阵列的加权子空间数据融合直接定位方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103713277A (zh) * | 2013-12-19 | 2014-04-09 | 中国航天科工集团八五一一研究所 | 一种基于位置信息场的辐射源定位算法 |
CN103857035A (zh) * | 2012-12-06 | 2014-06-11 | 中国电信股份有限公司 | 三点立体基站定位方法与装置 |
CN107229033A (zh) * | 2017-05-16 | 2017-10-03 | 电子科技大学 | 基于高度维分段搜索的多目标到达时间差定位方法 |
CN107371129A (zh) * | 2017-08-29 | 2017-11-21 | 广州中海达定位技术有限公司 | 基于高度辅助修正的室内定位的tdoa定位方法 |
CN108513254A (zh) * | 2018-02-12 | 2018-09-07 | 广州盛之焰信息科技有限公司 | 一种室内三维定位方法 |
CN108738132A (zh) * | 2017-04-24 | 2018-11-02 | 温州市鹿城区中津先进科技研究院 | 一种基于tdoa的三基站移动通信定位方法 |
CN110209758A (zh) * | 2019-04-18 | 2019-09-06 | 同济大学 | 一种基于张量分解的文本增量降维方法 |
CN110412504A (zh) * | 2019-08-12 | 2019-11-05 | 电子科技大学 | 基于角度与时差信息的无源多站多目标关联与定位方法 |
CN110441732A (zh) * | 2019-08-15 | 2019-11-12 | 电子科技大学 | 一种陌生散射环境中的单站多目标定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6522296B2 (en) * | 2001-06-25 | 2003-02-18 | Harris Corporation | Method and system for calibrating wireless location systems |
US20140303929A1 (en) * | 2013-04-03 | 2014-10-09 | Umm Al-Qura University | Method to obtain accurate vertical component estimates in 3d positioning |
-
2020
- 2020-03-10 CN CN202010162061.3A patent/CN111308418B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103857035A (zh) * | 2012-12-06 | 2014-06-11 | 中国电信股份有限公司 | 三点立体基站定位方法与装置 |
CN103713277A (zh) * | 2013-12-19 | 2014-04-09 | 中国航天科工集团八五一一研究所 | 一种基于位置信息场的辐射源定位算法 |
CN108738132A (zh) * | 2017-04-24 | 2018-11-02 | 温州市鹿城区中津先进科技研究院 | 一种基于tdoa的三基站移动通信定位方法 |
CN107229033A (zh) * | 2017-05-16 | 2017-10-03 | 电子科技大学 | 基于高度维分段搜索的多目标到达时间差定位方法 |
CN107371129A (zh) * | 2017-08-29 | 2017-11-21 | 广州中海达定位技术有限公司 | 基于高度辅助修正的室内定位的tdoa定位方法 |
CN108513254A (zh) * | 2018-02-12 | 2018-09-07 | 广州盛之焰信息科技有限公司 | 一种室内三维定位方法 |
CN110209758A (zh) * | 2019-04-18 | 2019-09-06 | 同济大学 | 一种基于张量分解的文本增量降维方法 |
CN110412504A (zh) * | 2019-08-12 | 2019-11-05 | 电子科技大学 | 基于角度与时差信息的无源多站多目标关联与定位方法 |
CN110441732A (zh) * | 2019-08-15 | 2019-11-12 | 电子科技大学 | 一种陌生散射环境中的单站多目标定位方法 |
Non-Patent Citations (5)
Title |
---|
Single-Station TDOA Passive Location Under Multipath Propagation conditions;Luxi Zhang等;《IOP Conf.Series: Materials Science and Engineering》;20191231;第1-6页 * |
Some Applications of Tensor Algebra to Estimation Theory;Samuel Picton Drake等;《2008 3rd international symposium on wireless pervasive computing》;20080702;第106-110页 * |
利用三站TDOA及辅助高度信息的空间目标二维定位方法及误差分析;杨林等;《电子学报》;19981231;第26卷(第12期);第71-74页 * |
无源时差定位模糊区与定位无解区研究;王婷;《现代雷达》;20171130;第39卷(第11期);第57-62页 * |
无源雷达的二维定位精度分析和图形描述;陆捷等;《火力与指挥控制》;20110228;第36卷(第2期);第135-137页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111308418A (zh) | 2020-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111308418B (zh) | 一种对高度未知的目标进行二维定位的稳健方法 | |
CN107817469B (zh) | 基于非视距环境下超宽频测距实现室内定位方法 | |
CN110058236A (zh) | 一种面向三维地表形变估计的InSAR和GNSS定权方法 | |
CN109917333B (zh) | 融合aoa观测量与tdoa观测量的无源定位方法 | |
CN107991676B (zh) | 星载单航过InSAR系统对流层误差校正方法 | |
US10656240B2 (en) | Hybrid TDOA closed form hyperbolic and spherical iteration geo-location technique | |
CN110132281B (zh) | 一种基于询问应答模式的水下高速目标高精度自主声学导航方法 | |
CN113342059B (zh) | 基于位置和速度误差的多无人机跟踪移动辐射源方法 | |
CN101865656B (zh) | 一种使用少数共面点精确定位多摄像头系统位姿的方法 | |
CN110738275A (zh) | 基于ut-phd的多传感器序贯融合跟踪方法 | |
CN110412634A (zh) | 基于参考站网的伪距差分定位方法及装置 | |
JP2014530353A (ja) | 受信機システムと関連付けられた量を推定する方法 | |
CN112822633B (zh) | 一种基于误差补偿的定位方法,系统,终端设备及可读存储介质 | |
CN111025354A (zh) | 基于单差分电离层加权模型的中长基线rtk定位方法 | |
CN107592654B (zh) | 一种基于压缩感知的同频多辐射源场强定位方法 | |
CN112068099B (zh) | 基于误差补偿的多辐射源快速定位测速方法和装置 | |
CN111398956B (zh) | 多基高比星载sar三维定位rd方程优化配权方法 | |
CN110907960B (zh) | 基于K-Means动态聚类分析的周跳探测方法及装置 | |
CN114742141A (zh) | 基于icp点云的多源信息数据融合研判方法 | |
CN113923590A (zh) | 一种锚节点位置不确定情况下的toa定位方法 | |
CN108120954B (zh) | 一种高精度toa变化率定位方法 | |
CN115308685B (zh) | 一种基于多基站通信伪距信息的外辐射源目标定位方法 | |
CN114428292B (zh) | 近地表速度模型的构建方法和存储介质 | |
CN116990850B (zh) | 一种uwb测距辅助gnss动对动基线模糊度固定方法 | |
Ivanov et al. | Comparative Analysis of Search Methods for the Source of a Magnetic Field of a Dipole Model |
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 | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 09 / F, 13 / F, building 683, zone 2, No. 5, Zhongguancun South Street, Haidian District, Beijing 100081 Patentee after: HUIZHONG XINGZHI TECHNOLOGY (BEIJING) CO.,LTD. Patentee after: University of Electronic Science and technology Address before: 100193209, building 1, East District, yard 10, northwest Wangdong Road, Haidian District, Beijing Patentee before: HUIZHONG XINGZHI TECHNOLOGY (BEIJING) CO.,LTD. Patentee before: University of Electronic Science and technology |