CN106154225A - 基于波达方向歧义消除的定位方法及装置 - Google Patents

基于波达方向歧义消除的定位方法及装置 Download PDF

Info

Publication number
CN106154225A
CN106154225A CN201610705725.XA CN201610705725A CN106154225A CN 106154225 A CN106154225 A CN 106154225A CN 201610705725 A CN201610705725 A CN 201610705725A CN 106154225 A CN106154225 A CN 106154225A
Authority
CN
China
Prior art keywords
sensor node
ray
combination
node
source node
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
Application number
CN201610705725.XA
Other languages
English (en)
Other versions
CN106154225B (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.)
Sichuan University
Original Assignee
Sichuan 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 Sichuan University filed Critical Sichuan University
Priority to CN201610705725.XA priority Critical patent/CN106154225B/zh
Publication of CN106154225A publication Critical patent/CN106154225A/zh
Application granted granted Critical
Publication of CN106154225B publication Critical patent/CN106154225B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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/06Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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/0205Details
    • G01S5/021Calibration, monitoring or correction

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)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明提供的基于波达方向歧义消除的定位方法及装置,通过对获取的各个传感器节点上的数据进行空间谱估计,得到源节点相对各个传感器节点的俯仰角及方位角的估计值。根据预设的代价函数和得到的俯仰角及方位角的估计值,确定由各个传感器节点的阵列中心指向源节点的射线组合,最后对该射线组合进行最小二乘法估计,得到源节点的位置信息。该方法及装置无需复杂的硬件设备作为基础,仅根据传感器节点上的信息融合来进行三维空间中波达方向估计的歧义消除,实现无线传感器网络中的目标准确定。

Description

