CN105959044B - 联合方法的层次码本结构设计方法 - Google Patents
联合方法的层次码本结构设计方法 Download PDFInfo
- Publication number
- CN105959044B CN105959044B CN201610251727.6A CN201610251727A CN105959044B CN 105959044 B CN105959044 B CN 105959044B CN 201610251727 A CN201610251727 A CN 201610251727A CN 105959044 B CN105959044 B CN 105959044B
- Authority
- CN
- China
- Prior art keywords
- layer
- codebook
- search
- code word
- codeword
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 122
- 238000013461 design Methods 0.000 title claims abstract description 31
- 238000010276 construction Methods 0.000 title abstract 2
- 238000004891 communication Methods 0.000 claims abstract description 28
- 239000013598 vector Substances 0.000 claims abstract description 27
- 238000003491 array Methods 0.000 claims description 12
- 238000010845 search algorithm Methods 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 abstract description 26
- 238000000151 deposition Methods 0.000 abstract 1
- 230000014509 gene expression Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 6
- 238000012549 training Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 241000287196 Asthenes Species 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
- H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
- H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
- H04B7/0617—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Transmission System (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提供了一种联合方法的层次码本结构设计方法,用于无线通信技术。本发明设计的码本满足两个设计原则:同一层中的全体码字需要能够覆盖整个空间;码本的任意一层中任意一个码字的空间覆盖范围必须是下一层中数个码字的集合;同时还满足第k层的波瓣宽度为2/2k,并且最后一层的码本将采用信道指向向量的表达形式。在生成每层的码字时,先生成每层的第一个码字,进而旋转该码字,将其指向不同方向,就能够构建出该层的全部码字。采用本发明设计的层次码本进行搜索,搜索方法的收敛速率更快,搜索成功率更高,波束成型的增益更高,在任意阶段接收功率较大,发射功率利用率较高,在存在多径的情况下鲁棒性更好。
Description
技术领域
本发明属于无线通信技术领域,具体是指一种应用于子阵式和开关天线式联合方法的层次码本结构设计方法。
背景技术
毫米通信技术是具有非凡前景的下一代无线通信技术。因为毫米波段具有极其丰富的频谱资源,这使得毫米波通信技术可以提供比当前的通信技术高的多的信道容量。所以近年来,对于毫米波无线局域网络和毫米波移动通信网络的研究正如火如荼的展开。然而毫米波通信技术所面临的巨大技术难题是在毫米波段信号的传输衰减非常大。为了克服传输衰减,一般需要采用天线阵列的阵列增益加以弥补。由于毫米波波长很短,所以可以在单一设备中集中大量天线单元,这使得毫米波通信成为可能。
在毫米波领域,射频器件的极高能量开销和射频元件的高昂制造成本,使得几乎不可能在毫米波领域采用传统多输入多输出(Multiple Input Multiple Output,MIMO)中的数字波束成形技术。而采用单数据流的模拟波束成形技术,更为适用于毫米波波束成形。这种结构中,所有的天线都链接到同一个射频链路,每个天线幅度增益相同,相位增益可调。由于在天线数目非常多的毫米波通信领域,信道信息(Channel State Information,CSI)的估计极为耗时并且非常困难,所以业界一般采用以下两种方式来实现毫米波模拟波束成形,即迭代波束成形方法和波束切换的方法。前者是在训练过程中,固定一方的波束,然后另一端不断进行迭代,从而找最优波束。后者则是从预定义的码本中查找出最优的码字进行通信。
对于波束切换方法,穷举搜索是一种比较便利的搜索方法,只需要遍历整个码本查找出最优的码字即可,但是这种方法耗费的搜索时间非常长,特别是当天线阵列很大的时候。为了提高码本搜索效率,可以采用层次化码本搜索方法。将码本依据码字的覆盖范围的不同,划分为不同的层次,上一层的码字的覆盖范围是下一层数个相邻码字覆盖范围的集合。通过构建这样的码本可以大大提高码本搜索效率。公开号为104935366的中国专利申请在2015年9月23日公开了一种毫米波通信中的波束搜索方法,利用开关天线的方式设计了一种层次化码本设计结构,大大提高了码本搜索效率。Samsung公司的SridharRajasopal设计了一种基于子阵式的层次化码本设计(参考文献:Rajagopal S.Beambroadening for phased antenna arrays using multi-beam subarrays[C]//Communications(ICC),2012IEEE International Conference on.IEEE,2012:3637-3642.),有效的提高了毫米波通信的容量。但这种方法的设计实质上组成了四叉树,在码本搜索的时间复杂度太高,并且,该文献仅仅提出了码字设计方法并没有讨论码本搜索方法。
发明内容
本发明为了解决码本搜索的效率问题,以使得总发射功率更小、功率利用率更高、抗干扰性能更好,提供了一种基于联合方法的层次码本结构设计方法。
本发明提供的联合方法的层次码本结构设计方法,实现过程为:设发送端或接收端的天线单元为N个,码本被划分为log2(N)+1层,层序号k=0,1,...,log2(N),则:
(1.1)当k=log2(N)时,k层中的第n个码字其中,n=1,2,...,N;a(·)表示信道指向向量;
(1.2)当k=log2(N)-l时,l=1,2,...,log2(N),k层中的第n个码字w(k,n)根据下面过程产生:
(1.2.1)将w(k,1)划分为M个子阵,其中第m个子阵的天线权重向量fm为:
其中,NA表示开启的子阵数目;
(1.2.2)当l为偶数时,设置NA=M/2,当l为奇数时,设置NA=M;
(1.2.3)根据第一个天线单元w(k,1),获得k层的其他天线单元w(k,n)为:
其中,“ο”代表Hadamard算子;
(1.2.4)将所***字归一化。
基于本发明提供的层次码本结构,采用了二叉树的查找算法搜索最优的码字,实现本发明的联合波束搜索方法。在每一层的搜索中,比较当前父码字所覆盖的两个子码字所对应的接收信号的信噪比SNR,将SNR较大的子码字设为下一层搜索的父码字,然后进入下一层搜索。重复搜索过程,直到搜索到最后一层,在最后一层搜索得出的码字就是整个码本的最优码字,将最优码字用于数据通信。
本发明的优点与积极效果在于:
(1)采用本发明设计的层次码本进行联合搜索,在每次搜索过程中,所使用的天线单元数目要远远大于DEACT方法,从而本发明方法的信号发射和接收能量也要远大于DEACT方法,也就是说本发明的搜索方法能够比DEACT方法提供更大的增益。仿真实验表明,本发明提出的码本相比于子阵式和开关天线式码本具有显著的优点。
(2)采用本发明设计的层次码本进行搜索,搜索方法的收敛速率更快,搜索成功率更高;波束成型的增益更高,在任意阶段接收功率较大,发射功率利用率较高,在存在多径的情况下鲁棒性更好,而且***设计相对简单。
附图说明
图1是毫米波通信中的收发端均匀天线阵列模型示意图;
图2是层次码本覆盖关系示意图;
图3是DEACT方法的波束示意图;
图4是本发明的所生成的层次码字波束示意图;
图5是2-D平面天线阵***模型;
图6是整体传输模型中在LOS信道模型和NLOS信道模型情况下,本发明波束搜索方法和DEACT方法接收能量的对比示意图;
图7是单天线传输模型中在LOS信道模型和NLOS信道模型情况下,本发明波束搜索方法和DEACT方法接收能量的对比示意图;
图8是整体传输模型中本发明波束搜索方法、DEACT和Sparse方法在LOS信道模型下的搜索成功率对比示意图;
图9是整体传输模型中本发明波束搜索方法、DEACT和Sparse方法在NLOS信道模型下的搜索成功率对比示意图;
图10是单天线传输模型中本发明波束搜索方法、DEACT和Sparse方法在LOS信道模型下的搜索成功率对比示意图;
图11是单天线传输模型中本发明波束搜索方法、DEACT和Sparse方法在NLOS信道模型下的搜索成功率对比示意图。
具体实施方式
下面将结合附图和实施例对本发明的技术方案作进一步的详细说明。
不失一般性,本发明设毫米波通信模型中的收发端均匀天线阵列(UniformlyLinear Antenna,ULA)中的天线单元之间间距均为半波长,收发端天线数目分别为NR和NT。如图1所示,仅有一条射频链路(RF Chain)连接到ULA上。在发送端,每个天线单元都和移相器和放大器(PA)连接,而在接收端,每个天线和低噪声放大器(Low Noise Amplifier,LNA)连接。此外,需要注意的是,每个天线单元都可以被单独开启或者关闭,并且每个天线单元都具有相同的放大倍数。
设s为发送的信号,接收信号y表达式为:
其中,Ptot表示所有使用的天线的总发送功率。wT和wR分别表示发送和接收的天线权重向量(Antenna Weight Vector,AWV)。H表示信道矩阵,n表示高斯白噪声,噪声能量为N0。上角标H表示共轭转置。
W(N)表示N个天线权重向量集合,其表达式如下所示。
其中,υ是归一化因子,从而可以使每个向量的能量都归一。可以看到,W(N)中的每一个AWV的增益要么是υ,要么是0。βi表示第i个天线单元的幅度增益,θi表示第i个天线单元的相位。上角标j表示虚符号。此外定义wT∈W(NT)和wR∈W(NR)。
定义发送信噪比γtot=Ptot/N0,接收信噪比ηtot可以表示为:
此模型下的能量增益,即阵列增益Gtot为:
在另一方面,毫米波通信中的增益放大器的能力是有限的。单天线传输模型在现实中也是非常有必要存在的。也就是接收信号可表示为:
其中,Pper是每个天线的传输能量,NTact是wT中参与信号传输的天线数目。
单天线传输模型中的每根天线的发送信噪比SNR为:γper=Pper/N0,接收信噪比ηper为:
此模型的能量增益为:
需要注意的是,整体传输模型和单天线传输模型分别用于描述信号放大器(PowerAmplifier,PA)具有足够放大能力和放大能力有限两种情况。鲁棒性较好的码本设计并不需要针对两种不同情况进行优化。
毫米波信道具有高反射、低散射的特点,其信道中的多径效应主要是由于反射产生的,即毫米波具有方向性。不同点传输路径具有不同的发射角度(Angle of Departure,AoD)和到达角度(Angle of Arrival,AoA)。依据前述的***模型,毫米波信道矩阵可以表示为:
其中,λl表示第l条路径的信道系数,L表示多径(Multiple PropagationComponent,MPC)的总数。a(·)表示信道指向向量,Ωl和ψl分别表示第l条路径AoD和AoA的余弦,它们的取值范围为[-1,1]。为了方便,下面将Ωl和ψl分别称为AoD和AoA。信道参数λl为瑞利分布。a(·)的表达式为:
其中,N表示接收端或者发送端的天线阵列单元数目,Ω表示AoD或者AoA。从a(·)的表达式可以看出a(·)是周期函数,a(N,Ω)=a(N,Ω+2)。信道矩阵同样需要归一化,表达式为:
进行波束搜索可以用数学表达式表述为:
可以清楚的看出,如果收发端能够知道H,并且收发天线的增益不是固定而是可变的,那么可以利用矩阵的奇异值(Singular Value Decomposition,SVD)分解来求得最优的AWV。然而在毫米波通信当中,对每个天线面对的信道都进行信道估计极其耗费时间,此外天线阵列一般是固定的幅度,增益是恒定的。所以SVD方法几乎不可行。
依据公式(8),毫米波通信信道主要不确定因素在于AoD和AoA,所以在本发明讨论的信道模型中,问题可以简化为在多径中寻找一条最想的路径并将AWV指向这个路径方向。
从这个角度出发,最简单的解决问题的方法是将AoD和AoA在[-1,1]依照1/N步长进行采样,并利用得到的AWV测试出最大的MPC方向。可见这种方法的码本就是由信道指向向量组成的。尽管穷举搜索总能搜索到最优路径,但是随之而来的代价是搜索时间的增加。穷举搜索的时间复杂度为O(N2),这对于毫米波通信***而言太高的,所以目前分层搜索算法广为采用。不同的分层搜索算法性能主要依靠码本WT和WR的设计,这也是本发明关注的重点。
在介绍两种原则之前,首先引入两种定义。定义A(w,Ω)为w沿着角度为Ω的信道指向向量方向的增益,表示如下:
其中,w表示发送或接收的天线权重向量,N表示w的非零单元数目,[w]n是向量的第n个单元。
CV(w)定义为w的覆盖范围,数学定义为:
其中,ρ为定义w的覆盖范围的参数,取值范围为(0,1)。随着ρ的变大,w的覆盖范围就会随之变小。当时,w的覆盖宽度刚好是3dB带宽覆盖范围。不同的码本设计方法所设计的波束覆盖范围是不同的,也就是ρ是不同的。即使在同一码本中不同的码字的波束覆盖范围也是不同的。A(w,ω)表示向量w沿着角度为ω的信道指向向量方向增益,ω是角度的余弦值,表示求取增益最大值。
层次搜索简单来说就是分层搜索,码本中的各个码字依据各自波束覆盖范围的不同划分成不同的层级。设w(k,n)为码本中第k层码本中的第n个码字。
层次码本的两个设计原则为:
1.同一层中的全体码字需要能够覆盖整个空间,也即:其中k=0,1,...,K-1;其中Nk表示序号为k的一层码本中的码字数目,K表示码本的层数。
2.码本的任意一层中任意一个码字的空间覆盖范围必须是下一层中数个码字的集合。也就是:其中Ιk,n是k+1层码本中码字的序号集合。为了方便起见,称w(k,n)为父码字,相应的{w(k+1,m)|m∈Ιk,n}为子码字。
显而易见,原则一保证了码本的全覆盖特性,没有任何角度是此码本无法覆盖到的。原则二则在不同的码字之间构建了一种树形结构关系,这样就可以采用结构化的搜索方法。如果任何一个父码字包含了M个子码字,那么就可以利用天线覆盖关系设计M叉树搜索算法。
在搜索时,首先,发送端采用全向天线发送信号,接收端采用M叉树搜索方法搜索最优码字,这个过程需要logM(NR)次搜索。然后接收端采用最优码字向发送端发送信号,发送端采用相同的码本搜索算法,经过logM(NT)次搜索找到最优码字。在每次搜索中,都需要采用穷举搜索从M个子码字中挑选出最优的码字作为下次搜索的父码字。这样总的层次搜索算法的搜索次数Ts为:
Ts=M logM(NT)+M logM(NR) (14)
因为现实中的天线阵列中天线单元数目一般都是2的指数,因此本发明方法基于二叉树的搜索方法。本发明方法的基础是开关天线式(DEACT)方法,因此首先介绍基于二分查找树的DEACT方法,其层次码本的波束覆盖关系如图2所示。此处以参考文件(He T,XiaoZ.Suboptimal Beam Search Algorithm and Codebook Design for Millimeter-WaveCommunications[J].Mobile Networks and Applications,2015,20(1):86-97)的记载来说明DEACT方法。
DEACT方法中的码本被划分为log2(N)+1层,序号分别为0,1,...,log2(N)。而序号为k的码本层的码字总数为Nk=2k。这里的N对应接收端的天线单元数目NR或发送端的天线单元数目NT。那么有:
定义
也就是信道指向向量的波束宽度为2/N,指向方向为Ω。这就是天线数目为N时最窄的波束宽度了。这种码本的码字的参数ρ可以表达为:
当N很大时ρ≈0.64。
由于最后一层码本中有N个码字,并且这些码字的集合覆盖了[-1,1]的全部角度。所以每个码字的波瓣宽度应该是2/N,并且分别指向N个不同的角度方向。所以,最后一层码本的码字覆盖范围可以表示为:
同理,对于其他层的码本的覆盖范围可以表示为:
与前述对比可以得到w(k,n)的表达式为:
采用DEACT方法,k层中激活的天线数目为2k,其余的天线都被关闭。图3给出了天线单元数目为64时的波束覆盖方式。
对于DEACT方法,当码本层序号k比较小时,激活的天线数目非常少,这就限制了总的发送能量。因为在毫米波通信中每个天线的传输能量是有限的,如果激活的天线数目很小,那么总的传输能量就会变小。为了让总的发送能量尽可能的变大,本发明结合子阵式和开关天线式的方法,提出了一种联合方法的层次码本结构设计方法。
此外,希望能够设计的码本能够满足前述的两种码本设计准则,同时还需要满足条件为,k层的波瓣宽度为2/2k,并且最后一层的码本将采用信道指向向量的表达形式。由于每层码本中的码字具有同样的波束宽度,但是指向方向不同,所以,可以先设计每一层的第一个码字,进而旋转该码字,将其指向不同方向,就能够构建出该层的全部码字了。通过展宽该层码字的波瓣宽度并重复以上过程,就可以得到下一层的码字集合。
为了实现这种设计思路,首先引入波束转动公式及其推论。
波束转动的公式一:
其中ο代表Hadamard算子,N是w中天线单元的数目,ψ是任意角度,A+ψ是集合A中所有元素加了ψ之后得到的新的集合。
公式一说明,如果得到了第k层中指向某一方向的波瓣宽度为2/2k的天线权重向量时,可以利用公式一来得到指向其他方向的码字。
波束转动的推论一:假设序号为k的码本层的第一个天线单元为w(k,1),那么该层的其他天线单元表达式为:
证明:
由于
从而有|[w(k,n)]i|=|[w(k,1)]i|,且w(k,1),w(k,n)∈W(N)。由于根据公式一可以得到:
本发明通过波束展宽和前述的两个码本设计准则,就能构建整个层次码本。更直观的说,如果将整个天线阵列划分成不同的子阵,这些子阵分别指向不同的方向,那么这些子阵波瓣的联合就能够构成更宽的波束。
设将天线单元数目N阵列天线划分为M个子阵,每个子阵的天线数目为Ns,其中第m个子阵的天线权重向量为其中第m个子阵的第n个天线单元的权值向量为将AWV的增益可以改写为下式:
这说明w的增益可以表示为所有子阵增益的集合。
如果子阵的AWV可以表示为其中是子阵AWV的天线系数,θm表示第m个子阵的相位。那么第m个子阵的波束覆盖范围CV(fm)如下:
从而w的覆盖范围可以表达为:
也即原本的波束宽度被扩展了M2倍。
然而,上面并没有考虑到不同子阵之间的相互相干作用。设子阵表达式为那么w的增益表达式可以进一步表示为:
上角标H表示共轭转置。由于信道指向向量具有相互独立的特性,任意两个子阵的天线权重相邻之间的增益不会发生互相干扰。
此外,为了减少波束的波动效应,需要保证,相邻子阵波束覆盖的交叉区域的天线增益必须足够大。这一点可以通过调整子阵天线系数来实现。为了确定合理的子阵天线系数,可以通过求解下面的无约束最优化问题来解决:
其中,Δθ表示相邻两个子阵的相位差,f(Ns,Δθ)表示增益。
求解过程如下:
其中(a)用到了信道指向向量间相互独立的特性,(b)假设天线阵列单元数目Ns是偶数。从而可以得到结果为 表示整数集合。所以选择来减少波束的波动效应。
综上,通过采用子阵结构,可以得到波束覆盖范围为的码字,如果结合子阵和DEACT方法,那么就可以得到波束覆盖范围为的码字,进而可实现二分法查找。其中,NA表示开启的子阵数目。
如前述的两个码本设计准测,设计的层次码本结构特点应该是,第k层的码字的波束覆盖宽度为2/2k。本发明用于联合方法的层次码本结构的具体实施方法为:设发送端或接收端的天线单元为N个,码本被划分为log2(N)+1层,层序号k=0,1,...,log2(N),则:
(1.1)当k=log2(N)时,该层中的第n个码字其中,n=1,2,...,N;
(1.2)当k=log2(N)-l时,l=1,2,...,log2(N),根据下面过程产生k层中的第n个码字w(k,n);
(1.2.1)将w(k,1)划分为个子阵,每个子阵表达式为
第m个子阵的天线权重向量fm如下:
其中,NA表示开启的子阵数目;
(1.2.2)按照公式(28)确定子阵表达式,其中当l为偶数时,开启的子阵数目NA=M/2,当l为奇数时,开启的子阵数据NA=M;
(1.2.3)利用波束转动的推论一,生成本层码本的其他码字;
(1.2.4)将所***字归一化。
图4给出了当天线单元数目N=128时,本发明设计的层次码字的波束图。从图中可以看出,w(0,1)的覆盖范围是w(1,1)和w(1,2)覆盖范围的联合。
基于本发明所实现的层次码本结构,采用二叉树的查找算法搜索最优的码字,实现联合波束搜索方法。在每一层的搜索中,比较当前父码字所覆盖的两个子码字所对应的接收信号信噪比SNR,将SNR较大的子码字设为下次搜索的父码字,然后进入下一次搜索。重复搜索步骤,直到搜索到最后一层。在最后一层搜索得出的码字就是整个码本的最优码字,将最优码字用于数据通信。接收端和发送端均采用该方法进行最优码字搜索以进行数据通信。
将本发明的联合搜索方法和DEACT方法进行比较,可以看出,本发明方法在接收信号能量上具有明显优势。如果每个天线单元的能量是有限的,那么每次搜索过程中,本发明方法所使用的天线单元数目要远远大于DEACT方法,从而本发明的联合搜索方法的信号发射和接收能量也要远大于DEACT方法,也就是说本发明的联合搜索方法能够比DEACT方法提供更大的增益。
此外,在波速搜索过程结束后,收发双发都将从最后一层码本中找出各自最优的码字进行通信。由于最后一层码本的波束覆盖范围为2/N。这样的精度或许在某些情况下不能满足实际应用需求,为此,还可以专门设计一个波束覆盖范围小于2/N的码本,为更为苛刻的通信条件提供支持。
进一步地,本发明实现针对平面天线阵的层次码本结构设计及联合搜索。将本发明适用于均匀线天线(Uniformly Linear Antenna,ULAs)的层次码本结构设计及联合搜索方法扩展到均匀面天线(Uniformly Planar Antenna,UPAs),从而可实现将一维的层次码本结构设计及联合搜索方法扩展为二维的联合波束搜索方法。UPAs***模型如图5所示,窄带平坦衰落的UPAs的信道模型为:
其中at和ar分别是收发两端对应的信道指向向量,具体表达式为:
其中,λ为载波波长,d为UPAs中相邻天线单元的距离。Nx和Ny分别是UPAs沿X轴和Y轴方向的天线单元数目。到达UPAs的信号相对于X轴和Y轴的角度分别为θ,它们的余弦分别为Φ和Θ。ax和ay分别为沿X轴方向和Y轴方向的信道指向向量。
通过将UPAs信道指向向量分解为沿X轴和Y轴的ax和ay,可以将UPAs的码本及其搜索过程分解为X轴向和Y轴向的两个独立的ULA的码本和搜索过程,在X轴向和Y轴向两个方向分别采用一维的联合搜索方法,就得到了适用于UPAs的码本和相应的搜索方法,得到X轴向和Y轴向的码字表达式为:
其中,wx(kx,ix)和wy(ky,iy)分别是X轴向的位于kx层的第ix个码字和Y轴向的位于ky层的第iy个码字。wx(kx,ix)和wy(ky,iy)均采用本发明所述的层次码本结构设计方法生成。kx、ky分别为X轴向、Y轴向的码本的层序号。ix、iy为正整数。
搜索最优的二维层次码本的过程是:首先保持X轴向的码字表达式为wx(1,1),在Y轴向上采用二叉树的查找算法搜索最优的码字,设查找到Y轴向的最优码字然后保持Y轴向的码字表达式为不变,在X轴向上同样采用二叉树的查找算法搜索最优的码字,查找X轴向的最优码字最终可以得到UPAs***的最优码字如下:
图6是整体传输模型(指UPAs***模型)中在LOS信道模型和NLOS信道模型情况下,本发明波束搜索方法和DEACT方法接收能量的对比。将本发明波速搜索方法简记为BMW-SS。其中收发端天线单元数目NT=NR=64,信道模型中多径数目L=3,发射总功率为Ptot=1W,噪声功率为N0=10-4W。其中1到6步为发射端训练,7到12步为接收端训练。图6和图7中,横坐标表示搜索步数,纵坐标表示信号接收功率。
图7是单天线传输模型中在LOS信道模型和NLOS信道模型情况下,BMW-SS和DEACT方法接收能量的对比。其中收发端天线单元数目NT=NR=64,信道模型中多径数目L=3,发射总功率为Pper=1W,噪声功率为N0=10-4W。其中1到6步为发射端训练,7到12步为接收端训练。
图8是整体传输模型中BMW-SS,DEACT和Sparse方法在LOS信道模型下的搜索成功率。其中收发端天线单元数目NT=NR=64,信道模型中多径数目L=3,发射总功率为Ptot=1W,噪声功率为N0=10-4W。η为LOS径和NLOS径信号功率的比。图8~图11中,横坐标表示接收信号的信噪比,纵坐标表示成功率。
Sparse方法为参考文件(Alkhateeb,A.,Ayach,O.El,Leus,G.,&Jr,R.W.H.(2014).Channel Estimation and Hybrid Precoding for Millimeter Wave CellularSystems.IEEE Journal of Selected Topics in Signal Processing,4553(c),1–1.)记载的波束成形方法。
图9是整体传输模型中BMW-SS,DEACT和Sparse方法在NLOS信道模型下的搜索成功率。其中收发端天线单元数目NT=NR=64,信道模型中多径数目L=3,发射总功率为Ptot=1W,噪声功率为N0=10-4W。η为LOS径和NLOS径信号功率的比。
图10是单天线传输模型中BMW-SS,DEACT和Sparse方法在LOS信道模型下的搜索成功率。其中收发端天线单元数目NT=NR=64,信道模型中多径数目L=3,发射总功率为Pper=1W,噪声功率为N0=10-4W。η为LOS径和NLOS径信号功率的比。
图11是单天线传输模型中BMW-SS,DEACT和Sparse方法在NLOS信道模型下的搜索成功率。其中收发端天线单元数目NT=NR=64,信道模型中多径数目L=3,发射总功率为Pper=1W,噪声功率为N0=10-4W。η为LOS径和NLOS径信号功率的比。
通过实验性能测试表明,在天线阵总能量约束和单个天线能量约束条件下,本发明波束搜索方法的性能都要比DEACT方法更具有优势。具体体现在本发明方法具有更平坦的波束和更多的参与通信的天线数目。这两个优点使得本发明方法的搜索成功率和接收SNR要比其他方法、尤其是DEACT方法有显著提升,并且搜索的收敛速率更快。
Claims (3)
1.一种联合方法的层次码本结构设计方法,其特征在于,实现步骤如下:设发送端或接收端的天线单元为N个,码本被划分为log2(N)+1层,层序号k=0,1,...,log2(N),依据层次码本的两个设计原则:
(1)同一层中的全体码字需要能够覆盖整个空间,也即:其中k=0,1,...,K-1;其中Nk表示序号为k的一层码本中的码字数目,K表示码本的层数;
(2)码本的任意一层中任意一个码字的空间覆盖范围必须是下一层中数个码字的集合,也就是:其中Ik,n是k+1层码本中码字的序号集合;称w(k,n)为父码字,相应的{w(k+1,m)|m∈Ik,n}为子码字;
则:
(1.1)当k=log2(N)时,k层中的第n个码字其中,n=1,2,...,N;a(·)表示信道指向向量;
(1.2)当k=log2(N)-l时,l=1,2,...,log2(N),k层中的第n个码字w(k,n)根据下面过程产生:
(1.2.1)将w(k,1)划分为M个子阵,其中第m个子阵的天线权重向量fm为:
其中,NA表示开启的子阵数目;Ns表示天线阵列单元数目;
(1.2.2)当l为偶数时,设置NA=M/2,当l为奇数时,设置NA=M;
(1.2.3)根据第一个天线单元w(k,1),获得k层的其他天线单元w(k,n)为:
其中,代表Hadamard算子;
(1.2.4)将所***字归一化。
2.根据权利要求1所述的联合方法的层次码本结构设计方法,其特征在于,采用了二叉树的查找算法搜索最优的码字,实现联合波束搜索;在每一层的搜索中,比较当前父码字所覆盖的两个子码字所对应的接收信号的信噪比SNR,将SNR较大的子码字设为下一层搜索的父码字,然后进入下一层搜索;重复搜索过程,直到搜索到最后一层,在最后一层搜索得出的码字就是整个码本的最优码字,将最优码字用于数据通信。
3.根据权利要求1所述的联合方法的层次码本结构设计方法,其特征在于,应用于均匀面天线时,将信道指向向量分解为沿X轴和Y轴的分量ax和ay,在X轴向和Y轴向上分别采用所述的层次码本结构设计方法生成码字,设X轴向的位于kx层的第ix个码字为wx(kx,ix),Y轴向的位于ky层的第iy个码字为wy(ky,iy),kx、ky分别为X轴向、Y轴向的码本的层序号,则得到二维的码字如下:
在搜索最优的码字时,首先保持X轴向的码字wx(1,1)不变,采用了二叉树的查找算法搜索在Y轴向上的最优的码字,设查找到Y轴向的最优码字然后保持Y轴向的码字不变,在X轴向上采用了二叉树的查找算法搜索最优的码字,设查找X轴向的最优码字最终得到最优的码字为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610251727.6A CN105959044B (zh) | 2016-04-21 | 2016-04-21 | 联合方法的层次码本结构设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610251727.6A CN105959044B (zh) | 2016-04-21 | 2016-04-21 | 联合方法的层次码本结构设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105959044A CN105959044A (zh) | 2016-09-21 |
CN105959044B true CN105959044B (zh) | 2019-03-01 |
Family
ID=56917885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610251727.6A Active CN105959044B (zh) | 2016-04-21 | 2016-04-21 | 联合方法的层次码本结构设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105959044B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106355027B (zh) * | 2016-09-13 | 2020-06-02 | 佛山科学技术学院 | 一种分解大奇数的方法 |
CN108365876A (zh) * | 2017-01-26 | 2018-08-03 | ***通信有限公司研究院 | 一种数字波束成形码本生成方法及设备 |
CN107135023B (zh) * | 2017-04-28 | 2020-04-24 | 东南大学 | 用于毫米波通信***的三维训练码书设计方法及波束对准方法 |
CN107634787A (zh) * | 2017-08-22 | 2018-01-26 | 南京邮电大学 | 一种大规模mimo毫米波信道估计的方法 |
CN108462520B (zh) * | 2018-02-06 | 2021-05-14 | 重庆邮电大学 | 一种低分辨率APSs分层码本设计及搜索方法 |
CN109831236B (zh) * | 2018-11-13 | 2021-06-01 | 电子科技大学 | 一种基于蒙特卡洛树搜索辅助的波束选择方法 |
CN110518949B (zh) * | 2019-10-24 | 2020-01-31 | 江苏久鑫铜业有限公司 | 等功率约束下大规模天线训练码本的设计方法及训练方法 |
CN113825146B (zh) * | 2020-06-19 | 2023-01-06 | 华为技术有限公司 | 波束确定方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414864A (zh) * | 2008-12-08 | 2009-04-22 | 华为技术有限公司 | 多天线分层预编码的方法及装置 |
CN103067065A (zh) * | 2012-12-20 | 2013-04-24 | 清华大学 | 基于双层圆形天线阵的60GHz波束赋形方法和*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5551823B2 (ja) * | 2010-04-07 | 2014-07-16 | アルカテル−ルーセント | チャネル状態情報フィードバック |
-
2016
- 2016-04-21 CN CN201610251727.6A patent/CN105959044B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414864A (zh) * | 2008-12-08 | 2009-04-22 | 华为技术有限公司 | 多天线分层预编码的方法及装置 |
CN103067065A (zh) * | 2012-12-20 | 2013-04-24 | 清华大学 | 基于双层圆形天线阵的60GHz波束赋形方法和*** |
Non-Patent Citations (1)
Title |
---|
"毫米波通信中一种优化的波束搜索算法研究";刘凌;《自动化技术与应用》;20140529;全文 |
Also Published As
Publication number | Publication date |
---|---|
CN105959044A (zh) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105959044B (zh) | 联合方法的层次码本结构设计方法 | |
Perović et al. | Channel capacity optimization using reconfigurable intelligent surfaces in indoor mmWave environments | |
Cai et al. | Effect of wideband beam squint on codebook design in phased-array wireless systems | |
EP3577776B1 (en) | Transmitter and method for transmitting symbols over wireless communication channel | |
Alrabadi et al. | Beamforming via large and dense antenna arrays above a clutter | |
KR102228091B1 (ko) | 밀리미터파 기반 대규모 mimo 시스템의 하이브리드 빔포밍 장치 및 방법 | |
Priya et al. | Hybrid beamforming for massive MIMO using rectangular antenna array model in 5G wireless networks | |
Alluhaibi et al. | Capacity maximisation for hybrid digital-to-analog beamforming mm-wave systems | |
Li et al. | Energy-efficient Butler-matrix-based hybrid beamforming for multiuser mmWave MIMO system | |
Nguyen et al. | Beam squint effects in THz communications with UPA and ULA: Comparison and hybrid beamforming design | |
Abdallah et al. | Deep reinforcement learning based beamforming codebook design for RIS-aided mmWave systems | |
Yu et al. | Performance analysis and codebook design for mmWave beamforming system with beam squint | |
Zhang et al. | Cooperative multi-subarray beam training in millimeter wave communication systems | |
Torkzaban et al. | Codebook design for composite beamforming in next-generation mmwave systems | |
Chen et al. | Multi‐beam receive scheme for millimetre wave wireless communication system | |
Li et al. | Efficient analog beamforming with dynamic subarrays for mmWave MU-MISO systems | |
Xue et al. | An efficient beam-training scheme for the optimally designed subarray structure in mmWave LoS MIMO systems | |
Gao et al. | Wideband hybrid analog-digital beamforming massive MIMO systems based on Rotman lens | |
Almasi et al. | MmWave lens-based MIMO system for suppressing small-scale fading and shadowing | |
Mercy Sheeba et al. | Beamforming Techniques for Millimeter Wave Communications-A Survey | |
Zhang et al. | A closed-form PS-DFT codebook design for mmWave beam alignment | |
Lin et al. | Coordinated beamforming training for mmWave and sub-THz communications with antenna subarrays | |
Kareem et al. | Particle swarm optimization based beamforming in massive MIMO systems | |
Li et al. | Hierarchical Multi-Beam Search Based ChannelEstimation for Millimeter-Wave Massive MIMO Systems | |
Lin et al. | Antenna subarray partitioning with interference cancellation for multi-user indoor terahertz communications |
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 |