本发明专利技术属于无线通信技术领域，它特别涉及一种在无线多输入多输出(Multiple Input Multiple Output，MIMO)通信系统中利用时分双工系统的信道互易性进行迭代从而同时快速地获得多个波束成形奇异矢量的方法。本发明专利技术提出了一种基于LANCZOS算法的快速迭代天线训练方法，由于LANCZOS方法的收敛速度快于幂迭代方法，而且可以只用一个阶段的迭代就能获得多流波束成形的多对收发端波束成形奇异矢量，从而成倍地减少天线训练阶段的开销，快速地完成天线训练，本发明专利技术具有重大的实际意义。
【技术实现步骤摘要】
本专利技术属于无线通信
，它特别涉及一种在无线多输入多输出(Multiple Input Multiple Output，MIMO)通信系统中利用时分双工系统的信道互易性进行迭代从而同时快速地获得多个波束成形奇异矢量的方法
技术介绍
在MIMO系统中，根据接收端最大化信噪比准则，获得最优波束成形矩阵的方法是特征波束成形方法。在收发双方都已知信道状态信息(CSI)的情况下，最优的发送和接收波束成形矩阵可以通过对信道矩阵H进行SVD分解得到。具体原理叙述如下：假设MIMO系统的收发天线数目分别为NT、NR，信道矩阵可以进行SVD分解，表示为H＝UΛVH，其中，(·)H表示矩阵共轭转置，和分别是大小为NR×NR与NT×NT的酉矩阵，Λ是一个NR×NT对角阵，其对角元为按降序排列的H的奇异值(σ1,σ2,...σm),m＝min(NT,NR)。对于NS(NS≤m)维的波束成形，发送端与接收端波束成形矩阵分别采用H的右奇异矩阵V和左奇异矩阵U的前m列，即F＝[v1,v2,...,vm]，W＝[u1,u2,...,um]。假设发送符号为x＝[x1,x2,...,xm]T，接收符号为y＝[y1,y2,...,ym]T，噪声则可见，特征波束成形等效地将MIMO信道划分为m个并行独立的子信道，每个子信道都获得了最大化的信噪比。通常，接收端通过估计信道矩阵H并进行SVD分解来获得收发双方的波束成形矩阵，之后接收端将发送端的波束成形矩阵F反馈至发送端。这种直接估计和反馈的方法适用于天线数目较小的情况，而在天线数目较多的MIMO系统中，其计算复杂度和训练开销都变得无法承受。在时分双工(TDD)系统中，利用上行信道和下行信道的互易性，文献Yang Tang,Branka Vucetic,Yonghui Li.An Iterative Singular Vectors Estimation Scheme for Beamforming Transmission and Detection in MIMO Systems.IEEE Communications Letters,VOL.9,NO.6,June2005.提出了一种不用估计信道参数即可获得特征向量的迭代波束成形方法，即幂迭代方法。文献Pengfei Xia,Su-Khiong Yong,Jisung Oh and Chiu Ngo.Multi-Stage Iterative Antenna Training for Millimeter Wave Communications.IEEE Globecom Conference 2008.进一步将这种方法扩展到了多维的波束成形，即通过逐个阶段剥离的方式得到NS个波束成形矢量，也就是波束成形矩阵，每个阶段都要经历一轮幂迭代。在不增加发送功率的条件下提高MIMO系统容量的有效方法是采用预编码技术并行地发送多个数据流。具有NS个数据流的MIMO系统的波束成形需要NS对收发端波束成形矢量，这些波束成形矢量都通过天线训练获得。幂迭代天线训练算法一个阶段的迭代只能获得一对波束成形奇异矢量，如果要完成NS流的波束成形天线训练，必须经过NS个阶段的迭代。在天线数目较大和数据流较多的系统中，开销无疑是巨大的，天线训练过程会持续较长时间。众所周知，天线训练的一个基本假设是在一个合理的短时间内信道状态不变，如果训练时间持续过长，信道状态发生较大改变的概率就会大大增加，就会使这一假设失去作用，在实际应用中必然不能达到预期的效果。另外，除了第一阶段的迭代，以后的每一阶段的迭代都要进行零空间投影操作，该操作引入了额外的误差和开销。
技术实现思路
LANCZOS算法较幂迭代方法具有更快速的收敛特性，而且在一个阶段的迭代中可以同时得到多个特征向量，所以，对于多流的波束成形来说，可以用幂迭代方法单流波束成形的开销完成多流波束成形的天线训练，相当于将训练开销降低到幂迭代方法的，可以快速地完成多流的波束成形天线训练，优势十分明显。为了克服多流MIMO系统中幂迭代方法一轮迭代只能获得一个波束成形奇异矢量的缺陷，本专利技术提出了一种基于LANCZOS算法的快速迭代天线训练方法，由于LANCZOS方法的收敛速度快于幂迭代方法，而且可以只用一个阶段的迭代就能获得多流波束成形的多对收发端波束成形奇异矢量，从而成倍地减少天线训练阶段的开销，快速地完成天线训练，本专利技术具有重大的实际意义。本专利技术的原理：利用LANCZOS算法的快速收敛特性和在一个阶段的迭代中可以获得多个特征向量的性质，本专利技术设计了一种时分双工系统下的快速迭代波束成形方法，可以大大降低天线训练阶段的开销，快速地得到多流波束成形的多个波束成形矢量，此系统是用于硬件平台，每个天线元素只需调节相位而不需调整幅度。为了方便地描述本专利技术的内容，首先对本专利技术中所使用的概念和术语进行定义。空间稀疏性：无线信号由于较高的路径损耗和极差的散射性能，收发双方只由有限的几条电磁波传播路径相连接，和信道有关的信号计算问题可以方便地表达为稀疏重建问题。稀疏多径信道模型：稀疏多径信道可以建模为具有K路多径的几何模型其中，表示第i径的复信道增益，θi表示第i径的离开角，φi表示第i径的到达角，aT(φi)是发送端的天线阵列响应，aR(θi)是接收端的天线阵列响应,i＝1,2,...,K。所述天线阵列采用均匀线性阵列(ULAs)，则发送端的天线阵列响应可以表达成接收端的天线阵列响应可以表达成其中，λ是信号波长，d是天线阵元间距，一般取一种在毫米波预编码系统中的快速迭代波束成形方法，步骤如下：S1.定义接收端字典矩阵定义发送端字典矩阵其中，N表示接收端字典长度，M表示接收端字典长度；S2.进行角度量化码本的建立，接收端码本为发射端码本为其中为码本大小，NT为发送天线数目，NR为接收天线数目；S3.初始阶段处理，具体如下：S31.发送端生成一个NT×1向量[1,0,0,…,0]T，并放入码本中使用稀疏信号恢复算法估计得到初始向量r0，将r0的模值赋给变量β0：β0＝||r0||，定义零向量q0：q0＝0，定义空矩阵Q用于存储迭代过程中产生的向量：Q＝[]，其中，NT为接收天线数目；S32.分别定义此过程中的发送端和接收端的测量次数Nmr0,Nmt0；S33.在O个随机矩阵中选取发送端最优测量矩阵ΦT0和选取接收端的最优测量矩阵ΦR0，其中，所述O为不为零的自然数；S34.接收波束成形向量训练，具体如下：S34-1.计算初始向量q1：令S21所述空矩阵Q的第1列为q1，即Q＝[Q,q1]；发送方在连续发送Nmr0次向量q1至接收方，接收端接收过程中使用ΦR0的列作为波束成形加权合并向量，接收端得到其中，nR表示接收端的加性高斯白噪声向量，S34-2.接收端使用基于压缩感知的稀疏信号恢复算法计算出表示接收信号到达角在字典矩阵ARD中的位置的稀疏向量zR，其中，zR是一个N×1的列向量，N表示字典ARD的长度，zR中有K个非零元素，K＜＜N；S34-3.Hq≈ARDzR，所述Hf存储在NR×1向量g中，即g＝ARDzR，其中，信道矩阵接收端对向量g进行归一化，即并将返回至发送端；S35.发送波束成形向量训练，具体如下：S35-1.接收端发送向本文档来自技高网...
【技术保护点】
一种在毫米波预编码系统中的快速迭代波束成形方法，其特征在于，步骤如下：S1.定义接收端字典矩阵定义发送端字典矩阵其中，N表示接收端字典长度，M表示接收端字典长度；S2.进行角度量化码本的建立，接收端码本为发射端码本为其中为码本大小，NT为发送天线数目，NR为接收天线数目；S3.初始阶段处理，具体如下：S31.发送端生成一个NT×1向量[1,0,0,…,0]T，并放入码本中使用稀疏信号恢复算法估计得到初始向量r0，将r0的模值赋给变量β0：β0＝||r0||，定义零向量q0：q0＝0，定义空矩阵Q用于存储迭代过程中产生的向量：Q＝[]，其中，NT为接收天线数目；S32.分别定义此过程中的发送端和接收端的测量次数Nmr0,Nmt0；S33.在O个随机矩阵中选取发送端最优测量矩阵ΦT0和选取接收端的最优测量矩阵ΦR0，其中，所述O为不为零的自然数；S34.接收波束成形向量训练，具体如下：S34‑1.计算初始向量q1：令S21所述空矩阵Q的第1列为q1，即Q＝[Q,q1]；发送方在连续发送Nmr0次向量q1至接收方，接收端接收过程中使用ΦR0的列作为波束成形加权合并向量，接收端得到其中，nR表示接收端的加性高斯白噪声向量，S34‑2.接收端使用基于压缩感知的稀疏信号恢复算法计算出表示接收信号到达角在字典矩阵ARD中的位置的稀疏向量zR，其中，zR是一个N×1的列向量，N表示字典ARD的长度，zR中有K个非零元素，K＜＜N；S34‑3.Hq≈ARDzR，所述Hf存储在NR×1向量g中，即g＝ARDzR，其中，信道矩阵接收端对向量g进行归一化，即并将返回至发送端；S35.发送波束成形向量训练，具体如下：S35‑1.接收端发送向量至发送端，发送端在接收过程中使用ΦT0的列作为波束成形加权合并向量，因此发送端得到其中，nT表示第k次迭代发送端的加性高斯白噪声向量，S35‑2.发送端使用基于压缩感知的稀疏信号恢复算法计算出表示接收信号到达角在字典矩阵ATD中的位置的稀疏向量zT，其中，zT是一个M×1的列向量，M表示字典ATD的长度，zT中有K个非零元素，K＜＜M；S35‑3.所述存储在NT×1向量f中，即f＝ATDzT，对f取共轭，结果仍然存储在f中：f＝f*。S35‑4.发送方计算需要用于构建三对角矩阵的参数α1和向量r1：然后进行重正交化r1，重正交化方法如下：计算用于构建三对角矩阵的参数β1：β1＝||r1||。S4.迭代过程，定义迭代次数为NITER，定义迭代循环控制变量为k，初始化k＝1，具体如下：S41.分别定义此过程中的发送端和接收端的测量次数Nmr,Nmt；S42.找到S24‑2，S25‑2中K个非零元素的位置所对应的角度，并将K个角度放入个相位中找到与之最近的相位，最后得到相位将其转置即为发送端和接收端的测量矩阵。其中，个相位是将‑π～π进行个量化；S43.接收波束成形向量训练，具体如下：发送端发送向量至接收端，空矩阵Q＝[Q,qk+1],接收端接收过程中使用ΦR的列作为波束成形加权合并向量，因此接收端得到r＝ΦRH(Hqk+1+nR)，接收端使用最小二乘法计算出系数hR＝(ΦRAR)‑1r，并求得v＝ARhR，并将v放入码本中使用稀疏信号恢复算法进行估计，接收端对向量v进行归一化，即并将返回至发送端；S44.发送波束成形向量训练，具体如下：接收端发送向量至发送端，发送端接收过程中使用ΦT的列作为波束成形加权合并向量，因此接收端得到接收端使用最小二乘法计算出系数hT＝(ΦTAT)‑1t，并求得f'＝AThT，并将f'放入码本中使用稀疏信号恢复算法进行估计，对f'取共轭，结果仍然存储在f'中：f'＝f'*。S35.发送方计算需要用于构建三对角矩阵的参数αk+1和向量rk+1：然后进行重正交化rk+1，重正交化方法如下：计算用于构建三对角矩阵的参数βk+1：βk+1＝||rk+1||，令k＝k+1，返回S31继续迭代；S4.利用S35迭代过程中得到的αk+1和βk+1构建三对角矩阵S5.对S4所述T作特征值分解，并将特征值按降序排列，从大到小依次为λ1,λ2,...,λITER，相应的特征向量为v1,v2,...,vITER；S6.计算发送波束成形矩阵将发送方波束成形矩阵F的列fk(k＝1,2,...NS)依次发送至接收方，接收端接收过程中使用ΦR的列作为波束成形加权合并向量，因此接收端得到r＝ΦRH(Hfk+nR)，接收端使用最小二乘法计算出系数hR＝(ΦRAR)‑1r，并求得w＝ARhR，并将w放入码本中使用稀疏信号恢复算法进行估计，接收端对向量w进行归一化，即合并为接收端的波束成形矩阵S7.最后将迭代后得到的W,F输出即可。...
【技术特征摘要】
1.一种在毫米波预编码系统中的快速迭代波束成形方法，其特征在于，步骤如下：S1.定义接收端字典矩阵定义发送端字典矩阵其中，N表示接收端字典长度，M表示接收端字典长度；S2.进行角度量化码本的建立，接收端码本为发射端码本为其中为码本大小，NT为发送天线数目，NR为接收天线数目；S3.初始阶段处理，具体如下：S31.发送端生成一个NT×1向量[1,0,0,…,0]T，并放入码本中使用稀疏信号恢复算法估计得到初始向量r0，将r0的模值赋给变量β0：β0＝||r0||，定义零向量q0：q0＝0，定义空矩阵Q用于存储迭代过程中产生的向量：Q＝[]，其中，NT为接收天线数目；S32.分别定义此过程中的发送端和接收端的测量次数Nmr0,Nmt0；S33.在O个随机矩阵中选取发送端最优测量矩阵ΦT0和选取接收端的最优测量矩阵ΦR0，其中，所述O为不为零的自然数；S34.接收波束成形向量训练，具体如下：S34-1.计算初始向量q1：令S21所述空矩阵Q的第1列为q1，即Q＝[Q,q1]；发送方在连续发送Nmr0次向量q1至接收方，接收端接收过程中使用ΦR0的列作为波束成形加权合并向量，接收端得到其中，nR表示接收端的加性高斯白噪声向量，S34-2.接收端使用基于压缩感知的稀疏信号恢复算法计算出表示接收信号到达角在字典矩阵ARD中的位置的稀疏向量zR，其中，zR是一个N×1的列向量，N表示字典ARD的长度，zR中有K个非零元素，K＜＜N；S34-3.Hq≈ARDzR，所述Hf存储在NR×1向量g中，即g＝ARDzR，其中，信道矩阵接收端对向量g进行归一化，即并将返回至发送端；S35.发送波束成形向量训练，具体如下：S35-1.接收端发送向量至发送端，发送端在接收过程中使用ΦT0的列作为波束成形加权合并向量，因此发送端得到其中，nT表示第k次迭代发送端的加性高斯白噪声向量，S35-2.发送端使用基于压缩感知的稀疏信号恢复算法计算出表示接收信号到达角在字典矩阵ATD中的位置的稀疏向量zT，其中，zT是一个M×1的列向量，M表示字典ATD的长度，zT中有K个非零元素，K＜＜M；S35-3.所述存储在NT×1向量f中，即f＝ATDzT，对f取共轭，结果仍然存储在f中：f＝f*。S35-4.发送方计算需要用于构建三对角矩阵的参数α1和向量r1：然后进行重正交化r1，重正交化方法如下：计算用于构建三对角矩...
【专利技术属性】
技术研发人员：娄念念，成先涛，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。