CN107181511B - 一种毫米波mimo系统的混合预编码方法及系统 - Google Patents

一种毫米波mimo系统的混合预编码方法及系统 Download PDF

Info

Publication number
CN107181511B
CN107181511B CN201710422927.8A CN201710422927A CN107181511B CN 107181511 B CN107181511 B CN 107181511B CN 201710422927 A CN201710422927 A CN 201710422927A CN 107181511 B CN107181511 B CN 107181511B
Authority
CN
China
Prior art keywords
precoding
initial
digital
coding
analog
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
CN201710422927.8A
Other languages
English (en)
Other versions
CN107181511A (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.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
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 Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201710422927.8A priority Critical patent/CN107181511B/zh
Publication of CN107181511A publication Critical patent/CN107181511A/zh
Application granted granted Critical
Publication of CN107181511B publication Critical patent/CN107181511B/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
    • 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

Landscapes

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

Abstract

本发明涉及无线通信技术领域，特别涉及一种毫米波MIMO系统的混合预编码方法及系统。所述毫米波MIMO系统的混合预编码方法包括：步骤a：建立毫米波MIMO系统的混合预编码问题模型；步骤b：给定初始模拟预编码，基于所述初始模拟预编码，按照最小均方原则得到数字预编码；步骤c：给定初始数字预编码，根据所述初始数字预编码对混合预编码问题模型进行优化，并采用搜索法得到模拟预编码；步骤d：迭代执行步骤b和步骤c，直到达到混合预编码迭代次数值，并输出所述数字预编码和模拟预编码。本发明降低了求导操作带来的高计算复杂度，同时大大减少了运算时间，直接搜索算法能够灵活调整搜索步长，避免基误匹配现象，提高成形波束的增益。

Description

一种毫米波MIMO系统的混合预编码方法及系统
技术领域
本发明涉及无线通信技术领域，特别涉及一种毫米波MIMO系统的混合预编码方法及系统。
背景技术
毫米波通信自由路损的急剧增加，早期关于毫米波通信的研究大多集中于室内场景。然而，得益于毫米等级的波长，大规模的毫米波通信系统天线阵列可以密集排列在较小体积的印制电路板上，因此基于大规模天线阵列的波束成形技术被广泛应用于实际的毫米波通信系统以对抗剧烈的路径损耗。然而，全数字的波束成形技术要求为每一根天线配置单独的数字收发通道，在使用大规模天线阵列的毫米波通信系统中，其维度可能会使得数字信号处理的复杂度和功率消耗难以接受，因此全数字的波束成形技术并不适合实际的毫米波通信系统。针对这一问题，一种模拟和数字器件相结合的波束成形方案在毫米波通信系统中受到了广泛的关注。在此方案中，发射端的基带源信号首先被送入数字基带处理器，然后数字基带处理器的输出由若干条数字发射通道送入由相位转换器组成的模拟单元，最终转换成基带发射信号。接收端的流程与发射端的流程是一一对应的。由于相位转换器只能改变输入信号的相位，并不能调整输入信号的幅度，因此被称为模拟波束成形器。另一方面，数字收发通道不再与阵列天线对应相连，相比于全数字的波束成形方案其数量可以大幅降低。然而，混合波束成形设计问题属于非线性优化难题，用一般算法不可能在较短时间内找到最优解。因此，现有研究将混合模拟数字波束成形设计问题建模为常模限制的矩阵分解优化问题求其近似最优解。
对现有技术进行检索发现，Omar El Ayach等在2014年IEEE Transactions onWireless Communications上发表的Spatially sparse precoding in millimeter waveMIMO systems(毫米波多入多出系统中空间稀疏预编码)中基于毫米波波束空间的稀疏性，将波束成形设计问题转化为二维搜索问题，提出了一种搜索范围基于二维收发天线阵列响应候选集的正交匹配追踪方法成形波束。M.Kim and Y.Lee等在2015年的IEEETrans.Veh.Technol上发表的MSE-based hybrid RF/baseband processing formillimeter wave communication systems in MIMO interference channels、Rusu C,Mendez-Rial R,Gonzalez-Prelcicy N and Robert W.Heath.Jr等在2015年的IEEEInternational Conference on Communications上发表的Low complexity hybridsparse precoding and combining in millimeter wave MIMO systems，以及Yun-YuehLee,Ching-Hung Wang and Yuan-Hao Huang等在2014年的IEEE Transactions on SignalProcessing上发表的A hybrid RF/baseband precoding processor based on parallel-index-selection matrix-inversion-bypass simultaneous orthogonal matchingpursuit for millimeter wave MIMO systems中分别优化了二维收发天线阵列响应候选集、用简单的相关计算设计替代复杂的匹配追踪迭代以及简化矩阵逆操作。Chiao-En Chen在2015年IEEE Wireless Communications Letters上发表的An iterative hybridtransceiver design algorithm for millimeter wave MIMO systems(毫米波多入多出系统迭代的混合收发机设计算法)中，基于求导迭代的一维下山单纯形法实现局部搜索。Sohrabi Foad等在2016年IEEE Journal of Selected Topics in Signal Processing上发表的Hybrid Digital and Analog Beamforming Design for Large-Scale AntennaArrays(大规模天线阵列的数字和模拟波束成形设计)中，在收发机配置大规模天线阵列的假设条件下，近似得到直接求解混合波束形成的方法。Xianghao Yu等在2016年的IEEEJournal of Selected Topics in Signal Processing上发表的Alternatingminimization algorithms for hybrid precoding in millimeter wave MIMO systems(基于交替最小化算法的毫米波多入多出系统混合预编码)中基于导频序列，采用黎曼优化技术实现了一种混合波束成形设计方法并基于数字波束成形器相互正交的限制，提出一种快速混合波束形成算法。
上述中，正交匹配追踪本质上是一种基于穷尽搜索的贪婪算法，会造成不能容忍的训练时延、能耗和计算复杂度。此外，基于正交匹配追踪技术的波束成形方法及其相应的简化策略在实际使用中，将严重受限于搜索候选集合的设计，不符合成形波束参数连续取值的事实，极有可能造成基误匹配现象，使得形成的天线阵列波束增益不高，严重降低了波束成形技术的性能表现。而基于导数复杂计算寻找下降方向及步长的下山单纯形及黎曼优化方法都会造成不可避免的高计算复杂度；特别地，基于每次更新一个波束形成矩阵元素的一维局部搜索方法加剧了系统的时延代价。基于大规模天线阵列假设的方法并不适用于小规模系统，而且方式实施中需要反复对矩阵求逆，给系统带来较高的运算复杂度。基于数字波束成形矢量相互正交假设前提的快速算法因为隐性的附加限制会带来不可避免的性能表现损失。
发明内容
本发明提供了一种毫米波MIMO系统的混合预编码方法及系统，旨在至少在一定程度上解决现有技术中的上述技术问题之一。
为了解决上述问题，本发明提供了如下技术方案：
一种毫米波MIMO系统的混合预编码方法，包括以下步骤：
步骤a：建立毫米波MIMO系统的混合预编码问题模型；
步骤b：给定初始模拟预编码，基于所述初始模拟预编码，按照最小均方原则得到数字预编码；
步骤c：给定初始数字预编码，根据所述初始数字预编码对混合预编码问题模型进行优化，并采用搜索法得到模拟预编码；
步骤d：迭代执行步骤b和步骤c，直到达到混合预编码迭代次数值，并输出所述数字预编码和模拟预编码。
本发明实施例采取的技术方案还包括：在所述步骤a中，所述混合预编码问题模型为：
Figure BDA0001315461150000051
Figure BDA0001315461150000052
上述公式中，FRF和FBB分别为模拟预编码和数字预编码，Ns为数据流的数目，Fopt＝[u1,u2,...,uNs]是最优全数字预编码，[u1,u2,...,uNs]是毫米波信道矩阵H的Ns个最大的右特征向量，|·|F表示Frobineous范数，[A]i,j表示矩阵A的第(i,j)个元素，|·|表示取模操作。
本发明实施例采取的技术方案还包括：在所述步骤b中，所述按照最小均方原则得到数字预编码具体为：
给定初始模拟预编码FRF，所述混合预编码问题模型简化为：
Figure BDA0001315461150000053
Figure BDA0001315461150000054
得到：
Figure BDA0001315461150000055
上述公式中，
Figure BDA0001315461150000056
本发明实施例采取的技术方案还包括：在所述步骤c中，所述根据初始数字预编码对混合预编码问题模型进行优化具体为：给定初始数字预编码FBB，则所述混合预编码问题模型改写为：
Figure BDA0001315461150000061
Figure BDA0001315461150000062
将FRF分割为
Figure BDA0001315461150000063
个列向量，
Figure BDA0001315461150000064
为发射段天线阵元数目，FBB分割为
Figure BDA0001315461150000065
个行向量以及Fopt分割为Nt个行向量，则上述公式改写为：
Figure BDA0001315461150000066
Figure BDA0001315461150000067
将
Figure BDA0001315461150000068
带入上述公式，则所述混合预编码问题模型重新建模为：
Figure BDA0001315461150000069
上述公式中，
Figure BDA00013154611500000610
Arg([FRF]n,l)∈[0,2π)，
Figure BDA00013154611500000611
是由矩阵中所有元素的相位标量组成的矩阵，[·]n,:表示矩阵的第n行向量。
本发明实施例采取的技术方案还包括：在所述步骤c中，所述采用搜索法得到模拟预编码具体包括：
步骤c1：将模拟预编码矩阵分割为向量；
步骤c2：设定模拟预编码矩阵初值，并迭代计算搜索初值s(1)；
步骤c3：针对每个向量执行探测移动，直到没有可行的移动方向，并将所述搜索初值更新为
Figure BDA0001315461150000071
其中，k为一轮探测移动中执行的迭代次数；
步骤c4：对每个向量执行移动方向更新；
步骤c5：判断更新后的每个向量是否还有可行的移动方向，或相邻两次搜索初值是否满足
Figure BDA0001315461150000072
如果更新后的向量还有可行的移动方向，或相邻两次搜索初值不满足
Figure BDA0001315461150000073
则交替执行步骤c3和步骤c4，直到更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure BDA0001315461150000074
如果更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure BDA0001315461150000075
执行步骤c6；
步骤c6：输出模拟预编码：
Figure BDA0001315461150000076
本发明实施例采取的另一技术方案为：一种毫米波MIMO系统的混合预编码系统，包括：
模型建立模块：用于建立毫米波MIMO系统的混合预编码问题模型；
数字预编码计算模块：用于给定初始模拟预编码，基于所述初始模拟预编码，按照最小均方原则得到数字预编码；
模拟预编码计算模块：用于给定初始数字预编码，根据所述初始数字预编码对混合预编码问题模型进行优化，并采用搜索法得到模拟预编码；
迭代判断模块：用于判断混合预编码迭代次数是否达到设定的迭代次数值，如果没有达到设定的迭代次数值，通过数字预编码计算模块和模拟预编码计算模块迭代计算数字预编码和模拟预编码，直到达到混合预编码迭代次数值，并输出所述数字预编码和模拟预编码。
本发明实施例采取的技术方案还包括：所述混合预编码问题模型为：
Figure BDA0001315461150000081
Figure BDA0001315461150000082
上述公式中，FRF和FBB分别为模拟预编码和数字预编码，Ns为数据流的数目，Fopt＝[u1,u2,...,uNs]是最优全数字预编码，[u1,u2,...,uNs]是毫米波信道矩阵H的Ns个最大的右特征向量，|·|F表示Frobineous范数，[A]i,j表示矩阵A的第(i,j)个元素，|·|表示取模操作。
本发明实施例采取的技术方案还包括：所述数字预编码计算模块按照最小均方原则得到数字预编码具体为：
给定初始模拟预编码FRF，所述混合预编码问题模型简化为：
Figure BDA0001315461150000091
Figure BDA0001315461150000092
得到：
Figure BDA0001315461150000093
上述公式中，
Figure BDA0001315461150000094
本发明实施例采取的技术方案还包括：所述模拟预编码计算模块根据初始数字预编码对混合预编码问题模型进行优化具体为：给定初始数字预编码FBB，则所述混合预编码问题模型改写为：
Figure BDA0001315461150000095
Figure BDA0001315461150000096
将FRF分割为
Figure BDA0001315461150000097
个列向量，
Figure BDA0001315461150000098
为发射段天线阵元数目，FBB分割为
Figure BDA0001315461150000099
个行向量以及Fopt分割为Nt个行向量，则上述公式改写为：
Figure BDA00013154611500000910
Figure BDA00013154611500000911
将
Figure BDA00013154611500000912
带入上述公式，则所述混合预编码问题模型重新建模为：
Figure BDA00013154611500000913
上述公式中，
Figure BDA0001315461150000101
Arg([FRF]n,l)∈[0,2π)，
Figure BDA0001315461150000102
是由矩阵中所有元素的相位标量组成的矩阵，[·]n,:表示矩阵的第n行向量。
本发明实施例采取的技术方案还包括：所述模拟预编码计算模块包括：
向量分割单元：用于将模拟预编码矩阵分割为向量；
初值设定单元：用于设定模拟预编码矩阵初值，并迭代计算搜索初值s(1)；
探测移动执行单元：用于针对每个向量执行探测移动，直到没有可行的移动方向，并将所述搜索初值更新为
Figure BDA0001315461150000103
其中，k为一轮探测移动中执行的迭代次数；
移动方向更新单元：用于对每个向量执行移动方向更新；
迭代判断单元：用于判断更新后的每个向量是否还有可行的移动方向，或相邻两次搜索初值是否满足
Figure BDA0001315461150000104
如果更新后的向量还有可行的移动方向，或相邻两次搜索初值不满足
Figure BDA0001315461150000105
则通过探测移动执行单元和移动方向更新单元交替执行，直到更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure BDA0001315461150000106
如果更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure BDA0001315461150000107
执行步骤c6；
模拟预编码输出单元：用于输出模拟预编码：
Figure BDA0001315461150000108
相对于现有技术，本发明实施例产生的有益效果在于：本发明实施例的毫米波MIMO系统的混合预编码方法及系统采用非导数矩阵分解技术为毫米波MIMO系统建立一种超高天线阵列增益的快速波束形成方法。非导数技术降低了求导操作带来的高计算复杂度，同时大大减少了运算时间。此外，直接搜索算法能够灵活调整搜索步长，避免基误匹配现象，提高成形波束的增益。同时，本发明通过利用矩阵元素之间的关系给出合理的搜索初值，大大减少了搜索算法达到收敛所需的迭代次数，大幅降低了算法复杂度和缩短了执行时间。
附图说明
图1是本发明实施例的毫米波MIMO系统的混合预编码方法的流程图；
图2是本发明实施例的采用Rosenbrock搜索法得到待求模拟预编码的方法的流程图；
图3是本发明实施例的毫米波MIMO系统的混合预编码系统的结构示意图；
图4(a)至图4(c)是本发明实施例的仿真结果对比示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
本发明实施例的毫米波MIMO系统的混合预编码方法及系统针对现有技术存在的问题，通过将配置常幅度移相器的毫米波MIMO((多输入多输出))系统混合波束成形问题建模为常模受限制的矩阵分解问题，基于依次交替方法的思想提供一种基于Rosenbrock搜索算法的非导数矩阵分解方法。一方面，给定初始模拟预编码，则根据最小均方原则得到最终的数字预编码；另一方面，给定初始数字预编码，设定待求模拟预编码的经验初值和算法参数，则采用Rosenbrock搜索算法得到最终的模拟预编码：首先将待求模拟预编码矩阵分割为向量，然后给出待求模拟预编码矩阵初值，再针对每个向量执行探测移动和移动方向更新，模拟预编码和数字预编码交替执行，直到达到预设的最大迭代次数，并输出最终的模拟预编码和数字预编码。
具体地，请参阅图1，是本发明实施例的毫米波MIMO系统的混合预编码方法的流程图。本发明实施例的毫米波MIMO系统的混合预编码方法包括以下步骤：
步骤100：建立毫米波MIMO系统的混合预编码问题模型；
在步骤100中，假设毫米波MIMO系统采用最优接收机，则混合预编码问题模型P建模为：
Figure BDA0001315461150000131
Figure BDA0001315461150000132
公式(1)中，FRF和FBB分别为模拟预编码和数字预编码，Ns为数据流的数目，
Figure BDA0001315461150000133
是最优全数字预编码，[u1,u2,...,uNs]是毫米波信道矩阵H的Ns个最大的右特征向量，|·|F表示Frobineous范数，[A]i,j表示矩阵A的第(i,j)个元素，|·|表示取模操作。
步骤200：给定初始模拟预编码，并基于初始模拟预编码，按照最小均方原则得到待求数字预编码；
在步骤200中，待求数字预编码由以下过程推导而得：
给定初始模拟预编码FRF，则公式(1)可以简化为：
Figure BDA0001315461150000134
Figure BDA0001315461150000135
得到：
Figure BDA0001315461150000136
公式(3)中，
Figure BDA0001315461150000137
步骤300：给定初始数字预编码，根据初始数字预编码对混合预编码问题模型进行重建和优化，并采用Rosenbrock搜索法得到待求模拟预编码；
如图2所示，采用Rosenbrock搜索法得到待求模拟预编码的方式具体包括以下步骤：
步骤301：将待求模拟预编码矩阵分割为向量；
在步骤301中，混合预编码问题模型简化具体为：
给定初始数字预编码FBB，则公式(1)可改写成：
Figure BDA0001315461150000141
Figure BDA0001315461150000142
将FRF分割为
Figure BDA0001315461150000143
个列向量，(
Figure BDA0001315461150000144
为发射段天线阵元数目)，FBB分割为
Figure BDA0001315461150000145
个行向量以及Fopt分割为Nt个行向量，则公式(4)中的混合预编码问题模型可以改写为：
Figure BDA0001315461150000146
Figure BDA0001315461150000147
不失一般性，考虑
Figure BDA0001315461150000148
将其带入公式(5)，则优化的混合预编码问题模型可以重新建模为：
Figure BDA0001315461150000149
公式(6)中，
Figure BDA00013154611500001410
Arg([FRF]n,l)∈[0,2π)，
Figure BDA00013154611500001411
是由矩阵中所有元素的相位标量组成的矩阵，[·]n,:表示矩阵的第n行向量。
步骤302：设定模拟预编码矩阵初值；
在步骤302中，Rosenbrock算法的初值设计过程如下：
不失一般性，令
Figure BDA0001315461150000151
和
Figure BDA0001315461150000152
α为实数，并带入
Figure BDA0001315461150000153
得到模拟预编码矩阵各元素之间的关系：
Figure BDA0001315461150000154
因此，设定模拟预编码矩阵的初值：
Figure BDA0001315461150000155
根据公式(7)中原则更新模拟预编码的所有元素，即可迭代出搜索初值s(1)。
步骤303：针对每个向量执行探测移动，直到没有可行的移动方向，并将搜索初值s(1)更新为
Figure BDA0001315461150000156
其中，k为一轮探测移动中执行的迭代次数；
在步骤303中，探测移动沿着
Figure BDA0001315461150000157
个正交方向依次执行，更新搜索初值s(1)，具体过程如下：
1)给定搜索初值
Figure BDA0001315461150000158
n＝1,...,Nt，放大因子μ＞1，收缩因子ν∈(-1,0)，初始正交搜索方向
Figure BDA0001315461150000159
Figure BDA0001315461150000161
初始移动步长：
Figure BDA0001315461150000162
其中
Figure BDA0001315461150000163
表示正实数。
2)首先沿d1方向探测移动，如果
Figure BDA0001315461150000164
则该移动方向被选中，令
Figure BDA0001315461150000165
Figure BDA0001315461150000166
并且更新探测步长：ξ1＝μξ1；如果
Figure BDA0001315461150000167
则该移动方向落选，令
Figure BDA0001315461150000168
Figure BDA0001315461150000169
探测步长更新为：ξ1＝νξ1。执行完沿d1方向的探测移动后，再依次执行
Figure BDA00013154611500001610
方向的探测移动。当前所有一轮探测移动结束后得到
Figure BDA00013154611500001611
然后将搜索初值更新为：
Figure BDA00013154611500001612
并开始新一轮沿当前可行的移动方向执行探测移动，直到所有移动方向都落选，接着将搜索初值更新为：
Figure BDA00013154611500001613
步骤304：对每个向量执行移动方向更新；
在步骤304中，当所有移动方向都落选时，更新移动方向，具体过程如下：
一轮k次迭代的探测移动完成后，可得：
Figure BDA0001315461150000171
公式(8)中，λi表示沿di方向累积的移动步长。进一步，公式(8)可以变换为：
Figure BDA0001315461150000172
可得：P＝s(k+1)-s(k)是最速下降方向。更新移动方向适应充分考虑该方向，因此，新的正交搜索方向定义为：
Figure BDA0001315461150000173
进一步采用施密特正交化过程得到正交搜索方向：
Figure BDA0001315461150000174
并进行归一化，得到正交搜索方向为：
Figure BDA0001315461150000175
然后，令
Figure BDA0001315461150000176
沿更新后的移动方向开始下一轮探测移动。
步骤305：判断更新后的向量是否还有可行的移动方向，或相邻两次搜索初值是否不低于门限值；如果更新后的向量还有可行的移动方向，或相邻两次搜索初值不低于门限值，则交替执行步骤303和步骤304，直到更新后的向量没有可行的移动方向，或相邻两次搜索初值低于门限值；如果更新后的向量没有可行的移动方向，或相邻两次搜索初值低于门限值，执行步骤306；
在步骤305中，探测移动与移动方向更新操作依次交替执行，直到更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure BDA0001315461150000181
步骤306：输出待求模拟预编码：
Figure BDA0001315461150000182
步骤400：判断混合预编码迭代次数是否达到设定的迭代次数值，如果没有达到设定的迭代次数值，依次交替执行步骤200和步骤300，直到达到设定的迭代次数值，并得到最终的FBB和
Figure BDA0001315461150000183
作为输出的数字预编码和模拟预编码；如果达到设定的迭代次数值，执行步骤500。
步骤500：评价本发明实施例的混合预编码在系统频谱效率方面的性能表现。
在步骤500中，采用的频谱效率评价标准定义为：
Figure BDA0001315461150000184
公式(12)中，ρ和
Figure BDA0001315461150000185
分别为信号和噪声功率。
请参阅图3，是本发明实施例的毫米波MIMO系统的混合预编码系统的结构示意图。本发明实施例的毫米波MIMO系统的混合预编码系统包括模型建立模块、数字预编码计算模块、模拟预编码计算模块、迭代判断模块和性能评价模块。
模型建立模块：用于建立毫米波MIMO系统的混合预编码问题模型；其中，假设毫米波MIMO系统采用最优接收机，则混合预编码问题模型P建模为：
Figure BDA0001315461150000191
Figure BDA0001315461150000192
公式(1)中，FRF和FBB分别为模拟预编码和数字预编码，Ns为数据流的数目，
Figure BDA0001315461150000193
是最优全数字预编码，
Figure BDA0001315461150000194
是毫米波信道矩阵H的Ns个最大的右特征向量，|·|F表示Frobineous范数，[A]i,j表示矩阵A的第(i,j)个元素，|·|表示取模操作。
数字预编码计算模块：用于给定初始模拟预编码，并基于初始模拟预编码，按照最小均方原则得到待求数字预编码；其中，待求数字预编码由以下过程推导而得：
给定初始模拟预编码FRF，则公式(1)可以简化为：
Figure BDA0001315461150000195
Figure BDA0001315461150000196
得到：
Figure BDA0001315461150000201
公式(3)中，
Figure BDA0001315461150000202
模拟预编码计算模块：用于给定初始数字预编码，根据初始数字预编码对混合预编码问题模型进行重建和优化，并采用Rosenbrock搜索法得到待求模拟预编码；具体地，模拟预编码计算模块包括向量分割单元、初值设定单元、探测移动执行单元、移动方向更新单元、迭代判断单元和模拟预编码输出单元；
向量分割单元：用于将待求模拟预编码矩阵分割为向量；混合预编码问题模型简化具体为：
给定初始数字预编码FBB，则公式(1)可改写成：
Figure BDA0001315461150000203
Figure BDA0001315461150000204
将FRF分割为
Figure BDA0001315461150000205
个列向量，(
Figure BDA0001315461150000206
为发射段天线阵元数目)，FBB分割为
Figure BDA0001315461150000207
个行向量以及Fopt分割为Nt个行向量，则公式(4)中的混合预编码问题模型可以改写为：
Figure BDA0001315461150000208
Figure BDA0001315461150000209
不失一般性，考虑
Figure BDA0001315461150000211
将其带入公式(5)，则优化的混合预编码问题模型可以重新建模为：
Figure BDA0001315461150000212
公式(6)中，
Figure BDA0001315461150000213
Arg([FRF]n,l)∈[0,2π)，
Figure BDA0001315461150000214
是由矩阵中所有元素的相位标量组成的矩阵，[·]n,:表示矩阵的第n行向量。
初值设定单元：用于设定模拟预编码矩阵初值；其中，Rosenbrock算法的初值设计过程如下：
不失一般性，令
Figure BDA0001315461150000215
和
Figure BDA0001315461150000216
α为实数，并带入
Figure BDA0001315461150000217
得到模拟预编码矩阵各元素之间的关系：
Figure BDA0001315461150000218
因此，设定模拟预编码矩阵的初值：
Figure BDA0001315461150000219
根据公式(7)中原则更新模拟预编码的所有元素，即可迭代出搜索初值s(1)。
探测移动执行单元：用于针对每个向量执行探测移动，直到没有可行的移动方向，并将搜索初值s(1)更新为
Figure BDA0001315461150000221
其中，k为一轮探测移动中执行的迭代次数；探测移动沿着
Figure BDA0001315461150000222
个正交方向依次执行，更新搜索初值s(1)，具体过程如下：
1)给定搜索初值
Figure BDA0001315461150000223
n＝1,...,Nt，放大因子μ＞1，收缩因子ν∈(-1,0)，初始正交搜索方向
Figure BDA0001315461150000224
Figure BDA0001315461150000225
初始移动步长：
Figure BDA0001315461150000226
其中
Figure BDA0001315461150000227
表示正实数。
2)首先沿d1方向探测移动，如果
Figure BDA0001315461150000228
则该移动方向被选中，令
Figure BDA0001315461150000229
Figure BDA00013154611500002210
并且更新探测步长：ξ1＝μξ1；如果
Figure BDA00013154611500002211
则该移动方向落选，令
Figure BDA00013154611500002212
Figure BDA00013154611500002213
探测步长更新为：ξ1＝νξ1。执行完沿d1方向的探测移动后，再依次执行
Figure BDA0001315461150000231
方向的探测移动。当前所有一轮探测移动结束后得到
Figure BDA0001315461150000232
然后将搜索初值更新为：
Figure BDA0001315461150000233
并开始新一轮沿当前可行的移动方向执行探测移动，直到所有移动方向都落选，接着将搜索初值更新为：
Figure BDA0001315461150000234
移动方向更新单元：用于对每个向量执行移动方向更新；具体过程如下：
一轮k次迭代的探测移动完成后，可得：
Figure BDA0001315461150000235
公式(8)中，λi表示沿di方向累积的移动步长。进一步，公式(8)可以变换为：
Figure BDA0001315461150000236
可得：P＝s(k+1)-s(k)是最速下降方向。更新移动方向适应充分考虑该方向，因此，新的正交搜索方向定义为：
Figure BDA0001315461150000237
进一步采用施密特正交化过程得到正交搜索方向：
Figure BDA0001315461150000238
并进行归一化，得到正交搜索方向为：
Figure BDA0001315461150000241
然后，令
Figure BDA0001315461150000242
沿更新后的移动方向开始下一轮探测移动。
迭代判断单元：用于判断更新后的向量是否还有可行的移动方向，或相邻两次搜索初值是否不低于门限值；如果更新后的向量还有可行的移动方向，或相邻两次搜索初值不低于门限值，则通过探测移动执行单元和移动方向更新单元交替执行，直到更新后的向量没有可行的移动方向，或或相邻两次搜索初值满足
Figure BDA0001315461150000243
如果更新后的向量没有可行的移动方向，或相邻两次搜索初值低于门限值，则通过模拟预编码输出单元输出待求模拟预编码；
模拟预编码输出单元：用于输出待求模拟预编码：
Figure BDA0001315461150000244
迭代判断模块：用于判断混合预编码迭代次数是否达到设定的迭代次数值，如果没有达到设定的迭代次数值，通过数字预编码计算模块和模拟预编码计算模块迭代计算数字预编码和模拟预编码，直到达到设定的迭代次数值，并得到最终的FBB和
Figure BDA0001315461150000245
Figure BDA0001315461150000246
作为输出的数字预编码和模拟预编码；如果达到设定的迭代次数值，通过性能评价模块评价本发明实施例的混合预编码在系统频谱效率方面的性能表现。
性能评价模块：用于评价本发明实施例的混合预编码在系统频谱效率方面的性能表现。其中，采用的频谱效率评价标准定义为：
Figure BDA0001315461150000251
公式(12)中，ρ和
Figure BDA0001315461150000252
分别为信号和噪声功率。
请参阅图4(a)至图4(c)，为本发明实施例的仿真结果对比示意图。经在MATLAB平台进行仿真验证，假设数据流的数目Ns与射频链的数目NRF满足Ns≤NRF≤2Ns，角度扩散满足均值7.5°的拉普拉斯分布，信噪比SNR定义为
Figure BDA0001315461150000253
从仿真结果可以得出结论，本发明实施例的频谱效率表现随着信噪比的增加而增加，随着角度扩散增加而缓慢减小，随着射频链的数目增加而增加，随着数据流数目的增加而增加。相较于现有主要相关算法，本发明表现出频谱效率性能的优越性。
本发明实施例的毫米波MIMO系统的混合预编码方法及系统采用非导数矩阵分解技术为毫米波MIMO系统建立一种超高天线阵列增益的快速波束形成方法。非导数技术降低了求导操作带来的高计算复杂度，同时大大减少了运算时间。此外，直接搜索算法能够灵活调整搜索步长，避免基误匹配现象，提高成形波束的增益。同时，本发明通过利用矩阵元素之间的关系给出合理的搜索初值，大大减少了搜索算法达到收敛所需的迭代次数，大幅降低了算法复杂度和缩短了执行时间。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (6)

