CN105158730B - 基于mds子空间第4和第5特征向量的tdoa定位方法 - Google Patents

基于mds子空间第4和第5特征向量的tdoa定位方法 Download PDF

Info

Publication number
CN105158730B
CN105158730B CN201510500766.0A CN201510500766A CN105158730B CN 105158730 B CN105158730 B CN 105158730B CN 201510500766 A CN201510500766 A CN 201510500766A CN 105158730 B CN105158730 B CN 105158730B
Authority
CN
China
Prior art keywords
vector
characteristic vector
sensor
mds
position coordinates
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.)
Expired - Fee Related
Application number
CN201510500766.0A
Other languages
English (en)
Other versions
CN105158730A (zh
Inventor
蒋武扬
徐昌庆
裴凌
郁文贤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201510500766.0A priority Critical patent/CN105158730B/zh
Publication of CN105158730A publication Critical patent/CN105158730A/zh
Application granted granted Critical
Publication of CN105158730B publication Critical patent/CN105158730B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Complex Calculations (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

一种基于MDS子空间第4和第5特征向量的TDOA定位方法，通过采集分布在平面内的传感器的位置坐标，测量信号源到达各传感器与到达参考传感器的到达时间差，计算对应的到达距离差；然后基于到达距离差生成标量乘积矩阵并进行子空间分析，对得到的特征值依绝对值降序排列，提取其中的第4和第5特征向量并作线性叠加；最后通过线性叠加后的向量作为组合系数，对位置坐标矩阵的列向量进行线性组合得到信号源位置坐标估计值；本发明运用多维标度法，通过测量得到的到达时间差在较小误差范围内估计信号源的位置坐标。

Description

基于MDS子空间第4和第5特征向量的TDOA定位方法
技术领域
本发明涉及的是一种信息技术领域的技术，具体是一种基于MDS子空间第4和第5特征向量的到达时间差(Time Difference Of Arrival，TDOA)定位方法。
背景技术
在雷达、声纳、移动通信、多媒体、无线传感器网络等应用领域中，常常面临一个重要问题，即依据TDOA信息，对一个信号源进行定位。所谓的到达时间差是指，由信号源发出信号，由分布在空间中、位置已知、而且时间相互同步的传感器接收该信号，并测量信号到达各个传感器的时间，由此计算得到信号源所发出的信号到达各个传感器的时间与到达参考传感器的时间之差，即到达时间差。
He‐Wen Wei等在《Multidimensional scaling analysis for passive movingtarget localization with TDOA and FDOA measurements》(IEEE Transactions onSignal Processing,vol.58,no.3,2010,pp.1677‐1688)中提出一种多维标度(Multidimensional scaling,MDS)分析方法，根据在无噪情况下信号子空间和噪声子空间之间的正交性，推导出基于MDS的到达时间差标量乘积矩阵，用以估计移动目标的位置和速度。
经过对现有技术的检索发现，中国专利文献号CN104090261A，公告日2014.10.8，公开了一种TDOA定位系统中采用距离建模的定位方法，自身位置未知的目标接收两个固定定位站发射的信号，测量这两个发射信号到达目标的时间差，计算出目标到固定定位站的距离差，用距离差进行建模，从而得到目标的当前位置。但该技术需要知道目标到基准站的距离，在已知量较少的情况下难以根据所述的方程求得目标的具体坐标。
发明内容
本发明针对现有技术存在的上述不足，提出一种基于MDS子空间第4和第5特征向量的TDOA定位方法，通过到达时间差计算相应的到达距离差，生成标量乘积矩阵并进行子空间分析，提取依绝对值降序排列的第4和第5特征值所对应的特征向量，作线性叠加后作为组合系数对位置坐标矩阵的列向量进行线性组合，得到信号源位置坐标的估计值。
本发明是通过以下技术方案实现的：
本发明通过采集分布在平面内的传感器的位置坐标，测量信号源到达各传感器与到达参考传感器的到达时间差，计算对应的到达距离差；然后基于到达距离差生成标量乘积矩阵并进行子空间分析，对得到的特征值依绝对值降序排列，提取其中的第4和第5特征向量并作线性叠加；最后通过线性叠加后的向量作为组合系数，对位置坐标矩阵的列向量进行线性组合得到信号源位置坐标估计值。
所述的子空间分析是指：生成修正的标量乘积矩阵并进行特征值分解，得到特征向量。
所述的传感器的数量M≥5。
技术效果
与现有技术相比，本发明通过在多维标度法的框架下的到达时间差计算到达距离差，无需知道信号源与参考传感器的距离，便可在较小的误差范围内得到信号源的位置坐标估计值，算法简单，准确率高。
附图说明
图1为本发明流程图。
具体实施方式
下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
实施例1
如图1所示，本实施例包括以下步骤：
步骤1、采集分布在平面上的传感器的位置坐标um＝[xm,ym]T(m＝1,...,M)，指定第1个传感器为参考传感器，测量信号源u0到达各传感器与到达参考传感器的到达时间差并根据信号传播速度c计算对应的到达距离差
所述的传感器数量M＝8。
设8个传感器的位置坐标分别为： 和
所述的信号源u0的真实位置假定为该位置待求。
测量信号源u0到达各传感器与到达参考传感器的时间差，即到达时间差分别为：和该到达时间差的测量误差的方差假设为信号传播速度c归一化为1。
计算信号源u0到达第m(m＝2,...,8)个传感器与到达参考传感器的距离差，即到达距离差且
将到达时间差和信号传播速度c代入计算，得到的到达距离差分别为： 和
根据到达时间差的测量误差方差和信号传播速度c，计算得到到达距离差的误差方差
步骤2、基于到达距离差和及其误差方差2σ2生成修正的标量乘积矩阵B1，然后进行子空间分析，对得到的特征值依绝对值降序排列，提取第4特征向量和第5特征向量并作线叠加，最后将线性叠加后得到的向量作为组合系数，对位置坐标矩阵的列向量作线性组合，得到信号源u0位置坐标估计值，具体包括：
步骤2.1)基于到达距离差和生成标量乘积矩阵(为实数域上的8阶方阵)，其第l行、第n列元素为
将各传感器的位置坐标um和对应的到达距离差代入，计算得到所述的标量乘积矩阵
步骤2.2)进行子空间分析，对得到的特征值依绝对值降序排列，提取第4特征向量和第5特征向量并作线叠加，得到叠加后的向量，具体包括：
步骤2.2.1)生成修正的标量乘积矩阵
所述的修正的标量乘积矩阵其中：为标量乘积矩阵，2σ2为(m＝2,...,8)的误差方差，I8为8×8单位矩阵，18为元素全部为1的8维列向量，T表示矩阵的转置。
将所述的标量乘积矩阵代入计算，得到修正的标量乘积矩阵
步骤2.2.2)对B1-σ2I8(2σ2为到达距离差的误差方差，I8为8×8单位矩阵)作特征值分解，得到m个特征值和特征向量，并对特征值依绝对值降序排列。
所述的特征值分解是指：B1-σ2I8＝[v1,...,v8]diag(s1,...,s8)[v1,...,v8]T，其中：是两两正交、而且模均为1的向量，diag(s1,...,s8)表示对角元为s1,...,s8的对角矩阵，T表示矩阵的转置。对s1,...,s8依绝对值进行降序排列，即|s1|≥...≥|s8|。
对B1-σ2I8作特征值分解得到： 和s1＝-1.1334×104、s2＝8.5322×103、s3＝2.8357×103、s4＝-66.9765、s5＝31.8773、s6＝5.7825×10-13、s7＝-5.1324×10-13和s8＝-1.1186×10-13。
步骤2.2.3)提取步骤2.2.2)中的第4特征向量v4和第5特征向量v5并作线叠加，得到叠加后的向量。
所述的线性叠加是指：计算得到其中：18为元素全部为1的8维列向量，T表示矩阵的转置；以k4和k5作为叠加系数，对第4特征向量v4和第5特征向量v5作线性叠加，得到叠加后的向量v＝k4v4+k5v5。
将步骤2.2.2)中的第4特征向量v4和第5特征向量v5代入，得到：a4＝0.7259，a5＝0.5155，
步骤2.3)以叠加后的向量v作为组合系数，对位置坐标矩阵的列向量进行线性组合，得到信号源u0位置坐标估计值
代入各值计算，得到信号源u0位置坐标的估计值
本实施例与现有技术相比，能够大幅度减小真实位置与上述估计值之间的误差。

