CN105259533A - Three-stage arrival time difference positioning method based on multidimensional scaling sub space analysis - Google Patents
Three-stage arrival time difference positioning method based on multidimensional scaling sub space analysis Download PDFInfo
- Publication number
- CN105259533A CN105259533A CN201510710744.7A CN201510710744A CN105259533A CN 105259533 A CN105259533 A CN 105259533A CN 201510710744 A CN201510710744 A CN 201510710744A CN 105259533 A CN105259533 A CN 105259533A
- Authority
- CN
- China
- Prior art keywords
- calculate
- arrival
- estimation point
- vector
- stage
- 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
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000004458 analytical method Methods 0.000 title claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 239000013598 vector Substances 0.000 claims abstract description 30
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 6
- 230000004807 localization Effects 0.000 claims description 6
- 230000000717 retained effect Effects 0.000 claims 2
- 230000017105 transposition Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000005755 formation reaction Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
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/0252—Radio frequency fingerprinting
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
一种无线网络和移动计算领域的基于多维标度子空间分析的三阶段到达时间差定位方法,包括三个阶段,第一阶段:获知平面上传感器位置坐标,到达时间差,到达时间差的测量误差的方差和信号传播速度,计算修正的测量标量乘积矩阵,通过特征值分解,计算得到第一估计点和一个距离估计值,并计算得到一个方向向量;第二阶段:在第一估计点的基础上,利用距离估计值和方向向量,计算得到第二估计点对应的参数;第三阶段:以第二估计点所对应的参数为初始值,通过二分法求根过程,得到第三估计点所对应的参数,进一步计算第三估计点,得到坐标的最终估计值,由此确定信号源的位置。本发明具有适用范围广,定位精度和准确性高的特点。
A three-stage time difference of arrival positioning method based on multidimensional scale subspace analysis in the field of wireless networks and mobile computing, including three stages, the first stage: obtaining the coordinates of the sensor position on the plane, the time difference of arrival, and the variance of the measurement error of the time difference of arrival and signal propagation speed, calculate the corrected measurement scalar product matrix, through eigenvalue decomposition, calculate the first estimated point and a distance estimated value, and calculate a direction vector; the second stage: on the basis of the first estimated point, Using the estimated distance value and direction vector, calculate the parameters corresponding to the second estimated point; the third stage: use the parameters corresponding to the second estimated point as the initial value, and obtain the corresponding parameters, and further calculate the third estimated point to obtain the final estimated value of the coordinates, thereby determining the position of the signal source. The invention has the characteristics of wide application range and high positioning precision and accuracy.
Description
技术领域technical field
本发明涉及的是无线网络和移动计算领域,具体是一种基于多维标度子空间分析的三阶段到达时间差定位方法。The invention relates to the field of wireless network and mobile computing, in particular to a three-stage time difference of arrival positioning method based on multi-dimensional scale subspace analysis.
背景技术Background technique
在雷达、声纳、移动通信、多媒体、无线传感器网络等应用领域中,常常面临一个重要问题,即依据到达时间差信息,对一个信号源进行定位。所谓到达时间差是指:由信号源发出信号,由分布在空间中、位置已知、而且时间相互同步的传感器接收该信号,并测量信号到达其他传感器的时间,由此计算得到信号源所发出的信号到达各个传感器的时间与到达第一个传感器的时间之差。In the application fields of radar, sonar, mobile communication, multimedia, wireless sensor network, etc., an important problem is often faced, that is, to locate a signal source according to the time difference of arrival information. The so-called time difference of arrival refers to: a signal is sent by a signal source, and the signal is received by a sensor distributed in space, whose position is known and time-synchronized with each other, and the time when the signal arrives at other sensors is measured, and the signal source is calculated from this. The difference between the time a signal arrives at each sensor and the time it arrives at the first sensor.
经过对现有技术的检索发现,He‐WenWei,RongPeng,QunWan,Zhang‐XinChen,&Shang‐FuYe在期刊IEEETransactionsonSignalProcessing的2010年3月第58卷3期发表的论文“MultidimensionalscalinganalysisforpassivemovingtargetlocalizationwithTDOAandFDOAmeasurements”,提出了利用多维标度法分析无源传感器阵列的到达时间差和到达频率差得到信号源位置的方法。但是该方法中当传感器阵列成为线状或近似于线状时,所需要做求逆运算的矩阵将成为秩亏矩阵或高度病态矩阵,强行对其求逆将产生巨大的定位误差。After searching the existing technologies, it was found that He‐WenWei, RongPeng, QunWan, Zhang‐XinChen, & Shang‐FuYe published the paper "Multidimensional scaling analysis for passive moving target localization with TDOAandFDOAmeasurements" in the journal IEEE Transactions on Signal Processing, Volume 58, Issue 3, March 2010. A method of analyzing the arrival time difference and arrival frequency difference of the passive sensor array to obtain the position of the signal source. However, in this method, when the sensor array becomes linear or nearly linear, the matrix that needs to be inverted will become a rank deficient matrix or a highly ill-conditioned matrix, and forcibly inverting it will cause a huge positioning error.
中国专利文献号CN103648164A,公开(公告)日2014.03.19,公开了一种基于到达时间差和Gossip算法的无线传感器网络分布式定位方法。其特点是,一、锚节点获取自身位置坐标;二、实现分布式时间同步;三、锚节点随机唤醒监测未知节点;四、唤醒锚节点保存接收信号时刻和本地坐标;五、所有锚节点是否全部完成信号监测和数据保存;六、锚节点j收到其所有M个相邻锚节点的数据;七、获取锚节点j对于未知节点未知的初始估计值;八、所有锚节点获取未知节点位置初始估计值;九、运行Gossip算法随机选择相邻锚节点交换定位数据;十、算法终止。但是该方法中不同锚节点采集到的未知节点位置估计值有可能不完全一致,如何从所获得的位置信息中计算出最大可能性的位置估计值在申请中并未提及。Chinese Patent Document No. CN103648164A, publication (announcement) date 2014.03.19, discloses a wireless sensor network distributed positioning method based on time difference of arrival and Gossip algorithm. Its characteristics are: 1. The anchor node obtains its own position coordinates; 2. Realizes distributed time synchronization; 3. The anchor node randomly wakes up and monitors unknown nodes; 4. Wakes up the anchor node to save the received signal time and local coordinates; Complete signal monitoring and data storage; 6. Anchor node j receives the data of all M adjacent anchor nodes; 7. Obtain the initial estimated value of anchor node j’s unknown unknown node; 8. All anchor nodes obtain the location of unknown nodes Initial estimated value; 9. Run the Gossip algorithm to randomly select adjacent anchor nodes to exchange positioning data; 10. Terminate the algorithm. However, in this method, the unknown node position estimates collected by different anchor nodes may not be completely consistent, and how to calculate the most likely position estimate from the obtained position information is not mentioned in the application.
发明内容Contents of the invention
本发明针对现有技术存在的上述不足,提出了一种基于多维标度子空间分析的三阶段到达时间差定位方法,通过对测量标量乘积矩阵作子空间分析,分三阶段对信号源位置进行估算,得到信号源较精确的位置。Aiming at the above-mentioned deficiencies in the prior art, the present invention proposes a three-stage time difference of arrival positioning method based on multi-dimensional scale subspace analysis, and estimates the position of the signal source in three stages by performing subspace analysis on the measurement scalar product matrix , to obtain a more accurate position of the signal source.
本发明是通过以下技术方案实现的,The present invention is achieved through the following technical solutions,
本发明包括以下步骤:The present invention comprises the following steps:
步骤1、根据平面上传感器位置坐标、到达时间差(即信号源所发出的信号到达各个传感器的时间与到达第一个传感器的时间之差)得到测量标量乘积矩阵,并通过所述的到达时间差的测量误差的方差对测量标量乘积矩阵进行修正;通过对修正后的测量标量乘积矩阵进行特征值分解处理,计算得到第一估计点,并计算出一个方向向量,和一个距离估计值,具体包括以下步骤:Step 1. Obtain the measurement scalar product matrix according to the position coordinates of the sensors on the plane and the time difference of arrival (i.e. the time difference between the time when the signal sent by the signal source reaches each sensor and the time when the signal arrives at the first sensor), and obtain the measured scalar product matrix through the time difference of arrival time The variance of the measurement error corrects the measurement scalar product matrix; by performing eigenvalue decomposition processing on the corrected measurement scalar product matrix, the first estimated point is calculated, and a direction vector and a distance estimate are calculated, specifically including the following step:
步骤1.1)、生成测量标量乘积矩阵所述的平面上传感器位置坐标为um=[xm,ym]T,m=1,...,M,其中:M表示传感器数量且数量大于等于5个,um表示第m个传感器的位置坐标,xm表示第m个传感器的x轴坐标,ym表示第m个传感器的y轴坐标;Step 1.1), generate measurement scalar product matrix The sensor position coordinates on the plane are u m =[x m ,y m ] T , m=1,...,M, where: M represents the number of sensors and the number is greater than or equal to 5, and u m represents the mth The position coordinates of the sensor, x m represents the x-axis coordinate of the mth sensor, and y m represents the y-axis coordinate of the mth sensor;
根据到达时间差与信号传播速度得到到达距离差:当m=1时当m=2,...,M时其中:表示已测量到的信号源u0到各个传感器um的到达时间与信号源u0到第1个传感器u1的到达时间之差,c表示信号传播速度;According to the arrival time difference and the signal propagation speed, the arrival distance difference is obtained: when m=1 When m=2,...,M in: Indicates the difference between the measured arrival time from signal source u 0 to each sensor u m and the arrival time from signal source u 0 to the first sensor u 1 , c indicates the signal propagation speed;
当m=2,...,M时,根据的测量误差方差与信号传播速度得到的误差方差其中:表示的测量误差方差,c表示信号传播速度;When m=2,...,M, according to The variance of the measurement error and the signal propagation speed get error variance of in: express The measurement error variance of , c represents the signal propagation speed;
所述的测量标量乘积矩阵为其中:矩阵的第i行、第j列元素为The measured scalar product matrix is where: matrix The elements of row i and column j of are
步骤1.2)、对测量标量乘积矩阵进行修正:得到修正后的测量标量乘积矩阵其中:IM表示M×M单位矩阵,1M表示元素全部为1的M维列向量;Step 1.2), correcting the measured scalar product matrix: obtain the corrected measured scalar product matrix Where: I M represents the M×M unit matrix, and 1 M represents the M-dimensional column vector whose elements are all 1;
步骤1.3)、计算第一估计点和距离估计值,具体包括以下步骤:Step 1.3), calculating the first estimated point and the estimated distance, specifically include the following steps:
步骤1.3.1)、对B1-σ2IM作特征值分解:Step 1.3.1), perform eigenvalue decomposition on B 1 -σ 2 I M :
B1-σ2IM=[v1,...,vM]diag(s1,...,sM)[v1,...,vM]T,其中:v1,...,vM∈RM是两两正交、而且模均为1的向量,特征值矩阵diag(s1,...,sM)表示对角元为s1,...,sM的对角矩阵,对特征值s1,...,sM依绝对值进行降序排列,即|s1|≥...≥|sM|;B 1 -σ 2 I M =[v 1 ,...,v M ]diag(s 1 ,...,s M )[v 1 ,...,v M ] T , where: v 1 ,. ..,v M ∈ R M is a pairwise orthogonal vector with a modulus of 1. The eigenvalue matrix diag(s 1 ,...,s M ) indicates that the diagonal elements are s 1 ,...,s The diagonal matrix of M , the eigenvalues s 1 ,...,s M are arranged in descending order of absolute value, that is, |s 1 |≥...≥|s M |;
步骤1.3.2)、由特征向量的线性组合得到系数向量:v1,...,vM为RM的一组标准正交基,故RM中的系数向量v可以表示为v1,...,vM的线性组合v=k1v1+...+k5v5,当m=1,...,5时,基vm的组合系数而当m≥6时,基vm的组合系数k6,...,kM均为0;Step 1.3.2), the coefficient vector is obtained by the linear combination of the eigenvectors: v 1 ,...,v M is a set of orthonormal basis of RM , so the coefficient vector v in RM can be expressed as v 1 , ..., the linear combination of v M v=k 1 v 1 +...+k 5 v 5 , when m=1,...,5, the combination coefficient of base v m And when m≥6, the combination coefficients k 6 ,...,k M of base v m are all 0;
步骤1.3.3)、以系数向量v作为加权系数,对位置坐标矩阵
步骤1.4)、计算方向向量p:当则否则
步骤2、在第一估计点的基础上,利用所述的距离估计值和方向向量p,计算得到第二估计点对应的中间参数t2:Step 2, on the basis of the first estimated point, using the estimated distance and the direction vector p, calculate the intermediate parameter t 2 corresponding to the second estimated point:
定义
步骤3、以第二估计点对应的参数为初始值,通过二分法求根过程,得到第三估计点对应的参数,进一步计算第三估计点,得到坐标的最终估计值,由此确定信号源的位置,具体步骤包括:Step 3. Using the parameter corresponding to the second estimated point as the initial value, obtain the parameter corresponding to the third estimated point through the dichotomy root-finding process, further calculate the third estimated point, and obtain the final estimated value of the coordinate, thereby determining the signal source location, the specific steps include:
步骤3.1)、设定第二估计点u(2)对应的参数t2为初始值:定义fmin=∞,n=1, Step 3.1), setting the parameter t2 corresponding to the second estimation point u (2) as the initial value: define f min =∞, n=1,
步骤3.2)、二分法求根g(t)=0,找出第三估计点u(3)对应的参数t3,将u(t)的分量记为xu,yu,将p的分量记为lx,ly,即:u(t)=[xu,yu]T,p=[lx,ly]T,对m=1,...,M,令wm=[xu-xm,yu-ym]T,dm=||wm||,
步骤3.2.1)、计算g(t2):u(t2)=u(1)+t2p,当g(t2)>0,则转步骤3.2.2,否则转步骤3.2.6;Step 3.2.1), calculate g(t 2 ): u(t 2 )=u (1) +t 2 p, when g(t 2 )>0, go to step 3.2.2, otherwise go to step 3.2.6 ;
步骤3.2.2)、计算g(t2-4σ):u(t2-4σ)=u(1)+(t2-4σ)p,当g(t2-4σ)≤0,则取tleft=t2-4σ,tright=t2,转步骤3.2.10,否则转步骤3.2.3;Step 3.2.2), calculate g(t 2 -4σ): u(t 2 -4σ)=u (1) +(t 2 -4σ)p, when g(t 2 -4σ)≤0, then take t left =t 2 -4σ, t right =t 2 , go to step 3.2.10, otherwise go to step 3.2.3;
步骤3.2.3)、计算g(t2-8σ):u(t2-8σ)=u(1)+(t2-8σ)p,计算g(t2-8σ),当g(t2-8σ)≤0,则取tleft=t2-8σ,tright=t2-4σ,转步骤3.2.10,否则转步骤3.2.4;Step 3.2.3), calculate g(t 2 -8σ): u(t 2 -8σ)=u (1) +(t 2 -8σ)p, calculate g(t 2 -8σ), when g(t 2 -8σ)≤0, then take t left =t 2 -8σ, t right =t 2 -4σ, go to step 3.2.10, otherwise go to step 3.2.4;
步骤3.2.4)、计算g(t2-12σ):u(t2-12σ)=u(1)+(t2-12σ)p,当g(t2-12σ)≤0,则取tleft=t2-12σ,tright=t2-8σ,转步骤3.2.10,否则转步骤3.2.5;Step 3.2.4), calculate g(t 2 -12σ): u(t 2 -12σ)=u (1) +(t 2 -12σ)p, when g(t 2 -12σ)≤0, take t left =t 2 -12σ, t right =t 2 -8σ, go to step 3.2.10, otherwise go to step 3.2.5;
步骤3.2.5)、计算g(t2-16σ):u(t2-16σ)=u(1)+(t2-16σ)p,当g(t2-16σ)≤0,则取tleft=t2-16σ,tright=t2-12σ,转步骤3.2.10,否则取t3=t2-16σ,转步骤3.3;Step 3.2.5), calculate g(t 2 -16σ): u(t 2 -16σ)=u (1) +(t 2 -16σ)p, when g(t 2 -16σ)≤0, take t left =t 2 -16σ, t right =t 2 -12σ, go to step 3.2.10, otherwise take t 3 =t 2 -16σ, go to step 3.3;
步骤3.2.6)、计算g(t2+4σ):u(t2+4σ)=u(1)+(t2+4σ)p,当g(t2+4σ)>0,则取tleft=t2,tright=t2+4σ,转步骤3.2.10,否则转步骤3.2.7;Step 3.2.6), calculate g(t 2 +4σ): u(t 2 +4σ)=u (1) +(t 2 +4σ)p, when g(t 2 +4σ)>0, then take t left =t 2 , t right =t 2 +4σ, go to step 3.2.10, otherwise go to step 3.2.7;
步骤3.2.7)、计算g(t2+8σ):u(t2+8σ)=u(1)+(t2+8σ)p,当g(t2+8σ)>0,则取tleft=t2+4σ,tright=t2+8σ,转步骤3.2.10,否则转步骤3.2.8;Step 3.2.7), calculate g(t 2 +8σ): u(t 2 +8σ)=u (1) +(t 2 +8σ)p, when g(t 2 +8σ)>0, then take t left =t 2 +4σ, t right =t 2 +8σ, go to step 3.2.10, otherwise go to step 3.2.8;
步骤3.2.8)、计算g(t2+12σ):u(t2+12σ)=u(1)+(t2+12σ)p,当g(t2+12σ)>0,则取tleft=t2+8σ,tright=t2+12σ,转步骤3.2.10,否则转步骤3.2.9;Step 3.2.8), calculate g(t 2 +12σ): u(t 2 +12σ)=u (1) +(t 2 +12σ)p, when g(t 2 +12σ)>0, then take t left =t 2 +8σ, t right =t 2 +12σ, go to step 3.2.10, otherwise go to step 3.2.9;
步骤3.2.9)、计算g(t2+16σ):u(t2+16σ)=u(1)+(t2+16σ)p,当g(t2+16σ)>0,则取tleft=t2+12σ,tright=t2+16σ,转步骤3.2.10,否则取t3=t2+16σ,转步骤3.3;Step 3.2.9), calculate g(t 2 +16σ): u(t 2 +16σ)=u (1) +(t 2 +16σ)p, when g(t 2 +16σ)>0, then take t left =t 2 +12σ, t right =t 2 +16σ, go to step 3.2.10, otherwise take t 3 =t 2 +16σ, go to step 3.3;
步骤3.2.10)、取计算g(tmiddle):u(tmiddle)=u(1)+tmiddle·p,当|g(tmiddle)|≤10-3,则取t3=tmiddle,转步骤3.3,否则转步骤3.2.11;Step 3.2.10), take Calculate g(t middle ): u(t middle )=u (1) +t middle ·p, when |g(t middle )|≤10 -3 , then take t 3 =t middle , go to step 3.3, otherwise go to Step 3.2.11;
步骤3.2.11)、当g(tmiddle)>0,则取tright=tmiddle,否则取tleft=tmiddle,当trightt-tleft≤10-2,则取转步骤3.3,否则转步骤3.2.10;Step 3.2.11), when g(t middle )>0, take t right =t middle , otherwise take t left =t middle , when t rightt -t left t≤10 -2 , take Go to step 3.3, otherwise go to step 3.2.10;
步骤3.3)、计算f(t3)和第三估计点u(3):u(3)=u(1)+t3p,当f(t3)≤fmin,则取fmin=f(t3),转步骤3.4;Step 3.3), calculate f(t 3 ) and the third estimated point u (3) : u (3) =u (1) +t 3 p, when f(t 3 )≤f min , then take f min = f(t 3 ), go to step 3.4;
步骤3.4)、对第二估计点u(2)所对应的参数t2设定新的初始值:取n=n+1,当有定义,则取转步骤3.2.1;否则保留步骤3.3中的fmin及其对应的转至步骤3.5;Step 3.4), set a new initial value for the parameter t2 corresponding to the second estimated point u (2) : take n=n+1, when is defined, take Go to step 3.2.1; otherwise keep f min and its corresponding f in step 3.3 Go to step 3.5;
步骤3.5)、决定信号源位置坐标的最终估计值:将步骤3.4保留的作为信号源位置坐标的最终估计值。Step 3.5), determine the final estimated value of the position coordinates of the signal source: keep the step 3.4 as the final estimate of the signal source location coordinates.
技术效果technical effect
与现有技术相比,本发明通过在多维标度法框架下的到达时间差计算到达距离差,适用于正常传感器阵型、近似线状传感器阵型、线状传感器阵型等各种传感器阵型,而且在传感器采集数据并不完全相容的情况下,仍然能够计算出均方误差最小的信号源位置估计值,提高了定位精度和准确性。Compared with the prior art, the present invention calculates the arrival distance difference through the time difference of arrival under the framework of the multidimensional scaling method, and is applicable to various sensor formations such as normal sensor formations, approximate linear sensor formations, and linear sensor formations. When the collected data are not completely compatible, the estimated value of the signal source position with the smallest mean square error can still be calculated, which improves the positioning accuracy and accuracy.
附图说明Description of drawings
图1为本发明流程图。Fig. 1 is the flow chart of the present invention.
具体实施方式detailed description
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。The embodiments of the present invention are described in detail below. This embodiment is implemented on the premise of the technical solution of the present invention, and detailed implementation methods and specific operating procedures are provided, but the protection scope of the present invention is not limited to the following implementation example.
实施例1Example 1
已知一个平面上8个传感器装置,这8个传感器装置的位置坐标分别为:It is known that there are 8 sensor devices on a plane, and the position coordinates of these 8 sensor devices are:
基于多维标度子空间分析的三阶段到达时间差定位方法包含以下三个阶段:The three-stage TDOA localization method based on multidimensional scaled subspace analysis includes the following three stages:
所述的第一阶段包括以下步骤:The first phase described includes the following steps:
步骤1.1)、生成测量标量乘积矩阵:Step 1.1), generate measurement scalar product matrix:
已知信号传播速度为c,归一化为c=1。已测量到信号源u0到各个传感器um的到达时间与信号源u0到第1个传感器u1的到达时间之差为:
已知当m=2,...,8时的测量误差的方差同为则的误差方差2σ2=2×0.32;It is known that when m=2,...,8 The variance of the measurement error is the same as but The error variance 2σ 2 =2×0.3 2 ;
生成测量标量乘积矩阵得到:
步骤1.2)、对测量标量乘积矩阵进行修正:
步骤1.3)、计算第一估计点和一个距离估计值包括以下步骤:Step 1.3), calculating the first estimated point and a distance estimate comprises the following steps:
步骤1.3.1)、对B1-σ2IM作特征值分解得到:
步骤1.3.2)、由特征向量的线性组合得到一个系数向量:
步骤1.3.3)、以系数向量v作为加权系数,对位置坐标矩阵
步骤1.4)、计算方向向量p:因为
所述的第二阶段利用方向向量p和距离估计值计算出第二估计点u(2)所对应的参数t2:a1=-64.7006,a0=50.7965,Δ=4.1354×103;因Δ>0,故 The second stage uses the direction vector p and the distance estimate Calculate the parameter t 2 corresponding to the second estimated point u (2) : a 1 =-64.7006, a 0 =50.7965, Δ=4.1354×10 3 ; since Δ>0, so
所述的第三阶段包括:The third phase described includes:
步骤3.1)、以第二估计点u(2)所对应的参数t2为初始值:取fmin=∞,n=1,则
步骤3.2)、根据参数t2的初始值,利用二分法求根g(t)=0,找出第三估计点u(3)所对应的参数t3:Step 3.2), according to the initial value of the parameter t2 , use the dichotomy to find the root g(t)=0, and find out the corresponding parameter t3 of the third estimated point u (3) :
步骤3.2.1)、计算
步骤3.2.6)、计算
步骤3.2.7)、计算
步骤3.2.8)、计算
步骤3.2.9)、计算
步骤3.3)、计算第三估计点
步骤3.4)、对第二估计点u(2)所对应的参数t2设定新的初始值:取n=n+1=2,
步骤3.2.1)、计算
步骤3.2.2)、计算
步骤3.2.10)、计算|g(tmiddle)|:tmiddle=-0.2063,
步骤3.2.11)、计算trightt-tleft:因为g(tmiddle)≤0,故取tleft=-0.2063,因为trightt-tleft>10-2,故转步骤3.2.10;Step 3.2.11), calculate t rightt -t left : since g(t middle )≤0, take t left t= -0.2063 , and since t rightt -t left >10 -2 , go to step 3.2.10;
步骤3.2.10)、计算|g(tmiddle)|:tmiddle=0.0937,
步骤3.2.11)、计算trightt-tleft:因为g(tmiddle)>0,故取tright=0.0937,因为trightt-tleft>10-2,故转步骤3.2.10;Step 3.2.11), calculate t right -t left : since g(t middle )>0, take t right = 0.0937 , and because t right -t left >10 -2 , go to step 3.2.10;
步骤3.2.10)、计算tmiddle=-0.0563,
步骤3.2.11)、计算trightt-tleft:因为g(tmiddle)≤0,故取tleft=-0.0563,因为trightt-tleft>10-2,故转步骤3.2.10;Step 3.2.11), calculate t rightt -t left : because g(t middle )≤0, take t left t= -0.0563 , and because t rightt -t left >10 -2 , go to step 3.2.10;
步骤3.2.10)、计算|g(tmiddle)|:tmiddle=0.0187,
步骤3.2.11)、计算trightt-tleft:因为g(tmiddle)≤0,故取tleft=0.0187,因为trightt-tleft>10-2,故转步骤3.2.10;Step 3.2.11), calculate t rightt -t left : since g(t middle )≤0, take t left =0.0187, and since t rightt -t left >10 -2 , go to step 3.2.10;
步骤3.2.10)、计算|g(tmiddle)|:tmiddle=0.0562,
步骤3.2.11)、计算trightt-tleft:因为g(tmiddle)>0,故取tright=0.0562,因为trightt-tleft>10-2,故转步骤3.2.10;Step 3.2.11), calculate t right -t left : since g(t middle )>0, take t right = 0.0562 , and because t right -t left >10 -2 , go to step 3.2.10;
步骤3.2.10)、计算|g(tmiddle)|:tmiddle=0.0375,
步骤3.2.11)、计算trightt-tleft:因为g(tmiddle)≤0,故取tleft=0.0375,则trightt-tleft>10-2,故转步骤3.2.10;Step 3.2.11), calculate t rightt -t left : since g(t middle )≤0, take t left =0.0375, then t rightt -t left >10 -2 , so go to step 3.2.10;
步骤3.2.10)、计算|g(tmiddle)|:tmiddle=0.0469,
步骤3.2.11)、计算trightt-tleft:因为g(tmiddle)>0,故取tright=0.0469,则trightt-tleft≤10-2,故取t3=0.0422,转步骤3.3;Step 3.2.11), calculate t rightt -t left : since g(t middle )>0, take t right =0.0469, then t rightt -t left ≤10 -2 , so take t 3 =0.0422, go to step 3.3 ;
步骤3.3)、计算第三估计点
步骤3.4)、以第二估计点u(2)所对应的参数t2为初始值:取n=n+1=3,因为无的定义,故保留步骤3.3中的fmin=1.0242×104及其对应的
步骤3.5)、决定信号源位置坐标的最终估计值:fmin=1.0242×104时,
经本实施例方式计算获得的信号源位置坐标的最终估计值
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510710744.7A CN105259533B (en) | 2015-10-28 | 2015-10-28 | The three stage reaching time-difference localization methods based on multidimensional scaling subspace analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510710744.7A CN105259533B (en) | 2015-10-28 | 2015-10-28 | The three stage reaching time-difference localization methods based on multidimensional scaling subspace analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105259533A true CN105259533A (en) | 2016-01-20 |
CN105259533B CN105259533B (en) | 2017-07-18 |
Family
ID=55099294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510710744.7A Expired - Fee Related CN105259533B (en) | 2015-10-28 | 2015-10-28 | The three stage reaching time-difference localization methods based on multidimensional scaling subspace analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105259533B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105891776A (en) * | 2016-04-06 | 2016-08-24 | 上海交通大学 | Direct method time difference of arrival positioning method based on multidimensional scaling (MDS) model |
CN108279411A (en) * | 2018-02-01 | 2018-07-13 | 电子科技大学 | A kind of passive MIMO time difference positioning methods based on MDS |
CN111551897A (en) * | 2020-04-25 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | TDOA localization method based on weighted multidimensional scaling and polynomial root finding in the presence of sensor location prior observation errors |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009010832A1 (en) * | 2007-07-18 | 2009-01-22 | Bang & Olufsen A/S | Loudspeaker position estimation |
CN101354435A (en) * | 2008-09-05 | 2009-01-28 | 清华大学 | Self-localization method of sensor network nodes based on distance order relationship |
US20110140969A1 (en) * | 2006-04-19 | 2011-06-16 | Mustafa Ergen | Method And System For Hybrid Positioning Using Partial Distance Information |
CN104038901A (en) * | 2014-05-30 | 2014-09-10 | 中南大学 | Indoor positioning method for reducing fingerprint data acquisition workload |
-
2015
- 2015-10-28 CN CN201510710744.7A patent/CN105259533B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110140969A1 (en) * | 2006-04-19 | 2011-06-16 | Mustafa Ergen | Method And System For Hybrid Positioning Using Partial Distance Information |
WO2009010832A1 (en) * | 2007-07-18 | 2009-01-22 | Bang & Olufsen A/S | Loudspeaker position estimation |
CN101354435A (en) * | 2008-09-05 | 2009-01-28 | 清华大学 | Self-localization method of sensor network nodes based on distance order relationship |
CN104038901A (en) * | 2014-05-30 | 2014-09-10 | 中南大学 | Indoor positioning method for reducing fingerprint data acquisition workload |
Non-Patent Citations (3)
Title |
---|
HE-WEN WEI等: "multidimensional scaling analysis for passive moving target localization with TDOA and FDOA measurement", 《IEEE TRANSACTIONS ON SIGNLA PROCESSING》 * |
彭鑫: "无线传感器网络中基于多维标度的节点定位算法", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
王琳等: "一种多维尺度分析到达时间差定位算法", 《导航定位学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105891776A (en) * | 2016-04-06 | 2016-08-24 | 上海交通大学 | Direct method time difference of arrival positioning method based on multidimensional scaling (MDS) model |
CN105891776B (en) * | 2016-04-06 | 2018-06-12 | 上海交通大学 | Direct method reaching time-difference localization method based on MDS models |
CN108279411A (en) * | 2018-02-01 | 2018-07-13 | 电子科技大学 | A kind of passive MIMO time difference positioning methods based on MDS |
CN108279411B (en) * | 2018-02-01 | 2020-04-14 | 电子科技大学 | A Passive MIMO Time Difference Location Method Based on MDS |
CN111551897A (en) * | 2020-04-25 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | TDOA localization method based on weighted multidimensional scaling and polynomial root finding in the presence of sensor location prior observation errors |
CN111551897B (en) * | 2020-04-25 | 2021-01-22 | 中国人民解放军战略支援部队信息工程大学 | TDOA localization method based on weighted multidimensional scaling and polynomial root finding under sensor position error |
Also Published As
Publication number | Publication date |
---|---|
CN105259533B (en) | 2017-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105866735B (en) | The reaching time-difference iteration localization method of amendment cost function based on MDS models | |
CN105954712B (en) | The direct localization method of the multiple target of associated wireless electric signal complex envelope and carrier phase information | |
CN106793087B (en) | Array antenna indoor positioning method based on AOA and PDOA | |
CN104853435B (en) | A kind of indoor orientation method based on probability and device | |
CN101221238B (en) | Dynamic deviation estimation method based on gauss average value mobile registration | |
CN105158730B (en) | TDOA localization methods based on MDS subspaces the 4th and the 5th characteristic vector | |
CN106405533B (en) | Radar target combined synchronization and localization method based on constraint weighted least-squares | |
CN104619020A (en) | WIFI Indoor Positioning Method Based on RSSI and TOA Ranging | |
CN104038901B (en) | Indoor positioning method for reducing fingerprint data acquisition workload | |
CN103369466B (en) | A kind of map match assists indoor orientation method | |
CN107124762B (en) | Wireless positioning method for efficiently eliminating non-line-of-sight errors | |
CN109195110B (en) | Indoor localization method based on hierarchical clustering technology and online extreme learning machine | |
CN104684081A (en) | Node localization algorithm for selecting anchor nodes based on distance clustering in wireless sensor networks | |
CN108051779A (en) | A kind of positioning node preferred method towards TDOA | |
CN106353720A (en) | Multi-station continuous positioning model based on TDOA/GROA (time different of arrival/gain ratio of arrival) | |
CN111551897A (en) | TDOA localization method based on weighted multidimensional scaling and polynomial root finding in the presence of sensor location prior observation errors | |
CN104581945B (en) | The WLAN indoor orientation methods of semi-supervised APC clustering algorithms based on distance restraint | |
CN102200573A (en) | Method for determining incoming wave direction of near-field target signal | |
CN105259533A (en) | Three-stage arrival time difference positioning method based on multidimensional scaling sub space analysis | |
CN110673196A (en) | Time difference positioning method based on multidimensional calibration and polynomial root finding | |
CN113923590B (en) | A TOA positioning method under the condition of uncertain anchor node position | |
CN105891776B (en) | Direct method reaching time-difference localization method based on MDS models | |
CN114325581A (en) | Elliptical target positioning method with clock synchronization error | |
CN104635206B (en) | A kind of method and device of wireless location | |
WO2017049914A1 (en) | Terminal positioning method, apparatus, and system |
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: 20191028 |