基于波达方向歧义消除的定位方法及装置
技术领域
本发明涉及无线传感器网络定位技术领域,具体而言,涉及一种基于波达方向歧义消除的定位方法及装置。
背景技术
近年来,无线传感器网络被广泛地应用于民用及军用场景中,特别是物联网,其中有诸多应用场景都需要获取目标的位置信息。所以作为无线传感器网络的一种主要应用,对未知节点或者自然环境中的目标进行定位引起了广泛地关注。传统的定位算法可以根据测量手段而分为几类,包括全球定位系统(GPS),基于接受信号强度(RSS)的方法,基于到达时间(TOA)或者到达时间差(TDOA)的方法,基于波达方向(DOA)的方法,以及综合上述手段的方法。
GPS定位方法需要昂贵的硬件设备作为基础,并且在封闭的环境如室内或隧道会失效。基于RSS的方法需要已知环境的信道传输模型,这一点对其应用范围产生极大限制。而基于TOA/TDOA的方法则需要网络中各节点之间苛刻的时间同步,同样限制其只适用于部分特殊的应用场景中。与上述方法不同的是,基于DOA的定位方法由于其简单、节能、以及较好的鲁棒性而可以更好地适用于无线传感器网络应用场景中。
但是,由于两个相同的关于平面阵列对称的信号会产生完全相同的空间频率响应即阵列流形,所以平面阵列上的DOA估计本质上会产生歧义,导致其定位失准。如何消除DOA估计中的歧义是业界目前亟待解决的难题。
发明内容
有鉴于此,本发明的目的在于提供一种基于波达方向歧义消除的定位方法及装置。
本发明较佳实施例提供一种基于波达方向歧义消除的定位方法,该方法包括:
获取每一个预设传感器节点所接收到的源节点发射的信号,其中每个所述传感器节点包括预设数目元素的平面天线阵列;
通过空间谱估计算法分别计算每个所述传感器节点接收到的信号的MUSIC谱中两个谱峰所对应的角度参量,所述角度参量包括该传感器节点相对所述源节点的方位角及俯仰角;
根据计算得到的各个角度参量,分别确定由该角度参量所对应的传感器节点的平面天线阵列中心指向所述源节点的射线;
对所有所述射线按照其对应的传感器节点进行排列组合,并根据预设代价函数计算每种组合中各个射线间的紧密程度;及
对具有最小代价函数值的射线组合进行最小二乘法估计,得到所述源节点的位置坐标。
本发明另一较佳实施例提供一种基于波达方向歧义消除的定位装置,该装置包括:
数据获取模块,用于获取每一个预设传感器节点所接收到的源节点发射的信号,其中每个所述传感器节点包括预设数目元素的平面天线阵列;
角度参量计算模块,用于通过空间谱估计算法分别计算每个所述传感器节点接收到的信号的MUSIC谱中两个谱峰所对应的角度参量,所述角度参量包括该传感器节点相对所述源节点的方位角及俯仰角;
射线确定模块,用于根据计算得到的各个角度参量,分别确定由该角度参量所对应的传感器节点的平面天线阵列中心指向所述源节点的射线;
代价函数计算模块,用于对所有所述射线按照其对应的传感器节点进行排列组合,并根据预设代价函数计算每种组合中各个射线间的紧密程度;及
位置坐标估计模块,用于对具有最小代价函数值的射线组合进行最小二乘法估计,得到所述源节点的位置坐标。
本发明提供的基于波达方向歧义消除的定位方法及装置,无需复杂的硬件设备作为基础,仅根据传感器节点上的信息融合来进行三维空间中波达方向估计的歧义消除,实现无线传感器网络中的目标准确定位。另外,通过增加每个传感器节点中的天线阵列的阵元数目以提高阵列的灵敏度,或者增大传感器节点的数目均可以有效降低本实施例提供的定位方法的误差,进一步确保实现目标的精准定位。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1是本发明较佳实施例提供的一种无线传感器网络源节点定位系统的示意性结构框图;
图2是本发明较佳实施例提供的一种基于波达方向歧义消除的定位方法的流程图;
图3是本发明较佳实施例提供的一种无线传感器网络源节点定位系统中源节点与传感器节点的相对位置分布示意图;
图4是本发明较佳实施例提供的一种基于波达方向歧义消除的定位装置的功能模块框图;
图5是本发明较佳实施例提供的模拟仿真实验中场景一的传感器节点与源节点的相对位置分布示意图;
图6为使用图5所示的L型天线阵列进行100次蒙特卡洛实验得到的平均绝对误差与信噪比之间的曲线关系;
图7为使用图5所示的L型天线阵列进行100次蒙特卡洛实验得到的平均相对误差与信噪比之间的曲线关系;
图8为使用模拟仿真实验中场景二的4×4方型天线阵列进行100次蒙特卡洛实验得到的平均绝对误差与信噪比之间的曲线关系;
图9为使用模拟仿真实验中场景二的4×4方型天线阵列进行100次蒙特卡洛实验得到的平均相对误差与信噪比之间的曲线关系;
图10为模拟仿真实验的场景三中进行100次蒙特卡洛实验得到的平均绝对误差与传感器节点个数之间的曲线关系;
图11为模拟仿真实验的场景三中进行100次蒙特卡洛实验得到的平均相对误差与传感器节点个数之间的曲线关系;
图12为模拟仿真实验的场景三中进行100次蒙特卡洛实验得到的平均计算时间与传感器节点个数之间的曲线关系。
附图标记:
无线传感器网络源节点定位系统 100
存储器 200
处理器 300
传感器节点 400
基于波达方向歧义消除的定位装置 500
数据获取模块 510
角度参量计算模块 520
射线确定模块 530
代价函数计算模块 540
位置坐标估计模块 550
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,是本发明较佳实施例提供的一种无线传感器网络源节点定位系统100的示意性结构框图。该系统包括存储器200、处理器300、多个传感器节点400以及基于波达方向歧义消除的定位装置500。
所述存储器200与所述处理器300之间直接或间接地电性连接,以实现数据的传输或交互。例如,两者通过通讯总线或信号线实现电性连接。所述多个传感器节点400分别与所述处理器300无线通信连接,以将感测的数据发送至处理器300进行处理。每一个所述传感器节点400包括预设数目元素的平面天线阵列。所述基于波达方向歧义消除的定位装置500包括至少一个可以软件或固件(firmware)的形式存储于所述存储器200中的软件功能模块。所述处理器300用于执行存储器200中存储的可执行模块,例如所述基于波达方向歧义消除的定位装置500包括的软件功能模块或计算机程序。
其中,所述存储器200可以是,但不限于,随机存取存储器、只读存储器、可编程只读存储器、可擦除只读存储器、电可擦除只读存储器等。其中,存储器200用于存储程序,所述处理器300在接收到执行指令后,执行所述程序,下述本发明任一实施例揭示的流过程定义的方法可以应用于处理器300中,或者由处理器300实现。
所述处理器300可能是一种集成电路芯片,具有信号的处理能力。上述的处理器300可以是通用处理器,包括中央处理器、网络处理器等;还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件,可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。
如图2所示,是本发明较佳实施例提供的一种基于波达方向歧义消除的定位方法的流程图,应用于图1所示的无线传感器网络源节点定位系统100中。所应说明的是,本发明提供的方法不以图2及以下所示的具体顺序为限制。下面将对图2所示的具体步骤进行详细描述。
步骤S101,获取每一个预设传感器节点400所接收到的源节点发射的信号,其中每个所述传感器节点400包括预设数目元素的平面天线阵列。
如图3所示,本实施例中,不失一般性地,假设所述无线传感器网络源节点定位系统100包括L个传感器节点400,该L个传感器节点400设置于三维空间中,每个传感器节点400由包含M个元素的平面天线阵列组成。不失一般性地,假设每个平面天线阵列均平行于x-o-y平面。为方便描述,将第l个传感器节点400的第m个天线记作pm,l=[xm,l,ym,l,zm,l]T,其中,m=1,2,...,M,l=1,2,...,L。则第l个传感器节点400的天线阵列中心为其中l=1,2,..L.,。
位于ps=[xs,ys,zs]T的源节点发射出频率为f的窄带信号s(t)。假设对于每个传感器节点400而言,源信号s(t)均为远场信号。如图3所示,源信号对于第l个传感器节点400的俯仰角为θl∈[0,π],方位角为φl∈[0,2π)。如此,源信号s(t)对于第l个传感器节点400的方向向量可以表示为:
ul=[sinθlcosφl,sinθlcosφl,cosφl]T
到达第l个传感器节点400的第m个天线的信号相对于原始信号的时延为其中c是信号在介质中的传播速度。所以第l个传感器节点400接收到的信号可以表示为:
xl(t)=a(φll)s(t)+nl(t)
其中λ是信号的波长,nl(t)是空间协方差为的高斯白噪声,IM表示秩为M的单位矩阵。
对第l个传感器节点400接收到的信号以预设采样频率进行采样,得到采样数据:
Xl=a(φll)S+Nl
其中,Xl=[xl(Ts),...,xl(NTs)],S=[s(Ts),...,s(NTs)],Nl=[nl(Ts),...,nl(NTs)]。Ts表示采样周期,N表示样本数。
由此,本实施例中,获取每一个预设传感器节点400所接收到的源节点发射的信号即是获取Xl,l=1,2,...,L。
步骤S103,通过空间谱估计算法分别计算每个所述传感器节点400接收到的信号的MUSIC谱中两个谱峰所对应的角度参量,所述角度参量包括该传感器节点400相对所述源节点的方位角及俯仰角。
本实施例中,通过MUSIC谱估计算法估计所述源信号相对每个传感器节点400的俯仰角及方位角。具体地,通过采集到的数据Xl,可以得到第l个传感器节点400所接收到的数据的协方差矩阵为:
R l = 1 N X l X l H
其中(·)H表示矩阵的共轭转置。对矩阵Rl做特征值分解得到:
R l = U l Σ l U l H = e 1 , l v 1 , l v 1 , l H + U l ( n u l l ) Σ l ( n u l l ) [ U l ( n u l l ) ] H
其中em,l表示特征向量,并且假设e1,l≥e2,l≥...≥eM,l,vm,l是其对应的特征向量。Σl是对角元素为e1,l,e2,l,...,eM,l的对角矩阵,Ul=[v1,l,...,vm,l]。是对角元素为e2,l,...,eM,l的对角矩阵,表示信号的噪声子空间。
由此,第l个传感器节点400的MUSIC谱为
p l ( φ , θ ) = { [ a ( φ l , θ l ) ] H U l ( n u l l ) [ U l ( n u l l ) ] H a ( φ l , θ l ) } - 1
则通过MUSIC谱估计算法得到的源信号相对于第l个传感器节点400的方位角及俯仰角的估计值为:
( φ ^ l , θ ^ l ) = m a x φ , θ { p l ( φ , θ ) } , l = 1 , 2 , ... , L
众所周知,平面天线阵列的波达方向估计会产生歧义,即在无噪声环境下,第l个传感器节点400的MUSIC谱会在(φll)及(φl,π-θl)分别出现谱峰,且它们关于平面z=zl对称。这是由于所以a(φl,π-θl)及a(φll)都在Rl的信号子空间中,并且都正交于噪声子空间
因此,计算得到的每个传感器节点400的MUSIC谱中两个谱峰分别对应的角度参量估计值为其中为第l个传感器节点400的MUSIC谱中两个谱峰分别对应的方位角,为第l个传感器节点400的MUSIC谱中两个谱峰分别对应的俯仰角,l=1,2,3,...,L。
步骤S105,根据计算得到的各个角度参量,分别确定由该角度参量所对应的传感器节点400的平面天线阵列中心指向所述源节点的射线。
本实施例中,根据计算得到的角度参量估计值分别得到所述源节点的波达方向估计的方向向量其中计算所述方向向量的公式为:
u ^ l ( 0 ) = [ sin θ ^ l ( 0 ) cosφ l ( 0 ) , sin θ ^ l ( 0 ) cos φ ^ l ( 0 ) , cos θ ^ l ( 0 ) ] T u ^ l ( 1 ) = [ sin θ ^ l ( 1 ) cos φ ^ l ( 1 ) , sin θ ^ l ( 1 ) cos φ ^ l ( 1 ) , cos θ ^ l ( 1 ) ] T , l = 1 , 2 , 3 , ... , L .
则由第l个传感器节点400的平面天线阵列中心pl及所述方向向量确定的两条指向所述源节点的射线为
对于空间中的L个传感器节点400,共可以确定2L条射线,表示为其中其中有L条射线共同指向所述源节点,另外L条射线是由波达方向歧义所产生的。
步骤S107,对所有所述射线按照其对应的传感器节点400进行排列组合,并根据预设代价函数计算每种组合中各个射线间的紧密程度。
本实施例中,对2L条射线按照其对应的传感器节点400进行排列组合得到2L种射线组合,其中每种组合包括L条射线,该L条射线分别对应不同的传感器节点400。歧义消除的目标就是找出指向源节点的L条射线,由于每个传感器节点400上的两条射线中只有一条指向源节点,所以问题等价于从2L种射线组合中找出一种正确的L条射线的组合方式。
对每种射线组合分别进行长度为L位的二进制编码,得到对应的编码向量b(m),m=1,2,...,2L,其中编码向量b(m)中的第l位表示对应第l个传感器节点400的两条射线的其中一个,即:
为了表示每种组合中L条射线之间的紧密程度,定义代价函数:
c ( b ( m ) ) = ∑ i = 1 L = 1 ∑ j = i + 1 L d i , j ( raay ) ( b ( m ) )
其中表示射线上任意一点与射线上任意一点之间的最短距离,表示为,
d i , j ( r a y ) ( b ( m ) ) = min | | ( p i + k i u ^ i [ b ( m ) ] i ) - ( p j + k j u ^ j [ b ( m ) ] j ) | | , i ≠ j , k i ≥ 0 , k j ≥ 0
上式中,约束条件ki≥0,kj≥0保证了是两条射线上任意两点的距离。
所以,指向所述源节点的L条射线的组合可以通过最小化代价函数得到,即:
b opt = min b m { c ( b ( m ) ) }
最终,指向源节点的L条射线可以表示为
步骤S109,对具有最小代价函数值的射线组合进行最小二乘法估计,得到所述源节点的位置坐标。
本实施例中,在消除了波达方向估计的歧义后,根据得到的L条指向源节点的射线的组合,采用最小二乘法来进行源节点定位。
对于空间中任意一点p=[x,y,z]T,到射线的投影为则点p到该射线的距离为
本实施例中,对所有的传感器节点400取相同的权重,则源节点位置的最小二乘估计为:
p ^ s = ( L I - V ^ V ^ T ) - 1 A W
其中,I为单位矩阵,W=[1,...,1]T,W长度为L,
如图4所示,是本发明较佳实施例提供的所述基于波达方向歧义消除的定位装置500的功能模块框图。该基于波达方向歧义消除的定位装置500包括数据获取模块510、角度参量计算模块520、射线确定模块530、代价函数计算模块540及位置坐标估计模块550。下面将对图4所示的功能模块进行详细描述。
所述数据获取模块510,用于获取每一个预设传感器节点400所接收到的源节点发射的信号,其中每个所述传感器节点400包括预设数目元素的平面天线阵列。具体地,该数据获取模块510可用于执行图2所示的步骤S101,具体的操作方法可参照步骤S101的详细描述。
所述角度参量计算模块520,用于通过空间谱估计算法分别计算每个所述传感器节点400接收到的信号的MUSIC谱中两个谱峰所对应的角度参量,所述角度参量包括该传感器节点400相对所述源节点的方位角及俯仰角。具体地,该角度参量计算模块520可用于执行图2所示的步骤S103,具体的操作方法可参照步骤S103的详细描述。
所述射线确定模块530,用于根据计算得到的各个角度参量,分别确定由该角度参量所对应的传感器节点400的平面天线阵列中心指向所述源节点的射线。具体地,该射线确定模块530可用于执行图2所示的步骤S105,具体的操作方法可参照步骤S105的详细描述。
所述代价函数计算模块540,用于对所有所述射线按照其对应的传感器节点400进行排列组合,并根据预设代价函数计算每种组合中各个射线间的紧密程度。具体地,该代价函数计算模块540可用于执行图2所示的步骤S107,具体的操作方法可参照步骤S107的详细描述。
所述位置坐标估计模块550,用于对具有最小代价函数值的射线组合进行最小二乘法估计,得到所述源节点的位置坐标。具体地,该位置坐标估计模块550可用于执行图2所示的步骤S109,具体的操作方法可参照步骤S109的详细描述。
下面将对本发明提供的方法选取几种不同的场景进行模拟仿真。在所有的仿真场景中,均由单一信号源向周围环境发射频率为f=4GHz,波长为λ=0.075m的单频信号s(t)=sin(ωt)。信号传播过程中,存在均值为0方差为1的高斯白噪声。采样频率为fs=4GHz,采样周期为N=800。
为了衡量算法的表现,定义源节点位置估计的绝对误差为源节点位置估计的相对误差为
1)场景一
如图5所示,三维空间中设有L=4个传感器节点400,分别为p1=[25.0125,25.0125,25]T,p2=[80.0125,4.0125,60]T p3=[37.0125,75.0125,27]T,p4=[98.0125,25.0125,5]T。每个传感器节点400包括3个元素的L型天线阵列,其对应的坐标分别为p1,1=[25,25,25]T,p2,1=[25.0375,25,25]T,p3,1=[25,25.0375,25]T,p1,2=[80,4,60]T,p2,2=[80.0375,4,60]T,p3,2=[80,4.0375,60]T,p1,3=[37,75,27]T,p2,3=[37.0735,75,27]T,p3,3=[37,75.0735,27]T,p1,4=[98,25,5]T p2,4=[98.0735,25,5]T,p3,4=[98,25.0735,5]T。源节点位于ps=[50,50,50]T。图中“source node”表示源节点,“sensor node”表示传感器节点。
图6表示使用L型天线阵列进行100次蒙特卡洛实验中平均绝对误差re与信噪比(SNR)之间的曲线关系。图7表示使用L型天线阵列进行100次蒙特卡洛实验中平均相对误差与信噪比(SNR)之间的曲线关系。
2)场景二
源节点及所有传感器节点400的空间位置与场景一相同,每个传感器节点400都是由一个包含16个元素的4×4方型天线阵列组成,且天线阵列内部沿着坐标轴的间隔距离为
图8表示使用4×4方型天线阵列进行100次蒙特卡洛实验中平均绝对误差re与信噪比(SNR)之间的曲线关系。图9表示使用4×4方型天线阵列进行100次蒙特卡洛实验中平均相对误差与信噪比(SNR)之间的曲线关系。
3)场景三
假设一个100×100×100的三维立体空间,源节点位于正中心的位置[50,50,50],共有L个传感器节点400均匀分布在整个空间中,每个传感器节点400都是由一个和场景一相同的L型天线阵列组成。
图10表示进行100次蒙特卡洛实验中平均绝对误差re与传感器节点400个数L之间的关系。图11表示进行100次蒙特卡洛实验中平均相对误差与传感器节点400个数L之间的关系。图12表示进行100次蒙特卡洛实验中平均计算时间与传感器节点400个数L之间的关系。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
需要说明的是,本文中术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求所述的保护范围为准。

