CN107340512B - 一种基于子阵划分的近远场混合源被动定位方法 - Google Patents
一种基于子阵划分的近远场混合源被动定位方法 Download PDFInfo
- Publication number
- CN107340512B CN107340512B CN201710511143.2A CN201710511143A CN107340512B CN 107340512 B CN107340512 B CN 107340512B CN 201710511143 A CN201710511143 A CN 201710511143A CN 107340512 B CN107340512 B CN 107340512B
- Authority
- CN
- China
- Prior art keywords
- sub
- array
- source
- subarray
- matrix
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims abstract description 41
- 230000003595 spectral effect Effects 0.000 claims abstract description 17
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 28
- 238000009499 grossing Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 28
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 6
- 238000003491 array Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000004807 localization Effects 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241000408529 Libra Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
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
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/14—Systems for determining distance or velocity not using reflection or reradiation using ultrasonic, sonic, or infrasonic waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于子阵划分的近远场混合源被动定位方法,属于阵列信号处理领域。奔放吗应用对称均匀线阵传感器阵列接收目标信号;通过恰当选择传感器构造子阵,对每个子阵进行协方差处理,选取合适的传感器数据,构造只包含角度信息的矩阵,对矩阵进行特征值分解,获得相应的噪声子空间,利用一维MUSIC谱峰搜素获得真实角度对应的子角度;再根据空间几何关系求出真实角度信息,从而进一步利用MUSIC谱峰搜索求出距离信息。本发明避免了四阶累积量的使用,有效减低了算法的计算复杂度,缩短算法运行时间,避免了额外的参数匹配过程。
Description
技术领域
本发明属于阵列信号处理领域,具体涉及一种对多目标近远场混合源的被动定位方法。
背景技术
信源被动定位是阵列信号处理领域里一个重要的研究内容,根据定位目标与传感器阵列之间的距离,可以分为近场源定位和远场源定位。传统的被动定位主要为远场源定位,其信源的波前为平面波,但是当有些信源位于阵列孔径的菲涅尔(Fresnel)区(例如:近场源)时,波前的固有弯曲不能忽略,即远场平面波的假设不再成立,则需要用球面波来进行描述。在实际应用场景中,例如当使用麦克风阵列对说话人进行定位的时候,这时目标信号即可能处于阵列孔径的夫琅禾费区域,也可能处于菲涅尔区域,即观测信号由近场源和远场源共同组成。若将纯远场或者纯近场的定位方法运用到近远场混合源的情况中来,则会存在计算复杂度高,近远场难以分离及参数估计错误等问题。因此研究近远场混合源被动定位的参数估计算法,对促进麦克风,雷达,声呐等系统中的应用,具有十分重要的意义。
远场近似(Far-Field Approximation,FFA)法可认为是最早解决远近场混合源定位问题的一个途径。该算法将近场协方差矩阵作为远场协方差矩阵的有损模型,根据远场协方差矩阵的Toeplitz特性来构造FFA协方差矩阵,在此基础上利用远场MUSIC(空间频谱估计)技术进行参量估计。1995年,Lee等人探索了阵列观测数据的循环相关(二阶循环矩)特性,将该算法进一步扩展,并提出了适用于循环平稳信源的改进算法。然而,FFA算法及其改进形式均基于近场源距离远远大于阵列孔径的假设条件,这导致当近场源比较接近传感器阵列时,相应的定位性能明显下降。
2010年,梁军利等人提出了基于四阶累积量的两步MUSIC算法。该算法通过选择特定的传感器观测数据构造两个特殊的四阶累积量矩阵,使得第一个方向矩阵仅包含角度信息,而第二个方向矩阵同时包含角度和距离参量,应用一维MUSIC谱峰搜索获得远场源与近场源的方位角,并将得到的DOA信息代入二维搜索实现距离估计。分析该算法的实现过程,可知存在如下两个不足:(1)高维四阶累积量矩阵的构建导致其计算复杂度较高;(2)当远场源与近场源具有相近甚至相同的方位角时,第一个方向矩阵不再满足列满秩条件,导致信号子空间和噪声子空间难以正确区分,出现估计错误问题。
2012年,He等人提出了基于二阶统计量的斜投影算法。该算法在通过一维MUSIC谱峰搜索获得远场源方位角的基础上,将斜投影技术应用到阵列观测数据,实现了远场源和近场源的分离,避免了因角度模糊引起的估计错误问题,进一步利用均匀线阵的对称性估计出近场源方位角和距离。该算法的实施过程仅依赖于二阶统计量,具有计算复杂度较低的优势。然而,由于在估计近场源方位角时仅利用了协方差矩阵的交叉对角线信息,这导致相应的近场源定位精度较低。
2013年,王波等人探索了阵列孔径扩展技术,提出了四阶累积量与二阶统计量相结合的混合阶MUSIC算法,改进了定位参量估计的分辨率。然而与两步MUSIC算法类似,该算法依然存在计算复杂度高和估计错误两个问题。
2014年,姜佳佳等人提出了无需任何谱峰搜索的远近场混合源定位参量估计新算法,但该算法本质上是近场ESPRIT-Like算法和远场求根MUSIC算法的直接结合,且远近场混合源的分离是在定位参量估计之后实现的。
发明内容
本发明的发明目的在于:针对上述存在的问题,提供一种基于子阵划分的近远场混合源被动定位方法。
步骤1:通过对称均匀线性传感器阵列接收目标信号,得到各个阵元的接收信号;
步骤2:将对称均匀线性传感器阵列进行子阵划分,得到第一子阵、第二子阵,其中第一子阵为:去除所述对称均匀线性传感器阵列的最右端阵元的子阵;第二子阵为:去除所述对称均匀线性传感器阵列的最左端阵元的子阵;
步骤3:基于阵元的接收信号分别计算第一、二子阵的协方差矩阵;
取第一子阵的协方差矩阵的反对角线元素,按照从右上角到左下角的顺序,得到向量r1;
取第二子阵的协方差矩阵的反对角线元素,按照从右上角到左下角的顺序,得到向量r2;
步骤4:分别对向量r1、r2进行空间平滑处理:
将向量ri划分为L个重叠的子向量,每个子向量包括(2M-L+1)个元素,其中L≥2,2M+1对应对称均匀线性传感器阵列的阵元数,下标i∈{1,2};
将每个子向量与自己的共轭转置矩阵相乘,得到每个子向量的协方差矩阵,由L个子向量的协方差矩阵的均值得到协方差矩阵Ri;
步骤6:确定K个目标信源的真实方位角:
步骤7:由步骤6得到的近场源的真实方位角得到对应的导向矢量,并通过一维MUSIC谱峰搜索法得到近场源的距离估计,完成对近远场混合源中的近场源的定位。
综上所述,由于采用了上述技术方案,本发明的有益效果是:运用二阶统计量实现远近场混合源定位,避免了四届累积量的使用,有效降低了算法的计算复杂度,缩短算法运行时间;同时,在定位近场源时,将获得的近场源方位角带入到导向矢量中去,进行一维music谱峰搜索,避免了二维谱峰搜索。
附图说明
图1是本发明使用的对称均匀线阵结构图示意图;
图2估计角度的均方根误差随信噪比变化的曲线;
图3估计距离的均方根误差随信噪比变化的曲线;
图4估计角度的均方根误差随快拍数变化的曲线;
图5估计距离的均方根误差随快拍数变化的曲线。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式和附图,对本发明作进一步地详细描述。
通过本发明的方法实现K个目标信源(简称信源)中的近场源的定位处理如下:
S1:通过对称均匀线性传感器阵列接收目标信号,得到各个阵元的接收信号xm(t),其中m为阵元标识符。
参见图1,对称均匀线性传感器阵列共包括2M+1个阵元,各阵元的间距为d。将中心阵元作为参考阵元,并将其阵元标识符设置为0,其左边的阵元标识符为负,右边的阵元标识符为正,则2M+1个阵元的标识符依次为:-M,-M+1,…,-1,0,1,…,M-1,M。
用sk(t)表示远场源或近场源包络,nm(t)表示阵元的噪声,τmk表示信源k(k=1,2,……,K)从参考阵元到第m个阵元的时延差,则接收信号xm(t)可表示为:其中,e表示自然底数,j为虚数单位。
S2:将对称均匀线性传感器阵列划分为两个子阵,每个子阵包括2M个阵元,其中子阵一为阵列的前2M个阵元构成,子阵二为阵列的2M个阵元构成。子阵一、子阵二的接收数据可以表示为:
X1(t)=[x-M(t),...,x0(t),...,xM-1(t)]T
X2(t)=[x-M+1(t),...,x0(t),...,xM(t)]T
符号(·)T表示矩阵转置。
S3:对每个子阵的接收数据进行协方差运算,其计算表达式为:其中i∈{1,2},m、n表示阵元标识符,符号(·)H表示矩阵的共轭转置。
取各矩阵的协方差矩阵的反对角线元素按照从右上角到左下角的顺序,得到向量ri,即:
r1=[r1(-M,M-1),r1(-M+1,M-2),…,r1(M-1,-M)];
r2=[r2(-M+1,M),r1(-M+2,M-1),…,r1(M,-M+1)];
从而将传统的二维估计转换为一维估计。
S4:由于向量ri存在空间系数的相关性,不能直接利用MUSIC谱峰搜索法,因此对其进行空间平滑处理,从而获得相应的噪声子空间,以便于利用MUSIC谱峰搜索法。
首先将向量ri分成L(L≥2)个重叠子空间ri,l,每个子向量ri,l包括(2M-L+1)个元素;
最后,对各子阵对应的协方差矩阵Ri进行特征值分解,由于零特征值所对应的特征向量仅包含噪声成分,所有的仅包含噪声成分的特征向量可组成噪声子空间。因此,选择所有的零特征值,将与之相对应的特征向量组成噪声子空间
S5:通过一维MUSIC谱峰搜索实现远近场混合源方位角在各个子阵的子角度(即各子阵的各信源的方位角θi,k)的同时估计:
依据MUSIC方法的基本原理,当将真实方位角的子角度代入下式时,函数P将出现最大值。因此远场源和近场方位角的子角度可通过寻找P的谱峰得到;
S6:获取信源的真实方位角。
因子阵一、子阵二的排列结构,可将对应的视为第k个信源相对于阵元-d/2和d/2位置的方位角(即),用rk,-d/2和rk,d/2表示第k个信源相对于阵元-d/2和d/2位置的距离,则基于几何关系可建立如下的等式:
rk,-d/2sinθk,-d/2-rk,d/2sinθk,d/2=d
rk,-d/2cosθk,-d/2-rk,d/2cosθk,d/2=0
S7:由步骤S6得到的远场源和近场源的真实方位角,将其值代入到导向矢量中利用一维MUSIC谱峰搜索,实现近场源距离估计;
依据MUSIC方法的基本原理,当将真实的近场源距离代入下式时,函数P将出现最大值;因此近场源的距离估计值可通过寻找P的谱峰得到。
下面通过仿真实验数据分析本发明提出的定位方法的定位性能及计算有效性,仿真实验1和仿真实验2所采用的仿真软件为MATLAB软件。
仿真实验1:该实验用以分析本发明提出的定位方法估计远近场混合源定位参量的性能。对称均匀线阵的传感器数目为9个(即阵元数为9),定位参量分别为(θ1,r1)=(-20°,2.7λ)和θ2=10°,当信噪比从0分贝递增到25分贝时,500次蒙特卡罗实验的仿真结果如图2、3所示。分析该仿真结果可知,本发明所提出的定位方法在估计近场源方位角时,相应的均方根误差随信噪比增加而平稳地变小,其方位角和距离的均方根误差小于现有的斜投影算法(BSOS算法)和两步MUSIC算法,即估计性能略高于这两种算法。
仿真实验2:该实验用以对比评价本发明所提出的定位方法的均方根误差。对称均匀线阵的传感器数目为9个,信噪比为10分贝,角度搜索步长和距离搜索步长分别为0.1和0.01,其他仿真条件与实验1相同。当采样点数从Ts=50变化到Ts=1400时,本发明所提出的定位方法和BSOS算法,两步MUSIC算法进行比较如图4、5所示。分析该实验结果可知,本发明的均方根误差随采样点数的递增,其结果都略小于以上两种方法,因此可能认为更具有实用性。
Claims (1)
1.一种基于子阵划分的近远场混合源被动定位方法,其特征在于,包括下列步骤:
步骤1:通过对称均匀线性传感器阵列接收目标信号,得到各个阵元的接收信号;
步骤2:将对称均匀线性传感器阵列进行子阵划分,得到第一子阵、第二子阵,其中第一子阵为:去除所述对称均匀线性传感器阵列的最右端阵元的子阵;第二子阵为:去除所述对称均匀线性传感器阵列的最左端阵元的子阵;
步骤3:基于阵元的接收信号分别计算第一、二子阵的协方差矩阵;
取第一子阵的协方差矩阵的反对角线元素,按照从右上角到左下角的顺序,得到向量r1;
取第二子阵的协方差矩阵的反对角线元素,按照从右上角到左下角的顺序,得到向量r2;
步骤4:分别对向量r1、r2进行空间平滑处理:
将向量ri划分为L个重叠的子向量,每个子向量包括(2M-L+1)个元素,其中L≥2,2M+1对应对称均匀线性传感器阵列的阵元数,下标i∈{1,2};
将每个子向量与自己的共轭转置矩阵相乘,得到每个子向量的协方差矩阵,由L个子向量的协方差矩阵的均值得到协方差矩阵Ri;
步骤6:确定K个目标信源的真实方位角:
步骤7:由步骤6得到的近场源的真实方位角得到对应的导向矢量,并通过一维MUSIC谱峰搜索法得到近场源的距离估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710511143.2A CN107340512B (zh) | 2017-06-29 | 2017-06-29 | 一种基于子阵划分的近远场混合源被动定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710511143.2A CN107340512B (zh) | 2017-06-29 | 2017-06-29 | 一种基于子阵划分的近远场混合源被动定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107340512A CN107340512A (zh) | 2017-11-10 |
CN107340512B true CN107340512B (zh) | 2020-02-18 |
Family
ID=60217954
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710511143.2A Expired - Fee Related CN107340512B (zh) | 2017-06-29 | 2017-06-29 | 一种基于子阵划分的近远场混合源被动定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107340512B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414985A (zh) * | 2018-02-13 | 2018-08-17 | 西北工业大学 | 基于三阶循环累积量的近场定位方法 |
CN108919178A (zh) * | 2018-08-06 | 2018-11-30 | 电子科技大学 | 一种基于对称嵌套阵列的混合场信源定位方法 |
CN109085531B (zh) * | 2018-08-27 | 2023-04-07 | 西安电子科技大学 | 基于神经网络的近场源到达角估计方法 |
CN109406627B (zh) * | 2018-09-11 | 2021-07-20 | 温州大学苍南研究院 | 一种基于旋转不变子空间的复合材料板冲击定位方法 |
CN110376547B (zh) * | 2019-07-20 | 2022-11-22 | 中国人民解放军国防科技大学 | 基于二阶统计量的近场源定位方法 |
CN111257822B (zh) * | 2020-03-05 | 2022-12-30 | 西北工业大学 | 基于近场稀疏阵列的准平稳信号参数估计方法 |
CN111830460B (zh) * | 2020-07-23 | 2022-11-01 | 哈尔滨工业大学(威海) | 基于序贯music的doa估计方法 |
CN113030850B (zh) * | 2021-03-04 | 2022-11-11 | 西北工业大学 | 基于复数域表征与学习的近场源智能定位方法 |
CN113702899B (zh) * | 2021-08-03 | 2023-09-29 | 哈尔滨工程大学 | 一种基于相位分数低阶矩的协方差差分的传播算法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175989A (zh) * | 2010-12-30 | 2011-09-07 | 电子科技大学 | 一种非相干分布式信号二维波达角的测定方法 |
CN105589056A (zh) * | 2015-12-15 | 2016-05-18 | 吉林大学 | 一种多目标远近场混合源定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150002446A1 (en) * | 2013-06-28 | 2015-01-01 | Synaptics Incorporated | Wireless communication enabling capacitive imaging sensor assembly |
-
2017
- 2017-06-29 CN CN201710511143.2A patent/CN107340512B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175989A (zh) * | 2010-12-30 | 2011-09-07 | 电子科技大学 | 一种非相干分布式信号二维波达角的测定方法 |
CN105589056A (zh) * | 2015-12-15 | 2016-05-18 | 吉林大学 | 一种多目标远近场混合源定位方法 |
Non-Patent Citations (2)
Title |
---|
Localization of mixed far-field and near-field sources under unknown;Jian Xie等;《Digital Signal Processing》;20161231;全文 * |
分布式信源低复杂度参数估计算法研究;郑植;《中国博士学位论文全文数据库信息科技辑》;20111215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107340512A (zh) | 2017-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107340512B (zh) | 一种基于子阵划分的近远场混合源被动定位方法 | |
CN106054123B (zh) | 一种稀疏l阵及其二维doa估计方法 | |
CN110927659B (zh) | 互耦条件下任意阵列流形doa估计与互耦校准方法及系统 | |
CN109444810B (zh) | 一种非负稀疏贝叶斯学习框架下的互质阵列非网格doa估计方法 | |
CN105589056B (zh) | 一种多目标远近场混合源定位方法 | |
CN106093921B (zh) | 基于稀疏分解理论的声矢量阵宽带测向方法 | |
CN107315162B (zh) | 基于内插变换和波束形成的远场相干信号doa估计方法 | |
CN110045323B (zh) | 一种基于矩阵填充的互质阵稳健自适应波束形成算法 | |
CN112379327A (zh) | 一种基于秩损估计的二维doa估计与互耦校正方法 | |
CN113673317B (zh) | 基于原子范数最小化可降维的二维离格doa估计方法 | |
CN109696657B (zh) | 一种基于矢量水听器的相干声源定位方法 | |
CN109557504B (zh) | 一种近场窄带信号源的定位方法 | |
CN113567913A (zh) | 基于迭代重加权可降维的二维平面doa估计方法 | |
CN111693947A (zh) | 基于互质阵列doa估计的改进music方法 | |
CN111352063A (zh) | 一种均匀面阵中基于多项式求根的二维测向估计方法 | |
CN110244258A (zh) | 用于双平行阵二维测向中扩展doa矩阵方法 | |
CN107450046B (zh) | 低仰角多径环境下的波达角估计方法 | |
CN110895325A (zh) | 基于增强四元数多重信号分类的到达角估计方法 | |
CN109507634B (zh) | 一种任意传感器阵列下的基于传播算子的盲远场信号波达方向估计方法 | |
CN116699511A (zh) | 一种多频点信号波达方向估计方法、系统、设备及介质 | |
CN115980721A (zh) | 一种无误差协方差矩阵分离的阵列自校正方法 | |
CN106844886B (zh) | 基于主分量分析的目标波达方向获取方法 | |
CN114184999B (zh) | 一种互耦小孔径阵列的生成式模型处理方法 | |
CN114966523A (zh) | 大规模mimo测向中基于幂迭代算法的doa估计方法 | |
Eskandari et al. | A novel solution for root-MUSIC with reduced complexity |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Zheng Zhi Inventor after: Sun Jie Inventor after: Zheng Yan Inventor after: Kong Yangyang Inventor before: Sun Jie Inventor before: Zheng Zhi Inventor before: Zheng Yan Inventor before: Kong Yangyang |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200218 |
|
CF01 | Termination of patent right due to non-payment of annual fee |