Claims (5)

1.一种基于MDS子空间第4和第5特征向量的定位方法，其特征在于，通过采集分布在平面内的传感器的位置坐标，测量信号源到达各传感器与到达参考传感器的到达时间差，计算对应的到达距离差；然后基于到达距离差生成标量乘积矩阵并进行子空间分析，对得到的特征值依绝对值降序排列，提取其中的第4和第5特征向量并作线性叠加；最后通过线性叠加后的向量与位置坐标矩阵的列向量进行线性组合得到信号源位置坐标估计值；
所述的子空间分析具体包括：
步骤1、生成修正的标量乘积矩阵即：其中：为标量乘积矩阵，2σ2为m＝2,...,M的误差方差，IM为M×M单位矩阵，1M为元素全部为1的M维列向量，T表示矩阵的转置；
步骤2、对B1-σ2I8作特征值分解，得到特征向量，即：B1-σ2IM＝[v1,...,vM]diag(s1,...,sM)[v1,...,vM]T，其中：是两两正交、而且模均为1的向量，diag(s1,...,sM)表示对角元为s1,...,sM的对角矩阵；对s1,...,s8依绝对值进行降序排列，即|s1|≥...≥|sM|。
2.根据权利要求1所述的基于MDS子空间第4和第5特征向量的定位方法方法，其特征是，所述的传感器的数量为至少5个。
3.根据权利要求1所述的基于MDS子空间第4和第5特征向量的定位方法方法，其特征是，所述的基于到达距离差生成标量乘积矩阵是指：基于到达距离差和m＝2,...,M，生成标量乘积矩阵其中：为实数域上的M阶方阵，其第l行、第n列元素为l,n＝1,2,...M，M为传感器的数量。
4.根据权利要求1所述的基于MDS子空间第4和第5特征向量的定位方法方法，其特征是，所述的提取其中的第4和第5特征向量并作线性叠加是指：计算得到以k4和k5作为叠加系数，对所述的第4特征向量v4和第5特征向量v5作线性叠加，得到叠加后的向量v＝k4v4+k5v5。
5.根据权利要求4所述的基于MDS子空间第4和第5特征向量的定位方法方法，其特征是，所述的信号源位置坐标估计值，通过以下方式得到：以叠加后的向量v作为组合系数与位置坐标矩阵的列向量进行线性组合，得到信号源u0位置坐标估计值
CN201510500766.0A 2015-08-14 2015-08-14 基于mds子空间第4和第5特征向量的tdoa定位方法 Expired - Fee Related CN105158730B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510500766.0A CN105158730B (zh) 2015-08-14 2015-08-14 基于mds子空间第4和第5特征向量的tdoa定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510500766.0A CN105158730B (zh) 2015-08-14 2015-08-14 基于mds子空间第4和第5特征向量的tdoa定位方法