Claims (10)

1.一种基于波达方向歧义消除的定位方法,其特征在于,该方法包括:
获取每一个预设传感器节点所接收到的源节点发射的信号,其中每个所述传感器节点包括预设数目元素的平面天线阵列;
通过空间谱估计算法分别计算每个所述传感器节点接收到的信号的MUSIC谱中两个谱峰所对应的角度参量,所述角度参量包括该传感器节点相对所述源节点的方位角及俯仰角;
根据计算得到的各个角度参量,分别确定由该角度参量所对应的传感器节点的平面天线阵列中心指向所述源节点的射线;
对所有所述射线按照其对应的传感器节点进行排列组合,并根据预设代价函数计算每种组合中各个射线间的紧密程度;及
对具有最小代价函数值的射线组合进行最小二乘法估计,得到所述源节点的位置坐标。
2.根据权利要求1所述的基于波达方向歧义消除的定位方法,其特征在于,所述通过空间谱估计算法分别计算每个所述传感器节点接收到的信号的MUSIC谱中两个谱峰所对应的角度参量的步骤包括:
分别对每个所述传感器节点接收到的信号进行采样,计算采样后的数据的协方差矩阵,并对所述协方差矩阵进行特征值分解得到相互正交的噪声子空间及信号子空间;
根据得到的所述相互正交的噪声子空间及信号子空间,对应构建所述传感器节点的MUSIC谱;
计算得到每个传感器节点的MUSIC谱中两个谱峰分别对应的角度参量其中为第l个传感器节点的MUSIC谱中两个谱峰分别对应的方位角,为第l个传感器节点的MUSIC谱中两个谱峰分别对应的俯仰角,l=1,2,3,...,L,L表示传感器节点的总数目。
3.根据权利要求2所述的基于波达方向歧义消除的定位方法,其特征在于,所述根据计算得到的各个角度参量,分别确定由该角度参量所对应的传感器节点的平面天线阵列中心指向所述源节点的射线的步骤包括:
根据所述角度参量分别对应计算得到所述源节点的波达方向估计的方向向量
由第l个传感器节点的平面天线阵列中心pl及所述方向向量确定两条指向所述源节点的射线其中,计算所述方向向量的公式为:
u ^ l ( 0 ) = [ sin θ ^ l ( 0 ) cos φ ^ l ( 0 ) , sin θ ^ l ( 0 ) cos φ ^ l ( 0 ) , cos θ ^ l ( 0 ) ] T u ^ l ( 1 ) = [ sin θ ^ l ( 1 ) cos φ ^ l ( 1 ) , sin θ ^ l ( 1 ) cos φ ^ l ( 1 ) , cos θ ^ l ( 1 ) ] T , l = 1 , 2 , 3 , ... , L .
4.根据权利要求3所述的基于波达方向歧义消除的定位方法,其特征在于,所述对所有所述射线按照其对应的传感器节点进行排列组合,并根据预设代价函数计算每种组合中各个射线间的紧密程度的步骤包括:
对2L条射线按照其对应的传感器节点进行排列组合得到2L种射线组合,其中每种组合包括L条射线,该L条射线分别对应不同的传感器节点;
对每种射线组合分别进行长度为L位的二进制编码,得到对应的编码向量b(m),m=1,2,...,2L,其中编码向量b(m)中的第l位表示对应第l个传感器节点的两条射线的其中一个,即:
其中m=1,2,...,2L,l=1,2,...,L;
根据代价函数计算每组射线组合中L条射线的紧密程度,其中表示射线上任意一点与射线上任意一点之间的最短距离,表示为,
d i , j ( r a y ) ( b ( m ) ) = min | | ( p i + k i u ^ i [ b ( m ) ] i ) - ( p j + k j u ^ j [ b ( m ) ] j ) | | , i ≠ j , k i ≥ 0 , k j ≥ 0 ;
获取具有最小代价函数值的射线组合bopt,其中opt∈{1,2,...,2L}。
5.根据权利要求4所述的基于波达方向歧义消除的定位方法,其特征在于,所述对具有最小代价函数值的射线组合进行最小二乘法估计,得到所述源节点的位置坐标的步骤包括:
根据计算式:
p ^ s = ( L I - V ^ V ^ T ) - 1 A W
对射线组合bopt进行最小二乘法估计,得到所述源节点的位置坐标其中,I为单位矩阵,W=[1,...,1]T,W长度为L,以及
6.一种基于波达方向歧义消除的定位装置,其特征在于,该装置包括:
数据获取模块,用于获取每一个预设传感器节点所接收到的源节点发射的信号,其中每个所述传感器节点包括预设数目元素的平面天线阵列;
角度参量计算模块,用于通过空间谱估计算法分别计算每个所述传感器节点接收到的信号的MUSIC谱中两个谱峰所对应的角度参量,所述角度参量包括该传感器节点相对所述源节点的方位角及俯仰角;
射线确定模块,用于根据计算得到的各个角度参量,分别确定由该角度参量所对应的传感器节点的平面天线阵列中心指向所述源节点的射线;
代价函数计算模块,用于对所有所述射线按照其对应的传感器节点进行排列组合,并根据预设代价函数计算每种组合中各个射线间的紧密程度;及
位置坐标估计模块,用于对具有最小代价函数值的射线组合进行最小二乘法估计,得到所述源节点的位置坐标。
7.根据权利要求6所述的基于波达方向歧义消除的定位装置,其特征在于,所述角度参量计算模块通过空间谱估计算法分别计算每个所述传感器节点接收到的信号的MUSIC谱中两个谱峰所对应的角度参量的实现方式包括:
分别对每个所述传感器节点接收到的信号进行采样,计算采样后的数据的协方差矩阵,并对所述协方差矩阵进行特征值分解得到相互正交的噪声子空间及信号子空间;
根据得到的所述相互正交的噪声子空间及信号子空间,对应构建所述传感器节点的MUSIC谱;
计算得到每个传感器节点的MUSIC谱中两个谱峰分别对应的角度参量其中为第l个传感器节点的MUSIC谱中两个谱峰分别对应的方位角,为第l个传感器节点的MUSIC谱中两个谱峰分别对应的俯仰角,l=1,2,3,...,L,L表示传感器节点的总数目。
8.根据权利要求7所述的基于波达方向歧义消除的定位装置,其特征在于,所述射线确定模块根据计算得到的各个角度参量,分别确定由该角度参量所对应的传感器节点的平面天线阵列中心指向所述源节点的射线的方式包括:
根据所述角度参量分别对应得到所述源节点的波达方向估计的方向向量
由第l个传感器节点的平面天线阵列中心pl及所述方向向量确定两条指向所述源节点的射线其中计算所述方向向量的公式为:
u ^ l ( 0 ) = [ sin θ ^ l ( 0 ) cos φ ^ l ( 0 ) , sin θ ^ l ( 0 ) cos φ ^ l ( 0 ) , cos θ ^ l ( 0 ) ] T u ^ l ( 1 ) = [ sin θ ^ l ( 1 ) cos φ ^ l ( 1 ) , sin θ ^ l ( 1 ) cos φ ^ l ( 1 ) , cos θ ^ l ( 1 ) ] T , l = 1 , 2 , 3 , ... , L .
9.根据权利要求8所述的基于波达方向歧义消除的定位装置,其特征在于,所述代价函数计算模块对所有所述射线按照其对应的传感器节点进行排列组合,并根据预设代价函数计算每种组合中各个射线间的紧密程度的实现方式包括:
对2L条射线按照其对应的传感器节点进行排列组合得到2L种射线组合,其中每种组合包括L条射线,该L条射线分别对应不同的传感器节点;
对每种射线组合分别进行长度为L位的二进制编码,得到对应的编码向量b(m),m=1,2,...,2L,其中编码向量b(m)中的第l位表示对应第l个传感器节点的两条射线的其中一个,即:
其中m=1,2,...,2L,l=1,2,...,L;
根据代价函数计算每组射线组合中L条射线的紧密程度,其中表示射线上任意一点与射线上任意一点之间的最短距离,表示为:
d i , j ( r a y ) ( b ( m ) ) = min | | ( p i + k i u ^ i [ b ( m ) ] i ) - ( p j + k j u ^ j [ b ( m ) ] j ) | | , i ≠ j , k i ≥ 0 , k j ≥ 0 ;
获取具有最小代价函数值的射线组合bopt,其中opt∈{1,2,...,2L}。
10.根据权利要求9所述的基于波达方向歧义消除的定位装置,其特征在于,所述位置坐标估计模块对具有最小代价函数值的射线组合进行最小二乘法估计,得到所述源节点的位置坐标的实现方式包括:
根据计算式:
p ^ s = ( L I - V ^ V ^ T ) - 1 A W
对射线组合bopt进行最小二乘法估计,得到所述源节点的位置坐标其中,I为单位矩阵,W=[1,...,1]T,W长度为L,以及
CN201610705725.XA 2016-08-22 2016-08-22 基于波达方向歧义消除的定位方法及装置 Active CN106154225B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610705725.XA CN106154225B (zh) 2016-08-22 2016-08-22 基于波达方向歧义消除的定位方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610705725.XA CN106154225B (zh) 2016-08-22 2016-08-22 基于波达方向歧义消除的定位方法及装置

