CN115397015A - 一种距离约束下联合aoa和rss的多源协同定位方法 - Google Patents
一种距离约束下联合aoa和rss的多源协同定位方法 Download PDFInfo
- Publication number
- CN115397015A CN115397015A CN202211014219.8A CN202211014219A CN115397015A CN 115397015 A CN115397015 A CN 115397015A CN 202211014219 A CN202211014219 A CN 202211014219A CN 115397015 A CN115397015 A CN 115397015A
- Authority
- CN
- China
- Prior art keywords
- radiation source
- distance
- weighted
- constraint
- nth
- 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 39
- 230000005855 radiation Effects 0.000 claims abstract description 147
- 238000005457 optimization Methods 0.000 claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 9
- 230000003068 static effect Effects 0.000 claims abstract description 7
- 230000014509 gene expression Effects 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000004807 localization Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001010 compromised effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
- H04W64/006—Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination
-
- 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/0257—Hybrid positioning
-
- 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/0295—Proximity-based methods, e.g. position inferred from reception of particular signals
- G01S5/02955—Proximity-based methods, e.g. position inferred from reception of particular signals by computing a weighted average of the positions of the signal transmitters
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种距离约束下联合AOA和RSS的多源协同定位方法,首先利用多个静止观测站同时获得多个静止辐射源的AOA观测量和RSS观测量,并引入辅助变量将联合AOA和RSS的观测方程伪线性处理,由此建立加权最小二乘估计准则;然后利用每个辐射源与观测站的距离等式约束、多辐射源之间的距离不等式约束,并结合加权最小二乘估计准则构建定位优化模型;最后利用半定松弛方法对此模型进行数值优化,进而得到辐射源位置向量的估计值。由于多辐射源之间距离不等式约束的存在,本发明提出的方法可以进一步提高对多个静止辐射源的定位精度。
Description
技术领域
本发明涉及无线信号定位技术领域,尤其涉及一种距离约束下联合AOA和RSS的多源协同定位方法。
背景技术
众所周知,无线信号定位技术广泛应用于通信、雷达、目标监测、导航遥测、地震勘测、射电天文、紧急救助、安全管理等领域,其在工业生产和军事应用中都发挥着重要作用。对目标进行定位(即位置参数估计)可以使用雷达、激光、声纳等有源设备来完成,该类技术称为有源定位技术,它具有全天候、高精度等优点。然而,有源定位系统通常需要依靠发射大功率电磁信号来实现,因此极易暴露自己的位置,容易被对方发现,从而遭到对方电子干扰的影响,导致定位性能急剧恶化,甚至会危及系统自身的安全性和可靠性。
目标定位还可以利用目标(主动)辐射或者(被动)散射的无线电信号来实现,该类技术称为无源定位技术,它是指在观测站不主动发射电磁信号的情况下,通过接收目标辐射或者散射的无线电信号来估计目标位置参数。与有源定位系统相比,无源定位系统具有不主动发射电磁信号、生存能力强、侦察作用距离远等优点,从而得到国内外学者的广泛关注和深入研究。无源定位系统根据观测站数目可以划分为单站无源定位系统和多站无源定位系统两大类,其中多站定位系统能够综合利用多个观测站的无源探测信息及空间差异性,有助于及时准确地获取定位信息。
辐射源定位的基本过程就是从电磁信号中提取出与其位置和速度有关的参数(也称为定位观测量),然后再利用这些参数解算出辐射源位置向量和速度向量。用于辐射源定位的观测量主要涉及到空、时、频、能量等多域参数,每个域的定位参数都有其特定的应用场景。在实际应用中,为了提高对辐射源的定位精度,可以考虑联合多域观测量进行定位。TOA、TDOA、FOA、FDOA观测量进行辐射源定位需要精确的同步,而小型智能天线及新型阵列技术的发展使AOA观测量应用越来越多、RSS定位具有低复杂性及易实现性等优点。
实际定位时,多辐射源定位更具有应用性。在多辐射定位场景中可以考虑通过多源协同定位来提升每个辐射源的定位精度。但是,现有的多辐射源定位中,学者们更多关注多观测站与多辐射源间的数据关联,忽视了多辐射源间的相关信息,如飞行编队、无人机等,导致定位精度不理想。
发明内容
本发明针对现有的多辐射源定位中,忽视了多辐射源间的相关信息,导致定位精度不理想的问题,提出一种距离约束下联合AOA和RSS的多源协同定位方法,可以进一步提高对多个静止辐射源的定位精度。
为了实现上述目的,本发明采用以下技术方案:
首先,利用多个静止观测站同时获得多个静止辐射源的AOA观测量(包括方位角观测量、仰角观测量)和RSS观测量(即接收信号强度观测量),并引入辅助变量将联合AOA和RSS的观测方程伪线性处理,由此建立加权最小二乘估计准则。然后,利用每个辐射源与观测站的距离等式约束、多辐射源之间的距离不等式约束,并结合加权最小二乘估计准则构建定位优化模型。最后,利用半定松弛方法对此模型进行数值优化,进而得到辐射源位置向量的估计值。由于多辐射源之间距离不等式约束的存在,本发明提出的方法可以进一步提高对多个静止辐射源的定位精度。
本发明的一种距离约束下联合AOA和RSS的多源协同定位方法,具体包括:
步骤1:利用空间中放置的M个静止观测站同时获得第n(1≤n≤N)个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角观测量。
步骤2:利用辅助变量依次对第n(1≤n≤N)个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角观测量进行伪线性处理,并建立第n(1≤n≤N)个辐射源信号的加权最小二乘估计准则,其中un表示第n个辐射源的位置向量,dn表示第n个辐射源的距离向量。
步骤3:确定第n(1≤n≤N)个辐射源信号的加权权重。
步骤4:利用N个辐射源的加权最小二乘估计准则和第n(1≤n≤N)个辐射源信号到达第m个观测站的距离等式约束,获得N个辐射源的加权最小二乘解关系式。
步骤5:确定两两辐射源之间的距离不等式约束。
步骤6:利用加权最小二乘估计准则、每个辐射源与观测站的距离等式约束、多辐射源之间的距离不等式约束构建定位优化模型。
步骤7:利用半定松弛方法对步骤6中的定位优化模型进行数值优化,从而确定辐射源位置向量的估计值。
式中分别表示第n个辐射源的三维空间坐标;分别表示第m个观测站的三维空间坐标;分别表示第n个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角的真实值;α表示路径损耗因子,其取值范围通常为2~5;P0为发射功率,其与天线增益和辐射源的能量有关,这里假设其为已知量;dnm为第n个辐射源和第m个观测站之间的距离;εp,nm、εθ,nm、εβ,nm分别表示第n个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角的观测误差,其相互独立且均服从均值为零的高斯分布,方差分别为
进一步地,所述步骤2中,需要依次对接收信号强度Pnm、方位角θnm、仰角观测量βnm进行伪线性处理,并建立第n(1≤n≤N)个辐射源信号的加权最小二乘估计准则。
首先根据接收信号强度(RSS)观测方程,将其表示为辐射源与观测站之间的距离:
在小噪声水平下,将上式等式右边使用一阶Taylor级数展开并忽略高阶项可得:
则建立加权最小二乘估计准则有:
Ap,nvn=ζp,n(1≤n≤N)
式中Ap,n=[OM×3 IM -λ]∈RM×(M+4),IM表示M×M的单位阵,λ=[λn1 λn2 … λnM]T;vn为辅助变量,dn=[dn1 dn2 … dnM]T;ζp,n=[ζp,n1 ζp,n2 … ζp,nM]T,ζp,nm=(λnmln10/10α)εp,nm。
将上式展开并忽略高阶项可得:
则建立加权最小二乘估计准则有:
Aθ,nvn=ζθ,n(1≤n≤N)
式中
将上式展开并忽略高阶项可得:
则建立加权最小二乘估计准则有:
Aβ,nvn=ζβ,n(1≤n≤N)
式中
最后将接收信号强度、方位角、仰角伪线性处理得到的加权最小二乘估计准则合并获得第n(1≤n≤N)个辐射源信号的加权最小二乘估计准则:
Anvn=ζn(1≤n≤N)
进一步地,所述步骤3中,需要根据第n(1≤n≤N)个辐射源信号的加权最小二乘估计准则获得权重Wn,相应的表达式为:
进一步地,所述步骤4中,首先需要将每个辐射源un观测方程进行线性转换并合并,获得N个辐射源的加权最小二乘估计准则:
Av=ζ
然后添加第n(1≤n≤N)个辐射源信号到达第m个观测站的距离作为等式约束条件:
dnm=||un-sm||=||Bnmv||(1≤m≤M,1≤n≤N)
最后得到N个辐射源的加权最小二乘解关系式:
s.t.dnm=||un-sm||=||Bnmv||(1≤m≤M,1≤n≤N)
式中Bnm=[O3×(n-1)(M+4)I3 O3×M-sm O3×(N-n)(M+4)],O3×M表示3×M的0矩阵,W=E(ζζT)=blkdiag{W1,W2,…,WN}为加权权重。
进一步地,所述步骤5中,两两辐射源之间的距离范围不等式约束条件可以表示为:
由于第n个辐射源un用v表示有:
un=Jnv=[O3×(n-1)(M+4) I3 O3×(M+1) O3×(N-n)(M+4)]v
式中Jn=[O3×(n-1)(M+4) I3 O3×(M+1) O3×(N-n)(M+4)],则距离不等式约束可以重新表示为:
不等式两边平方且应用aTRa=tr(aaTR)化简可得:
进一步地,所述步骤6中,利用加权最小二乘估计准则、每个辐射源与观测站的距离等式约束、多辐射源之间的距离不等式约束构建定位优化模型,如下式所示:
s.t.dnm=||un-sm||=||Bnmv||
(1≤m≤M,1≤n≤N)
进一步地,所述步骤7中,利用半定松弛方法对步骤6中的定位优化模型进行数值优化,从而确定辐射源位置向量的估计值。由于定位优化模型转换为非凸形式,这里需要定义V=vvT,依据aTRa=tr(aaTR),则步骤6中的定位优化模型可以松弛为:
s.t.tr(FnmV)=V(n-1)(M+4)+3+m,(n-1)(M+4)+3+m
V≥OM+4
VM+4,M+4=1
rank(V)=1
1≤m≤M,1≤n≤N
由于rank(V)=1的约束,上述模型仍是非凸形式。解决此非凸形式的方式是通过去除rank(V)=1的约束来进行凸松弛,从而产生下式所示的凸定位优化模型:
s.t.tr(FnmV)=V(n-1)(M+4)+3+m,(n-1)(M+4)+3+m
V≥OM+4
VM+4,M+4=1
1≤m≤M,1≤n≤N
通过使用内点法可以得到上式凸优化问题的解。从V的定义中可以提取出第n个辐射源位置向量为un=V(n-1)(M+4)+1:(n-1)(M+4)+3,M+4。
与现有技术相比,本发明具有的有益效果:
本发明针对现有的多辐射源定位中,忽视了多辐射源间的相关信息,导致定位精度不理想的问题,提出一种距离约束下联合AOA和RSS的多源协同定位方法。首先,利用多个静止观测站同时获得多个静止辐射源的AOA观测量(包括方位角观测量、仰角观测量)和RSS观测量(即接收信号强度观测量),并引入辅助变量将联合AOA和RSS的观测方程伪线性处理,由此建立加权最小二乘估计准则。然后,利用每个辐射源与观测站的距离等式约束、多辐射源之间的距离不等式约束,并结合加权最小二乘估计准则构建定位优化模型。最后,利用半定松弛方法对此模型进行数值优化,进而得到辐射源位置向量的估计值。本发明利用多(群)辐射源之间存在的距离范围不等式约束信息,增加定位信息对多辐射源进行定位,可以进一步提高对多个静止辐射源的定位精度。
附图说明
图1为本发明实施例一种距离约束下联合AOA和RSS的多源协同定位方法的基本流程图;
图2为AOA观测量示意图;
图3为辐射源间距离约束示意图;
图4为辐射源1位置估计RMSE随标准差δθ的变化曲线;
图5为辐射源2位置估计RMSE随标准差δθ的变化曲线;
图6为辐射源1位置估计RMSE随标准差δβ的变化曲线;
图7为辐射源2位置估计RMSE随标准差δβ的变化曲线;
图8为辐射源1位置估计RMSE随标准差δp的变化曲线;
图9为辐射源2位置估计RMSE随标准差δp的变化曲线。
具体实施方式
下面结合附图和具体的实施例对本发明做进一步的解释说明:
如图1所示,本发明的一种距离约束下联合AOA和RSS的多源协同定位方法,包括:
步骤1:利用空间中放置的M个静止观测站同时获得第n(1≤n≤N)个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角观测量。
步骤2:利用辅助变量依次对第n(1≤n≤N)个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角观测量进行伪线性处理,并建立第n(1≤n≤N)个辐射源信号的加权最小二乘估计准则,其中un表示第n个辐射源的位置向量,dn表示第n个辐射源的距离向量。
步骤3:确定第n(1≤n≤N)个辐射源信号的加权权重。
步骤4:利用N个辐射源的加权最小二乘估计准则和第n(1≤n≤N)个辐射源信号到达第m个观测站的距离等式约束,获得N个辐射源的加权最小二乘解关系式。
步骤5:确定两两辐射源之间的距离不等式约束。
步骤6:利用加权最小二乘估计准则、每个辐射源与观测站的距离等式约束、多辐射源之间的距离不等式约束构建定位优化模型。
步骤7:利用半定松弛方法对步骤6中的定位优化模型进行数值优化,从而确定辐射源位置向量的估计值。
式中分别表示第n个辐射源的三维空间坐标;分别表示第m个观测站的三维空间坐标;分别表示第n个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角的真实值;α表示路径损耗因子,其取值范围通常为2~5;P0为发射功率,其与天线增益和辐射源的能量有关,这里假设其为已知量;dnm为第n个辐射源和第m个观测站之间的距离;εp,nm、εθ,nm、εβ,nm分别表示第n个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角的观测误差,其相互独立且均服从均值为零的高斯分布,方差分别为
AOA观测量示意图如图2所示。
进一步地,所述步骤2中,需要依次对接收信号强度Pnm、方位角θnm、仰角观测量βnm进行伪线性处理,并建立第n(1≤n≤N)个辐射源信号的加权最小二乘估计准则。
首先根据接收信号强度(RSS)观测方程,将其表示为辐射源与观测站之间的距离:
在小噪声水平下,将上式等式右边使用一阶Taylor级数展开并忽略高阶项可得:
则建立加权最小二乘估计准则有:
Ap,nvn=ζp,n(1≤n≤N)
式中Ap,n=[OM×3 IM -λ]∈RM×(M+4),IM表示M×M的单位阵,λ=[λn1 λn2 … λnM]T;vn为辅助变量,dn=[dn1 dn2 … dnM]T;ζp,n=[ζp,n1 ζp,n2 … ζp,nM]T,ζp,nm=(λnmln10/10α)εp,nm。
将上式展开并忽略高阶项可得:
则建立加权最小二乘估计准则有:
Aθ,nvn=ζθ,n(1≤n≤N)
式中
将上式展开并忽略高阶项可得:
Aβ,nvn=ζβ,n(1≤n≤N)
式中
最后将接收信号强度、方位角、仰角伪线性处理得到的加权最小二乘估计准则合并获得第n(1≤n≤N)个辐射源信号的加权最小二乘估计准则:
Anvn=ζn(1≤n≤N)
值得说明的是,λnm、ζp,n、Ap,n、ζp,nm、bθ,nm、Aθ,n、ζθ,n、ζθ,nm、bβ,nm、Aβ,n、ζβ,n、ζβ,nm、An、ζn为表述方便而设置的中间参量,无具体物理含义。
进一步地,所述步骤3中,需要根据第n(1≤n≤N)个辐射源信号的加权最小二乘估计准则获得权重Wn,相应的表达式为:
进一步地,所述步骤4中,首先需要将每个辐射源un观测方程进行线性转换并合并,获得N个辐射源的加权最小二乘估计准则:
Av=ζ
然后添加第n(1≤n≤N)个辐射源信号到达第m个观测站的距离作为等式约束条件:
dnm=||un-sm||=||Bnmv||(1≤m≤M,1≤n≤N)
最后得到N个辐射源的加权最小二乘解关系式:
s.t.dnm=||un-sm||=||Bnmv||(1≤m≤M,1≤n≤N)
式中Bnm=[O3×(n-1)(M+4) I3 O3×M-sm O3×(N-n)(M+4)],O3×M表示3×M的0矩阵,W=E(ζζT)=blkdiag{W1,W2,…,WN}为加权权重。
值得说明的是,A、v、Bnm为表述方便而设置的中间参量,无具体物理含义。
进一步地,所述步骤5中,两两辐射源之间的距离范围不等式约束条件可以表示为:
由于第n个辐射源un用v表示有:
un=Jnv=[O3×(n-1)(M+4) I3 O3×(M+1) O3×(N-n)(M+4)]v
则距离不等式约束可以重新表示为:
不等式两边平方且应用aTRa=tr(aaTR)化简可得:
其中tr()为迹函数。
值得说明的是,aTRa=tr(aaTR)为一个固定表示式,a、R没有具体含义,仅为显示该表达式的具体形式,而本实施例中推导时直接套用该表达式。
辐射源间距离约束示例图如图3所述。
进一步地,所述步骤6中,利用加权最小二乘估计准则、每个辐射源与观测站的距离等式约束、多辐射源之间的距离不等式约束构建定位优化模型,如下式所示:
s.t.dnm=||un-sm||=||Bnmv||
(1≤m≤M,1≤n≤N)
进一步地,所述步骤7中,利用半定松弛方法对步骤6中的定位优化模型进行数值优化,从而确定辐射源位置向量的估计值。由于定位优化模型转换为非凸形式,这里需要定义V=vvT,依据aTRa=tr(aaTR),则步骤6中的定位优化模型可以松弛为:
s.t.tr(FnmV)=V(n-1)(M+4)+3+m,(n-1)(M+4)+3+m
V≥OM+4
VM+4,M+4=1
rank(V)=1
1≤m≤M,1≤n≤N
由于rank(V)=1的约束,上述模型仍是非凸形式。解决此非凸形式的方式是通过去除rank(V)=1的约束来进行凸松弛,从而产生下式所示的凸定位优化模型:
s.t.tr(FnmV)=V(n-1)(M+4)+3+m,(n-1)(M+4)+3+m
V≥OM+4
VM+4,M+4=1
1≤m≤M,1≤n≤N
通过使用内点法可以得到上式凸优化问题的解。从V的定义中可以提取出第n个辐射源位置向量为un=V(n-1)(M+4)+1:(n-1)(M+4)+3,M+4。
为验证本发明效果,进行如下具体示例:
使用6个观测站对2个辐射源进行定位仿真,观测站位置分布如表1所示。
表1观测站的位置坐标(单位:m)
假设辐射源的位置分别为u1=[31.2 52.4 80.2]T(m)、u2=[33.2 54.4 71.2]T(m),此时,两个辐射源之间的距离大约为9.4(m)。假设距离范围约束为d1,2=9.6(m),发射功率P0为-45dBm,路径损耗因子α为4,仿真进行10000次独立的蒙特卡洛实验。下面将本专利公开的定位方法与RLS-SDP定位方法(未利用距离约束)进行比较。
首先将RSS估计误差标准差δp和仰角估计误差标准差δβ分别设置为1.0dB和2.0°,图4和图5分别给出了辐射源1和辐射源2的位置估计RMSE随方位角误差标准差δθ的变化曲线。其次将RSS估计误差标准差δp和方位角误差标准差δθ分别设置为1.0dB和2.0°,图6和图7分别给出了辐射源1和辐射源2的位置估计RMSE随仰角估计误差标准差δβ的变化曲线。最后将方位角误差标准差δθ和仰角估计误差标准差δβ均设置为2.0°,图8和图9分别给出了辐射源1和辐射源2的位置估计RMSE随RSS估计误差标准差δp的变化曲线。从图4~图9中可以看出,相比于未利用距离约束的RLS-SDP定位方法,本专利公开的定位方法具有更高的定位精度,因为本专利公开的定位方法有效地利用了多(群)辐射源之间的距离范围不等式约束信息,并进行了多(群)辐射源协同定位,因而进一步提高对多(群)辐射源的定位精度。
以上所示仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种距离约束下联合AOA和RSS的多源协同定位方法,其特征在于,包括:
步骤1:利用空间中放置的M个静止观测站同时获得第n个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角观测量,其中1≤n≤N,1≤m≤M;
步骤2:利用辅助变量依次对第n个辐射源信号到达第m个观测站的接收信号强度、方位角、仰角观测量进行伪线性处理,并建立第n个辐射源信号的加权最小二乘估计准则,其中un表示第n个辐射源的位置向量,dn表示第n个辐射源的距离向量;
步骤3:确定第n个辐射源信号的加权权重;
步骤4:利用N个辐射源的加权最小二乘估计准则和第n个辐射源信号到达第m个观测站的距离等式约束,获得N个辐射源的加权最小二乘解关系式;
步骤5:确定两两辐射源之间的距离不等式约束;
步骤6:利用加权最小二乘估计准则、每个辐射源与观测站的距离等式约束、多辐射源之间的距离不等式约束构建定位优化模型;
步骤7:利用半定松弛方法对步骤6中的定位优化模型进行数值优化,从而确定辐射源位置向量的估计值。
2.根据权利要求1所述的一种距离约束下联合AOA和RSS的多源协同定位方法,其特征在于,所述步骤1中,当第n个辐射源的位置向量为第m个观测站的位置向量为时,第n个辐射源信号到达第m个观测站的接收信号强度Pnm、方位角θnm、仰角观测量βnm相应的表达式分别为:
3.根据权利要求2所述的一种距离约束下联合AOA和RSS的多源协同定位方法,其特征在于,所述步骤2包括:
首先根据接收信号强度RSS观测方程,将其表示为辐射源与观测站之间的距离:
将上式等式右边使用一阶Taylor级数展开并忽略高阶项可得:
建立加权最小二乘估计准则:
Ap,nvn=ζp,n;1≤n≤N
式中Ap,n=[OM×3 IM -λ]∈RM×(M+4),IM表示M×M的单位阵,λ=[λn1 λn2…λnM]T;vn为辅助变量,dn=[dn1 dn2…dnM]T;ζp,n=[ζp,n1 ζp,n2…ζp,nM]T,ζp,nm=(λnmln10/10α)εp,nm;
将上式展开并忽略高阶项可得:
建立加权最小二乘估计准:
Aθ,nvn=ζθ,n;1≤n≤N
式中
将上式展开并忽略高阶项可得:
建立加权最小二乘估计准则:
Aβ,nvn=ζβ,n;1≤n≤N
式中
最后将接收信号强度、方位角、仰角伪线性处理得到的加权最小二乘估计准则合并获得第n个辐射源信号的加权最小二乘估计准则:
Anvn=ζn;1≤n≤N
5.根据权利要求4所述的一种距离约束下联合AOA和RSS的多源协同定位方法,其特征在于,所述步骤4包括:
首先将每个辐射源观测方程进行线性转换并合并,获得N个辐射源的加权最小二乘估计准则:
Av=ζ
然后添加第n个辐射源信号到达第m个观测站的距离dnm作为等式约束条件:
dnm=||un-sm||=||Bnmv||;1≤m≤M,1≤n≤N
最后得到N个辐射源的加权最小二乘解关系式:
s.t.dnm=||un-sm||=||Bnmv||;1≤m≤M,1≤n≤N
式中Bnm=[O3×(n-1)(M+4) I3 O3×M -sm O3×(N-n)(M+4)],O3×M表示3×M的0矩阵;W=E(ζζT)=blkdiag{W1,W2,…,WN}为加权权重。
8.根据权利要求7所述的一种距离约束下联合AOA和RSS的多源协同定位方法,其特征在于,所述步骤7中,利用半定松弛方法对步骤6中的定位优化模型进行数值优化,从而确定辐射源位置向量的估计值,包括:
定义V=vvT,依据aTRa=tr(aaTR),则步骤6中的定位优化模型可以松弛为:
s.t.tr(FnmV)=V(n-1)(M+4)+3+m,(n-1)(M+4)+3+m
VM+4,M+4=1
rank(V)=1
1≤m≤M,1≤n≤N
通过去除rank(V)=1的约束来进行凸松弛,从而产生下式所示的凸定位优化模型:
s.t.tr(FnmV)=V(n-1)(M+4)+3+m,(n-1)(M+4)+3+m
VM+4,M+4=1
1≤m≤M,1≤n≤N
通过使用内点法得到上式凸优化问题的解,从V的定义中提取出第n个辐射源位置向量为un=V(n-1)(M+4)+1:(n-1)(M+4)+3,M+4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211014219.8A CN115397015B (zh) | 2022-08-23 | 2022-08-23 | 一种距离约束下联合aoa和rss的多源协同定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211014219.8A CN115397015B (zh) | 2022-08-23 | 2022-08-23 | 一种距离约束下联合aoa和rss的多源协同定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115397015A true CN115397015A (zh) | 2022-11-25 |
CN115397015B CN115397015B (zh) | 2024-05-10 |
Family
ID=84121111
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211014219.8A Active CN115397015B (zh) | 2022-08-23 | 2022-08-23 | 一种距离约束下联合aoa和rss的多源协同定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115397015B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109342993A (zh) * | 2018-09-11 | 2019-02-15 | 宁波大学 | 基于RSS-AoA混合测量的无线传感器网络目标定位方法 |
CN110662163A (zh) * | 2019-08-23 | 2020-01-07 | 宁波大学 | 基于rss和aoa的三维无线传感网络协作定位方法 |
US20200205011A1 (en) * | 2018-12-20 | 2020-06-25 | Slavisa Tomic | Apparatus and method for 3-D network localization |
KR102248758B1 (ko) * | 2020-05-18 | 2021-05-07 | 국방과학연구소 | 무선 센서 네트워크에서 Approximated Weighted Least Square (WLS)를 사용한 RSS/AOA 측정치 혼합 표적 위치 추정 기법 |
CN112986906A (zh) * | 2021-02-23 | 2021-06-18 | 中国计量大学上虞高等研究院有限公司 | 一种半正定规划的rss-toa联合定位方法 |
CN113030853A (zh) * | 2021-03-07 | 2021-06-25 | 中国人民解放军陆军工程大学 | 基于rss和aoa联合测量的多辐射源无源定位方法 |
-
2022
- 2022-08-23 CN CN202211014219.8A patent/CN115397015B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109342993A (zh) * | 2018-09-11 | 2019-02-15 | 宁波大学 | 基于RSS-AoA混合测量的无线传感器网络目标定位方法 |
US20200205011A1 (en) * | 2018-12-20 | 2020-06-25 | Slavisa Tomic | Apparatus and method for 3-D network localization |
CN110662163A (zh) * | 2019-08-23 | 2020-01-07 | 宁波大学 | 基于rss和aoa的三维无线传感网络协作定位方法 |
KR102248758B1 (ko) * | 2020-05-18 | 2021-05-07 | 국방과학연구소 | 무선 센서 네트워크에서 Approximated Weighted Least Square (WLS)를 사용한 RSS/AOA 측정치 혼합 표적 위치 추정 기법 |
CN112986906A (zh) * | 2021-02-23 | 2021-06-18 | 中国计量大学上虞高等研究院有限公司 | 一种半正定规划的rss-toa联合定位方法 |
CN113030853A (zh) * | 2021-03-07 | 2021-06-25 | 中国人民解放军陆军工程大学 | 基于rss和aoa联合测量的多辐射源无源定位方法 |
Non-Patent Citations (3)
Title |
---|
HENGNIAN QI; LUFENG MO; XIAOPING WU: "SDP Relaxation Methods for RSS/AOA-Based Localization in Sensor Networks", IEEE ACCESS, vol. 8, 18 March 2020 (2020-03-18), XP011780408, DOI: 10.1109/ACCESS.2020.2981639 * |
SHENGMING CHANG; YOUMING LI; XINJIE YANG; HUI WANG; WENFEI HU FACULTY OF ELECTRICAL ENGINEERING AND COMPUTER SCIENCE, NINGBO UNIVE: "A Novel Localization Method Based on RSS-AOA Combined Measurements by Using Polarized Identity", IEEE SENSORS JOURNAL, vol. 19, no. 4, 15 February 2019 (2019-02-15), XP011706543, DOI: 10.1109/JSEN.2018.2881892 * |
张瑞玲;: "一种基于接收信号强度和到达角的三维组合定位方法", 舰船电子对抗, no. 02, 25 April 2019 (2019-04-25) * |
Also Published As
Publication number | Publication date |
---|---|
CN115397015B (zh) | 2024-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107132505B (zh) | 直达与非直达混合场景中的多目标直接定位方法 | |
CN111308530B (zh) | 一种基于二维波达方向的短波多站和单星协同直接定位方法 | |
CN109975749B (zh) | 一种校正源存在条件下的短波单站直接定位方法 | |
CN108693511B (zh) | 时分复用mimo雷达的运动目标角度计算方法 | |
CN103869347A (zh) | 一种北斗卫星导航定位系统 | |
CN103412286A (zh) | 基于mimo雷达的发射极化优化doa估计方法 | |
CN105487072A (zh) | 一种基于t2/r的时差与多普勒频移联合定位方法及系统 | |
Li et al. | Review on positioning technology of wireless sensor networks | |
CN111199281B (zh) | 基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法 | |
CN116148829A (zh) | 一种基于非合作外辐射源的多基雷达定位方法 | |
Meles et al. | Drone localization based on 3D-AoA signal measurements | |
CN113419219B (zh) | 基于空域特征认知的外辐射源雷达同频干扰级联相消方法 | |
Nijsure et al. | WSN-UAV monitoring system with collaborative beamforming and ADS-B based multilateration | |
CN115397015B (zh) | 一种距离约束下联合aoa和rss的多源协同定位方法 | |
CN109991564B (zh) | 基于神经网络的短波单站定位结果纠偏方法 | |
CN110208741B (zh) | 一种基于多圆阵测相的超视距单目标直接定位方法 | |
CN107219498B (zh) | 基于mml-pmht的多基站单频网无源协同定位方法 | |
CN111079929B (zh) | 基于神经计算的短波单站多目标地理坐标快速估计方法 | |
Zhao et al. | Active phased array radar-based 2D beamspace MUSIC channel estimation for an integrated radar and communication system | |
CN108646238A (zh) | 一种基于副瓣对消系数映射的干扰源跟踪方法 | |
CN104459680A (zh) | Mimo雷达估计目标方向的快速方法 | |
Chen et al. | Compressed sensing-based DOA and DOD estimation in bistatic co-prime MIMO arrays | |
Wu et al. | A simple and effective iterative algorithm for elliptic localization | |
Cao et al. | Direct self-position awareness based on array-sensing multiple source data fitting | |
CN117202343B (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 |