1.一种毫米波MIMO系统的混合预编码方法，其特征在于，包括以下步骤：
步骤a：建立毫米波MIMO系统的混合预编码问题模型；
步骤b：给定初始模拟预编码，基于所述初始模拟预编码，按照最小均方原则得到数字预编码；
步骤c：给定初始数字预编码，根据所述初始数字预编码对混合预编码问题模型进行优化，并采用搜索法得到模拟预编码；
步骤d：迭代执行步骤b和步骤c，直到达到混合预编码迭代次数值，并输出所述数字预编码和模拟预编码；
在所述步骤c中，根据所述初始数字预编码对混合预编码问题模型进行优化具体为：给定初始数字预编码FBB，则所述混合预编码问题模型改写为：
Figure FDA0002473461910000011
Figure FDA0002473461910000012
将FRF分割为
Figure FDA0002473461910000013
个列向量，
Figure FDA0002473461910000014
为发射段天线阵元数目，FBB分割为
Figure FDA0002473461910000015
个行向量以及Fopt分割为Nt个行向量，则上述公式改写为：
Figure FDA0002473461910000016
Figure FDA0002473461910000017
将
Figure FDA0002473461910000021
带入上述公式，则所述混合预编码问题模型重新建模为：
Figure FDA0002473461910000022
上述公式中，
Figure FDA0002473461910000023
Figure FDA0002473461910000024
是由矩阵中所有元素的相位标量组成的矩阵，[·]n,:表示矩阵的第n行向量；
在所述步骤c中，所述采用搜索法得到模拟预编码具体包括：
步骤c1：将模拟预编码矩阵分割为向量；
步骤c2：设定模拟预编码矩阵初值，并迭代计算搜索初值s(1)；
步骤c3：针对每个向量执行探测移动，直到没有可行的移动方向，并将所述搜索初值更新为
Figure FDA0002473461910000025
其中，k为一轮探测移动中执行的迭代次数；
步骤c4：对每个向量执行移动方向更新；
步骤c5：判断更新后的每个向量是否还有可行的移动方向，或相邻两次搜索初值是否满足
Figure FDA0002473461910000026
如果更新后的向量还有可行的移动方向，或相邻两次搜索初值不满足
Figure FDA0002473461910000027
则交替执行步骤c3和步骤c4，直到更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure FDA0002473461910000028
如果更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure FDA0002473461910000029
执行步骤c6；
步骤c6：输出模拟预编码：
Figure FDA0002473461910000031
2.根据权利要求1所述的毫米波MIMO系统的混合预编码方法，其特征在于，在所述步骤a中，所述混合预编码问题模型为：
P:
Figure FDA0002473461910000032
Figure FDA0002473461910000033
上述公式中，FRF和FBB分别为模拟预编码和数字预编码，Ns为数据流的数目，
Figure FDA0002473461910000034
是最优全数字预编码，
Figure FDA0002473461910000035
是毫米波信道矩阵H的Ns个最大的右特征向量，|·|F表示Frobineous范数，[A]i,j表示矩阵A的第(i,j)个元素，|·|表示取模操作。
3.根据权利要求2所述的毫米波MIMO系统的混合预编码方法，其特征在于，在所述步骤b中，所述按照最小均方原则得到数字预编码具体为：
给定初始模拟预编码FRF，所述混合预编码问题模型简化为：
Figure FDA0002473461910000036
Figure FDA0002473461910000037
得到：
Figure FDA0002473461910000038
上述公式中，
Figure FDA0002473461910000039
4.一种毫米波MIMO系统的混合预编码系统，其特征在于，包括：
模型建立模块：用于建立毫米波MIMO系统的混合预编码问题模型；
数字预编码计算模块：用于给定初始模拟预编码，基于所述初始模拟预编码，按照最小均方原则得到数字预编码；
模拟预编码计算模块：用于给定初始数字预编码，根据所述初始数字预编码对混合预编码问题模型进行优化，并采用搜索法得到模拟预编码；
迭代判断模块：用于判断混合预编码迭代次数是否达到设定的迭代次数值，如果没有达到设定的迭代次数值，通过数字预编码计算模块和模拟预编码计算模块迭代计算数字预编码和模拟预编码，直到达到混合预编码迭代次数值，并输出所述数字预编码和模拟预编码；
所述模拟预编码计算模块根据初始数字预编码对混合预编码问题模型进行优化具体为：给定初始数字预编码FBB，则所述混合预编码问题模型改写为：
Figure FDA0002473461910000041
Figure FDA0002473461910000042
将FRF分割为
Figure FDA0002473461910000043
个列向量，
Figure FDA0002473461910000044
为发射段天线阵元数目，FBB分割为
Figure FDA0002473461910000045
个行向量以及Fopt分割为Nt个行向量，则上述公式改写为：
Figure FDA0002473461910000051
Figure FDA0002473461910000052
将
Figure FDA0002473461910000053
带入上述公式，则所述混合预编码问题模型重新建模为：
Figure FDA0002473461910000054
上述公式中，
Figure FDA0002473461910000055
Figure FDA0002473461910000056
是由矩阵中所有元素的相位标量组成的矩阵，[·]n,:表示矩阵的第n行向量；
所述模拟预编码计算模块包括：
向量分割单元：用于将模拟预编码矩阵分割为向量；
初值设定单元：用于设定模拟预编码矩阵初值，并迭代计算搜索初值s(1)；
探测移动执行单元：用于针对每个向量执行探测移动，直到没有可行的移动方向，并将所述搜索初值更新为
Figure FDA0002473461910000057
其中，k为一轮探测移动中执行的迭代次数；
移动方向更新单元：用于对每个向量执行移动方向更新；
迭代判断单元：用于判断更新后的每个向量是否还有可行的移动方向，或相邻两次搜索初值是否满足
Figure FDA0002473461910000058
如果更新后的向量还有可行的移动方向，或相邻两次搜索初值不满足
Figure FDA0002473461910000059
则通过探测移动执行单元和移动方向更新单元交替执行，直到更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure FDA0002473461910000061
模拟预编码输出单元：用于当更新后的向量没有可行的移动方向，或相邻两次搜索初值满足
Figure FDA0002473461910000062
时，输出模拟预编码：
Figure FDA0002473461910000063
5.根据权利要求4所述的毫米波MIMO系统的混合预编码系统，其特征在于，所述混合预编码问题模型为：
P:
Figure FDA0002473461910000064
Figure FDA0002473461910000065
上述公式中，FRF和FBB分别为模拟预编码和数字预编码，Ns为数据流的数目，
Figure FDA0002473461910000066
是最优全数字预编码，
Figure FDA0002473461910000067
是毫米波信道矩阵H的Ns个最大的右特征向量，|·|F表示Frobineous范数，[A]i,j表示矩阵A的第(i,j)个元素，|·|表示取模操作。
6.根据权利要求5所述的毫米波MIMO系统的混合预编码系统，其特征在于，所述数字预编码计算模块按照最小均方原则得到数字预编码具体为：
给定初始模拟预编码FRF，所述混合预编码问题模型简化为：
Figure FDA0002473461910000068
Figure FDA0002473461910000069
得到：
Figure FDA0002473461910000071
上述公式中，
Figure FDA0002473461910000072
CN201710422927.8A 2017-06-07 2017-06-07 一种毫米波mimo系统的混合预编码方法及系统 Active CN107181511B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710422927.8A CN107181511B (zh) 2017-06-07 2017-06-07 一种毫米波mimo系统的混合预编码方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710422927.8A CN107181511B (zh) 2017-06-07 2017-06-07 一种毫米波mimo系统的混合预编码方法及系统

