CN105262550A - 一种高阶累积量来波方向快速估计的方法 - Google Patents
一种高阶累积量来波方向快速估计的方法 Download PDFInfo
- Publication number
- CN105262550A CN105262550A CN201510603287.1A CN201510603287A CN105262550A CN 105262550 A CN105262550 A CN 105262550A CN 201510603287 A CN201510603287 A CN 201510603287A CN 105262550 A CN105262550 A CN 105262550A
- Authority
- CN
- China
- Prior art keywords
- matrix
- theta
- signal
- centerdot
- doa
- 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
Landscapes
- Radio Transmission System (AREA)
Abstract
本发明公开了一种高阶累积量来波方向快速估计的方法,旨在公开一种解决系统估计容量和阵列估计精度,而且在高斯色噪声环境下,可以实现相关信号的DOA估计,且阶数越高,估计性能越好的来波方向快速估计方法。本发明首先利用首先构造2q阶累积量矩阵,并对构造的2q阶累积量矩阵进行旋转变换,然后对累积量矩阵取均值,得到空间平滑处理后的2q阶累积量矩阵,再进行特征值分解,得到信号子空间和噪声子空间,利用谱函数求出最优解。上述改进优化的DOA算法不依赖于信号的具体特征,适合于高斯白噪声和高斯色噪声环境。该方法能降低计算复杂度,大幅度缩短了60GHz无线通信系统中DOA计算耗时,适用于60GHz无线通信系统。
Description
技术领域
本发明涉及通信技术领域,尤其是涉及一种应用于60GHz无线通信系统的来波方向快速估计的方法。
背景技术
近年来移动通信事业在全球范围内得到了迅猛的发展,对带宽的要求越来越高。60GHz无线通信技术由于其超高速的数据传输能力、成本低、抗干扰性和安全性高、具有国际通用性和免许可特性的特点,将成为第4代无线通信技术的代表。然而它面临的最大挑战是过高的路径损失。6OGHz频段的路径损耗决定了其必须采用波束形成技术和来波方向快速估计来改善链路质量。来波方向(DOA)估计又可称为来波方向估计或来波角度估计,是阵列信号处理的一个重要分支。DOA估计的基本问题就是如何利用天线阵列的输出,从背景噪声中检测出同时处于空间某一区域内的多个期望用户的位置信息,即各个信号到达阵列天线的入射角。DOA估计算法性能评价标准主要包括:估计偏差、估计误差、分辨力、成功概率、计算量及给定阵元数目情况下算法可估计的信号数目。
DOA是阵列信号处理的一项重要内容,其应用涉及雷达、通信、声纳、地震、勘探、射电天文以及生物医学工程等众多军事及民用经济领域。DOA估计是利用一组按一定方式布置在空间不同位置的传感器对空间信源在时域和空域同时进行采样,再由传感器阵列采样数据的分析处理来实现对空间信源的方位向估计。在DOA估计中,分辨力,精度,稳健性,速度,可分辨目标数的多少等,决定了方位向估计性能的好坏。近几十年,国内外均开展了有关DOA估计技术的理论和应用研究,使得DOA估计理论得到迅速发展。但是,在实际应用中,由于受到物理器件水平和成本的影响,阵列孔径不可能无限增大,DOA估计的分辨力始终受到阵列孔径的限制,最大可分辨目标数受到阵元数的限制。
中华人民共和国国家知识产权局于2007年04月25日公开了公告号为CN1953348A的专利文献,名称是多扇区智能天线系统的来波方向估计方法,其包括以下步骤:(1)利用谱估计算法对需要进行来波方向估计的用户进行初始来波方向估计,得到估计的来波功率谱最大值对应的来波角度;(2)根据所述来波角度判断用户是否位于目标扇区的边缘,若是,则分别计算用户在和目标扇区相邻的两个扇区内的功率;(3)若所述来波角度所在区域靠近功率值大的相邻扇区,则所述来波角度即为最终的估计结果,否则,搜索来波角度所在区域的相反区域,得到所述相反区域内来波功率谱最大值对应的来波角度即为最终的估计结果。此方案仍然收到阵元数的限制,在分辨率和分辨度上都不够理想。
发明内容
本发明主要是解决现有技术所存在的对信号入射角相差较小的信号分辨率低、对入射信号低信噪比下分辨率不高、对入射信号的分辨度不高等的技术问题,提供一种应用于60GHz无线通信系统的DOA的快速计算方法,能在改善无线接收器、发送器性能的基础上降低波束形成算法的计算复杂度,提高分辨率,从而在保证高速通讯信号品质的同时,大幅度缩短了DOA计算耗时。
该DOA改进算法,不仅具有空间平滑算法去强相干信号的效果,同时弥补其在低信噪比下分辨率不高和对信号的分辨度不高的不足之处,另外,该算法也比空间平滑算法的算法复杂度更小,大大减少了整个算法的运算量,而且不牺牲天线的孔径,比较实用,实验证明了该算法的有效性。
本发明针对上述技术问题主要是通过下述技术方案得以解决的:一种改进的DOA估计的方法,包括以下步骤:
S01、天线阵接收原始信号;
S02、构造2q阶累积量矩阵,并对构造的2q阶累积量矩阵进行旋转变换;
S03、对累积量矩阵取均值,得到空间平滑处理后的2q阶累积量矩阵;
S04、进行特征值分解,得到信号子空间和噪声子空间
S05、构造空间谱函数;
S06、判断空间谱函数是否收敛,如果不收敛,则对权向量进行迭代;如果收敛,则进入步骤S07;
S07、求出空间谱函数谱峰θ域的最优解。
q为初始值为1的正整数,每迭代一次增大1。
步骤S02及后续步骤的具体过程如下:
1、根据下式(1)构造2q阶累积量矩阵C2q,x;
其中:
B(θ)是广义方向矢量矩阵,B(θ)=[b(θ1),b(θ2),…,b(θk)];
C2q,s(l)是信号s(t)的2q阶累积量矩阵,
是阵列的输入信号的协处理信号,用于信
号预测和部分信号的补充;
H为矩阵共轭转置符号;
2、根据下式(2)得到旋转后的2q阶累积量矩阵C2q,x';
其中:
J是一个Mq×Mq的置换阵;
矩阵B(θ)经过旋转后变成JB*(θ),JB*(θ)=B(θ)Φ-(M-1),
矩阵经过旋转后变成
Φ为旋转矩阵,表示为
H为矩阵共轭转置符号;
3、对式(1)与式(2)的累积量矩阵取均值,得到空间平滑处理后的2q阶累积量矩阵为
其中
4、对2q阶累积量矩阵进行平滑处理后,再对C(l)进行特征值分解,得到信号子空间US和噪声子空间Un。构造空间谱函数,即
其中:
ξ为权向量系数;α(θ)为阵列方向矩阵;
Γ是校正函数,Γ=ρejsinθ+Ωejtgθ,ρ为信号矢量系数,Ω为噪声矢量系数;
m=0,1,2,…,M-1;k=0,1,2,…K-1
j为虚数单位,m为矩阵的行编号,也可以看作是对天线元件的编号;
k为矩阵的列编号,也可以看作是对波束的编号,K为波束数量;floor函数表示取小于或者等于括号内指定表达式的最大整数;mod表示取余运算函数,mod(X,Y)为X除以Y所得的余数;
判断空间谱函数是否收敛,如果不收敛,则对权向量进行迭代,q增大1以后重复步骤S02到S05;如果收敛,则从式(5)求出谱函数的θ域的最优解,所得到θ即为入射信号DOA的估计。
本发明提供的60GHz无线通信系统改进的DOA快速计算方法,所述DOA快速计算方法应用在超高带宽60GHz无线接收器的控制程序中。
在60GHz无线通信系统的实际应用中,由于外部干扰阻挡的存在,如发送和接收之间瞬间被障碍物阻挡等等,天线接收到的信号中可能存在某种比较少见、一般情况不考虑的干扰信号,以及其他无法预知的有效信号迅速衰减,通过本方法的处协理信号权矢量中的校正函数的综合的补偿和调整作用,可以有效消除这些干扰。
本发明提供的60GHz无线通信系统改进的DOA快速计算方法,其基本思想就是首先构造2q阶累积量矩阵,并对构造的2q阶累积量矩阵进行旋转变换,然后对累积量矩阵取均值,得到空间平滑处理后的2q阶累积量矩阵,再进行特征值分解,得到信号子空间和噪声子空间。利用谱函数求出最优解。上述改进优化的DOA算法不依赖于信号的具体特征,适合于高斯白噪声和高斯色噪声环境。进行权矢量的迭代,从而降低计算复杂度,大幅度缩短了60GHz无线通信系统中DOA计算耗时。
本发明带来的实质性效果是,算法有着更好的特性,收敛速度较快。算法具有以下优点:可抑制高斯白噪声和高斯色噪声;可增加虚拟阵元,实现阵列孔径扩展,提高可估计信号数目;不依赖信号的具体结构。不会因为导频信号占用信道资源而影响效率。
附图说明
图1是本发明的一种计算过程流程图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例:本实施例的一种DOA快速估计的方法,如图1所示,包括以下步骤:
信号预处理阶段:
根据下式(1)构造2q阶累积量矩阵C2q,x;
其中:
B(θ)是广义方向矢量矩阵,B(θ)=[b(θ1),b(θ2),…,b(θk)];
C2q,s(l)是信号s(t)的2q阶累积量矩阵,
是阵列的输入信号的协处理信号,用于信号预测和部分信号的补充;H为矩阵共轭转置符号;
(2)根据下式(2)得到旋转后的2q阶累积量矩阵C2q,x';
其中:
J是一个Mq×Mq的置换阵;
矩阵B(θ)经过旋转后变成JB*(θ),JB*(θ)=B(θ)Φ-(M-1),
矩阵经过旋转后变成
Φ为旋转矩阵,表示为
H为矩阵共轭转置符号。
(3)信号处理及后续阶段:
A.对式(1)与式(2)的累积量矩阵取均值,得到空间平滑处理后的2q阶累积量矩阵为
B.对2q阶累积量矩阵进行平滑处理后,再对C(l)进行特征值分解,得到信号子空间VS和噪声子空间VU。
C.构造空间谱函数,即
其中:
ξ为权向量系数;α(θ)为阵列方向矩阵;
Γ是校正函数,Γ=ρejsinθ+Ωejtgθ,ρ为信号矢量系数,Ω为噪声矢量系数;
m=0,1,2,…,M-1;k=0,1,2,…K-1
j为虚数单位,m为矩阵的行编号,也可以看作是对天线元件的编号;k为矩阵的列编号,也可以看作是对波束的编号,K为波束数量;floor函数表示取小于或者等于括号内指定表达式的最大整数;mod表示取余运算函数,mod(X,Y)为X除以Y所得的余数;
判断空间谱函数是否收敛,如果不收敛,则对权向量进行迭代,q增大1以后重复步骤S02到S05;如果收敛,则从式(5)求出谱函数的θ域的最优解,所得到θ即为入射信号DOA的估计。
按照本发明提供的60GHz无线通信系统DOA快速计算方法,所述DOA快速计算方法还包括对计算出来的信号的误码率采取抑制。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
尽管本文较多地使用了信号子空间、时隙等术语,但并不排除使用其它术语的可能性。使用这些术语仅仅是为了更方便地描述和解释本发明的本质;把它们解释成任何一种附加的限制都是与本发明精神相违背的。
Claims (5)
1.一种高阶累积量来波方向快速估计的方法,其特征在于,包括以下步骤:
S01、天线阵接收原始信号;
S02、构造2q阶累积量矩阵,并对构造的2q阶累积量矩阵进行旋转变换;
S03、对累积量矩阵取均值,得到空间平滑处理后的2q阶累积量矩阵;
S04、进行特征值分解,得到信号子空间和噪声子空间;
S05、构造空间谱函数;
S06、判断空间谱函数是否收敛,如果不收敛,则对权向量进行迭代;如果收敛,则进入步骤S07;
S07、求出空间谱函数谱峰θ域的最优解;
q为初始值为1的正整数,每迭代一次增大1。
2.根据权利要求1所述的一种高阶累积量来波方向快速估计的方法,其特征在于,步骤S02具体为:
A、根据下式(1)构造2q阶累积量矩阵C2q,x;
其中:
B(θ)是广义方向矢量矩阵,B(θ)=[b(θ1),b(θ2),…,b(θk)];
C2q,s(l)是信号s(t)的2q阶累积量矩阵,
是阵列的输入信号的协处理信号,用于信号预测和部分信号的补充;H为矩阵共轭转置符号;
B、根据下式(2)得到旋转后的2q阶累积量矩阵C2q,x';
其中:
J是一个Mq×Mq的置换阵;
矩阵B(θ)经过旋转后变成JB*(θ),JB*(θ)=B(θ)Φ-(M-1),
矩阵经过旋转后变成
Φ为旋转矩阵,表示为
H为矩阵共轭转置符号。
3.根据权利要求2所述的一种高阶累积量来波方向快速估计的方法,其特征在于,步骤S03具体为:对式(1)与式(2)的累积量矩阵取均值,得到空间平滑处理后的2q阶累积量矩阵为
其中
4.根据权利要求3所述的一种高阶累积量来波方向快速估计的方法,其特征在于,步骤S04和步骤S05具体为:
对2q阶累积量矩阵进行平滑处理后,再对C(l)进行特征值分解,得到信号子空间VS和噪声子空间VU;
构造空间谱函数,即
其中:ξ为权向量系数;α(θ)为阵列方向矩阵;Γ是校正函数,Γ=ρejsinθ+Ωejtgθ,ρ为信号矢量系数,Ω为噪声矢量系数;
j为虚数单位,m为矩阵的行编号,也可以看作是对天线元件的编号;k为矩阵的列编号,也可以看作是对波束的编号,K为波束数量;floor函数表示取小于或者等于括号内指定表达式的最大整数;mod表示取余运算函数,mod(X,Y)为X除以Y所得的余数。
5.根据权利要求1所述的一种高阶累积量来波方向快速估计的方法,其特征在于,权向量迭代具体为,q增大1,返回步骤S02。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510603287.1A CN105262550B (zh) | 2015-09-21 | 2015-09-21 | 一种高阶累积量来波方向快速估计的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510603287.1A CN105262550B (zh) | 2015-09-21 | 2015-09-21 | 一种高阶累积量来波方向快速估计的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105262550A true CN105262550A (zh) | 2016-01-20 |
CN105262550B CN105262550B (zh) | 2017-12-22 |
Family
ID=55102080
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510603287.1A Expired - Fee Related CN105262550B (zh) | 2015-09-21 | 2015-09-21 | 一种高阶累积量来波方向快速估计的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105262550B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106019215A (zh) * | 2016-07-15 | 2016-10-12 | 西安电子科技大学 | 基于四阶累量的嵌套阵列波达方向角估计方法 |
CN106226729A (zh) * | 2016-07-15 | 2016-12-14 | 西安电子科技大学 | 基于四阶累量的互质阵列波达方向角估计方法 |
CN108845191A (zh) * | 2018-06-20 | 2018-11-20 | 中国电子科技集团公司第四十研究所 | 一种高斯色噪声通用产生装置及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060292994A1 (en) * | 2005-06-27 | 2006-12-28 | Nec Corporation | Transmitter-receiver of mobile communication system |
CN101369014A (zh) * | 2008-08-08 | 2009-02-18 | 西安电子科技大学 | 应用于多输入多输出雷达的双边约束自适应波束形成方法 |
CN103546221A (zh) * | 2013-10-25 | 2014-01-29 | 东南大学 | 一种宽带相干信号波达角估计方法 |
US20140328487A1 (en) * | 2013-05-02 | 2014-11-06 | Sony Corporation | Sound signal processing apparatus, sound signal processing method, and program |
-
2015
- 2015-09-21 CN CN201510603287.1A patent/CN105262550B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060292994A1 (en) * | 2005-06-27 | 2006-12-28 | Nec Corporation | Transmitter-receiver of mobile communication system |
CN101369014A (zh) * | 2008-08-08 | 2009-02-18 | 西安电子科技大学 | 应用于多输入多输出雷达的双边约束自适应波束形成方法 |
US20140328487A1 (en) * | 2013-05-02 | 2014-11-06 | Sony Corporation | Sound signal processing apparatus, sound signal processing method, and program |
CN103546221A (zh) * | 2013-10-25 | 2014-01-29 | 东南大学 | 一种宽带相干信号波达角估计方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106019215A (zh) * | 2016-07-15 | 2016-10-12 | 西安电子科技大学 | 基于四阶累量的嵌套阵列波达方向角估计方法 |
CN106226729A (zh) * | 2016-07-15 | 2016-12-14 | 西安电子科技大学 | 基于四阶累量的互质阵列波达方向角估计方法 |
CN106019215B (zh) * | 2016-07-15 | 2018-06-05 | 西安电子科技大学 | 基于四阶累量的嵌套阵列波达方向角估计方法 |
CN108845191A (zh) * | 2018-06-20 | 2018-11-20 | 中国电子科技集团公司第四十研究所 | 一种高斯色噪声通用产生装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105262550B (zh) | 2017-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408278A (zh) | 一种基于干扰噪声协方差矩阵估计的稳健波束形成方法 | |
CN102830387B (zh) | 一种基于数据预处理的协方差矩阵正交化波束形成方法 | |
Liu et al. | Joint estimation of DOA and TDOA of multiple reflections in mobile communications | |
Wei et al. | Multiple IRS-assisted wideband dual-function radar-communication | |
Balamurugan et al. | DOA tracking for seamless connectivity in beamformed IoT-based drones | |
CN102401888B (zh) | 一种电磁矢量传感器阵列耦合误差的自校正方法 | |
CN103760527B (zh) | 单基地mimo雷达相干源波达方向估计方法 | |
CN103605117B (zh) | 基于干扰抑制的相控阵信号失真实时校正方法 | |
CN102801455A (zh) | 波束码本生成方法、波束搜索方法及相关装置 | |
CN104035069B (zh) | 基于部分矫正对称均匀线阵的窄带近场信号源定位方法 | |
CN114900400A (zh) | 一种基于智能反射面辅助物联网的联合稀疏信道估计方法 | |
KR20160012284A (ko) | 송신신호 파형의 직교성을 이용한 재머 신호 억제 및 원신호 도래각 추정 방법 및 장치 | |
CN103728601A (zh) | 雷达信号运动干扰空域-极化域联合稳健滤波方法 | |
CN105262550A (zh) | 一种高阶累积量来波方向快速估计的方法 | |
CN106443569A (zh) | 一种基于导向矢量修正的稳健自适应波束形成方法 | |
CN107290732A (zh) | 一种量子大爆炸的单基地mimo雷达测向方法 | |
Yang et al. | Robust adaptive beamforming of coherent signals in the presence of the unknown mutual coupling | |
CN105306117A (zh) | 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 | |
Ghaseminajm et al. | RIS-aided mobile localization error bounds under hardware impairments | |
CN105068090A (zh) | 一种利用单快拍数据进行gnss天线阵干扰抑制的方法 | |
CN107390188A (zh) | 一种用于mimo雷达主瓣抗干扰的发射初相优化方法 | |
CN106599551A (zh) | 一种用于阵列天线足球机器人的快速自适应波束形成算法 | |
CN105262530B (zh) | 一种来波方向快速估计的方法 | |
CN115334524A (zh) | 一种基于全向智能超表面的通信和雷达目标检测方法 | |
CN103792509B (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 | ||
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: 20171222 Termination date: 20180921 |