CN110174643A - 一种无需噪声功率信息的基于到达时间差的定位方法 - Google Patents
一种无需噪声功率信息的基于到达时间差的定位方法 Download PDFInfo
- Publication number
- CN110174643A CN110174643A CN201910406753.5A CN201910406753A CN110174643A CN 110174643 A CN110174643 A CN 110174643A CN 201910406753 A CN201910406753 A CN 201910406753A CN 110174643 A CN110174643 A CN 110174643A
- Authority
- CN
- China
- Prior art keywords
- difference
- noise
- equation
- arrival rate
- sensor
- 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
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/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/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/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)
- Electromagnetism (AREA)
- Probability & Statistics with Applications (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种无需噪声功率信息的基于到达时间差的定位方法,通过引入冗余变量的方式将到达时间差的定位方程及到达频率差的定位方程转化为到达时间差的矩阵方程及到达频率差的矩阵方程,并根据到达时间差测量噪声的协方差矩阵以及到达频率差测量噪声的协方差矩阵提取并分离测量噪声功率分量,再运用半正定规划方法精确求解待定位目标的位置和速度,从而可以在测量噪声功率未知的条件下精确估计待定位目标的位置和速度。
Description
技术领域
本发明涉及无线移动通信技术领域,尤其是一种无需噪声功率信息的基于到达时间差的定位方法。
背景技术
无源定位是指定位方自身不发射信号,而通过接收目标发射或反射的信号来对其进行定位的技术。由于无源定位***具有探测距离远、抗干扰能力强和隐蔽性好等优点,其在雷达、声呐、电子对抗、无线通信和无线传感器网络等领域都得到了广泛的应用,同时也是信号处理领域的重要研究内容之一。
无源定位的主要原理是先由目标发出的信号中提取定位参数,再由这些定位参数完成目标定位。在实际应用中,无源定位需要的定位参数主要有信号强度(ReceivedSignal Strength,RSS)、到达角度(Angle Of Arrival,AOA)、到达时间(Time Of Arrival,AOA)、到达时间差(Time Difference Of Arrival,TDOA)、到达频率差(FrequencyDifference Of Arrival,FDOA)等。其中,基于到达时间差(TDOA)的和到达频率差(FDOA)的定位具有测量精度高、抗干扰能力强、设备简单、无需与目标时间同步等优势,在实际应用中得到了普遍的关注。因此,研究基于到达时间差的定位技术有着十分重要的现实意义。
目前常用的基于到达时频差的定位算法主要有泰勒级数展开算法、两步加权最小二乘算法、限制加权最小二乘算法、半正定规划算法等算法。其中,两步加权最小二乘算法计算简单且具有闭式解,但当背景噪声的强度较高时,其定位精度非常不理想;泰勒级数展开算法和限制加权最最小二乘算法求解过程需要迭代,因此其对于初值的选取有一定要求且无法保证能收敛到全局最优解。半正定规划算法具有最好的定位精度,但其计算负担较重。另外,这些算法的实现依赖准确的噪声功率估计,当噪声功率无法得到准确的估计或者其估计误差较大时,这些方法的定位性能会严重恶化。由于在实际中经常面对噪声功率无法得到准确估计的情形,因此这个缺陷严重制约了这些基于到达时频差算法的实际应用。
目前在基于到达时频差的定位研究中,不依赖噪声功率先验信息的定位研究还很少,在已知的算法中,只有总体最小二乘法和偏差修正法不需要知道噪声的功率先验信息,但是这两种方法的定位精度均比较低,而且这两种方法均存在严重的“门限效应”:当噪声强度达到某个电平时,算法的定位误差会急剧增加。因此,研究在功率信息未知的条件下仍能保持良好的定位精度的基于到达时间差的定位算法的非常必要的。
发明内容
本发明的目的在于提供一种无需噪声功率信息的基于到达时间差的定位方法,以提高现有的无源定位方法的精度。
为了达到上述目的,本发明提供了一种无需噪声功率信息的基于到达时间差的定位方法,包括:
建立到达时间差的定位方程及到达频率差的定位方程,并得到到达时间差测量噪声的协方差矩阵以及到达频率差测量噪声的协方差矩阵;
对所述到达时间差的定位方程进行线性化以形成到达时间差的矩阵方程;
将所述到达时间差的矩阵方程转换为第一加权最小二乘问题,并利用所述到达时间差测量噪声的协方差矩阵分离噪声功率以形成第二加权最小二乘问题;
利用半正定规划算法求解所述第二加权最小二乘问题,以得到估计目标位置;
将所述估计目标位置带入所述到达频率差的定位方程中,并对所述到达频率差的定位方程进行线性化以形成到达频率差的矩阵方程;
将所述到达频率差的矩阵方程转换为第三加权最小二乘问题,并利用所述到达频率差测量噪声的协方差矩阵分离噪声功率以形成第四加权最小二乘问题;
利用半正定规划算法求解所述第四加权最小二乘问题,以得到估计目标速度。
可选的,设定待定位目标是三维的,所述待定位目标的位置为u=[x0,y0,z0],速度为接收M个所述待定位目标发射的信号的传感器的位置为si=[xi,yi,zi]T,速度为则所述到达时间差的定位方程为:
其中,其中,c为信号的传播速度,ni1为到达时间差信号的测量噪声,ti1代表所述待定位目标发出的信号到达第1个传感器和第i个传感器的时间差,代表所述待定位目标到第一个传感器的距离与第i个传感器的距离的真实距离差,ri1代表由测量得到的到达时间计算得出的含噪声的距离差,ri代表所述待定位目标到第i个传感器之间的距离,r1代表所述待定位目标到第1个传感器之间的距离;
所述到达频率差的定位方程为:
其中,代表所述待定位目标到第i个传感器之间的距离差分率,fi1代表所述待定位目标发出的信号到达第1个传感器和第i个传感器的频率差,代表所述待定位目标到第1个传感器的距离差分率和第i个传感器的距离差分率之间的真实差值,代表到达各传感器之间的测量多普勒频率差计算得到的含噪声的差值,为到达频率差信号的测量噪声。
可选的,令n=[n21,...,nM1]T,则所述到达时间差测量噪声的协方差矩阵为:
所述到达频率差测量噪声的协方差矩阵为:
其中,n21,...,nM1为ni1的分量信息,为的分量信息,Rt和Rf分别为所述到达时间差测量噪声的协方差矩阵及所述到达频率差测量噪声的协方差矩阵的结构,和分别为到达时间差测量噪声的功率及到达频率差测量噪声的功率。
可选的,对所述到达时间差的定位方程进行线性化以形成到达时间差的矩阵方程的步骤包括:
将所述到达时间差的定位方程略去c,并在等式两边平方后除以以变形为如下公式:
其中i=2,...,M,将r1视作与u0无关的冗余变量,并共同构成关于位置的待求向量θt=[u,r1],则所述到达时间差的矩阵方程为εt=ht-Gtθt,
其中,
εt=Btn为到达时间差的测量噪声部分,Bt为基于到达时间差的定位方程组矩阵化后各噪声分量的权值矩阵,且Bt满足如下公式:
可选的,所述第一加权最小二乘问题为:
其中,Wt为第一权值矩阵,且
所述第二加权最小二乘问题为:
其中,Wt_new为第二权值矩阵,且
可选的,利用半正定规划算法求解所述第二加权最小二乘问题,以得到估计目标位置的步骤包括:
令则所述第二加权最小二乘问题为:
其中,
由等式r1=||u-s1||构建约束,并将松弛为以得到第一半正定规划问题:
求解所述第一半正定规划问题并取θt的前三项为所述估计目标位置。
可选的,对所述到达频率差的定位方程进行线性化以形成到达频率差的矩阵方程的步骤包括:
将所述估计目标位置带入所述到达频率差的定位方程中,并将所述到达频率差的定位方程略去c,并在等式两边平方后除以以变形为如下公式:
其中i=2,...,M,将视作与无关的冗余变量,并共同构成关于速度的待求向量则所述到达频率差的矩阵方程为εf=hf-Gfθf,
其中,
为到达频率差的测量噪声部分,Bf为基于到达频率差的定位方程组矩阵化后各噪声分量的权值矩阵,且Bf满足如下公式:
可选的,所述第三加权最小二乘问题为:
其中,Wf为第三权值矩阵,且
所述第二加权最小二乘问题为:
其中,Wf_new为第四权值矩阵,且
可选的,利用半正定规划算法求解所述第四加权最小二乘问题,以得到估计目标速度的步骤包括:
令则所述第四加权最小二乘问题为:
其中,
由等式构建约束,并将松弛为以得到第二半正定规划问题:
求解所述第一半正定规划问题并取θf的前三项为所述估计目标速度。
在本发明提供的无需噪声功率信息的基于到达时间差的定位方法中,通过引入冗余变量的方式将到达时间差的定位方程及到达频率差的定位方程转化为到达时间差的矩阵方程及到达频率差的矩阵方程,并根据到达时间差测量噪声的协方差矩阵以及到达频率差测量噪声的协方差矩阵提取并分离测量噪声功率分量,再运用半正定规划方法精确求解待定位目标的位置和速度,从而可以在测量噪声功率未知的条件下精确估计待定位目标的位置和速度。
附图说明
图1为本发明实施例提供的定位方法的流程图;
图2为本发明实施例提供的近场源位置的均方根误差在高斯白噪声下随噪声强度变化关系示意图;
图3为本发明实施例提供的近场源的速度的均方根误差在高斯白噪声下随噪声强度变化关系示意图;
图4为本发明实施例提供的远场源的位置的均方根误差在高斯白噪声下随噪声强度变化关系示意图;
图5为本发明实施例提供的远场源的速度的均方根误差在高斯白噪声下随噪声强度变化关系示意图。
具体实施方式
下面将结合示意图对本发明的具体实施方式进行更详细的描述。根据下列描述和权利要求书,本发明的优点和特征将更清楚。需说明的是,附图均采用非常简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施例的目的。
如图1所示,本实施例提供了一种无需噪声功率信息的基于到达时间差的定位方法,包括如下步骤:
步骤1:设置M个接收待定位目标发射的信号的传感器
假设待定位目标是三维的,其位置和速度分别为u=[x0,y0,z0]和M个接收目标发射信号的传感器的位置和速度可以分别表示为si=[xi,yi,zi]T和例如传感器数设置5个,5个传感器位置的三维坐标分别为[300,100,150]、[400,150,100]、[300,500,200]、[350,200,100]和[-100,-100,-100],单位为米,速度的三维分量分别为[30,-20,20]、[-30,10,20]、[10,-20,10]、[10,20,30]和[-20,10,10],单位为米每秒,本实施例分别定位一个近场信号源和一个远场信号源,近场信号源位置的三维坐标为[600,650,550],单位为米,速度的三维分量为[-20,15,40],单位为米每秒;远场信号源位置的三维坐标为[2000,2500,3000],单位为米,速度的三维分量为[-20,15,40],单位为米每秒,传感器上的噪声为高斯白噪声。
选择第一个传感器为参考传感器,则第i个传感器和第一个传感器之间的理论到达时间差为
其中c代表信号的传播速度,ri=||u-si||代表信号源(远场信号源或近场信号源)到第i个传感器之间的距离。
考虑到噪声影响,第i个传感器和第一个传感器之间的真实到达时间差为:
其中,ni1是到达时间差测量噪声,ti1代表所述待定位目标发出的信号到达第1个传感器和第i个传感器的时间差,代表所述待定位目标到第一个传感器的距离与第i个传感器的距离的真实距离差,ri1代表由测量得到的到达时间计算得出的含噪声的距离差,ri代表所述待定位目标到第i个传感器之间的距离,r1代表所述待定位目标到第1个传感器之间的距离;
待定位目标和第i个传感器之间的速度变化率可表示为:
因此i个传感器和第一个传感器之间的理论到达频率差为
同样考虑噪声影响,第i个传感器和第一个传感器之间的真实到达频率差为:
其中,代表所述待定位目标到第i个传感器之间的距离差分率,fi1代表所述待定位目标发出的信号到达第1个传感器和第i个传感器的频率差,代表所述待定位目标到第1个传感器的距离差分率和第i个传感器的距离差分率之间的真实差值,代表到达各传感器之间的测量多普勒频率差计算得到的含噪声的差值,为到达频率差信号的测量噪声。由于ti1和ri1、fi1和之间仅存在常数c倍的差别,因此为表示简洁,在后文中均使用ri1和代表到达时间差和到达频率差。将到达时间差信号的测量噪声及到达频率差信号的测量噪声分量用向量的形式表示为n=[n21,...,nM1]T和则到达时间差测量噪声的协方差矩阵为到达频率差测量噪声的协方差矩阵为Rt和Rf分别为所述到达时间差测量噪声的协方差矩阵及所述到达频率差测量噪声的协方差矩阵的结构,和分别为到达时间差测量噪声的功率及到达频率差测量噪声的功率。
步骤2:线性化到达时间差的定位方程:
将到达时间差的定位方程移项得ri1-ni1+r1=ri,将该式两边平方并忽略二阶误差项后,可将到达时间差的定位方程变形为:
由于r1=||u-s1||,因此该非线性关系使得上述方程是非线性的,为了将到达时间差的定位方程线性化,将r1视作与u无关的冗余变量,并定义关于位置的待求向量θt=[u,r1],则到达时间差的定位方程可以用矩阵的形式表示为:
εt=ht-Gtθt
其中
εt=Btnt代表到达时间差的测量噪声部分,Bt为基于到达时间差的定位方程组矩阵化后各噪声分量的权值矩阵,且Bt满足如下公式:
步骤3:利用到达时间差的矩阵方程构成加权最小二乘问题:
依据最小二乘估计的思想,到达时间差的矩阵方程的解可以由求解如下最小二乘问题的极小值得到:
其中,Wt为第一权值矩阵,且
由于所以最小二乘问题也可以表示为:
可见,到达时间差测量噪声功率在最小二乘问题中对θt的结果没有影响,因此可以将其忽略构建新的最小二乘问题:
其中,Wt_new为第二权值矩阵,且
步骤4:运用半正定规划算法求解待定位目标的估计目标位置:
定义则新构建的加权最小二乘问题可以等价地表示为:
其中,
利用已知的u和r1之间的关系r1=||u-s1||,将其两边平方后我们可以得到:
θ(4)2=(u1-s1)T(u1-s1)
其可以等价地写为如下形式:
再利用半正定松弛的思想,将非凸的约束条件松弛为得到以下半正定规划问题:
该问题是凸优化问题,直接使用成熟的凸优化求解工具CVX进行求解获得向量θt,取该向量的前三个元素作为所述估计目标位置的三维坐标。
步骤5:线性化到达频率差的定位方程:
在得到估计目标位置u以后,对到达频率差的定位方程求导,并利用可将到达频率差的定位方程变形为如下公式:
将视作与无关的冗余变量,并共同构成关于速度的待求向量可将到达频率差的定位方程转换为矩阵的形式:
εf=hf-Gfθf
其中
εf=Bfnf代表到达频率差的测量噪声部分,Bf为基于到达频率差的定位方程组矩阵化后各噪声分量的权值矩阵,且Bf满足如下公式:
步骤6:利用到达频率差的矩阵方程构成加权最小二乘问题:
由于在步骤4中,预估目标位置u已经估计得出,因此步骤5中矩阵Bf中的元素ri可以由u计算得出。因此到达频率差的矩阵方程可以构成如下加权最小二乘问题:
其中,Wf为第三权值矩阵,且
再利用到达频率差测量噪声的协方差矩阵的结构,可以将到达频率差测量噪声功率分离构建新的最小二乘问题:
其中,Wf_new为第四权值矩阵,且
步骤7:运用半正定规划算法求解预估目标速度:
定义则新构建的加权最小二乘问题可以等价地表示为:
其中,
再由等式构建约束,并将松弛为可得半正定规划问题:
其中,
利用CVX最优化工具包直接求解该问题并取θf的前三项为待定位目标的预估目标速度。
为验证算法的性能,本实施例设计两组仿真实验。仿真分别定位一个近场信号源和一个远场信号源,信号源的位置和速度如上所述。仿真中采取的噪声设置为典型的高斯白噪声,且各传感器受到的噪声影响相等,测量噪声的协方差矩阵的结构Rt和Rf可以表示为主对角线元素为1,其余元素为0.5的(M-1)×(M-1)维矩阵。在仿真实验中,将本实施例中的定位方法与2种其它的方法的性能进行了对比,其余的方法分别为:总体最小二乘算法和偏差修正算法。
两组实验的随机实验次数均为5000,半正定规划的求解工具为CVX工具箱中的SeDuMi求解器。四组实验的结果分别如图2、3、4、5所示。由图2、3、4、5可知,在噪声干扰较强的条件下,本发明提供的定位方法较其它的定位方法得到的定位精度更高,且估计性能不会随着背景噪声的增强出现大幅恶化。
综上,在本发明实施例提供的定位方法中,通过引入冗余变量的方式将到达时间差的定位方程及到达频率差的定位方程转化为到达时间差的矩阵方程及到达频率差的矩阵方程,并根据到达时间差测量噪声的协方差矩阵以及到达频率差测量噪声的协方差矩阵提取并分离测量噪声功率分量,再运用半正定规划方法精确求解待定位目标的位置和速度,从而可以在测量噪声功率未知的条件下精确估计待定位目标的位置和速度。
上述仅为本发明的优选实施例而已,并不对本发明起到任何限制作用。任何所属技术领域的技术人员,在不脱离本发明的技术方案的范围内,对本发明揭露的技术方案和技术内容做任何形式的等同替换或修改等变动,均属未脱离本发明的技术方案的内容,仍属于本发明的保护范围之内。
Claims (9)
1.一种无需噪声功率信息的基于到达时间差的定位方法,其特征在于,包括:
建立到达时间差的定位方程及到达频率差的定位方程,并得到到达时间差测量噪声的协方差矩阵以及到达频率差测量噪声的协方差矩阵;
对所述到达时间差的定位方程进行线性化以形成到达时间差的矩阵方程;
将所述到达时间差的矩阵方程转换为第一加权最小二乘问题,并利用所述到达时间差测量噪声的协方差矩阵分离噪声功率以形成第二加权最小二乘问题;
利用半正定规划算法求解所述第二加权最小二乘问题,以得到估计目标位置;
将所述估计目标位置带入所述到达频率差的定位方程中,并对所述到达频率差的定位方程进行线性化以形成到达频率差的矩阵方程;
将所述到达频率差的矩阵方程转换为第三加权最小二乘问题,并利用所述到达频率差测量噪声的协方差矩阵分离噪声功率以形成第四加权最小二乘问题;
利用半正定规划算法求解所述第四加权最小二乘问题,以得到估计目标速度。
2.如权利要求1所述的定位方法,其特征在于,设定待定位目标是三维的,所述待定位目标的位置为u=[x0,y0,z0],速度为接收M个所述待定位目标发射的信号的传感器的位置为si=[xi,yi,zi]T,速度为则所述到达时间差的定位方程为:
其中,c为信号的传播速度,ni1为到达时间差信号的测量噪声,ti1代表所述待定位目标发出的信号到达第1个传感器和第i个传感器的时间差,代表所述待定位目标到第一个传感器的距离与第i个传感器的距离的真实距离差,ri1代表由测量得到的到达时间计算得出的含噪声的距离差,ri代表所述待定位目标到第i个传感器之间的距离,r1代表所述待定位目标到第1个传感器之间的距离;
所述到达频率差的定位方程为:
其中,代表所述待定位目标到第i个传感器之间的距离差分率,fi1代表所述待定位目标发出的信号到达第1个传感器和第i个传感器的频率差,代表所述待定位目标到第1个传感器的距离差分率和第i个传感器的距离差分率之间的真实差值,代表到达各传感器之间的测量多普勒频率差计算得到的含噪声的差值,为到达频率差信号的测量噪声。
3.如权利要求2所述的定位方法,其特征在于,令n=[n21,...,nM1]T,则所述到达时间差测量噪声的协方差矩阵为:
所述到达频率差测量噪声的协方差矩阵为:
其中,n21,...,nM1为ni1的分量信息,为的分量信息,Rt和Rf分别为所述到达时间差测量噪声的协方差矩阵及所述到达频率差测量噪声的协方差矩阵的结构,和分别为到达时间差测量噪声的功率及到达频率差测量噪声的功率。
4.如权利要求3所述的定位方法,其特征在于,对所述到达时间差的定位方程进行线性化以形成到达时间差的矩阵方程的步骤包括:
将所述到达时间差的定位方程略去c,并在等式两边平方后除以以变形为如下公式:
其中i=2,...,M,将r1视作与u0无关的冗余变量,并共同构成关于位置的待求向量θt=[u,r1],则所述到达时间差的矩阵方程为εt=ht-Gtθt,
其中,
εt=Btn为到达时间差的测量噪声部分,Bt为基于到达时间差的定位方程组矩阵化后各噪声分量的权值矩阵,且Bt满足如下公式:
n=[n21,…,nM1]T。
5.如权利要求4所述的定位方法,其特征在于,所述第一加权最小二乘问题为:
其中,Wt为第一权值矩阵,且
所述第二加权最小二乘问题为:
其中,Wt_new为第二权值矩阵,且
6.如权利要求5所述的定位方法,其特征在于,利用半正定规划算法求解所述第二加权最小二乘问题,以得到估计目标位置的步骤包括:
令则所述第二加权最小二乘问题为:
其中,
由等式r1=||u-s1||构建约束,并将松弛为以得到第一半正定规划问题:
求解所述第一半正定规划问题并取θt的前三项为所述估计目标位置。
7.如权利要求6所述的定位方法,其特征在于,对所述到达频率差的定位方程进行线性化以形成到达频率差的矩阵方程的步骤包括:
将所述估计目标位置带入所述到达频率差的定位方程中,并将所述到达频率差的定位方程略去c,并在等式两边平方后除以以变形为如下公式:
其中i=2,...,M,将视作与无关的冗余变量,并共同构成关于速度的待求向量则所述到达频率差的矩阵方程为εf=hf-Gfθf,
其中,
为到达频率差的测量噪声部分,Bf为基于到达频率差的定位方程组矩阵化后各噪声分量的权值矩阵,且Bf满足如下公式:
8.如权利要求7所述的定位方法,其特征在于,所述第三加权最小二乘问题为:
其中,Wf为第三权值矩阵,且
所述第二加权最小二乘问题为:
其中,Wf_new为第四权值矩阵,且
9.如权利要求8所述的定位方法,其特征在于,利用半正定规划算法求解所述第四加权最小二乘问题,以得到估计目标速度的步骤包括:
令则所述第四加权最小二乘问题为:
其中,
由等式构建约束,并将松弛为以得到第二半正定规划问题:
求解所述第一半正定规划问题并取θf的前三项为所述估计目标速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910406753.5A CN110174643B (zh) | 2019-05-16 | 2019-05-16 | 一种无需噪声功率信息的基于到达时间差的定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910406753.5A CN110174643B (zh) | 2019-05-16 | 2019-05-16 | 一种无需噪声功率信息的基于到达时间差的定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110174643A true CN110174643A (zh) | 2019-08-27 |
CN110174643B CN110174643B (zh) | 2021-01-05 |
Family
ID=67691366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910406753.5A Expired - Fee Related CN110174643B (zh) | 2019-05-16 | 2019-05-16 | 一种无需噪声功率信息的基于到达时间差的定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110174643B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110568406A (zh) * | 2019-09-03 | 2019-12-13 | 浙江万里学院 | 一种能量衰减因子未知条件下基于声能的定位方法 |
CN111239718A (zh) * | 2020-01-17 | 2020-06-05 | 电子科技大学 | 一种基于单星辐射源的多基站目标定位方法 |
CN111505584A (zh) * | 2020-05-25 | 2020-08-07 | 中国科学技术大学 | 基于相对速度和多普勒效应的音频信号到达时间修正方法 |
CN111551897A (zh) * | 2020-04-25 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | 传感器位置先验观测误差存在下基于加权多维标度和多项式求根的tdoa定位方法 |
CN112540342A (zh) * | 2020-11-19 | 2021-03-23 | 安徽大学 | 一种移动信号源定位方法 |
CN112986906A (zh) * | 2021-02-23 | 2021-06-18 | 中国计量大学上虞高等研究院有限公司 | 一种半正定规划的rss-toa联合定位方法 |
CN113084801A (zh) * | 2021-03-30 | 2021-07-09 | 深圳市人工智能与机器人研究院 | 基于半正定规划优化的多机器人初始位姿相对定位方法 |
CN114690116A (zh) * | 2022-03-25 | 2022-07-01 | 西安交通大学 | 基于tdoa与fdoa量测的无源目标定位方法及*** |
CN114910864A (zh) * | 2022-06-14 | 2022-08-16 | 中国人民解放军战略支援部队信息工程大学 | 一种信号传播速度未知且存在信号频率漂移的多平台多普勒定位方法 |
CN115508774A (zh) * | 2022-10-12 | 2022-12-23 | 中国电子科技集团公司信息科学研究院 | 基于两步加权最小二乘的时差定位方法、装置和存储介质 |
CN116017281A (zh) * | 2022-12-30 | 2023-04-25 | 深圳市中承科技有限公司 | 一种基于超宽带通信技术的室内定位方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040233105A1 (en) * | 2001-05-04 | 2004-11-25 | Lockheed Martin Corporation | System and method for central association and tracking in passive coherent location applications |
EP2368133B1 (fr) * | 2008-12-23 | 2012-09-05 | Thales | Procede de localisation par estimation multi-voies des tdoa et fdoa des multi-trajets d'une source avec ou sans aoa |
CN105759241A (zh) * | 2015-12-18 | 2016-07-13 | 中国航天科工集团八五研究所 | 基于时差频差的直接定位方法 |
US20170079007A1 (en) * | 2013-03-15 | 2017-03-16 | DGS Global Systems, Inc. | Systems, methods, and devices for electronic spectrum management with remote access to data in a virtual computing network |
CN107526073A (zh) * | 2017-08-22 | 2017-12-29 | 哈尔滨工程大学 | 一种运动多站无源时差频差联合定位方法 |
CN109581281A (zh) * | 2018-11-08 | 2019-04-05 | 宁波大学 | 基于到达时间差和到达频率差的移动目标定位方法 |
CN109633592A (zh) * | 2019-01-18 | 2019-04-16 | 杭州电子科技大学 | 运动观测站误差下外辐射源雷达时差与频差协同定位方法 |
CN109633581A (zh) * | 2018-12-26 | 2019-04-16 | 杭州电子科技大学 | 基于外辐射源tdoa/fdoa误差校正下的定位方法 |
-
2019
- 2019-05-16 CN CN201910406753.5A patent/CN110174643B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040233105A1 (en) * | 2001-05-04 | 2004-11-25 | Lockheed Martin Corporation | System and method for central association and tracking in passive coherent location applications |
EP2368133B1 (fr) * | 2008-12-23 | 2012-09-05 | Thales | Procede de localisation par estimation multi-voies des tdoa et fdoa des multi-trajets d'une source avec ou sans aoa |
US20170079007A1 (en) * | 2013-03-15 | 2017-03-16 | DGS Global Systems, Inc. | Systems, methods, and devices for electronic spectrum management with remote access to data in a virtual computing network |
CN105759241A (zh) * | 2015-12-18 | 2016-07-13 | 中国航天科工集团八五研究所 | 基于时差频差的直接定位方法 |
CN107526073A (zh) * | 2017-08-22 | 2017-12-29 | 哈尔滨工程大学 | 一种运动多站无源时差频差联合定位方法 |
CN109581281A (zh) * | 2018-11-08 | 2019-04-05 | 宁波大学 | 基于到达时间差和到达频率差的移动目标定位方法 |
CN109633581A (zh) * | 2018-12-26 | 2019-04-16 | 杭州电子科技大学 | 基于外辐射源tdoa/fdoa误差校正下的定位方法 |
CN109633592A (zh) * | 2019-01-18 | 2019-04-16 | 杭州电子科技大学 | 运动观测站误差下外辐射源雷达时差与频差协同定位方法 |
Non-Patent Citations (3)
Title |
---|
RUI-RUI LIU 等: ""Semidefinite programming for NLOS localization U sing TDOA and FDOA measurements"", 《2017 3RD IEEE INTERNATIONAL CONFERENCE ON COMPUTER AND COMMUNICATIONS (ICCC)》 * |
YUNLONG WANG 等: ""An Efficient Semidefinite Relaxation Algorithm for Moving Source Localization Using TDOA and FDOA Measurements"", 《IEEE COMMUNICATIONS LETTERS》 * |
张杰 等: ""无线传感器网络中基于TDOA/FDOA的增强半正定松弛定位算法研究"", 《传感技术学报》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110568406A (zh) * | 2019-09-03 | 2019-12-13 | 浙江万里学院 | 一种能量衰减因子未知条件下基于声能的定位方法 |
CN110568406B (zh) * | 2019-09-03 | 2021-07-20 | 浙江万里学院 | 一种能量衰减因子未知条件下基于声能的定位方法 |
CN111239718A (zh) * | 2020-01-17 | 2020-06-05 | 电子科技大学 | 一种基于单星辐射源的多基站目标定位方法 |
CN111551897B (zh) * | 2020-04-25 | 2021-01-22 | 中国人民解放军战略支援部队信息工程大学 | 传感器位置误差下基于加权多维标度和多项式求根的tdoa定位方法 |
CN111551897A (zh) * | 2020-04-25 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | 传感器位置先验观测误差存在下基于加权多维标度和多项式求根的tdoa定位方法 |
CN111505584A (zh) * | 2020-05-25 | 2020-08-07 | 中国科学技术大学 | 基于相对速度和多普勒效应的音频信号到达时间修正方法 |
CN111505584B (zh) * | 2020-05-25 | 2022-07-15 | 中国科学技术大学 | 基于相对速度和多普勒效应的音频信号到达时间修正方法 |
CN112540342B (zh) * | 2020-11-19 | 2024-01-19 | 安徽大学 | 一种移动信号源定位方法 |
CN112540342A (zh) * | 2020-11-19 | 2021-03-23 | 安徽大学 | 一种移动信号源定位方法 |
CN112986906A (zh) * | 2021-02-23 | 2021-06-18 | 中国计量大学上虞高等研究院有限公司 | 一种半正定规划的rss-toa联合定位方法 |
CN112986906B (zh) * | 2021-02-23 | 2022-05-03 | 中国计量大学上虞高等研究院有限公司 | 一种半正定规划的rss-toa联合定位方法 |
CN113084801A (zh) * | 2021-03-30 | 2021-07-09 | 深圳市人工智能与机器人研究院 | 基于半正定规划优化的多机器人初始位姿相对定位方法 |
CN114690116A (zh) * | 2022-03-25 | 2022-07-01 | 西安交通大学 | 基于tdoa与fdoa量测的无源目标定位方法及*** |
CN114910864B (zh) * | 2022-06-14 | 2023-08-15 | 中国人民解放军战略支援部队信息工程大学 | 一种信号传播速度未知且存在信号频率漂移的多平台多普勒定位方法 |
CN114910864A (zh) * | 2022-06-14 | 2022-08-16 | 中国人民解放军战略支援部队信息工程大学 | 一种信号传播速度未知且存在信号频率漂移的多平台多普勒定位方法 |
CN115508774A (zh) * | 2022-10-12 | 2022-12-23 | 中国电子科技集团公司信息科学研究院 | 基于两步加权最小二乘的时差定位方法、装置和存储介质 |
CN116017281A (zh) * | 2022-12-30 | 2023-04-25 | 深圳市中承科技有限公司 | 一种基于超宽带通信技术的室内定位方法 |
CN116017281B (zh) * | 2022-12-30 | 2023-10-24 | 深圳市中承科技有限公司 | 一种基于超宽带通信技术的室内定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110174643B (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110174643A (zh) | 一种无需噪声功率信息的基于到达时间差的定位方法 | |
CN103969622B (zh) | 一种基于多运动接收站的时差定位方法 | |
CN104316903B (zh) | 一种三站时差定位性能试验评估方法 | |
CN105072581B (zh) | 一种基于路径衰减系数建库的室内定位方法 | |
CN106793087B (zh) | 一种基于aoa和pdoa的阵列天线室内定位方法 | |
CN106019217A (zh) | 基于aoa的二维无线传感器网络半定规划定位方法 | |
CN103596267A (zh) | 一种基于欧氏距离的指纹图匹配方法 | |
CN112986906B (zh) | 一种半正定规划的rss-toa联合定位方法 | |
CN107690184A (zh) | 联合tdoa‑aoa无线传感器网络半定规划定位方法 | |
CN105898865A (zh) | 非线性非高斯条件下基于ekf和pf的协同定位方法 | |
CN106842118A (zh) | 用于分布式多点定位监视***的时差定位方法 | |
CN105353351A (zh) | 一种基于多信标到达时间差改进型定位方法 | |
CN109151714A (zh) | 一种三维鲁棒估计定位方法 | |
CN110515037A (zh) | 非视距环境下能时频多域联合的被动定位方法 | |
Meissner et al. | Analysis of position-related information in measured UWB indoor channels | |
CN106249197A (zh) | 一种多点定位***中接收机位置误差的自校准方法 | |
CN105572639B (zh) | 一种室内超声波差分定位方法 | |
Zwirello et al. | Localization in industrial halls via ultra-wideband signals | |
CN110657806A (zh) | 一种基于CKF、chan解算和Savitzky-Golay平滑滤波的位置解算方法 | |
CN111157943A (zh) | 在异步网络中基于toa的传感器位置误差抑制方法 | |
CN107613458B (zh) | 一种tdoa条件下最优联合时间同步与定位的定位方法 | |
Retscher et al. | Indoor positioning using Wi-Fi lateration—Comparison of two common range conversion models with two novel differential approaches | |
CN110426672A (zh) | 一种基于区间分析的双迭代定位***及其应用 | |
Arai et al. | Color radiomap interpolation for efficient fingerprint wifi-based indoor location estimation | |
CN110365420A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210105 Termination date: 20210516 |
|
CF01 | Termination of patent right due to non-payment of annual fee |