CN103558581A - 一种多径环境下的到来角估计算法 - Google Patents
一种多径环境下的到来角估计算法 Download PDFInfo
- Publication number
- CN103558581A CN103558581A CN201310340255.8A CN201310340255A CN103558581A CN 103558581 A CN103558581 A CN 103558581A CN 201310340255 A CN201310340255 A CN 201310340255A CN 103558581 A CN103558581 A CN 103558581A
- Authority
- CN
- China
- Prior art keywords
- signal
- arrival angle
- algorithm
- sigma
- rank
- 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
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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/46—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mobile Radio Communication Systems (AREA)
- Noise Elimination (AREA)
Abstract
本发明涉及一种多径环境下的到来角估计算法,属于无线通信技术领域。本发明针对移动通信中存在的多径问题,基于多径信号的衰落特性,提出了将时间平均与合成空间谱方法结合的算法,该方法解决了多径环境下信号源相干、天线数相对少时的到来角估计问题,不存在阵列孔径的损失;同时也有效的解决了多径环境下的低信噪比、来波方位间隔较近时的到来角估计问题。
Description
技术领域
本发明涉及一种多径环境下的到来角估计算法,属于无线通信技术领域。
背景技术
为了降低移动通信***中的多址干扰、降低发射功率和提高***容量,智能天线技术已成为移动通信领域的研究热点,天线到来角(DOA,Direction ofArrival)的研究一直是天线技术研究的热点之一。到来角是指同一天线接收到的各路信号到达接收端时的方向不同而形成的角。作为反映用户空间位置的重要参量在智能天线中扮演着非常重要的角色。到来角估计,是信号处理的一个重要研究方向,是将多个天线设置在空间的不同位置组成天线阵列,并利用这一阵列对空间信号场进行接收和处理,提取阵列所接收的信号及其特征信息,其应用涉及雷达、通信、地震、声呐以及生物医学工程等众多领域。
然而,在移动通信环境中,同一移动用户信号经过各种反射体、折射体形成多径信号,这种经反射、折射的多径信号和直达波信号矢量相加后被天线接收,从而引起接收端合成接收波振幅和相位的变化。由于这两种信号是相干的,将直达波从多径信号中分离出来比较困难,对到来角估计引起较大误差。
由于相干多径的存在,造成普通的阵列超分辨到来角估计算法的失效,如经典的多重信号分类(MUSIC,Multiple Signal Classification)算法、旋转不变子空间法(ESPRIT,Estimation of Signal Parameters via Rotational InvarianceTechniques)等。对此,很多学者提出了空间平滑算法及其改进算法,此类算法能有效的对相干信号进行解相干,但信源的去相干往往是以牺牲阵列孔径为代价的,可估计的最大到来角数为2M/3(M为阵元数),而移动通信中的多径数往往远远大于此。
发明内容
本发明针对以上问题的提出,而研制一种多径环境下的到来角估计算法。
本发明的技术方案如下:
一种多径环境下的到来角估计算法如下:
在多径环境下,假设M元天线阵列,障碍物相对信号源较近,Q个独立的窄带信号源S(k)(Q个用户),经过P条多径(共P个波前)形成P个远场信号辐射到阵面,其中Q<M。
阵列接收信号向量可表述为:
X(k)=Au(k)+n(k),k=1,2,L,K, (1)
其中,阵列导向矢量矩阵A=[A1 A2…AQ],Aq=[a(θq1) a(θq2)…a(θqpq)],q=1,2,…,Q。a(θqp)是第q个信源的第P条多径的方向向量。u(k)是P×1的由S(k)经P条多径产生的波前向量;n(k)是噪声向量;K是快拍次数;θ为信号的到来角。
对于均匀线阵,ωqp=(2πd/λ)sin(θqp),d是阵元间距,λ是载波波长。
因此,X(k)也可以表示为:
X(k)=ACs(k)+n(k), (2)
矩阵C的形式为
式中0是由若干个0组成的零向量,cq是第q个独立信号源的pq条多径的衰落因子(复幅度)αqp构成的列向量。显然有,p1+p2+…+pQ=P。
在典型的移动通信环境中,相对于移动用户与基站之间的大距离,移动用户的移动是缓慢的。因此,通常认为多径的到来角是慢变的,但是αqp是快变的,特别是衰落因子的相位可以认为是在[0,2π]内均匀分布的随机数。考虑接收一个移动的移动用户的m个时间间隔的信号,在这些时间间隔,到来角的所有波前认为是不变的。在每个时间间隔,接收信号的协方差矩阵为:
由于多径的影响,取m个时间间隔的协方差矩阵,做平均,有:
其中,
由于各个时间间隔之间是快变的,所以又因为rank(A)=P,所以,平均后的协方差矩阵RT-X也满秩,即rank(RT-X)=P。然后再结合合成空间谱算法(SSM,Synthetic SpatialSpectrum Method)提出本发明的基于时间平均的合成空间谱算法,即经过时间平均后的合成空间谱算法(T-SSM,Time averaging-Synthetic Spatial SpectrumMethod)。该算法提高了对信号子空间、噪声子空间和主特征值中的信息的利用,结合了信号子空间处理稳健性高和噪声子空间处理估计精度高的优点,实现了高分辨率和很强的稳健性,低信噪比时,算法依然有效。阵列输出功率与到来角的关系为:
其中,λj,qj分别为平均后的RT-X的第j个特征值和对应的特征向量。a(θ)是信号方向向量,H是转置的数学符号。通过找PT-SSM(θ)的最大值,得到到来角的估计值。
本发明的有益效果:
1、本发明解决了多径环境下信号源相干、天线数相对少时的到来角估计问题,不存在阵列孔径的损失,可估计的最大到来角数为M-1。
2、本发明与传统的方法相比,有效的解决了多径环境下的低信噪比、来波方位间隔较近时的到来角估计问题。
3、仿真建模分析表明,本发明具有很大的工程应用价值,对于国家军事和民用等均有积极的意义。
附图说明
图1:到来角估计的***框图。
图2:到来角估计的***示意图。
图3(a):传统到来角估计算法的简化框图。
图3(b):基于时间平均的到来角估计算法的简化框图。
图4:信噪比为15dB时,MUSIC算法、SSM算法、T-MUSIC算法(时间平均后的MUSIC算法),T-SSM算法的性能比较。
图5:信噪比为-5dB时,T-MUSIC算法与T-SSM算法的性能比较。
具体实施方式
下面结合附图对本发明进一步描述。
本发明的一个实施例的实施方式如下:
如图1所示:本发明涉及的整个***由三部分组成:空间信号入射、空间阵列接收和参数估计。我们的目标是利用一些方法从复杂的环境中估计出信号的到来角度。图2是本发明***的整体示意图:同一移动用户信号s(k),经过各种反射体、折射体后形成含有多径信号的入射信号u(k),入射到接收阵列上。这种经反射、折射的多径信号和直达波信号矢量相加后被天线接收,从而引起接收端合成接收波振幅和相位的变化。由于这两种信号是相干的,将直达波从多径信号中分离出来比较困难,对到来角估计引起较大误差。本发明提出的方法有效地解决了此问题。
为了深入说明与理解本发明的实施方式,我们首先比较了传统的到来角估计算法,其原理框图如图3(a)所示:将阵列接收信号X求取协方差,得到协方差矩阵RX,然后对协方差矩阵进行矩阵分解,得到特征值和对应的特征向量。按特征值从大到小的顺序,把与信号个数相等的大特征值和对应的特征向量看做信号子空间,把剩下的小特征值和对应的特征向量看做噪声子空间。
传统的到来角估计算法中,接收到的数据矢量的协方差矩阵可以表示为
RX=E[XXH]=AE{ssH}AH+E{nnH}, (7)
或 RX=ARSAH+RN, (8)
式中,RS为信号的协方差矩阵,RN为噪声的协方差矩阵,对RX进行特征分解,得到M个特征值,将这些特征值按从大到小的顺序排列有特征值λi的特征向量为qi,其中大特征值对应的特征向量构成信号子空间(主特征子空间)US,小特征值对应的特征向量构成噪声子空间(非主特征子空间)UN。
MUSIC类算法是利用信号子空间(主特征子空间)US和噪声子空间(非主特征子空间)UN的正交性,通过在所有可能的阵列方向向量中搜寻那些与非主特征向量张成的空间相垂直的向量,从而通过MUSIC空间谱的峰值就可以估计出多个入射信号的到来角。这些峰值可以由下式表示:
其中,a(θ)是方向向量,qi为第i个特征值对应的特征向量,P+1≤i≤M,H是共轭转置的数学符号。通过找PMUSIC(θ)的最大值,得到到来角的估计值。
但是,在多径环境中,信号源由于多径效应,使得入射波前为强相关信号,破坏了信号子空间与噪声子空间的正交性,导致接收阵列的协方差矩阵奇异,以致MUSIC算法、SSM算法等子空间类算法均失效。
本发明提出的多径环境下的到来角估计算法的实现如下,如图3(b)所示:在典型的移动通信环境中,相对于移动用户与基站之间的大距离,移动用户的移动是缓慢的。因此,通常认为多径的到来角是慢变的,但是αqp是快变的,特别是衰落因子的相位可以认为是在[0,2π]内均匀分布的随机数。考虑接收一个移动的移动用户的m个时间间隔的信号,在这些时间间隔,到来角的所有波前认为是不变的。对m个时间间隔的信号的协方差矩阵求取平均得到RT-X,然后再对RT-X进行矩阵分解,得到特征值和对应的特征向量。
具体的算法:在T-SSM算法中,由于多径的影响,所以取m个时间间隔的协方差矩阵,做平均,有:
其中,
由于各个时间间隔之间是快变的,所以又因为 rank(A)=P,所以,平均后的协方差矩阵RT-X也满秩,即rank(RT-X)=P。接下来对RT-X进行矩阵分解,得到特征值和对应的特征向量。阵列输出功率与到来角的关系为:
通过找PT-SSM(θ)的最大值,得到到来角的估计值。其中,λj,qj分别为平均后的RT-X的第j个特征值和对应的特征向量。a(θ)是方向向量,H是转置的数学符号。该方法提高了对信号子空间、噪声子空间和主特征值中的信息的利用,结合了信号子空间处理稳健性高和噪声子空间处理估计精度高的优点,实现了高分辨率和很强的稳健性。低信噪比时,算法依然有效。
我们对传统方法与本发明方法的仿真进行了比较。如图4所示:信噪比为15dB,入射波前为-30°,-40°,0°,45°时,MUSIC算法、SSM算法、T-MUSIC算法(时间平均后的MUSIC算法)与T-SSM算法的性能比较,从图4中可以看出,MUSIC算法、SSM算法均失效,T-MUSIC算法与T-SSM能够有效的估计出到来角,谱峰“尖锐”且平坦性好。为了验证低信噪比、来波方位间隔较近时的算法性能,我们在入射波前角度分别为-30°,-33.5°,0°,45°,信噪比为-5dB时,在高斯白噪声的条件下对T-MUSIC算法与T-SSM算法的性能作了比较,如图5所示:从图5的结果可知,多径环境下,入射波前为强相关信号且小信噪比时,T-MUSIC算法失效,然而T-SSM算法却可以有效的估计出到来角,且谱峰“尖锐”。
针对移动通信中存在的多径问题,我们提出的方法能有效的解决多径引起的相干问题,且不存在阵列孔径的损失,低信噪比、来波方位间隔较近时算法依然有效。仿真建模分析表明,该方法具有很大的工程应用价值。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
本发明涉及到的缩略语和关键术语:
MUSIC(Multiple Signal Classification)算法:多重信号分类算法。
ESPRIT(Estimation of Signal Parameters via Rotational Invariance Techniques)算法:旋转不变子空间法。
SSM(Synthetic Spatial Spectrum Method)算法:合成空间谱算法。
T-SSM(Time averaging-Synthetic Spatial Spectrum Method):时间平均后的SSM算法。
T-MUSIC算法(Time averaging-Multiple Signal Classification算法):时间平均后的MUSIC算法。
Claims (2)
2.如权利要求1所述的一种多径环境下的到来角估计算法,其特征在于:阵列输出功率与到来角的关系为:
其中,λj,qj分别为平均后的RT-X的第j个特征值和对应的特征向量,a(θ)是信号的方向向量,H是共轭转置符号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310340255.8A CN103558581B (zh) | 2013-08-06 | 2013-08-06 | 一种多径环境下的到来角估计算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310340255.8A CN103558581B (zh) | 2013-08-06 | 2013-08-06 | 一种多径环境下的到来角估计算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103558581A true CN103558581A (zh) | 2014-02-05 |
CN103558581B CN103558581B (zh) | 2016-05-04 |
Family
ID=50012883
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310340255.8A Expired - Fee Related CN103558581B (zh) | 2013-08-06 | 2013-08-06 | 一种多径环境下的到来角估计算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103558581B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104598732A (zh) * | 2015-01-19 | 2015-05-06 | 西安交通大学 | 一种相干与非相干混合信号的半参数个数估计方法 |
CN109471065A (zh) * | 2018-09-28 | 2019-03-15 | 中国电子科技集团公司第三十六研究所 | 一种相干信号的测向方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1739221A (zh) * | 2002-11-22 | 2006-02-22 | 本·古里安大学 | 改进偏振信号源定位的智能天线*** |
-
2013
- 2013-08-06 CN CN201310340255.8A patent/CN103558581B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1739221A (zh) * | 2002-11-22 | 2006-02-22 | 本·古里安大学 | 改进偏振信号源定位的智能天线*** |
Non-Patent Citations (1)
Title |
---|
岳晋: "移动通信中自适应阵列信号处理关键技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104598732A (zh) * | 2015-01-19 | 2015-05-06 | 西安交通大学 | 一种相干与非相干混合信号的半参数个数估计方法 |
CN104598732B (zh) * | 2015-01-19 | 2018-03-02 | 西安交通大学 | 一种相干与非相干混合信号的半参数个数估计方法 |
CN109471065A (zh) * | 2018-09-28 | 2019-03-15 | 中国电子科技集团公司第三十六研究所 | 一种相干信号的测向方法 |
CN109471065B (zh) * | 2018-09-28 | 2021-05-25 | 中国电子科技集团公司第三十六研究所 | 一种相干信号的测向方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103558581B (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101325807B (zh) | 信号波达方向估计方法 | |
CN107167763B (zh) | 基于非圆特性的远近场混合信号波达方向估计方法 | |
CN102175989B (zh) | 一种非相干分布式信号二维波达角的测定方法 | |
CN106054123A (zh) | 一种稀疏l阵及其二维doa估计方法 | |
CN103353596A (zh) | 基于压缩感知的波束空间域米波雷达测高方法 | |
CN104698433A (zh) | 基于单快拍数据的相干信号doa估计方法 | |
CN104515969B (zh) | 一种基于六角形阵列的相干信号二维doa估计方法 | |
CN103353595A (zh) | 基于阵列内插压缩感知的米波雷达测高方法 | |
CN101431354B (zh) | 一种波达角估计方法 | |
CN102694588A (zh) | 一种基于共轭扩展的波达方向估计方法 | |
CN103323827B (zh) | 基于快速傅里叶变换的mimo雷达***角度估计方法 | |
Hu et al. | Hybrid-field channel estimation for extremely large-scale massive MIMO system | |
CN114095318B (zh) | 智能超表面辅助的混合构型毫米波通信***信道估计方法 | |
CN107332594A (zh) | 一种mimo车载通信***及其性能评估方法 | |
CN101908918B (zh) | 无线通信接收机中的一种波束合成方法 | |
CN101483280B (zh) | 一种稳健波束合成器权值求解方法 | |
CN109521393A (zh) | 一种基于信号子空间旋转特性的波达方向估计算法 | |
CN105306117A (zh) | 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 | |
CN104506256A (zh) | 一种mimo多天线***性能评估方法及多天线*** | |
CN103558581A (zh) | 一种多径环境下的到来角估计算法 | |
CN105068090A (zh) | 一种利用单快拍数据进行gnss天线阵干扰抑制的方法 | |
CN117075033A (zh) | 基于修正降维music的智能超表面辅助的用户定位方法 | |
CN101499816B (zh) | 一种来波到达角的估计方法及装置 | |
CN105510872A (zh) | 一种适用于大规模mimo***的二维波达角的测定方法 | |
Cui et al. | Direction-of-arrival estimation for coherent GPS signals based on oblique projection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: Wang Zhisen Inventor after: Shen Ying Inventor after: Tang Yiming Inventor after: Yu Zhiyuan Inventor before: Wang Zhisen Inventor before: Shen Ying |
|
COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: WANG ZHISEN SHEN YING TO: WANG ZHISEN SHEN YING TANG YIMING YU ZHIYUAN |
|
C14 | Grant of patent or utility model | ||
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: 20160504 Termination date: 20190806 |