Publications (2)

Publication Number Publication Date
CN107181511A CN107181511A (zh) 2017-09-19
CN107181511B true CN107181511B (zh) 2020-07-24

Family

ID=59835859

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710422927.8A Active CN107181511B (zh) 2017-06-07 2017-06-07 一种毫米波mimo系统的混合预编码方法及系统

Country Status (1)

Country Link
CN (1) CN107181511B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107566017B (zh) * 2017-10-24 2020-04-17 电子科技大学 用于毫米波mimo系统的预编码方法
CN107809275B (zh) * 2017-11-30 2020-06-23 杭州电子科技大学 一种基于毫米波mimo系统的有限反馈混合预编码方法
CN108173584B (zh) * 2017-12-13 2020-10-23 电子科技大学 一种fdd模式下数模混合低反馈量的多用户波束成形方法
CN108599825A (zh) * 2018-02-12 2018-09-28 大连理工大学 一种基于mimo-ofdm毫米波结构的混合编码方法
CN108494455B (zh) * 2018-02-27 2021-04-30 同济大学 采用单比特模数转换器mimo系统的混合预编码设计方法
CN108847874B (zh) * 2018-06-15 2020-12-18 电子科技大学 一种用于毫米波rsm-mimo系统的数模混合预编码方法
CN109302215B (zh) * 2018-09-18 2020-12-01 北京邮电大学 一种基于行向量优化的混合预编码方法
CN110071751B (zh) * 2019-03-19 2021-10-26 西安电子科技大学 模拟域波束赋形参数非恒模的部分连接混合波束赋形方法
CN110611526A (zh) * 2019-09-20 2019-12-24 中国科学院上海微系统与信息技术研究所 基于改进黎曼流形优化的毫米波混合模/数波束成形方法
CN111726144B (zh) * 2020-06-24 2021-04-16 中南大学 基于初值优化的混合预编码设计方法、装置、介质及设备
CN113131981B (zh) * 2021-03-23 2022-08-26 湖南大学 一种混合波束成形方法、装置及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104486044A (zh) * 2014-12-30 2015-04-01 北京航空航天大学 一种大规模mimo系统中的宽带模数混合预处理方法
WO2016168128A1 (en) * 2015-04-15 2016-10-20 Ping Liang Hybrid beamforming multi-antenna wireless systems
CN106160809A (zh) * 2015-04-10 2016-11-23 上海贝尔股份有限公司 多用户多天线系统的混合预编码方法及其装置
CN106253956A (zh) * 2016-08-24 2016-12-21 东南大学 基于码本的模数混合预编码方法
CN106330273A (zh) * 2015-07-03 2017-01-11 电信科学技术研究院 一种预编码方法及装置
CN106506051A (zh) * 2015-09-08 2017-03-15 上海贝尔股份有限公司 基于可重构天线的混合预编码的方法与设备
CN106571858A (zh) * 2016-11-03 2017-04-19 北京邮电大学 一种混合波束成形传输系统和方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9444534B2 (en) * 2012-02-06 2016-09-13 Samsung Electronics Co., Ltd. Apparatus and method for low complexity spatial division multiple access in a millimeter wave mobile communication system
US20150358063A1 (en) * 2014-06-09 2015-12-10 Collision Communications, Inc. Reduction Of Interference In A Wireless Network Through Spectrum Analysis And Transmit Beamforming

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104486044A (zh) * 2014-12-30 2015-04-01 北京航空航天大学 一种大规模mimo系统中的宽带模数混合预处理方法
CN106160809A (zh) * 2015-04-10 2016-11-23 上海贝尔股份有限公司 多用户多天线系统的混合预编码方法及其装置
WO2016168128A1 (en) * 2015-04-15 2016-10-20 Ping Liang Hybrid beamforming multi-antenna wireless systems
CN106330273A (zh) * 2015-07-03 2017-01-11 电信科学技术研究院 一种预编码方法及装置
CN106506051A (zh) * 2015-09-08 2017-03-15 上海贝尔股份有限公司 基于可重构天线的混合预编码的方法与设备
CN106253956A (zh) * 2016-08-24 2016-12-21 东南大学 基于码本的模数混合预编码方法
CN106571858A (zh) * 2016-11-03 2017-04-19 北京邮电大学 一种混合波束成形传输系统和方法

