本发明涉及无线通信领域，特别涉及大规模mimo系统中基于线性搜索的混合结构下到达角估计方法。
背景技术：
由于具有很高的频谱效率和空间谱估计精度，大规模mimo系统引起了学术界和工业界的广泛注意。同时，由于到达角应用范围非常广泛，包括物联网、方向调制、无人机网络、智能交通、无线传感网络等领域，而大规模mimo能够提供较高的空间谱估计精度，使得在到达角估计中应用大规模mimo能够大幅的提高定位精度，保证定位的准确性。
无线定位问题具有悠久的历史，可以追溯到无线通信的开端，而在未来，到达角估计问题将涉及物联网、方向调制系统、无人机、智能交通、无线传感网络等领域，从而具有广阔的市场需求。capon首先提出了以最大化信噪比为目的的功率最大似然估计算法。schmidt在80年代提出了基于特征值分解的高精度music算法。为减少music算法的复杂度，ren提出了一种能够得到到达角估计解析解的根值music算法。
在到达角估计中应用大规模mimo系统所带来的硬件资源损耗也就成为了不可忽视的问题，于是zhang首次提出混合波束成形结构以达到算法性能和成本的平衡。作为一项新兴的波束成形技术，混合波束成形在毫米波通信领域得到了国内外学者的广泛关注和研究。传统的数字波束成形方法需对天线阵列中的每副天线都配备一条专有的rf链路进行单独的数据处理。当天线阵列趋向于大规模发展，硬件成本代价将会巨幅增加。同时考虑到大规模天线阵列的高维度接收数据，传统的数字波束形成方法计算复杂度高，难以满足实际应用高实时性的要求。而模拟波束成形可仅采用一条rf链路处理接收信号，然而其性能往往难以与数字波束成形相比拟。混合波束成形通常采用远小于天线数的rf链路以降低系统开销，同时采用大量的相移器增加天线阵列的增益，能够实现系统性能与硬件成本的折衷，从而成为了5g毫米波通信系统的主要技术之一。
考虑到现有的移相器精度普遍较高，采用线性搜索改变移相器的相移，以使得接收来波信号功率最大是恰当的到达角估计方法。
技术实现要素：
发明目的：在大规模天线系统中，混合波束成形器能有效地实现系统性能与硬件成本的折衷。通过利用混合波束成形器的优势，本发明提供大规模mimo系统中基于线性搜索的混合结构下到达角估计方法，能有效提高定位精度，同时降低定位算法的复杂度。
技术方案：为实现上述目的，本发明采用的技术方案为：
一种基于线性搜索混合波束成形方法，具体过程包括：
(1)天线阵列划分子阵：
考虑由n个全向阵元组成的线性均匀天线阵列，且该天线阵列位于信号源的远场范围。将阵列均匀地分为k个子阵，子阵天线数目均为m，即n＝km。天线阵列的导向向量表示为：
第k个子阵的相移矢量为
其中，va,k,m表示第k个子阵中第m个阵元的相移因子。整个天线阵列所形成的波束方向图可表示为
其中，vk表示第k个子阵的权值。对于混合结构下的波束成形，将分别设计模拟波束成形矩阵va与数字波束成形矢量vd。
(2)模拟波束成形矩阵设计：
假设来波信号源为窄带信号，来波方向为θ0。设计模拟波束成形矩阵va，可使移相器对准接收功率最大的方向。假设每个子阵的相移矢量均相同，即，
其中，
θ为
则n×k维的模拟波束成形矩阵va为
(3)数字波束成形向量设计：
考虑到模拟波束成形矩阵设计中，每个子阵的相移均相同，于是，在数字波束成形端，我们对每一个输入进行不同的相移，即，数字波束成形向量设计为：
其中，
同样的，θ为
设来波信号为
其中，αk,m为第k个子阵的第m根天线对应的模拟波束成形矩阵的相移，τk,m是由来波信号方向决定的传播时延，
τ0是从发射器到阵列参考点的传播时延，c为光速，d是阵列距离。将k个子阵的输出写成向量的形式则可以得到，
其中
经过数字波束成形，输出信号变为
根据设计完成的模拟波束成形矩阵和数字波束成形向量，输出信号转化成
假设快拍数为l，我们可以得到平均接收功率
当θ在
有益效果：本发明提供的在大规模mimo系统中基于线性搜索的混合结构下到达角估计方法，具有如下优点：1、本方法对大规模天线系统能有效地实现系统性能与硬件成本的折衷；2.本方法可有效估计来波信号到达角；3.本方法使得模拟波束成形矩阵的调相个数从n减少到m，大大降低了算法的复杂度。
本发明附加的方面和优点将在下面的描述中部分给出，这些将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
图1示出了大规模mimo系统中基于线性搜索的混合结构下到达角估计方法的系统流程图。
图2示出了混合架构下传统模拟相位对齐与基于线性搜索的到达角估计方法在低信噪比时角度估计误差随着搜索步长的变化图。
图3示出了混合架构下传统模拟相位对齐与基于线性搜索的到达角估计方法复杂度比较图。
具体实施方式
下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明提供了一种大规模mimo系统中基于线性搜索的混合结构下到达角估计方法，本发明中，将模拟波束成形与数字波束成形相结合，利用线性搜索进行模拟波束成形和数字波束成形设计，从而将波束调向使得接收信号功率最大的方向，从而有效估计出来波信号到达角。本发明的混合波束成形算法能有效地实现系统性能与硬件成本的折衷。并且，可有效估计来波信号方向角，同时降低估计复杂度。
(1)天线阵列划分子阵：
考虑由n个全向阵元组成的线性均匀天线阵列，且该天线阵列位于信号源的远场范
当θ在
有益效果：本发明提供的在大规模mimo系统中基于线性搜索的混合结构下到达角估计方法，具有如下优点：1、本方法对大规模天线系统能有效地实现系统性能与硬件成本的折衷；2.本方法可有效估计来波信号到达角；3.本方法使得模拟波束成形矩阵的调相个数从n减少到m，大大降低了算法的复杂度。
本发明附加的方面和优点将在下面的描述中部分给出，这些将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
图1示出了大规模mimo系统中基于线性搜索的混合结构下到达角估计方法的系统流程图。
图2示出了混合架构下传统模拟相位对齐与基于线性搜索的到达角估计方法在低信噪比时角度估计误差随着搜索步长的变化图。
图3示出了混合架构下传统模拟相位对齐与基于线性搜索的到达角估计方法复杂度比较图。
具体实施方式
下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明提供了一种大规模mimo系统中基于线性搜索的混合结构下到达角估计方法，本发明中，将模拟波束成形与数字波束成形相结合，利用线性搜索进行模拟波束成形和数字波束成形设计，从而将波束调向使得接收信号功率最大的方向，从而有效估计出来波信号到达角。本发明的混合波束成形算法能有效地实现系统性能与硬件成本的折衷。并且，可有效估计来波信号方向角，同时降低估计复杂度。
(1)天线阵列划分子阵：
考虑由n个全向阵元组成的线性均匀天线阵列，且该天线阵列位于信号源的远场范围。将阵列均匀地分为k个子阵，子阵天线数目均为m，即n＝km。天线阵列的导向矢量表示为：
第k个子阵的相移矢量为
其中，va,k,m表示第k个子阵中第m个阵元的相移因子。整个天线阵列所形成的波束方向图可表示为
其中，vk表示第k个子阵的权值。对于混合结构下的波束成形，将分别设计模拟波束成形矩阵va与数字波束成形矢量vd。
(2)模拟波束成形矩阵设计：
假设来波信号源为窄带信号，来波方向为θ0。设计模拟波束成形矩阵va，可使移相器对准接收功率最大的方向。假设每个子阵的相移矢量均相同，即，
其中，
θ为
则n×k维的模拟波束成形矩阵为
(3)数字波束成形向量设计：
考虑到模拟波束成形矩阵设计中，每个子阵的相移均相同，于是，在数字波束成形端，我们对每一个输入进行不同的相移，即，数字波束成形向量设计为：
其中，
同样的，θ为
设来波信号为
其中，αk,m为第k个子阵的第m根天线对应的模拟波束成形矩阵的相移，τk,m是由来波信号方向决定的传播时延，
τ0是从发射器到阵列参考点的传播时延，c为光速，d是阵列距离。将k个子阵的输出写成矩阵的形式则可以得到，
其中
经过数字波束成形，输出信号变为
根据设计完成的模拟波束成形矩阵和数字波束成形向量，输出信号转化成
假设快拍数为l，我们可以得到平均接收功率
当θ在
图1示出了大规模mimo系统中基于线性搜索的混合结构下到达角估计方法的系统流程图。
图2反映了信号源的来波方向为41.165°，信噪比为0db，总的天线数为n＝128，子阵个数k＝16时，混合架构下传统模拟相位对齐与基于线性搜索的到达角估计方法角度估计误差随着搜索步长的变化图。从图中可以看出无论是模拟相位对齐还是本发明提出的基于线性搜索的到达角估计方法，角度估计误差均随着搜索步长的减小而缩小，并且本发明提出的基于线性搜索的算法估计精度优于传统的模拟相位对齐。
图3示出了总的天线数为n＝128，子阵个数k＝16时，混合架构下传统模拟相位对齐与基于线性搜索的到达角估计方法的复杂度比较。从图中可以看出，与传统模拟相位对齐相比，本发明提出的算法具有更低的复杂度。