CN107135023B - 用于毫米波通信***的三维训练码书设计方法及波束对准方法 - Google Patents

用于毫米波通信***的三维训练码书设计方法及波束对准方法 Download PDF

Info

Publication number
CN107135023B
CN107135023B CN201710293707.XA CN201710293707A CN107135023B CN 107135023 B CN107135023 B CN 107135023B CN 201710293707 A CN201710293707 A CN 201710293707A CN 107135023 B CN107135023 B CN 107135023B
Authority
CN
China
Prior art keywords
training
codebook
dimensional
dimensional training
node
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
CN201710293707.XA
Other languages
English (en)
Other versions
CN107135023A (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 CN201710293707.XA priority Critical patent/CN107135023B/zh
Publication of CN107135023A publication Critical patent/CN107135023A/zh
Application granted granted Critical
Publication of CN107135023B publication Critical patent/CN107135023B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/0619—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 using feedback from receiving side
    • H04B7/0621—Feedback content
    • H04B7/0632—Channel quality parameters, e.g. channel quality indicator [CQI]

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Transmission System (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种用于毫米波通信***的三维训练码书设计方法及波束对准方法，其中三维训练码书设计方法包括如下步骤：1、根据分辨率和波束空间
Figure DDA0001282519760000011
的范围建立深度为S、度为N的三维训练码书树结构；2、将
Figure DDA0001282519760000012
等分为N个矩形区域
Figure DDA0001282519760000013
记为集合
Figure DDA0001282519760000014
Figure DDA0001282519760000015
有波束集合
Figure DDA0001282519760000016
与之对应；树结构的根节点为
Figure DDA0001282519760000017
与
Figure DDA0001282519760000018
的组合；3、确定第s层第b个节点Cs,b，Cs,b为
Figure DDA0001282519760000019
与波束集合
Figure DDA00012825197600000110
的组合；依次确定2到S层中的各节点
Figure DDA00012825197600000111
4、求解每个节点对应波束集合包含的波束成形向量
Figure DDA00012825197600000112
得到三维训练码书，其中q＝1,…,N。该方法生成的三维训练码书可用于实现高精度的波束对准和信道估计，并且能够显著降低***的训练开销。

Description

用于毫米波通信***的三维训练码书设计方法及波束对准 方法
技术领域
本发明属于无线通信技术领域，具体涉及到收发端天线呈均匀面阵排布的毫米波通信***的一种三维训练码书设计方法。
背景技术
随着无线通信技术的不断发展，高速数据业务以及无处不在的接入需求正呈现出一种***式的增长。下一代5G移动通信技术，对容量、能耗和带宽的需求将越来越高。工作在30-300GHz相对空闲频段的毫米波通信技术，因其工作频段中包含的大量未授权的带宽被认为是下一代无线局域网和移动通信的关键技术之一，因为直接扩展频谱带宽对于提升***容量来说既简单又有效。已有研究表明，毫米波通信能实现10Gbps的峰值传输速率。
然而，相较于传统微波频段，毫米波传输面临着更大的路径损耗，因此通信距离与覆盖范围十分有限。好在毫米波信号的波长极短，可以将大量天线以较小尺寸封装，进而大规模阵列天线与数模混合波束成形技术相结合，其提供的阵列增益和空分复用增益能弥补***的部分衰减，提高***传输速率和传输质量。除此之外，为更好地获取阵列增益，在通信开始时对收发波束进行波束对准也很有必要，高精度的波束对准在建立可靠的毫米波通信链路、获取所需的传输数据以及扩大小区通信覆盖范围方面发挥着关键作用。通过估计包括到达角(AoA)、离去角(AoD)以及路径增益等相关参数，准确的波束对准还可用于提高毫米波***的信道估计性能。
在实际的毫米波通信***中，要实现准确的波束对准有一定的困难。首先，毫米波段的高频带意味着信道可能在较短时间内快速变化，波束对准需要在极短的信道相干时间内完成，因此穷尽的波束搜索算法在此并不适用。第二，为充分利用大型天线阵列的阵列增益，训练波束应当足够窄，这无疑会增加波束对准的复杂度，因此提出高效的波束码本设计方法以及波束搜索算法很有必要。为减少波束对准的训练开销，一种有效的方法是使用基于分层多分辨率训练码书的树搜索算法。分层训练码书通常由不同级别的子码书组成，在高阶层上，子码书中包含少量低分辨率的训练波束，覆盖预定的角度范围；在低阶层上，子码书中包含的训练波束增多，分辨率也有了一定的提高。
虽然分层搜索能够显著降低***的训练开销，但其性能在很大程度上取决于所使用的层次训练码书。针对层次码书设计方法的研究已有很多，但这些研究主要集中在均匀线性阵列(ULA)结构，而对均匀平面阵列(UPA)结构的研究非常少。为实现可变精度的三维波束覆盖，并且获得更大的波束成形增益，本发明提出了一种适用于毫米波通信***的三维训练码书设计方法，该方法将训练波束的设计问题描述成波束图案的优化问题，并提出了一种有效的迭代算法来求解此优化问题。在本发明中，训练波束图案的过冲以及主旁瓣的纹波被适当地约束，因此每个训练波束都具有相对平坦的幅度响应以及较为理想的过渡频带。本发明设计的三维训练码书即使在低信噪比条件下，始终可以在毫米波通信***中实现高精度的波束对准以及信道估计。
发明内容
发明目的：针对现有技术中存在的问题，本发明公开了一种用于毫米波通信***的三维训练码书设计方法，该方法生成的三维训练码书可用于实现高精度的波束对准和信道估计，并且能够显著降低***的训练开销。
技术方案：本发明一方面公开了一种用于毫米波通信***的三维训练码书设计方法，包括如下步骤：
(1)根据分辨率rs和波束空间B的范围建立深度为S、度为N的三维训练码书树结构；
设分辨率rs＝ve×va，波束空间B在e方向的宽度为We，在a方向的宽度为Wa，即B为We×Wa的矩形区域，树结构的度N＝Ne×Na，深度S和度N满足条件：
Figure GDA0002371391930000021
所述三维训练码书树结构中1到S-1层中每一个节点都有N个子节点；
(2)将波束空间B在e方向进行Ne等分，在a方向进行Na等分，即将B等分为N个矩形区域
Figure GDA0002371391930000022
N＝Ne×Na，记为集合
Figure GDA0002371391930000023
每个矩形区域
Figure GDA0002371391930000024
有波束成形向量
Figure GDA0002371391930000025
与之对应，i＝1,…,N，N个波束成形向量组成波束集合
Figure GDA0002371391930000026
即
Figure GDA0002371391930000027
所述三维训练码书树结构的根节点C1,1为矩形区域集合
Figure GDA0002371391930000031
与波束集合
Figure GDA0002371391930000032
的组合
Figure GDA0002371391930000033
(3)确定所述三维训练码书树结构的第s层第b个节点Cs,b，Cs,b为矩形区域集合
Figure GDA0002371391930000034
与波束集合
Figure GDA0002371391930000035
的组合
Figure GDA0002371391930000036
其中
Figure GDA0002371391930000037
与第s-1层的节点
Figure GDA0002371391930000038
具有如下关系：
Figure GDA0002371391930000039
其中s＝2,…,S，b＝1,…,Ns-1，
Figure GDA00023713919300000310
为向上取整运算，％为取余运算；
对
Figure GDA00023713919300000311
在e方向进行Ne等分，在a方向进行Na等分，得到N个子矩形区域
Figure GDA00023713919300000312
每个子矩形区域
Figure GDA00023713919300000313
有波束成形向量
Figure GDA00023713919300000314
与之对应，i＝1,…,N，N个波束成形向量组成波束集合
Figure GDA00023713919300000315
即
Figure GDA00023713919300000316
按上述步骤依次确定2到S层中的各节点
Figure GDA00023713919300000317
(4)求解所述三维训练码书树结构中每个节点包含的波束成形向量
Figure GDA00023713919300000318
得到三维训练码书，其中q＝1,…,N。
优选地，所述三维训练码书树结构中每个节点包含的波束成形向量
Figure GDA00023713919300000319
在主瓣和旁瓣具有恒定但不相同的幅度响应值，即
Figure GDA00023713919300000320
其中ψ＝(ψe,ψa)＝(sin(θe),sin(θa))是俯仰角θe与方位角θa的正弦值的组合，
Figure GDA00023713919300000321
为水平方向上的天线阵列响应矢量ah(2π/λdh sin(θe)cos(θa))与垂直方向上的天线阵列响应矢量av(2π/λdv sin(θa))的Kronecker乘积，
Figure GDA00023713919300000322
是波束成形向量
Figure GDA00023713919300000323
的幅度响应值；其中q＝1,…,N。
优选地，步骤(4)中求解波束成形向量
Figure GDA0002371391930000041
为求解优化问题：
Figure GDA0002371391930000042
其中ε为训练波束主旁瓣的纹波，是一个极小的正实数，
Figure GDA0002371391930000043
分别为主旁瓣对应区域。
优选地，对
Figure GDA0002371391930000044
进行采样离散化为离散波束主、旁瓣对应区域
Figure GDA0002371391930000045
Figure GDA0002371391930000046
区域。
为避免离散采样之后过冲现象的发生，引入约束条件
Figure GDA0002371391930000047
步骤(4)中求解波束成形向量
Figure GDA0002371391930000048
为求解优化问题：
Figure GDA0002371391930000049
其中ε为训练波束主旁瓣的纹波，是一个极小的正实数；Es为常数，取值随着s的增大而减小；
Figure GDA00023713919300000410
分别为离散波束主、旁瓣对应区域。
优选地，将约束条件松弛化，采用约束凹凸过程迭代算法求解优化问题，步骤如下：
(6.1)确定迭代初始值f0，构建如下约束优化问题：
Figure GDA00023713919300000411
其中
Figure GDA00023713919300000413
表示取括号内数值的实部，
Figure GDA00023713919300000412
分别表示离散化的主旁瓣对应区域；
(6.2)迭代求解步骤(6.1)中的优化问题，检查本次迭代得到的r值是否满足收敛标准，若满足，本次迭代的最优解fn即为最终解
Figure GDA00023713919300000511
若不满足，根据本次迭代的最优解fn，再次求解步骤(6.1)中的优化问题；
(6.3)输出最终解
Figure GDA00023713919300000512
获得所需的三维训练波束。
优选地，步骤(6.1)中初始值f0通过求解如下优化问题来确定：
Figure GDA0002371391930000051
本发明另一方面公开了一种波束对准方法，包括如下步骤：
(8.1)参数配置与初始化：收发端按照上述任一种三维训练码书设计方法设计三维训练码书；
设所设计的三维训练码书含有S个训练子码书
Figure GDA0002371391930000052
分别对应S个训练阶段；所设计的三维训练码书树结构的度为N；波束成形向量
Figure GDA0002371391930000053
中的下标值初始化为：s＝1，bf＝1；
(8.2)发射端连续使用第一层子码书中的低分辨率训练波束向量
Figure GDA0002371391930000054
来传输训练信号z，且每个训练波束重复训练Ms次以增加***的接收信噪比；
(8.3)接收端根据对应的接收信号
Figure GDA0002371391930000055
选择产生最高接收能量的波束成形向量，即
Figure GDA0002371391930000056
其中
Figure GDA0002371391930000057
p为发送功率，h为信道矩阵，hH为h的共轭转置矩阵，
Figure GDA0002371391930000058
为加性高斯噪声；
接收端将索引值
Figure GDA00023713919300000513
反馈给发射端，由此，发射端根据
Figure GDA0002371391930000059
及其对应的矩形区域
Figure GDA00023713919300000510
初步确定最强波束对应的离去角AoD所在区域；
(8.4)发射端根据索引值
Figure GDA0002371391930000063
在三维码书的下一层子码书中选择一组具有较高分辨率的波束成形矢量
Figure GDA0002371391930000061
来传输训练信号，从而进一步更精确地确定AoD所在区域；
其中
Figure GDA0002371391930000064
(8.5)重复步骤(8.3)～(8.4)，直到查找到满足最高分辨率要求的波束成形向量，该向量对应的矩形域即为AoD所在区域。
有益效果：本发明提供的适用于毫米波通信***的三维训练码书设计方法，可用于实现高精度的波束对准和信道估计，其有益效果在于：
(1)现有的研究主要针对天线呈ULA排布的毫米波***，本发明充分考虑了UPA结构，弥补了现有研究的不足。在UPA结构下，有限的二维空间内能够配置更多的传输天线，因此可以检测侦听更广的地理区域，实现可变精度的三维波束覆盖，并且获得更大的波束成形增益。
(2)本发明对训练波束图案中未采样点上可能发生的过冲现象进行了处理，同时还对主旁瓣的纹波进行了适当地约束，因此每个训练波束都具有相对平坦的幅度响应以及较为理想的过渡频带，有利于提高波束对准精度。
(3)本发明生成的三维训练码书适用于基于树搜索算法的波束对准算法，能显著降低***的训练开销，并可在基于UPA结构的毫米波通信***中实现高精度的波束对准和信道估计性能。
附图说明
图1为均匀平面阵列(UPA)的结构示意图；
图2为实施例中生成的三维训练码书的结构示意图；
图3为主旁瓣幅度波动示意图；
图4为实施例中生成的4个训练波束矢量
Figure GDA0002371391930000062
的幅度响应图；
图5为两种码书设计方法(BPSA、LSA)下的训练波束图；
图6为本发明提出的三维训练码书设计方法的具体流程图以及码书生成后用于波束训练的流程图；
图7为两种码书设计方法(BPSA、LSA)下，平均波束对准错误率随信噪比的变化趋势图。
具体实施方式
下面结合附图和具体实施方式，进一步阐明本发明。
如图1所示，本发明公开的三维训练码书设计方法适用于天线呈均匀面阵排布的毫米波***，以发射端为例，图中Nh、Nv分别表示水平方向和垂直方向上的天线数，因此总发射天线数NT＝NhNv，dh和dv分别为水平方向、垂直方向上相邻两根天线之间的间距。
波束空间B为一个矩形区域，设e方向与a方向为垂直的两个方向，B在e方向的宽度为We，在a方向的宽度为Wa，即B为We×Wa的矩形区域。
本实施例考虑发射端天线呈均匀面阵排布的毫米波通信***，其中Nh＝Nv＝32，相邻两根天线之间的间距dh＝dv＝3λ/8，λ为毫米波信号的波长。为简单起见，假设接收端只配备有单根天线。需要注意的是，虽然本实例只以发射端的三维训练码书设计为例，但这种方法同样适用于接收端的码书设计。
一种用于毫米波通信***的三维训练码书设计方法，包括如下步骤：
步骤1、根据分辨率rs和波束空间B的范围建立深度为S、度为N的三维训练码书树结构；
设分辨率rs＝ve×va；若发射波束的俯仰角θe、方位角θa分别位于角度范围
Figure GDA0002371391930000073
之内，且俯仰角与方位角是独立不相关变量，可以用垂直的两个方向e方向与a方向的变量来表示，则本实施例中波束空间B可定义为sin(θe)、sin(θa)对应范围的乘积，即B＝[-sin(θe),sin(θe)]×[-sin(θa),sin(θa)]＝[-0.64,0.64)]×[-0.64,0.64)]，其在(sin(θe),sin(θa))平面是一个矩形域，即We＝2sin(θe)，Wa＝2sin(θa)。
三维训练码书树结构的度N＝Ne×Na，深度S和度N满足条件：
Figure GDA0002371391930000071
Figure GDA0002371391930000072
即将波束空间B等分为NS个小矩形区域后，每个小矩形区域要小于等于分辨率。
所述三维训练码书树结构中1到S-1层中每一个节点都有N个子节点；
步骤2、将波束空间B在e方向进行Ne等分，在a方向进行Na等分，即将B等分为N个矩形区域
Figure GDA0002371391930000081
N＝Ne×Na，记为集合
Figure GDA0002371391930000082
每个矩形区域
Figure GDA0002371391930000083
有波束成形向量
Figure GDA0002371391930000084
与之对应，i＝1,…,N，N个波束成形向量组成波束集合
Figure GDA0002371391930000085
即
Figure GDA0002371391930000086
所述三维训练码书树结构的根节点C1,1为矩形区域集合
Figure GDA0002371391930000087
与波束集合
Figure GDA0002371391930000088
的组合
Figure GDA0002371391930000089
步骤3、确定所述三维训练码书树结构的第s层第b个节点Cs,b，Cs,b为矩形区域集合
Figure GDA00023713919300000810
与波束集合
Figure GDA00023713919300000811
的组合
Figure GDA00023713919300000812
其中
Figure GDA00023713919300000813
与第s-1层的节点
Figure GDA00023713919300000814
具有如下关系：
Figure GDA00023713919300000815
其中s＝2,…,S，b＝1,…,Ns-1，
Figure GDA00023713919300000816
为向上取整运算，％为取余运算；
对
Figure GDA00023713919300000817
在e方向进行Ne等分，在a方向进行Na等分，得到N个子矩形区域
Figure GDA00023713919300000818
每个子矩形区域
Figure GDA00023713919300000819
有波束成形向量
Figure GDA00023713919300000820
与之对应，i＝1,…,N，N个波束成形向量组成波束集合
Figure GDA00023713919300000821
即
Figure GDA00023713919300000822
且
Figure GDA00023713919300000823
对应；
按上述步骤依次确定2到S层中的各节点
Figure GDA00023713919300000824
步骤4、求解所述三维训练码书树结构中每个节点包含的波束成形向量
Figure GDA00023713919300000825
得到三维训练码书，其中q＝1,…,N。
三维训练码书(CodeBook)的树结构有S层，每一层为一个子码书(Sub-CodeBook)，对应一个训练阶段；第s层的子码书有Ns-1个节点，每个节点包括N个码字(CodeWord)，每个节点由波束集合
Figure GDA00023713919300000826
及其对应的矩形区域
Figure GDA0002371391930000091
来表示；
Figure GDA0002371391930000092
及其对应的矩形区域
Figure GDA0002371391930000093
的组合
Figure GDA0002371391930000094
为一个码字。
如图2所示，本实施例中生成的三维训练码书为一个四叉树结构，即三维训练码书树结构的度N＝Ne×Na＝2×2。由对应S个训练阶段的S个子码书构成；在用于第s个训练阶段的子码书
Figure GDA0002371391930000095
中，总共包含4s个具有相同主瓣宽度的波束成形向量，且这些向量构成4s-1个波束集合
Figure GDA0002371391930000096
其中第b个集合
Figure GDA0002371391930000097
包含波束成形向量
Figure GDA0002371391930000098
子码书
Figure GDA0002371391930000099
中的波束成形向量个数是上一层子码书
Figure GDA00023713919300000910
中训练波束个数的4倍。在第s个训练阶段，波束空间被分割成4s-1个相同大小的矩形域
Figure GDA00023713919300000911
其中第b个矩形域
Figure GDA00023713919300000912
与波束集合
Figure GDA00023713919300000913
相关，且
Figure GDA00023713919300000914
同时波束成形向量
Figure GDA00023713919300000915
分别与
Figure GDA00023713919300000916
一一对应；在下一个训练阶段，矩形域
Figure GDA00023713919300000917
又被等分成4个更小的子区域，j∈{LL,LR,RL,RR}，即
Figure GDA00023713919300000918
Figure GDA00023713919300000919
为提高毫米波通信的波束对准性能，所设计的波束成形向量
Figure GDA00023713919300000920
分别在主瓣和旁瓣上具有恒定但不相同的幅度响应值，即：
Figure GDA00023713919300000921
其中ψ＝(ψe,ψa)＝(sin(θe),sin(θa))是俯仰角θe与方位角θa的正弦值的组合，
Figure GDA00023713919300000922
为水平方向、垂直方向上的天线阵列响应矢量ah(2π/λdh sin(θe)cos(θa))及av(2π/λdv sin(θa))的Kronecker乘积，dh和dv分别为水平方向、垂直方向上相邻两根天线之间的间距，λ为毫米波信号的波长，
Figure GDA00023713919300000923
是波束向量
Figure GDA00023713919300000924
的幅度响应值。
若训练波束主旁瓣的纹波ε被严格限制为0，则难以成功实现三维码书的设计，为避免这种不可行性，同时为确保高精度的波束对准性能，如图3所示，可允许波束图案的主旁瓣存在小幅度波动，同时在纹波ε的值固定的情况下，主瓣的幅度响应值应尽可能大，因此训练波束的设计可以表示为如下优化问题：
Figure GDA0002371391930000101
其中ε是一个极小的正实数，
Figure GDA0002371391930000102
分别为主、旁瓣对应区域，
Figure GDA0002371391930000103
为矩形区域
Figure GDA0002371391930000104
内，
Figure GDA0002371391930000105
为
Figure GDA0002371391930000106
外。由于
Figure GDA0002371391930000107
连续不可数，因此必须对波束主旁瓣对应区域进行采样或离散化；又由于采样点数量有限，在未采样的区域可能存在过冲现象，为避免离散采样之后过冲现象的发生，引入约束条件||f||≤Es，Es的取值随着s的增大而减小，||·||为求向量的2范数，即优化问题为：
Figure GDA0002371391930000108
本实施例采用如下的有限散射信道模型：
Figure GDA0002371391930000109
其中L为信道总数，αl是第l条路径的复增益，β是平均路径损耗，
Figure GDA00023713919300001010
是天线阵列响应矢量，且ab(b∈{h,v})具有如下形式：
Figure GDA00023713919300001011
Figure GDA00023713919300001012
为Kronecker乘积运算符。
以三维训练码书中的一个训练码字中的波束成形向量
Figure GDA00023713919300001013
为例，该码字在离散波束区域
Figure GDA00023713919300001014
的幅度响应值较大，在其他区域
Figure GDA00023713919300001015
的幅度响应极小，该训练波束的设计问题是一个非凸优化问题，求解困难。通过将部分非凸约束条件松弛化，该非凸优化问题可被转换成凸优化问题，并可利用约束凹凸过程(CCCP)迭代算法求解。若fn表示第n次迭代得到的最优解，则下一次迭代的最优解fn+1可通过求解如下问题得到：
Figure GDA0002371391930000111
其中
Figure GDA0002371391930000112
表示取括号内数值的实部，
Figure GDA0002371391930000113
分别表示离散化的主旁瓣对应区域。
由于初始化对CCCP算法的收敛性能有较大影响，为得到一个良好的初始值，可求解如下优化问题：
Figure GDA0002371391930000114
此问题是一个二阶锥规划问题(SOCP)，可通过MATLAB仿真平台中的CVX工具箱求解。
在得到良好的初始值f0之后，基于CCCP算法的三维训练波束设计方法具体包括如下步骤：
(6.1)构建并求解如下约束优化问题：
Figure GDA0002371391930000115
其中
Figure GDA0002371391930000116
表示取括号内数值的实部，
Figure GDA0002371391930000117
分别表示离散化的主旁瓣对应区域；
(6.2)迭代求解步骤(6.1)中的优化问题，检查本次迭代得到的r值是否满足收敛标准，若满足，本次迭代的最优解fn即为最终解f＊；若不满足，根据本次迭代的最优解fn，再次求解步骤(6.1)中的优化问题；
(6.3)输出最终解f＊，获得所需的三维训练波束。
图4以三维码书中的第一层子码书为例，给出了利用本发明的设计方法生成的4个训练码字
Figure GDA0002371391930000121
的幅度响应图。可以发现，这几个码字的波束图案都具有相对平坦的幅度响应以及较为理想的过渡频带，而且在未采样点上也没有过冲现象发生。图5还将本发明与基于最小二乘(LS)的波束设计方法相比较，其中BPSA表示本发明提出的三维码书设计算法，LSA表示最小二乘波束设计算法。从图中可以看出，利用LS法生成的训练波束的旁瓣具有相对较大波动，此外，其过渡带的宽度也较大，这进一步体现了本发明提出的三维训练码书设计方法的优越性。
本发明生成的三维训练码书可用于实现基于树搜索算法的波束对准，检测出单径毫米波信道的最强波束，也可用于估计毫米波通信信道的相关参数，如离去角(AoD)和到达角(AoA)，具体包括如下步骤：
(8.1)参数配置与初始化：收发端按照权利要求1-7中任一项所述的方法设计三维训练码书；
设所设计的三维训练码书树结构S层，即含有S个训练子码书
Figure GDA0002371391930000122
分别对应S个训练阶段；所设计的三维训练码书树结构的度为N；波束成形向量
Figure GDA0002371391930000123
中的下标值初始化为：s＝1，bf＝1；
(8.2)发射端连续使用第一层子码书中的低分辨率训练波束向量
Figure GDA0002371391930000124
来传输训练信号z，且每个训练波束重复训练Ms次以增加***的接收信噪比；
(8.3)接收端根据对应的接收信号
Figure GDA0002371391930000125
选择产生最高接收能量的波束成形向量，即
Figure GDA0002371391930000126
其中
Figure GDA0002371391930000127
p为发送功率，h为信道矩阵，hH为h的共轭转置矩阵，
Figure GDA0002371391930000128
为加性高斯噪声；
接收端将索引值q＊反馈给发射端，由此，发射端根据
Figure GDA0002371391930000129
及其对应的矩形区域
Figure GDA0002371391930000131
初步确定最强波束对应的离去角AoD所在区域；
(8.4)发射端根据索引值q＊，在三维码书的下一层子码书中选择一组具有较高分辨率的波束成形矢量
Figure GDA0002371391930000132
来传输训练信号，从而进一步更精确地确定AoD所在区域；
其中bf+1＝bf*N-N+q＊；
对于本实施例中的四叉树结构，bf的更新遵循如下原则：
a)若q＊＝LL，则bf+1＝4bf-3；
b)若q＊＝LR，则bf+1＝4bf-2；
c)若q＊＝RL，则bf+1＝4bf-1；
d)若q＊＝RR，则bf+1＝4bf；
(8.5)重复步骤(8.3)～(8.4)，直到查找到满足最高分辨率要求的波束成形向量，该向量对应的矩形域即为AoD所在区域。
图6给出了本发明提出的三维训练码书设计方法(BPSA)的具体流程，除此之外，图6也给出了码书生成之后用于波束训练的具体过程。如图7所示，该图表示当毫米波通信信道的路径增益固定为1的情况下，利用LS法以及本发明方法生成的码书进行波束训练，产生的平均波束对准错误概率(BAER)随信噪比(SNR)的变化趋势。从图中可以看出，本发明方法生成的码本性能明显优于LSA码本，能够实现更精确的波束对准和信道估计性能。随着SNR的增大，两种码书对应的BAER随之减小，但BPSA码本对应的错误概率衰减明显，而LSA码本对应的错误概率衰减缓慢。

