一种来波方向快速估计的方法
【技术领域】
[0001] 本发明涉及通信技术领域，尤其是涉及一种应用于60GHz无线通信系统的来波方 向快速估计的方法。
【背景技术】
[0002] 近年来移动通信事业在全球范围内得到了迅猛的发展，对带宽的要求越来越高。 60GHz无线通信技术由于其超高速的数据传输能力、成本低、抗干扰性和安全性高、具有国 际通用性和免许可特性的特点，将成为第4代无线通信技术的代表。然而它面临的最大挑 战是过高的路径损失。60GHz频段的路径损耗决定了其必须采用波束形成技术和来波方向 快速估计来改善链路质量。来波方向（D0A)估计又可称为来波方向估计或来波角度估计， 是阵列信号处理的一个重要分支。D0A估计的基本问题就是如何利用天线阵列的输出，从背 景噪声中检测出同时处于空问某一区域内的多个期望用户的位置信息，即各个信号到达阵 列天线的入射角。D0A估计算法性能评价标准主要包括：估计偏差、估计误差、分辨力、成功 概率、计算量及给定阵元数目情况下算法可估计的信号数目。
[0003] D0A是阵列信号处理的一项重要内容，其应用涉及雷达、通信、声纳、地震、勘探、射 电天文以及生物医学工程等众多军事及民用经济领域。D0A估计是利用一组按一定方式布 置在空间不同位置的传感器对空间信源在时域和空域同时进行采样，再由传感器阵列采样 数据的分析处理来实现对空间信源的方位向估计。在D0A估计中，分辨力，精度，稳健性，速 度，可分辨目标数的多少等，决定了方位向估计性能的好坏。近几十年，国内外均开展了有 关D0A估计技术的理论和应用研究，使得D0A估计理论得到迅速发展。但是，在实际应用中， 由于受到物理器件水平和成本的影响，阵列孔径不可能无限增大，D0A估计的分辨力始终受 到阵列孔径的限制，最大可分辨目标数受到阵元数的限制。
[0004] 中华人民共和国国家知识产权局于2007年04月25日公开了公告号为 CN1953348A的专利文献，名称是多扇区智能天线系统的来波方向估计方法，其包括以下步 骤：（1)利用谱估计算法对需要进行来波方向估计的用户进行初始来波方向估计，得到估 计的来波功率谱最大值对应的来波角度；（2)根据所述来波角度判断用户是否位于目标扇 区的边缘，若是，则分别计算用户在和目标扇区相邻的两个扇区内的功率；（3)若所述来波 角度所在区域靠近功率值大的相邻扇区，则所述来波角度即为最终的估计结果，否则，搜索 来波角度所在区域的相反区域，得到所述相反区域内来波功率谱最大值对应的来波角度即 为最终的估计结果。此方案仍然收到阵元数的限制，在分辨率和分辨度上都不够理想。
【发明内容】
[0005] 本发明主要是解决现有技术所存在的对信号入射角相差较小的信号分辨率低、对 入射信号低信噪比下分辨率不高、对入射信号的分辨度不高等的技术问题，提供一种应用 于60GHz无线通信系统的D0A的快速计算方法，能在改善无线接收器、发送器性能的基础上 降低波束形成算法的计算复杂度，提高分辨率，从而在保证高速通讯信号品质的同时，大幅 度缩短了 DOA计算耗时。
[0006] 该D0A改进算法，不仅具有空间平滑算法去强相干信号的效果，同时弥补其在低 信噪比下分辨率不高和对信号的分辨度不高的不足之处，另外，该算法也比空间平滑算法 的算法复杂度更小，大大减少了整个算法的运算量，而且不牺牲天线的孔径，比较实用，实 验证明了该算法的有效性。
[0007] 本发明针对上述技术问题主要是通过下述技术方案得以解决的：一种来波方向快 速估计的方法，包括以下步骤：
[0008] S01、天线阵接收原始信号；
[0009] S02、对信号进行预处理；
[0010] S03、构造协方差矩阵；
[0011] S04、权向量迭代；
[0012] S05、构造空间谱函数；
[0013] S06、求出空间谱函数谱峰Θ域的最优解；
[0014] S07、判断空间谱函数是否收敛，如果收敛，则所得最优解即为入射信号的估计值； 如果不收敛，则返回步骤S04。
[0015] 假设有K个信号，其中前Kc个信号为相关信号，其余信号为独立信号且与前Kc个 信号相互独立，入射到均匀天线阵，该天线阵由Μ = 2N+1个间距为d的全向阵元组成，并将 阵元〇作为参考阵元，则第i个阵元的输出信号矢量可表示为
[0017] 式中，sk(t)为第k个信号的复包络，λ为中心波长，d为阵元间距，P 多径衰 落系数，ni(t)是加性高斯噪声。
[0018] 则阵列的输出信号矢量可表示为
[0019] X (t) = [x N (t), ···, X ! (t), x0 (t), Xi (t), ···, xN (t) ]T
[0020] = A ( Θ ) s (t) +n (t) (2)
[0021] 式中，A( Θ ) = [a ( Θ丄a ( Θ 2), ···，a ( Θ k)]为阵列方向矩阵，s(t)为阵列接收 信号矢量，n(t)为噪声矢量。其中，a (0k)对应的方向矢量为
[0023] 根据式（2)估计阵列输出信号矢量的协方差矩阵，即
[0025] 其中：
[0026] N为全向阵元总个数的一半，总全向阵元个数为2N+1 (其中包括作为参考阵元的0 阵元）；
[0027] x = [x(-N)，…，x(0),…，x(N)]是阵列的输出信号矢量；
是阵列的输出信号的协处理信号，用于信号预测和 部分信号的补充；Η为矩阵共辄转置符号；T为矩阵转置符号；
[0029] 构造新的协方差矩阵Rx，即
[0031] 式中：J为MXM的反向单位矩阵，巧为息的共辄。
[0033] 特征分解上面得到的协方差矩阵Rx，即
[0034] [U, S, V] = SVD (Rx)
[0035] 取￥"=1](:，0+1满来构造噪声子空间向量Vn。SVD表示奇异值分解，U(:，D+1:M) 表示矩阵U中的任何一行，其列从D+1到Μ。
[0036] 把满秩矩阵馬用低秩矩阵替代，构造出新的协方差矩阵，然后对新的协方差矩 阵特征分解，再次构造得到噪声子空间向量Vnn。
[0037] 计算得到的平均噪声子空间，设为％，即
[0039] 构造空间谱函数，即
[0041] 其中：
[0042] ξ为权向量系数；α ( Θ )为阵列方向矩阵；
[0043] Γ是校正函数，
Ρ为信号矢量系数，Ω为噪声矢量系数；
[0045] j为虚数单位，m为矩阵的行编号，也可以看作是对天线元件的编号；k为矩阵的列 编号，也可以看作是对波束的编号，K为波束数量；floor函数表示取小于或者等于括号内 指定表达式的最大整数;mod表示取余运算函数，mod (X，Y)为X除以Y所得的余数；
[0046] 从式（5)求出户(扔谱峰Θ域的最优解，在空间谱函数收敛的情况下所得到谱峰 Θ即为入射信号D0A的估计，如果不收敛则继续迭代计算。
[0047] 本发明提供的60GHz无线通信系统D0A快速计算方法，所述D0A快速计算方法应 用在超高带宽60GHz无线接收器的控制程序中。
[0048] 在60GHz无线通信系统的实际应用中，由于外部干扰阻挡的存在，如发送和接收 之间瞬间被障碍物阻挡等等，天线接收到的信号中可能存在某种比较少见、一般情况不考 虑的干扰信号，以及其他无法预知的有效信号迅速衰减，通过本方法的处协理信号火/)、权 矢量中的
校正函数的综合的补偿和调整作用，可以有效消除这些干 扰。
[0049]