Also Published As

Publication number Publication date
CN107181511A (zh) 2017-09-19

Similar Documents

Publication Publication Date Title
CN107181511B (zh) 一种毫米波mimo系统的混合预编码方法及系统
Gao et al. Machine learning inspired energy-efficient hybrid precoding for mmWave massive MIMO systems
Hu et al. Statistical CSI based design for intelligent reflecting surface assisted MISO systems
CN104935366B (zh) 毫米波通信中的波束搜索方法
Song et al. Codebook design for hybrid beamforming in millimeter wave systems
Hosoya et al. Multiple sector ID capture (MIDC): A novel beamforming technique for 60-GHz band multi-Gbps WLAN/PAN systems
CN101867402B (zh) 一种自适应天线选择的mimo系统及其应用方法
CN113315547A (zh) 智能反射面辅助多小区的鲁棒联合传输波束优化方法
CN108650200B (zh) 高低频混合组网系统的低频辅助信道估计方法
CN110011712B (zh) 一种面向毫米波大规模多输入多输出的混合预编码方法
KR102228091B1 (ko) 밀리미터파 기반 대규모 mimo 시스템의 하이브리드 빔포밍 장치 및 방법
CN110365388B (zh) 一种低复杂度毫米波多播波束成形方法
CN113612508B (zh) 基于机器学习的irs辅助毫米波通信波束成形设计方法
CN102801455A (zh) 波束码本生成方法、波束搜索方法及相关装置
CN114268350A (zh) 智能反射面辅助通信中非理想信道下的鲁棒波束赋形设计方法
KR102266761B1 (ko) 빔형성 방법 및 디바이스
CN110611526A (zh) 基于改进黎曼流形优化的毫米波混合模/数波束成形方法
Ren et al. Machine learning-based hybrid precoding with robust error for UAV mmWave massive MIMO
CN109495150A (zh) 一种毫米波通信码本设计方法
Alluhaibi et al. Hybrid digital-to-analog beamforming approaches to maximise the capacity of mm-wave systems
CN110719127B (zh) 一种具有恒模约束的毫米波mimo系统波束成形方法
Abbas et al. Millimeter wave communications over relay networks
CN115021779A (zh) Ris辅助的mimo系统的离散相移设计方法和装置
CN114866126A (zh) 智能反射面辅助毫米波系统的低开销信道估计方法
Wang et al. New Environment Adaptation with Few Shots for OFDM Receiver and mmWave Beamforming

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