Claims (2)

1.一种用于毫米波通信***的三维训练码书设计方法，其特征在于，包括如下步骤：
(1)根据分辨率rs和波束空间
Figure FDA0002371391920000011
的范围建立深度为S、度为N的三维训练码书树结构；
设分辨率rs＝ve×va，波束空间
Figure FDA0002371391920000012
为矩形区域，e方向与a方向为垂直的两个方向；
Figure FDA0002371391920000013
在e方向的宽度为We，在a方向的宽度为Wa，即
Figure FDA0002371391920000014
为We×Wa的矩形区域，树结构的度N＝Ne×Na，深度S和度N满足条件：
Figure FDA0002371391920000015
且
Figure FDA0002371391920000016
所述三维训练码书树结构中1到S-1层中每一个节点都有N个子节点；
(2)将波束空间
Figure FDA0002371391920000017
在e方向进行Ne等分，在a方向进行Na等分，即将
Figure FDA0002371391920000018
等分为N个矩形区域
Figure FDA0002371391920000019
N＝Ne×Na，记为集合
Figure FDA00023713919200000110
每个矩形区域
Figure FDA00023713919200000111
有波束成形向量
Figure FDA00023713919200000112
与之对应，i＝1，…，N，N个波束成形向量组成波束集合
Figure FDA00023713919200000113
即
Figure FDA00023713919200000114
所述三维训练码书树结构的根节点C1，1为矩形区域集合
Figure FDA00023713919200000115
与波束集合
Figure FDA00023713919200000116
的组合
Figure FDA00023713919200000117
(3)确定所述三维训练码书树结构的第s层第b个节点Cs，b，Cs，b为矩形区域集合
Figure FDA00023713919200000118
与波束集合
Figure FDA00023713919200000119
的组合
Figure FDA00023713919200000120
其中
Figure FDA00023713919200000121
与第s-1层的节点
Figure FDA00023713919200000122
具有如下关系：
Figure FDA00023713919200000123
其中s＝2，…，S，b＝1，…，Ns-1，
Figure FDA00023713919200000124
为向上取整运算，％为取余运算；
对
Figure FDA00023713919200000125
在c方向进行Ne等分，在a方向进行Na等分，得到N个子矩形区域
Figure FDA00023713919200000126
每个子矩形区域
Figure FDA00023713919200000127
有波束成形向量
Figure FDA00023713919200000128
与之对应，i＝1，…，N，N个波束成形向量组成波束集合
Figure FDA0002371391920000021
即
Figure FDA0002371391920000022
按上述步骤依次确定2到S层中的各节点
Figure FDA0002371391920000023
(4)求解所述三维训练码书树结构中每个节点包含的波束成形向量
Figure FDA0002371391920000024
得到三维训练码书，其中q＝1，…，N；所述三维训练码书树结构中每个节点包含的波束成形向量
Figure FDA0002371391920000025
在主瓣和旁瓣具有恒定但不相同的幅度响应值，即
Figure FDA0002371391920000026
其中ψ＝(ψe，ψa)＝(sin(θe)，sin(θa))是俯仰角θe与方位角θa的正弦值的组合，
Figure FDA0002371391920000027
为水平方向上的天线阵列响应矢量ah(2π/λdh sin(θe)cos(θa))与垂直方向上的天线阵列响应矢量av(2π/λdvsin(θa))的Kronecker乘积，
Figure FDA0002371391920000028
是波束成形向量
Figure FDA0002371391920000029
的幅度响应值；其中q＝1，...，N；dn和dv分别为水平方向、垂直方向上相邻两根天线之间的间距；
对主旁瓣对应区域
Figure FDA00023713919200000210
和
Figure FDA00023713919200000211
进行采样离散化为离散波束主、旁瓣对应区域
Figure FDA00023713919200000212
和
Figure FDA00023713919200000213
区域；
求解波束成形向量
Figure FDA00023713919200000214
为求解优化问题：
Figure FDA00023713919200000215
其中ε为训练波束主旁瓣的纹波，是一个极小的正实数；Es为常数，取值随着s的增大而减小；
Figure FDA00023713919200000216
和
Figure FDA00023713919200000217
分别为离散波束主、旁瓣对应区域；
将约束条件松弛化，采用约束凹凸过程迭代算法求解优化问题，步骤如下：
(6.1)确定迭代初始值f0：
初始值f0通过求解如下优化问题来确定：
Figure FDA0002371391920000031
构建如下约束优化问题：
Figure FDA0002371391920000032
其中
Figure FDA0002371391920000033
Figure FDA0002371391920000034
表示取括号内数值的实部，
Figure FDA0002371391920000035
分别表示离散化的主旁瓣对应区域；
(6.2)迭代求解步骤(6.1)中的优化问题，检查本次迭代得到的r值是否满足收敛标准，若满足，本次迭代的最优解fn即为最终解f＊；若不满足，根据本次迭代的最优解fn，再次求解步骤(6.1)中的优化问题；
(6.3)输出最终解f＊，获得所需的三维训练波束。
2.一种波束对准方法，其特征在于，包括如下步骤：
(8.1)参数配置与初始化：收发端按照权利要求1中所述的方法设计三维训练码书；
设所设计的三维训练码书含有S个训练子码书
Figure FDA0002371391920000036
分别对应S个训练阶段；所设计的三维训练码书树结构的度为N；波束成形向量
Figure FDA0002371391920000037
中的下标值初始化为：s＝1，bf＝1；
(8.2)发射端连续使用第一层子码书中的低分辨率训练波束向量
Figure FDA0002371391920000038
来传输训练信号z，且每个训练波束重复训练Ms次以增加***的接收信噪比；
(8.3)接收端根据对应的接收信号
Figure FDA0002371391920000039
选择产生最高接收能量的波束成形向量，即
Figure FDA0002371391920000041
其中
Figure FDA0002371391920000042
p为发送功率，h为信道矩阵，hH为h的共轭转置矩阵，
Figure FDA0002371391920000043
为加性高斯噪声；
接收端将索引值q＊反馈给发射端，由此，发射端根据
Figure FDA0002371391920000044
及其对应的矩形区域
Figure FDA0002371391920000045
初步确定最强波束对应的离去角AoD所在区域；
(8.4)发射端根据索引值q＊，在三维码书的下一层子码书中选择一组具有较高分辨率的波束成形矢量
Figure FDA0002371391920000046
来传输训练信号，从而进一步更精确地确定AoD所在区域；
其中bf+1＝bf*N-N+q＊；
(8.5)重复步骤(8.3)～(8.4)，直到查找到满足最高分辨率要求的波束成形向量，该向量对应的矩形域即为AoD所在区域。
CN201710293707.XA 2017-04-28 2017-04-28 用于毫米波通信***的三维训练码书设计方法及波束对准方法 Active CN107135023B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710293707.XA CN107135023B (zh) 2017-04-28 2017-04-28 用于毫米波通信***的三维训练码书设计方法及波束对准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710293707.XA CN107135023B (zh) 2017-04-28 2017-04-28 用于毫米波通信***的三维训练码书设计方法及波束对准方法