Publications (2)

Publication Number Publication Date
CN106154225A true CN106154225A (zh) 2016-11-23
CN106154225B CN106154225B (zh) 2018-05-04

Family

ID=57343169

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610705725.XA Active CN106154225B (zh) 2016-08-22 2016-08-22 基于波达方向歧义消除的定位方法及装置

Country Status (1)

Country Link
CN (1) CN106154225B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107436969A (zh) * 2017-07-03 2017-12-05 四川大学 一种基于遗传算法的三维多目标定位方法
CN110210384A (zh) * 2019-05-31 2019-09-06 北京科技大学 一种道路全局信息实时提取与表示系统
CN110716177A (zh) * 2019-10-22 2020-01-21 哈尔滨工程大学 分布式声学定位网络中多目标分类方法
CN112763974A (zh) * 2020-12-17 2021-05-07 深圳先进技术研究院 定位系统中传感器的部署方法、装置及介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101442823A (zh) * 2008-12-17 2009-05-27 西安交通大学 基于波达方向估计的wsn分布式节点定位方法
CN102520389A (zh) * 2011-11-23 2012-06-27 华中科技大学 一种二维空间谱谱峰搜索方法
CN104392114A (zh) * 2014-11-11 2015-03-04 西北大学 一种基于空时数据的高分辨目标方位估计方法
CN104614653A (zh) * 2015-02-09 2015-05-13 国家电网公司 基于天线阵列的局放检测装置的多局放点定位和分辨方法
CN104993860A (zh) * 2015-05-08 2015-10-21 西北工业大学 基于阵列冲激响应的多径信号波达角估计方法
US9279880B2 (en) * 2014-07-15 2016-03-08 Applied Signals Intelligence, Inc. Electrically small, range and angle-of-arrival RF sensor and estimation system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101442823A (zh) * 2008-12-17 2009-05-27 西安交通大学 基于波达方向估计的wsn分布式节点定位方法
CN102520389A (zh) * 2011-11-23 2012-06-27 华中科技大学 一种二维空间谱谱峰搜索方法
US9279880B2 (en) * 2014-07-15 2016-03-08 Applied Signals Intelligence, Inc. Electrically small, range and angle-of-arrival RF sensor and estimation system
CN104392114A (zh) * 2014-11-11 2015-03-04 西北大学 一种基于空时数据的高分辨目标方位估计方法
CN104614653A (zh) * 2015-02-09 2015-05-13 国家电网公司 基于天线阵列的局放检测装置的多局放点定位和分辨方法
CN104993860A (zh) * 2015-05-08 2015-10-21 西北工业大学 基于阵列冲激响应的多径信号波达角估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
徐宏宇 等: "《基于扩展卡尔曼滤波的GPS定位数据处理方法研究》", 《科学技术与工程》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107436969A (zh) * 2017-07-03 2017-12-05 四川大学 一种基于遗传算法的三维多目标定位方法
CN107436969B (zh) * 2017-07-03 2020-07-24 四川大学 一种基于遗传算法的三维多目标定位方法
CN110210384A (zh) * 2019-05-31 2019-09-06 北京科技大学 一种道路全局信息实时提取与表示系统
CN110716177A (zh) * 2019-10-22 2020-01-21 哈尔滨工程大学 分布式声学定位网络中多目标分类方法
CN110716177B (zh) * 2019-10-22 2022-09-09 哈尔滨工程大学 分布式声学定位网络中多目标分类方法
CN112763974A (zh) * 2020-12-17 2021-05-07 深圳先进技术研究院 定位系统中传感器的部署方法、装置及介质
CN112763974B (zh) * 2020-12-17 2024-04-12 深圳先进技术研究院 定位系统中传感器的部署方法、装置及介质