Publications (2)

Publication Number Publication Date
CN105158730A CN105158730A (zh) 2015-12-16
CN105158730B true CN105158730B (zh) 2017-07-18

Family

ID=54799648

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510500766.0A Expired - Fee Related CN105158730B (zh) 2015-08-14 2015-08-14 基于mds子空间第4和第5特征向量的tdoa定位方法

Country Status (1)

Country Link
CN (1) CN105158730B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105891776B (zh) * 2016-04-06 2018-06-12 上海交通大学 基于mds模型的直接法到达时间差定位方法
CN105866735B (zh) * 2016-04-06 2018-04-10 上海交通大学 基于mds模型的修正代价函数的到达时间差迭代定位方法
CN108279411B (zh) * 2018-02-01 2020-04-14 电子科技大学 一种基于mds的被动mimo时差定位方法
CN112346014B (zh) * 2020-09-23 2022-06-21 宁波大学 一种基于信号到达时间差的多基地声纳的定位方法
CN113835064B (zh) * 2021-08-13 2022-09-23 中国人民解放军战略支援部队信息工程大学 一种协同校正源观测信息的加权多维标度tdoa定位方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101014877A (zh) * 2003-08-14 2007-08-08 塞恩塞斯有限公司 使用tdoa分布型天线的目标定位法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
PT1654561E (pt) * 2003-08-14 2012-10-24 Saab Sensis Corp Localização de um alvo utilizando antena distribuída para tdoa
US7663547B2 (en) * 2007-04-13 2010-02-16 Glowlink Communications Technology, Inc. Determining a geolocation solution of an emitter on earth based on weighted least-squares estimation
CN102186242B (zh) * 2011-05-09 2014-04-16 江南大学 固定区域无线传感器网络移动节点定位方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101014877A (zh) * 2003-08-14 2007-08-08 塞恩塞斯有限公司 使用tdoa分布型天线的目标定位法