Publications (2)

Publication Number Publication Date
CN107135023A CN107135023A (zh) 2017-09-05
CN107135023B true CN107135023B (zh) 2020-04-24

Family

ID=59716151

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710293707.XA Active CN107135023B (zh) 2017-04-28 2017-04-28 用于毫米波通信***的三维训练码书设计方法及波束对准方法

Country Status (1)

Country Link
CN (1) CN107135023B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108809397B (zh) * 2018-06-27 2020-06-30 东南大学 多天线***中高功效数模混合波束成形方法、装置及设备
CN109379122B (zh) * 2018-11-08 2021-03-16 东南大学 一种毫米波通信多径信道动态波束训练方法
CN109474549B (zh) * 2018-12-04 2021-08-17 青岛矽昌通信技术有限公司 一种基于三维波束图案的三维信道估计方法
CN111698007B (zh) * 2019-03-15 2021-04-16 大唐移动通信设备有限公司 一种基于混合波束赋形架构的校准补偿方法及装置
CN111478724A (zh) * 2020-04-15 2020-07-31 南京航空航天大学 一种面向无人机毫米波平台的三维波束搜索方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104935366A (zh) * 2015-04-21 2015-09-23 北京航空航天大学 毫米波通信中的波束搜索方法
CN105610478A (zh) * 2016-02-15 2016-05-25 东南大学 毫米波mimo通信多子阵协作波束对准方法及装置
CN105790811A (zh) * 2016-04-15 2016-07-20 东南大学 一种高分辨率模拟波束快速训练方法及装置
CN105959044A (zh) * 2016-04-21 2016-09-21 北京航空航天大学 联合方法的层次码本结构设计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9100068B2 (en) * 2008-03-17 2015-08-04 Qualcomm, Incorporated Multi-resolution beamforming in MIMO systems

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104935366A (zh) * 2015-04-21 2015-09-23 北京航空航天大学 毫米波通信中的波束搜索方法
CN105610478A (zh) * 2016-02-15 2016-05-25 东南大学 毫米波mimo通信多子阵协作波束对准方法及装置
CN105790811A (zh) * 2016-04-15 2016-07-20 东南大学 一种高分辨率模拟波束快速训练方法及装置
CN105959044A (zh) * 2016-04-21 2016-09-21 北京航空航天大学 联合方法的层次码本结构设计方法

