[发明专利]一种在毫米波预编码系统中的快速迭代波束成形方法

摘要

本发明属于无线通信技术领域，它特别涉及一种在无线多输入多输出(Multiple Input Multiple Output，MIMO)通信系统中利用时分双工系统的信道互易性进行迭代从而同时快速地获得多个波束成形奇异矢量的方法。本发明提出了一种基于LANCZOS算法的快速迭代天线训练方法，由于LANCZOS方法的收敛速度快于幂迭代方法，而且可以只用一个阶段的迭代就能获得多流波束成形的多对收发端波束成形奇异矢量，从而成倍地减少天线训练阶段的开销，快速地完成天线训练，本发明具有重大的实际意义。

主权项

一种在毫米波预编码系统中的快速迭代波束成形方法，其特征在于，步骤如下：S1.定义接收端字典矩阵定义发送端字典矩阵其中，N表示接收端字典长度，M表示接收端字典长度；S2.进行角度量化码本的建立，接收端码本为发射端码本为其中为码本大小，N_T为发送天线数目，N_R为接收天线数目；S3.初始阶段处理，具体如下：S31.发送端生成一个N_T×1向量[1,0,0,…,0]^T，并放入码本中使用稀疏信号恢复算法估计得到初始向量r₀，将r₀的模值赋给变量β₀：β₀＝||r₀||，定义零向量q₀：q₀＝0，定义空矩阵Q用于存储迭代过程中产生的向量：Q＝[]，其中，N_T为接收天线数目；S32.分别定义此过程中的发送端和接收端的测量次数Nmr⁰,Nmt⁰；S33.在O个随机矩阵中选取发送端最优测量矩阵Φ_T⁰和选取接收端的最优测量矩阵Φ_R⁰，其中，所述O为不为零的自然数；S34.接收波束成形向量训练，具体如下：S34‑1.计算初始向量q₁：令S21所述空矩阵Q的第1列为q₁，即Q＝[Q,q₁]；发送方在连续发送Nmr⁰次向量q₁至接收方，接收端接收过程中使用Φ_R⁰的列作为波束成形加权合并向量，接收端得到其中，n_R表示接收端的加性高斯白噪声向量，S34‑2.接收端使用基于压缩感知的稀疏信号恢复算法计算出表示接收信号到达角在字典矩阵A_RD中的位置的稀疏向量z_R，其中，z_R是一个N×1的列向量，N表示字典A_RD的长度，z_R中有K个非零元素，K＜＜N；S34‑3.Hq≈A_RDz_R，所述Hf存储在N_R×1向量g中，即g＝A_RDz_R，其中，信道矩阵接收端对向量g进行归一化，即并将返回至发送端；S35.发送波束成形向量训练，具体如下：S35‑1.接收端发送向量至发送端，发送端在接收过程中使用Φ_T⁰的列作为波束成形加权合并向量，因此发送端得到其中，n_T表示第k次迭代发送端的加性高斯白噪声向量，S35‑2.发送端使用基于压缩感知的稀疏信号恢复算法计算出表示接收信号到达角在字典矩阵A_TD中的位置的稀疏向量z_T，其中，z_T是一个M×1的列向量，M表示字典A_TD的长度，z_T中有K个非零元素，K＜＜M；S35‑3.所述存储在N_T×1向量f中，即f＝A_TDz_T，对f取共轭，结果仍然存储在f中：f＝f^*。S35‑4.发送方计算需要用于构建三对角矩阵的参数α₁和向量r₁：然后进行重正交化r₁，重正交化方法如下：计算用于构建三对角矩阵的参数β₁：β₁＝||r₁||。S4.迭代过程，定义迭代次数为N_ITER，定义迭代循环控制变量为k，初始化k＝1，具体如下：S41.分别定义此过程中的发送端和接收端的测量次数Nmr,Nmt；S42.找到S24‑2，S25‑2中K个非零元素的位置所对应的角度，并将K个角度放入个相位中找到与之最近的相位，最后得到相位将其转置即为发送端和接收端的测量矩阵。其中，个相位是将‑π～π进行个量化；S43.接收波束成形向量训练，具体如下：发送端发送向量至接收端，空矩阵Q＝[Q,q_k+1],接收端接收过程中使用Φ_R的列作为波束成形加权合并向量，因此接收端得到r＝Φ_R^H(Hq_k+1+n_R)，接收端使用最小二乘法计算出系数h_R＝(Φ_RA_R)^‑1r，并求得v＝A_Rh_R，并将v放入码本中使用稀疏信号恢复算法进行估计，接收端对向量v进行归一化，即并将返回至发送端；S44.发送波束成形向量训练，具体如下：接收端发送向量至发送端，发送端接收过程中使用Φ_T的列作为波束成形加权合并向量，因此接收端得到接收端使用最小二乘法计算出系数h_T＝(Φ_TA_T)^‑1t，并求得f'＝A_Th_T，并将f'放入码本中使用稀疏信号恢复算法进行估计，对f'取共轭，结果仍然存储在f'中：f'＝f'^*。S35.发送方计算需要用于构建三对角矩阵的参数α_k+1和向量r_k+1：然后进行重正交化r_k+1，重正交化方法如下：计算用于构建三对角矩阵的参数β_k+1：β_k+1＝||r_k+1||，令k＝k+1，返回S31继续迭代；S4.利用S35迭代过程中得到的α_k+1和β_k+1构建三对角矩阵S5.对S4所述T作特征值分解，并将特征值按降序排列，从大到小依次为λ₁,λ₂,...,λ_ITER，相应的特征向量为v₁,v₂,...,v_ITER；S6.计算发送波束成形矩阵将发送方波束成形矩阵F的列f_k(k＝1,2,...N_S)依次发送至接收方，接收端接收过程中使用Φ_R的列作为波束成形加权合并向量，因此接收端得到r＝Φ_R^H(Hf_k+n_R)，接收端使用最小二乘法计算出系数h_R＝(Φ_RA_R)^‑1r，并求得w＝A_Rh_R，并将w放入码本中使用稀疏信号恢复算法进行估计，接收端对向量w进行归一化，即合并为接收端的波束成形矩阵S7.最后将迭代后得到的W,F输出即可。