Also Published As

Publication number Publication date
CN105158730A (zh) 2015-12-16

Similar Documents

Publication Publication Date Title
CN105158730B (zh) 基于mds子空间第4和第5特征向量的tdoa定位方法
Cheung et al. A multidimensional scaling framework for mobile location using time-of-arrival measurements
CN105866735B (zh) 基于mds模型的修正代价函数的到达时间差迭代定位方法
CN105954712B (zh) 联合无线电信号复包络和载波相位信息的多目标直接定位方法
CN103713277B (zh) 一种基于位置信息场的辐射源定位方法
Liu et al. Joint TDOA and AOA location algorithm
CN106405533B (zh) 基于约束加权最小二乘的雷达目标联合同步与定位方法
CN104519571B (zh) 一种基于rss的室内定位方法
CN108872932B (zh) 基于神经网络的超视距目标直接定位结果纠偏方法
CN107707324B (zh) 一种基于相位差和最大似然估计的声信号时延估计方法
CN106658713A (zh) 基于多参数估计的单基站移动用户定位方法
CN104038901A (zh) 一种减少指纹数据采集工作量的室内定位方法
CN106353720A (zh) 基于tdoa/groa的多站连续定位模型
Wei et al. A new NLOS error mitigation algorithm in location estimation
Yu et al. Practical constrained least-square algorithm for moving source location using TDOA and FDOA measurements
CN111551897B (zh) 传感器位置误差下基于加权多维标度和多项式求根的tdoa定位方法
Sun et al. 3D hybrid TOA-AOA source localization using an active and a passive station
CN113960532B (zh) 一种基于假想源的二次定位计算的微地震定位方法
CN105259533B (zh) 基于多维标度法子空间分析的三阶段到达时间差定位方法
Sun et al. Total least-squares solution of active target localization using TDOA and FDOA measurements in WSN
CN105891776B (zh) 基于mds模型的直接法到达时间差定位方法
CN105704676A (zh) 利用信号时间相关性提高指纹室内定位精度的方法
Lee et al. Analytic equivalence of iterated TOA and TDOA techniques under structured measurement characteristics
Li et al. Robust kernel-based machine learning localization using NLOS TOAs or TDOAs
Spiesberger Geometry of locating sounds from differences in travel time: Isodiachrons

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170718

Termination date: 20200814