Also Published As

Publication number Publication date
CN107135023A (zh) 2017-09-05

Similar Documents

Publication Publication Date Title
CN107135023B (zh) 用于毫米波通信***的三维训练码书设计方法及波束对准方法
Chen et al. Sum-rate maximization for intelligent reflecting surface based terahertz communication systems
EP3577776B1 (en) Transmitter and method for transmitting symbols over wireless communication channel
CN104935366B (zh) 毫米波通信中的波束搜索方法
Song et al. Codebook design for hybrid beamforming in millimeter wave systems
CN109946654B (zh) 雷达阵元-脉冲编码与处理方法
KR20230045455A (ko) 무선 통신 시스템에서 지능형 반사 평면을 기반한 데이터 통신을 위한 장치 및 방법
US11843434B2 (en) Millimeter wave coarse beamforming using outband sub-6GHz reconfigurable antennas
Lu et al. A tutorial on near-field XL-MIMO communications towards 6G
Wu et al. Two-stage 3D codebook design and beam training for millimeter-wave massive MIMO systems
CN107231180B (zh) 一种基于波束成形的动态方向调制方法
Brady et al. Multi-beam MIMO prototype for real-time multiuser communication at 28 GHz
KR20200026983A (ko) 빔형성 방법 및 디바이스
Ghermezcheshmeh et al. Parametric channel estimation for LoS dominated holographic massive MIMO systems
Singh et al. Fast beam training for RIS-assisted uplink communication
You et al. IP aided OMP based channel estimation for millimeter wave massive MIMO communication
Pulipati et al. Xilinx RF-SoC-based digital multi-beam array processors for 28/60 GHz wireless testbeds
Zou et al. Step-wisely refinement based beam searching scheme for 60 GHz communications
Xue et al. An efficient beam-training scheme for the optimally designed subarray structure in mmWave LoS MIMO systems
Mahfoudi et al. Joint range extension and localization for low‐power wide‐area network
Gao et al. Wideband hybrid analog-digital beamforming massive MIMO systems based on Rotman lens
Wang et al. Terahertz Near-Field Communications and Sensing
Hadji et al. Channel estimation based on low‐complexity hierarchical codebook design for millimeter‐wave MIMO systems
Song et al. High throughput line-of-sight MIMO systems for next generation backhaul applications
Zhang et al. Codebook design and beam training for intelligent omni-surface aided communications

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