Also Published As

Publication number Publication date
CN106154225B (zh) 2018-05-04

Similar Documents

Publication Publication Date Title
Azzouzi et al. New measurement results for the localization of uhf rfid transponders using an angle of arrival (aoa) approach
US10571544B2 (en) Direction finding using signal power
CN107015191B (zh) 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法
CN105676171B (zh) 单通道双基站超短波信号空间定位方法
CN106872935B (zh) 一种基于四元数的电磁矢量传感器阵列波达方向估计方法
CN103746757B (zh) 一种基于星载多波束天线的单星干扰源定位方法
CN106154225A (zh) 基于波达方向歧义消除的定位方法及装置
CN102540138B (zh) 一种多基线相位搜索式二维空间谱测向方法
CN109298388A (zh) 基于方位信息的超视距目标地理坐标直接估计方法
CN104330768B (zh) 一种基于声矢量传感器的机动声源方位估计方法
CN103018713A (zh) 基于导航数字多波束接收阵列天线的卫星跟踪测角方法
CN102841344A (zh) 一种少阵元近场宽带信号源参数估计方法
CN102520399B (zh) 基于电磁矢量阵列的米波雷达角度估计方法
CN106501770A (zh) 基于幅相误差阵列的远近场宽带混合源中近场源定位方法
CN102200573B (zh) 一种对近场目标信号来波方向进行测定的方法
CN110389316B (zh) 基于两步矩阵差分的近场和远场混合信源定位方法
CN104793177B (zh) 基于最小二乘法的麦克风阵列测向方法
CN107907853A (zh) 一种基于均匀圆阵差分相位的单分布源doa估计方法
CN108020812A (zh) 基于特殊三平行线阵结构的二维doa估计方法
CN104181513A (zh) 一种雷达天线阵元位置的校正方法
CN112130111A (zh) 一种大规模均匀十字阵列中单快拍二维doa估计方法
CN106874642A (zh) 基于旋转均匀圆阵的近场源参数估计解模糊方法
CN113835063A (zh) 一种无人机阵列幅相误差与信号doa联合估计方法
CN110376547B (zh) 基于二阶统计量的近场源定位方法
CN110376546B (zh) 基于协方差矩阵差分的远场和近场混合信源定位方法

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