CN109379122B - 一种毫米波通信多径信道动态波束训练方法 - Google Patents

一种毫米波通信多径信道动态波束训练方法 Download PDF

Info

Publication number
CN109379122B
CN109379122B CN201811324848.4A CN201811324848A CN109379122B CN 109379122 B CN109379122 B CN 109379122B CN 201811324848 A CN201811324848 A CN 201811324848A CN 109379122 B CN109379122 B CN 109379122B
Authority
CN
China
Prior art keywords
codeword
codebook
channel
receiving
layer
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
Application number
CN201811324848.4A
Other languages
English (en)
Other versions
CN109379122A (zh
Inventor
戚晨皓
陈康建
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201811324848.4A priority Critical patent/CN109379122B/zh
Publication of CN109379122A publication Critical patent/CN109379122A/zh
Application granted granted Critical
Publication of CN109379122B publication Critical patent/CN109379122B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0456Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613Diversity 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/0615Diversity 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/0617Diversity 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
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0868Hybrid systems, i.e. switching and combining
    • H04B7/088Hybrid systems, i.e. switching and combining using beam selection

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明公开了一种毫米波通信多径信道动态波束训练方法,本发明充分研究了毫米波通信多径信道的波束训练方法,设计了动态可调且具有闭合表达式的发送端和接收端的分层码本,通过移除之前估计的路径干扰的方法,将多径信道估计的问题转换序贯的多个单径信道估计问题。其中,移除已估计径的影响通过设计码字实现,使得多径信道的波束训练复杂度大幅降低,准确度大幅提高。

Description

一种毫米波通信多径信道动态波束训练方法
技术领域
本发明属于毫米波无线通信领域,涉及一种毫米波通信多径信道动态波束训练的方法。
背景技术
随着移动终端设备的流行,人们对于无线通信的需求量也越来越大。为满足人们日益增长的对数据流量的需求,毫米波通信(30GHz-300GHz)因其丰富的频谱资源和极高的传输速率得到了人们的广泛关注。
由于比现有微波频段通信更高的载频,毫米波通信在空间传播时有更大的路径损耗。另一方面,更高的载频使得毫米波通信具有更小的天线尺寸。这使得在有限的区域内可以被封装更大规模天线阵列,并利用天线阵列的增益弥补信道路径损耗。现行的微波频段通信一般为每个天线分配专用的射频链路。然而,由于毫米波通信一般应用大规模天线阵列,为每根天线分配专用的射频链路将会产生高昂的射频链路成本。为节省射频链路的成本,一种使用少量射频链路的混合预编码结构在毫米波通信中被广泛采用。在混合预编码结构中,每个射频链路通过天线数目的移相器链接到所有的天线上。
为了获得毫米波通信信道的信息,一种基于预先设定的码本的波束训练方式被广泛采用。为了进一步加快波束训练的速度,一种基于分等级码本的分层波束训练方式被提出。在分等级码本中,每个码字均是覆盖一定空间范围的波束,且上层码字的波束覆盖范围是底层相应的两个码字波束覆盖范围的叠加。波束训练通常寻找测试的码字中具有最大能量的码字。运用分等级码本获得单径信道的信道状态信息通常较为容易且表现较好。然而,将分等级波束训练应用于多径信道时,如果想要获得多个径的信息,需要对这些路径进行序贯搜索。在搜索新的路径的时候,前面路径的信息需要被移除。
如何移除之前估计的路径的影响时毫米波通信系统多径信道波束训练的重点。虽然多径信道在文献[1]中被考虑([1]Z.Xiao,T.He,P.Xia,and X.G.Xia,“Hierarchicalcodebook design for beamforming training in millimeter-wave communication,”IEEE Trans.Wireless Commun.,vol.15,no.5,pp.3380–3392,May.2016.),然而其方法只适合于寻找一条路径。将文献[1]的方法应用于获取多条路径的信息,将会有极大的可能产生错误的估计。作为对文献[1]方法的提升,文献[2]通过增加额外的训练,获取更精确的信道信息。虽然,比起文献[1],文献[2]提高了获取多径信道信息的准确度,其性能依然较低([2]Z.Xiao,H.Dong,L.Bai,P.Xia,and X.Xia,“Enhanced channel estimation andcodebook design for millimeter-wave communication,”IEEE Trans.Veh.Technol.,vol.67,no.10,pp.9393-9405,Oct.2018)。
文献[1]和[2]通过计算出已估计路径的参数(信道的衰落系数、发送角及到达角),然后从接收的信号中减去已估计路径的影响。然而,通过码本难以得到多径的精确参数。这使得在估计当前径的时候残留大量已估计路径的影响,极大地影响了信道信息的获取。
发明内容
发明目的:针对以上问题,本发明提出一种有效的多径信道动态波束训练的方法。该方法分成动态码本设计和动态波束训练两部分。动态码本设计主要是设计具有闭合表达式的可调码本。动态波束训练是根据已估计的路径,为每条路径设计适合的码本,从而避免了对信道衰落系数的估计,提高了获取毫米波多径信道的准确性。
技术方案:为实现本发明的目的,本发明所采用的技术方案是:一种毫米波通信多径信道动态波束训练方法,该方法包括如下步骤:
(1)设计毫米波通信接收端的信号接收模型;
(2)设计毫米波通信中的传输信道模型;
(3)设计动态可调且具有闭合表达式的发送端和接收端的分层码本;
(4)搜索分层码本,估计多径信道的第一条路径;
(5)基于前Lf条已估计的路径,对第Lf+1条路径进行估计,Lf≥1;
(6)预设Ld条路径,重复步骤(5)计算,直到所有的Ld条路径均已被找到。
进一步,步骤(1)中的信号接收模型设计如下:
Figure BDA0001858436850000021
其中,fBB、FRF、H、WRF和wBB分别表示数字预编码向量、模拟预编码矩阵、信道矩阵、模拟合并矩阵和数字合并向量;y、P、x和η表示接收信号、发送功率、发送信号和加性高斯白噪声向量,(·)H表示做共轭转置。
进一步,步骤(2)中,设计毫米波通信中的传输信道模型:
Figure BDA0001858436850000022
其中,L、λl
Figure BDA0001858436850000024
Nt和Nr分别表示路径数目、第l条路径的信道增益、信道接收角、信道发送角、发送天线数目和接收天线数目;α(N,Ω)表示信道引导向量,定义为:
Figure BDA0001858436850000023
其中,N为天线数目,Ω为信道发送角或接收角。
进一步,步骤(3)中,设计动态可调且具有闭合表达式的发送端和接收端的分层码本,方法如下:
(3.1)设发送端的天线数为Nt,当Nt设为2的指数次方时,码本具有如下特征:
①发送端的码本共有S+1层,S=log2Nt
②码本的第s层共有2s个码字,s=0,1,2,…S;码本的最底层共有Nt个码字,每个码字均是一个信道引导向量,且第i个码字表示为fi=α(Nt,-1+(2i-1)/Nt),i=1,2,…,Nt
③第s层每个码字覆盖的宽度为2/2s,最底层码字覆盖的宽度为2/Nt
④每个码字覆盖的宽度均可表示为数个最底层码字覆盖的宽度的集合;
(3.2)待设计码字假设为v=Vt(s,m),即待设计的码字是分层码本Vt的第s层的第m个码字,m=1,2,…,2s,则用于设计v的引导向量索引集合可以表示为
Figure BDA0001858436850000031
Figure BDA0001858436850000032
未功率归一化的码字
Figure BDA0001858436850000033
设计为:
Figure BDA0001858436850000034
其中,θi表示第i个引导向量fi的加权相位;
(3.3)步骤(3.2)中的θi采用如下方法来设计:
Figure BDA0001858436850000035
(3.4)根据步骤(3.3)中获得的
Figure BDA0001858436850000036
将步骤(3.2)中的
Figure BDA0001858436850000037
写为:
Figure BDA0001858436850000038
对设计的码字进行归一化即可得到最终设计的码字:
Figure BDA0001858436850000039
其中,|| ||2表示向量的
Figure BDA00018584368500000318
范数。
(3.5)对于发送端的码本Vt中所有码字Vt(s,m)均按照步骤(3.2)-(3.4)设计,即可得到发送端码本Vt,s=0,1,2,…S,m=1,2,…2s
(3.6)对于接收端,给定天线数目Nr,依照步骤(3.1)到步骤(3.5),即可得到接收端码本Vr
进一步,步骤(4)中,估计第一条路径的方法如下:
(4.1)在估计第一条径时,发送端和接收端分别使用步骤(3)中设计的分层码本Vt和Vr,记为
Figure BDA00018584368500000310
Figure BDA00018584368500000311
其中上标“1”表示用于估计第一条路径的码本,其第s层第m个码字分别表示为
Figure BDA00018584368500000312
Figure BDA00018584368500000313
s=0,1,2,…,S,m=0,1,2,…,2s,步骤(3)中用于设计发端码字和收端码字的引导向量的索引的集合分别表示为
Figure BDA00018584368500000314
Figure BDA00018584368500000315
Figure BDA00018584368500000316
Figure BDA00018584368500000317
代入步骤(3.4)的公式中即可得到
Figure BDA0001858436850000041
Figure BDA0001858436850000042
(4.2)将步骤(1)中的信号接收模型写为:
Figure BDA0001858436850000043
其中,
Figure BDA0001858436850000044
表示发送的码字,
Figure BDA0001858436850000045
表示接收的码字,给定搜索的起始层数S0,其中1≤S0≤S,分层波束训练从分层码本的第S0层开始,发送端依次发送
Figure BDA0001858436850000046
接收端依次使用
Figure BDA0001858436850000047
接收;在发送端使用第p个码字发送、接收端使用第q个码字接收时,
Figure BDA0001858436850000048
Figure BDA0001858436850000049
接收到的信号可表示为:
Figure BDA00018584368500000410
第S0层最佳的一个收发码字索引组合
Figure BDA00018584368500000411
可通过下式得到:
Figure BDA00018584368500000412
其中,
Figure BDA00018584368500000413
表示|y(p,q)|取得最大值时的p和q的取值;
Figure BDA00018584368500000414
Figure BDA00018584368500000415
分别表示发送端和接收端在S0层找到的最佳码字对的索引,即发送端和接收端找到的最佳码字分别为
Figure BDA00018584368500000416
Figure BDA00018584368500000417
(4.3)基于步骤(4.2),假设第s-1层(s≥S0+1)找到的收发码字索引组合为(ps-1,qs-1),对应于发送端码字
Figure BDA00018584368500000418
接收端码字
Figure BDA00018584368500000419
在进行第s层搜索时,发送端分别发送
Figure BDA00018584368500000420
Figure BDA00018584368500000421
接收端分别使用
Figure BDA00018584368500000422
Figure BDA00018584368500000423
接收,在发送端使用第p个码字发送、接收端使用第q个码字接收,p=1,2,q=1,2,接收到的信号可表示为:
Figure BDA00018584368500000424
计算(c,d)如下:
Figure BDA00018584368500000425
其中,
Figure BDA00018584368500000426
表示|y(p,q)|取得最大值时的p和q的取值,找到的第s层的最佳收发码字索引组合(ps,qs)=(2(ps-1-1)+c,2(qs-1-1)+d),即发送端和接收端在s层找到的最佳码字分别为
Figure BDA00018584368500000427
Figure BDA00018584368500000428
(4.4)重复执行步骤(4.3),直到s=S,假设第S层的最佳收发码字索引组合为(p1,q1),其中上标“1”表示第1条路径,将
Figure BDA00018584368500000429
Figure BDA00018584368500000430
分别作为信道的第一条路径的发送引导向量和接收引导向量的估计,将p1和q1分别存入集合T和R中:
T=T∪p1,R=R∪q1
其中,计算第一条路径之前,T和R均初始化为空集。
进一步,步骤(5)中,基于前Lf条已估计的路径,对第Lf+1条路径进行估计,方法如下:
(5.1)发送端和接收端使用的分层码本分别表示为
Figure BDA0001858436850000051
Figure BDA0001858436850000052
其第s层第m个码字分别表示为
Figure BDA0001858436850000053
Figure BDA0001858436850000054
s=0,1,2,…,S;m=0,1,2,…,2s,对应的引导向量索引集合分别表示为
Figure BDA0001858436850000055
Figure BDA0001858436850000056
计算
Figure BDA0001858436850000057
Figure BDA0001858436850000058
其中,A\B表示从集合A中删除在集合B中的元素,将
Figure BDA0001858436850000059
Figure BDA00018584368500000510
代入步骤(3.4)的公式中计算
Figure BDA00018584368500000511
Figure BDA00018584368500000512
从而将Lf条已估计的信道路径信息分别从
Figure BDA00018584368500000513
Figure BDA00018584368500000514
中移除,最终可得到
Figure BDA00018584368500000515
Figure BDA00018584368500000516
(5.2)估计第Lf+1条径时,发送端和接收端分别使用步骤(5.1)中获得的码本
Figure BDA00018584368500000517
Figure BDA00018584368500000518
执行步骤(4.2)至(4.4)的波束训练得到第Lf+1条径的发送角和接收角信息分别表示为
Figure BDA00018584368500000519
Figure BDA00018584368500000520
Figure BDA00018584368500000521
Figure BDA00018584368500000522
存入信道发送角和接收角信息存储集合T和R中:
Figure BDA00018584368500000523
有益效果:与现有技术相比,本发明的技术方案具有以下有益技术效果:
(1)本发明设计的动态分等级码本具有闭合表达式,可以快速地根据实际需求设计出合适的分等级码本。
(2)本发明应用动态的分等级码本进行动态波束训练,避免了复杂的计算,通过闭合表达式,可以快速地设计用于多径信道波束训练的码本。
(3)本发明设计的多径信道波束训练方案,通过码本更新移除已估计路径的影响,避免了对信道衰落系数的估计,在更少训练开销的前提下,估计的精度远高于现有方案。
附图说明
图1是本发明实施例使用的毫米波通信系统模型的示意图;
图2是本发明实施例使用的毫米波通信信道模型的示意图;
图3是本发明实施例毫米波通信波束训练使用的分等级码本的示意图;
图4是本发明实施例设计的动态分层码本中码字的示意图;
图5是本发明与文献[1]和[2]的方法获取信道信息准确度的比较。
具体实施方式
下面结合附图和实施例对本发明的技术方案做进一步的详细描述。
(1)如图1所示,本发明使用的毫米波通信系统模型描述如下:
信号在发送端经过数字预编码、射频链路、模拟预编码之后,通过天线阵列发送出去。信号在无线信道中传输之后到达接收端,接收端天线阵列收到的信号经过模拟合并、射频链路和数字合并之后得到最终接收信号。发送端和接收端的天线阵列均为间隔为半波长的均匀线性阵,分别有Nt和Nr根天线。则接收到的信号可以表示为:
Figure BDA0001858436850000061
其中,fBB、FRF、H、WRF和wBB分别表示数字预编码向量、模拟预编码矩阵、信道矩阵、模拟合并矩阵和数字合并向量;y、P、x和η表示接收信号、发送功率、发送信号和加性高斯白噪声向量,(·)H表示做共轭转置。
(2)如图2所示,本发明毫米波通信系统模型中的信道模型描述如下:
设发送端到传输端共有L条传输路径,每条传输路径用发送角,接收角和信道增益来表示。根据广泛应用的Saleh-Valenzuela(S-V)模型,这个端对端毫米波通信系统的信道一般被建模为:
Figure BDA0001858436850000062
其中,L、λl
Figure BDA0001858436850000063
Figure BDA0001858436850000064
分别表示路径数目、第l条路径的信道增益、信道的到达角(angle-of-arrival,AOA)和信道发送角(angle-of-departure,AOD)。α(N,Ω)表示信道引导向量,定义为:
Figure BDA0001858436850000065
其中,N为天线数目,Ω为信道AOA或者AOD。事实上,假设第l条路径的空间的发送角和到达角分别为
Figure BDA0001858436850000066
Figure BDA0001858436850000067
Figure BDA0001858436850000068
因此,我们可以得到
Figure BDA0001858436850000069
(3)本发明提出的用于估计步骤(2)中信道路径发送角和到达角的分层码本设计方案如下所述:
(3.1)分层码本的示意图如图3所示。其中发送天线数目Nt为16。其中Vt(s,m)表示分等级码本第s层的第m个码字,s=0,1,…,4,m=1,2,…,2s
(3.2)考虑步骤(3.1)所述码本的通用情况,设发送端的天线数为Nt,当Nt设为2的指数次方时,码本具有如下特征:
①发送端的码本共有S+1层,S由天线数目决定,S=log2Nt
②码本的第s层共有2s个码字,s=0,1,2,…S。
特别地,码本的最底层共有Nt个码字,每个码字均是一个信道引导向量,且第i个码字表示为fi=α(Nt,-1+(2i-1)/Nt),i=1,2,…,Nt
③码本中每个码字均覆盖一定的范围,第s层每个码字覆盖的宽度为2/2s。特别地,最底层码字覆盖的宽度为2/Nt
④每个上层码字覆盖的宽度均可表示为数个底层码字覆盖的宽度的集合。例如,在图3中,Vt(3,2)的覆盖范围可以表示为f3和f4覆盖范围的叠加,Vt(2,1)的覆盖范围可以表示为f1,f2,f3和f4覆盖范围的叠加。
(3.3)根据步骤(3.2)的分层码本特征④,分层码本中的任一码字的覆盖范围均可表示为一个或数个最底层码字的覆盖范围的叠加。待设计码字假设为v=Vt(s,m),即待设计的码字是分层码本Vt第s层的第m个码字,m=1,2,…,2s,则用于设计v的fi的索引集合可以表示为
Figure BDA0001858436850000071
例如,Nt=16,s=2,m=1时,Ψ2,1={i|1≤i≤4,i=1,2,…,Nt}={1,2,3,4},这与步骤(3.2)的分层码本特征④的例子相一致。未功率归一化的码字
Figure BDA0001858436850000072
设计为:
Figure BDA0001858436850000073
其中,θi表示第i个引导向量fi的加权相位,用于避免波束叠加过程中波束覆盖范围内低波束增益的出现,i=1,2,…,Nt
(3.4)步骤(3.3)中的θi采用如下方法来设计:对于任意的码字
Figure BDA0001858436850000074
其沿着任意方向Ω(Ω∈[-1,1])的波束增益可写为
Figure BDA0001858436850000075
其中,[v]n表示向量v的第n个元素。为了避免波束覆盖范围内低波束增益的出现,θi可以通过
Figure BDA0001858436850000076
来设计。由此可以得到
Figure BDA0001858436850000077
再设置
Figure BDA0001858436850000078
即可得到
Figure BDA0001858436850000079
(3.5)根据步骤(3.4)中获得的
Figure BDA00018584368500000710
可以将步骤(3.3)
Figure BDA00018584368500000711
写为
Figure BDA00018584368500000712
之后,对设计的码字进行归一化即可得到最终设计的码字:
Figure BDA0001858436850000081
(3.6)对于发送端的码本Vt中所有码字Vt(s,m)均按照步骤(3.3)-(3.5)设计,即可得到发送端码本Vt,s=0,1,2,…S,m=1,2,…2s
(3.7)对于接收端,给定天线数目Nr,依照步骤(3.1)到步骤(3.6),即可得到接收端码本Vr
(4)基于步骤(3)提出的分层码本设计方案,下面介绍估计多径信道方法。通常,系统只需要估计信道L条路径中的Ld条,Ld由并行传输的数据流数目决定。当运用波束训练估计信道的Ld条路径时,它们被序贯地估计。首先使用分层码本估计多径信道估计第一条路径,即信道增益最大的一条路径,方法如下所述:
(4.1)在估计第一条径时,发送端和接收端分别使用步骤(3)中设计的分层码本Vt和Vr,分别将他们重命名为
Figure BDA0001858436850000082
Figure BDA0001858436850000083
其中上标“1”表示用于估计第一条路径的码本。其第s层第m个码字分别表示为
Figure BDA0001858436850000084
Figure BDA0001858436850000085
s=0,1,2,…,S,m=0,1,2,…,2s。步骤(3)中用于设计发端码字和收端码字的引导向量的索引集合分别表示为
Figure BDA0001858436850000086
Figure BDA0001858436850000087
显然,将
Figure BDA0001858436850000088
Figure BDA0001858436850000089
代入步骤(3.5)的公式中即可得到
Figure BDA00018584368500000810
Figure BDA00018584368500000811
(4.2)步骤(1)中的系统模型可以写为:
Figure BDA00018584368500000812
其中
Figure BDA00018584368500000813
表示发送的码字,
Figure BDA00018584368500000814
表示接收的码字。给定搜索的起始层数S0,其中S≥S0≥1,分层波束训练从分层码本的第S0层开始,发送端依次发送
Figure BDA00018584368500000815
接收端依次使用
Figure BDA00018584368500000816
接收。在发送端使用第p个码字发送,接收端使用第q个码字接收时,
Figure BDA00018584368500000817
Figure BDA00018584368500000818
接收到的信号可表示为:
Figure BDA00018584368500000819
第S0层的最佳的一个收发码字索引组合
Figure BDA00018584368500000820
可通过下式得到:
Figure BDA00018584368500000821
其中,
Figure BDA00018584368500000822
表示|y(p,q)|取得最大值时的p和q的取值;
Figure BDA00018584368500000823
Figure BDA00018584368500000824
分别表示发送端和接收端在S0层找到的最佳码字的索引。具体地,发送端和接收端找到的最佳码字分别为
Figure BDA0001858436850000091
Figure BDA0001858436850000092
(4.3)基于步骤(4.2),假设第s-1层(s≥S0+1)找到的最佳码字对的次序为(ps-1,qs-1),即发送端和接收端在s-1层找到的最佳码字分别为
Figure BDA0001858436850000093
Figure BDA0001858436850000094
在进行第s层搜索时,发送端分别发送
Figure BDA0001858436850000095
Figure BDA0001858436850000096
接收端分别使用
Figure BDA0001858436850000097
Figure BDA0001858436850000098
接收,在发送端使用第p个码字发送、接收端使用第q个码字接收时,p=1,2,q=1,2,接收到的信号可表示为:
Figure BDA0001858436850000099
计算(c,d)如下:
Figure BDA00018584368500000910
其中,
Figure BDA00018584368500000911
表示|y(p,q)|取得最大值时的p和q的取值,找到的第s层的最佳收发码字索引组合(ps,qs)=(2(ps-1-1)+c,2(qs-1-1)+d),即发送端和接收端在s层找到的最佳码字分别为
Figure BDA00018584368500000912
Figure BDA00018584368500000913
(4.4)重复执行步骤(4.3),直到s=S,假设第S层的最佳收发码字索引组合为(p1,q1),其中上标“1”表示第1条路径,将
Figure BDA00018584368500000914
Figure BDA00018584368500000915
分别作为信道的第一条路径的发送引导向量和接收引导向量的估计,将p1和q1分别存入集合T和R中:
T=T∪p1,R=R∪q1
其中,T和R均初始化为空集。
(5)假设Lf条路径已经被估计,Lf≥1,那么还有Ld-Lf路径需要被估计,此时的接收信号可以变形为:
Figure BDA00018584368500000916
上式中,
Figure BDA00018584368500000917
为待估计项,wHη为噪声项,
Figure BDA00018584368500000918
为已估计项。
(5.1)为了对当前的第Lf+1条路径进行估计,前Lf条路径的影响需要被移除。移除前Lf条路径的影响,等效于最小化已估计项
Figure BDA00018584368500000919
显然,在系统和信道给定的情况下,上式中的变量只有w和v。因此,可以通过设计w和v来移除已估计路径的影响。这即是设计一个码本对已估计的Lf条路径进行抑制。
(5.2)发送端和接收端使用的分层码本分别表示为
Figure BDA0001858436850000101
Figure BDA0001858436850000102
其第s层第m个码字分别表示为
Figure BDA0001858436850000103
Figure BDA0001858436850000104
s=0,1,2,…,S;m=0,1,2,…,2s,对应的引导向量索引集合分别表示为
Figure BDA0001858436850000105
Figure BDA0001858436850000106
计算
Figure BDA0001858436850000107
Figure BDA0001858436850000108
其中,A\B表示从集合A中删除在集合B中的元素,将
Figure BDA0001858436850000109
Figure BDA00018584368500001010
代入步骤(3.4)的公式中计算
Figure BDA00018584368500001011
Figure BDA00018584368500001012
从而将Lf条已估计的信道路径信息分别从
Figure BDA00018584368500001013
Figure BDA00018584368500001014
中移除,最终可得到
Figure BDA00018584368500001015
Figure BDA00018584368500001016
(5.3)估计第Lf+1条径时,发送端和接收端分别使用步骤(5.2)中获得的码本
Figure BDA00018584368500001017
Figure BDA00018584368500001018
执行步骤(4.2)至(4.4)的波束训练得到第Lf+1条径的发送角和接收角信息分别表示为
Figure BDA00018584368500001019
Figure BDA00018584368500001020
Figure BDA00018584368500001021
Figure BDA00018584368500001022
存入信道发送角和接收角信息存储集合T和R中:
Figure BDA00018584368500001023
(6)剩余的路径重复步骤(5)计算,直到所有的Ld条路径均已被找到,每一条路径的发送角和接收角对应的引导向量次序分别存储在T和R中。
下面结合仿真条件与结果对本发明做进一步的描述:
(1)图4是Nt=16时,本发明实施例设计的码字
Figure BDA00018584368500001024
Figure BDA00018584368500001025
的示意图。首先,通过步骤(3.3)可以得到
Figure BDA00018584368500001026
然后根据步骤(3.5),可以得到
Figure BDA00018584368500001027
如图中点划线所示。设第一条估计路径的发送角对应码本
Figure BDA00018584368500001028
底层的第12个引导向量。将12从
Figure BDA00018584368500001029
中删除,于是可以得到
Figure BDA00018584368500001030
同样地,然后根据步骤(3.5),可以得到
Figure BDA00018584368500001031
如图中实线所示。由图中可以看出,相比于
Figure BDA00018584368500001032
沿着Ω=7/16方向的波束增益为0,这就实现了对已估计径的抑制。
(2)图5是Nt=Nr=32,L=Ld=3,S0=2时本发明实施例的波束训练方法和文献[1],[2]方法获取信道信息准确度的比较。其中,获取信息的准确度用搜索成功率衡量。对于一次训练,如果训练找到的引导向量和信道的Ld条路径完全匹配,那么我们就说这次训练是成功的。由图中可以看出,本发明所提出的方法,远远优于文献[1]和[2]的方法。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步地详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种毫米波通信多径信道动态波束训练方法,其特征在于,该方法包括如下步骤:
(1)设计毫米波通信接收端的信号接收模型;
(2)设计毫米波通信中的传输信道模型;
(3)设计动态可调且具有闭合表达式的发送端和接收端的分层码本;
(4)搜索分层码本,估计多径信道的第一条路径;
(5)基于前Lf条已估计的路径,对第Lf+1条路径进行估计,Lf≥1;
(6)预设Ld条路径,重复步骤(5)计算,直到所有的Ld条路径均已被找到;
步骤(1)中的信号接收模型设计如下:
Figure FDA0002822572810000011
其中,fBB、FRF、H、WRF和wBB分别表示数字预编码向量、模拟预编码矩阵、信道矩阵、模拟合并矩阵和数字合并向量;y、P、x和η表示接收信号、发送功率、发送信号和加性高斯白噪声向量,(·)H表示做共轭转置;
步骤(2)中,设计毫米波通信中的传输信道模型:
Figure FDA0002822572810000012
其中,L、λl
Figure FDA0002822572810000013
Nt和Nr分别表示路径数目、第l条路径的信道增益、信道接收角、信道发送角、发送天线数目和接收天线数目;α(N,Ω)表示信道引导向量,定义为:
Figure FDA0002822572810000014
其中,N为天线数目,Ω为信道发送角或接收角;
步骤(3)中,设计动态可调且具有闭合表达式的发送端和接收端的分层码本,方法如下:
(3.1)设发送端的天线数为Nt,当Nt设为2的指数次方时,码本具有如下特征:
①发送端的码本共有S+1层,S=log2Nt
②码本的第s层共有2s个码字,s=0,1,2,…S;码本的最底层共有Nt个码字,每个码字均是一个信道引导向量,且第i个码字表示为fi=α(Nt,-1+(2i-1)/Nt),i=1,2,…,Nt
③第s层每个码字覆盖的宽度为2/2s,最底层码字覆盖的宽度为2/Nt
④每个码字覆盖的宽度均可表示为多个最底层码字覆盖的宽度的集合;
(3.2)待设计码字假设为v=Vt(s,m),即待设计的码字是分层码本Vt的第s层的第m个码字,m=1,2,…,2s,则用于设计v的引导向量索引集合可以表示为
Figure FDA0002822572810000021
Figure FDA0002822572810000022
未功率归一化的码字
Figure FDA0002822572810000023
设计为:
Figure FDA0002822572810000024
其中,θi表示第i个引导向量fi的加权相位;
(3.3)步骤(3.2)中的θi采用如下方法来设计:
Figure FDA0002822572810000025
(3.4)根据步骤(3.3)中获得的
Figure FDA0002822572810000026
将步骤(3.2)中的
Figure FDA0002822572810000027
写为:
Figure FDA0002822572810000028
对设计的码字进行归一化即可得到最终设计的码字:
Figure FDA0002822572810000029
其中,|| ||2表示向量的l2范数;
(3.5)对于发送端的码本Vt中所有码字Vt(s,m)均按照步骤(3.2)-(3.4)设计,即可得到发送端码本Vt,s=0,1,2,…S,m=1,2,…2s
(3.6)对于接收端,给定天线数目Nr,依照步骤(3.1)到步骤(3.5),即可得到接收端码本Vr
步骤(4)中,估计第一条路径的方法如下:
(4.1)在估计第一条径时,发送端和接收端分别使用步骤(3)中设计的分层码本Vt和Vr,记为
Figure FDA00028225728100000210
Figure FDA00028225728100000211
其中上标“1”表示用于估计第一条路径的码本,其第s层第m个码字分别表示为
Figure FDA00028225728100000212
Figure FDA00028225728100000213
s=0,1,2,…,S,m=0,1,2,…,2s,步骤(3)中用于设计发端码字和收端码字的引导向量的索引的集合分别表示为
Figure FDA00028225728100000214
Figure FDA00028225728100000215
Figure FDA00028225728100000216
Figure FDA00028225728100000217
代入步骤(3.4)的公式中即可得到
Figure FDA00028225728100000218
Figure FDA00028225728100000219
(4.2)将步骤(1)中的信号接收模型写为:
Figure FDA00028225728100000220
其中,
Figure FDA00028225728100000221
表示发送的码字,
Figure FDA00028225728100000222
表示接收的码字,给定初始层数S0,其中S≥S0≥1,分层波束训练从分层码本的第S0层开始,发送端依次发送
Figure FDA0002822572810000031
接收端依次使用
Figure FDA0002822572810000032
接收;在发送端使用第p个码字发送、接收端使用第q个码字接收时,
Figure FDA0002822572810000033
Figure FDA0002822572810000034
接收到的信号可表示为:
Figure FDA0002822572810000035
第S0层最佳的一个收发码字索引组合
Figure FDA0002822572810000036
可通过下式得到:
Figure FDA0002822572810000037
其中,
Figure FDA0002822572810000038
表示|y(p,q)|取得最大值时的p和q的取值;
Figure FDA0002822572810000039
Figure FDA00028225728100000310
分别表示发送端和接收端在S0层找到的最佳码字对的索引,即发送端和接收端找到的最佳码字分别为
Figure FDA00028225728100000311
Figure FDA00028225728100000312
(4.3)基于步骤(4.2)假设第s-1层(s≥S0+1)找到的收发码字索引组合为(ps-1,qs-1),对应于发送端码字
Figure FDA00028225728100000313
接收端码字
Figure FDA00028225728100000314
在进行第s层搜索时,发送端分别发送
Figure FDA00028225728100000315
Figure FDA00028225728100000316
接收端分别使用
Figure FDA00028225728100000317
Figure FDA00028225728100000318
接收,在发送端使用第p个码字发送、接收端使用第q个码字接收,p=1,2,q=1,2,接收到的信号可表示为:
Figure FDA00028225728100000319
计算(c,d)如下:
Figure FDA00028225728100000320
其中,
Figure FDA00028225728100000321
表示|y(p,q)|取得最大值时的p和q的取值,找到的第s层的最佳收发码字索引组合(ps,qs)=(2(ps-1-1)+c,2(qs-1-1)+d),即发送端和接收端在s层找到的最佳码字分别为
Figure FDA00028225728100000322
Figure FDA00028225728100000323
(4.4)重复执行步骤(4.3),直到s=S,假设第S层的最佳收发码字索引组合为(p1,q1),其中上标“1”表示第1条路径,将
Figure FDA00028225728100000324
Figure FDA00028225728100000325
分别作为信道的第一条路径的发送引导向量和接收引导向量的估计,将p1和q1分别存入集合T和R中:
T=T∪p1,R=R∪q1
其中,计算第一条路径之前,T和R均初始化为空集。
2.根据权利要求1所述的一种毫米波通信多径信道动态波束训练方法,其特征在于,步骤(5)中,基于前Lf条已估计的路径,对第Lf+1条路径进行估计,方法如下:
(5.1)发送端和接收端使用的分层码本分别表示为
Figure FDA0002822572810000041
Figure FDA0002822572810000042
其第s层第m个码字分别表示为
Figure FDA0002822572810000043
Figure FDA0002822572810000044
s=0,1,2,…,S;m=0,1,2,…,2s,对应的引导向量索引集合分别表示为
Figure FDA0002822572810000045
Figure FDA0002822572810000046
计算
Figure FDA0002822572810000047
Figure FDA0002822572810000048
其中,A\B表示从集合A中删除在集合B中的元素,将
Figure FDA0002822572810000049
Figure FDA00028225728100000410
代入步骤(3.4)的公式中计算
Figure FDA00028225728100000411
Figure FDA00028225728100000412
从而将Lf条已估计的信道路径信息分别从
Figure FDA00028225728100000413
Figure FDA00028225728100000414
中移除,最终可得到
Figure FDA00028225728100000415
Figure FDA00028225728100000416
(5.2)估计第Lf+1条径时,发送端和接收端分别使用步骤(5.1)中获得的码本
Figure FDA00028225728100000417
Figure FDA00028225728100000418
执行步骤(4.2)至(4.4)的波束训练得到第Lf+1条径的发送角和接收角信息分别表示为
Figure FDA00028225728100000419
Figure FDA00028225728100000420
Figure FDA00028225728100000421
Figure FDA00028225728100000422
存入信道发送角和接收角信息存储集合T和R中:
Figure FDA00028225728100000423
CN201811324848.4A 2018-11-08 2018-11-08 一种毫米波通信多径信道动态波束训练方法 Active CN109379122B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811324848.4A CN109379122B (zh) 2018-11-08 2018-11-08 一种毫米波通信多径信道动态波束训练方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811324848.4A CN109379122B (zh) 2018-11-08 2018-11-08 一种毫米波通信多径信道动态波束训练方法

Publications (2)

Publication Number Publication Date
CN109379122A CN109379122A (zh) 2019-02-22
CN109379122B true CN109379122B (zh) 2021-03-16

Family

ID=65383843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811324848.4A Active CN109379122B (zh) 2018-11-08 2018-11-08 一种毫米波通信多径信道动态波束训练方法

Country Status (1)

Country Link
CN (1) CN109379122B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110233649B (zh) * 2019-06-09 2020-09-01 杭州电子科技大学 一种基于毫米波mimo系统的动态子阵列设计方法
CN110417444B (zh) * 2019-07-08 2020-08-04 东南大学 一种基于深度学习的毫米波信道波束训练方法
CN110365388B (zh) * 2019-07-31 2022-07-29 东南大学 一种低复杂度毫米波多播波束成形方法
CN110401476B (zh) * 2019-08-05 2022-07-08 东南大学 一种基于码本的毫米波通信多用户并行波束训练方法
CN112073106B (zh) * 2020-08-14 2022-04-22 清华大学 毫米波波束预测方法及装置、电子设备、可读存储介质
CN112803977B (zh) * 2021-01-07 2021-12-10 电子科技大学 波束偏移效应下毫米波通信系统的混合预编码方法
CN115150227B (zh) * 2021-03-30 2023-10-13 复旦大学 一种点对点多天线通信系统中的角度估计方法
CN113300747B (zh) * 2021-05-28 2022-06-17 东南大学 一种智能反射表面辅助的毫米波系统中的波束训练方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106027133A (zh) * 2016-05-20 2016-10-12 北京邮电大学 一种多径信道下的分级波束搜索方法
CN107135023A (zh) * 2017-04-28 2017-09-05 东南大学 用于毫米波通信系统的三维训练码书设计方法及波束对准方法
CN107453795A (zh) * 2017-08-21 2017-12-08 东南大学 多用户毫米波通信系统的波束分配方法及其装置和系统
CN107733479A (zh) * 2016-08-09 2018-02-23 北京信威通信技术股份有限公司 一种波束训练的方法及装置
CN108023618A (zh) * 2016-11-01 2018-05-11 华为技术有限公司 基于mimo天线的波束训练方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10122434B2 (en) * 2016-06-30 2018-11-06 Intel Corporation Apparatus, system and method of hybrid beamforming training

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106027133A (zh) * 2016-05-20 2016-10-12 北京邮电大学 一种多径信道下的分级波束搜索方法
CN107733479A (zh) * 2016-08-09 2018-02-23 北京信威通信技术股份有限公司 一种波束训练的方法及装置
CN108023618A (zh) * 2016-11-01 2018-05-11 华为技术有限公司 基于mimo天线的波束训练方法及装置
CN107135023A (zh) * 2017-04-28 2017-09-05 东南大学 用于毫米波通信系统的三维训练码书设计方法及波束对准方法
CN107453795A (zh) * 2017-08-21 2017-12-08 东南大学 多用户毫米波通信系统的波束分配方法及其装置和系统

Also Published As

Publication number Publication date
CN109379122A (zh) 2019-02-22

Similar Documents

Publication Publication Date Title
CN109379122B (zh) 一种毫米波通信多径信道动态波束训练方法
CN110401476B (zh) 一种基于码本的毫米波通信多用户并行波束训练方法
CN110519189B (zh) 高度移动场景下基于压缩感知的毫米波信道估计方法
EP3348001B1 (en) Systems and methods for beam selection for hybrid beamforming
JP6685623B2 (ja) 圧縮センシングに基づく効率的なスパースチャネル推定
CN102118825B (zh) 实现多点联合传输的方法、终端及系统
CN103067062B (zh) 多小区系统中基于干扰对齐的基站天线选择方法
CN108650200B (zh) 高低频混合组网系统的低频辅助信道估计方法
CN107566305B (zh) 一种低复杂度的毫米波系统信道估计方法
CN111917508A (zh) 基于多天线抗干扰通信模型及动态空间谱抗干扰方法
US8553797B2 (en) Channel information prediction system and channel information prediction method
CN110380997B (zh) 基于自适应压缩感知的毫米波信道估计方法
CN109495150B (zh) 一种毫米波通信码本设计方法
CN110870216B (zh) 一种波束成形方法及设备
US8462869B2 (en) Channel information prediction system and channel information prediction method
CN107743043A (zh) 多用户毫米波系统中基于带外空间信息的用户分组方法
CN102098705B (zh) 多维信道参数提取方法
CN111446999A (zh) 基于多臂强盗的位置辅助波束对准方法及其系统
CN108259397A (zh) 基于自适应正则化子空间追踪压缩感知算法的大规模mimo系统信道估计
CN106788631A (zh) 一种基于局部校准的大规模mimo互易性校准方法
CN104168046B (zh) 利用压缩感知的单端频域波束搜索方法
CN114448535B (zh) 基于fdd部分互易性的信道预测方法、网络侧设备及终端
CN104218984B (zh) 利用压缩感知的双端频域波束搜索方法
Odarchenko et al. Estimation of the communication range and bandwidth of UAV communication systems
CN105610476B (zh) 分离型数模混合收发通信系统中的模拟波束矢量优化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant