基于MDS子空间第4和第5特征向量的TDOA定位方法技术

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

【技术实现步骤摘要】

本专利技术涉及的是一种信息
的技术，具体是一种基于MDS子空间第4和第5特征向量的到达时间差(TimeDifferenceOfArrival，TDOA)定位方法。
技术介绍
在雷达、声纳、移动通信、多媒体、无线传感器网络等应用领域中，常常面临一个重要问题，即依据TDOA信息，对一个信号源进行定位。所谓的到达时间差是指，由信号源发出信号，由分布在空间中、位置已知、而且时间相互同步的传感器接收该信号，并测量信号到达各个传感器的时间，由此计算得到信号源所发出的信号到达各个传感器的时间与到达参考传感器的时间之差，即到达时间差。He‐WenWei等在《MultidimensionalscalinganalysisforpassivemovingtargetlocalizationwithTDOAandFDOAmeasurements》(IEEETransactionsonSignalProcessing,vol.58,no.3,2010,pp.1677‐1688)中提出一种多维标度(Multidimensionalscaling,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个传感器的位置坐标分别为：u1=11.899850.5957,u2=49.836469.9077,u3=95.974489.0903,]]>u4=34.038695.9291,u5=58.526854.7216,u6=22.381213.8624,u7=75.126714.9294]]>和u8=25.509525.7508.]]>所述的信号源u0的真实位置假定为u0=84.071725.4282,]]>该位置待求。测量信号源u0到达各传感器与到达参考传感器的时间差，即到达时间差分别为：t^21=-20.7060,t^31=-11.4547,t^41=-10.5034,t^51=-37.7748,t^61=-13.4115,t^71=-62.2653]]>和该到达时间差的测量误差的方差假设为信号传播速度c归一化为1。计算信号源u0到达第m(m＝2,...,8)个传感器与到达参考传感器的距离差，即到达距离差d^m1=c·t^m1(m=2,...,8),]]>且d^11=0.]]>将到达时间差和信号传播速度c代入计算，得到的到达距离差分别为：d^21=c·t^21=-20.7060,d^31=c·t^31=-11.4547,d^41=c·t^41=-10.5034,d^51=c·t^51=-37.7748,]]>d^61=c·t^61=-13.4115,d^71=c·t^71=-62.2653]]>和d^81=c·t^81=-18.3493.]]>根据到达时间差的测量误差方差和信号传播速度c，计算得到到达距离差的误差方差2σ2=c2·2σt2=2×0.32.]]>步骤2、基于到达距离差和及其误差方差2σ2生成修正的标量乘积矩阵B1，然后进行子空间分析，对得到的特征值依绝对值降序排列，提取第4特征向量和第5特征向量并作线叠加，最后将线性叠加后得到的向量作为组合系数，对位置坐标矩阵x1...x8y1...y8]]>的列向量作线性组合，得到信号源u0位置坐标估计值，具体包括：步骤2.1)基于到达距离差和生成标量乘积矩阵(为实数域上的8阶方阵)，其第l行、第n列元素为[B^]l,n=12(d^l1-d^n1)2-12[(xl-xn)2+(yl-yn)2](l,n=1,2,...8).]]>将各传感器的位置坐标um和对应的到达距离差代入，计算得到所述的标量乘积矩阵B^=103×0-0.6917-4.2096-1.2175-0.3821-0.6397-0.6964-0.2329-0.69170-1.20560.0237-0.0074-1.9208-0.9675-1.2680-4.2096-1.20560-1.7003-0.9454-5.5357-1.6764-4.4648-1.21750.02本文档来自技高网...

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

【技术特征摘要】
1.一种基于MDS子空间第4和第5特征向量的TDOA定位方法，其特征在于，通过采集
分布在平面内的传感器的位置坐标，测量信号源到达各传感器与到达参考传感器的到达时间差，
计算对应的到达距离差；然后基于到达距离差生成标量乘积矩阵并进行子空间分析，对得到的
特征值依绝对值降序排列，提取其中的第4和第5特征向量并作线性叠加；最后通过线性叠加
后的向量作为组合系数，对位置坐标矩阵的列向量进行线性组合得到信号源位置坐标估计值。
2.根据权利要求1所述的基于MDS子空间第4和第5特征向量的TDOA定位方法，其特
征是，所述的传感器的数量为至少5个。
3.根据权利要求1所述的基于MDS子空间第4和第5特征向量的TDOA定位方法，其特
征是，所述的基于到达距离差生成标量乘积矩阵是指：基于到达距离差和m＝2,...,M，生
成标量乘积矩阵其中：为实数域上的M阶方阵，其第l行、第n列元素为
[B^]l,n=12(d^l1-d^n1)2-12[(xl-xn)2+(yl-yn)2],]]>l,n＝1,2,...M，M为传感器的数量。
4.根据权利要求1所述的基于MDS子空间第4和第5特征向量的TDOA定位方法，其特
征是，所述的子空间分析是指：生成修正的标量乘积矩阵并进行特征值分解，得到特征向量。
5.根据权利要求1或4所述的基于MDS子空间第4和...

【专利技术属性】
技术研发人员：蒋武扬，徐昌庆，裴凌，郁文贤，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：上海;31