CN113839695A - Fdd大规模mimo和速率最优统计预编码方法及设备 - Google Patents

Fdd大规模mimo和速率最优统计预编码方法及设备 Download PDF

Info

Publication number
CN113839695A
CN113839695A CN202111087692.4A CN202111087692A CN113839695A CN 113839695 A CN113839695 A CN 113839695A CN 202111087692 A CN202111087692 A CN 202111087692A CN 113839695 A CN113839695 A CN 113839695A
Authority
CN
China
Prior art keywords
statistical
rate
precoding
matrix
traversal
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.)
Granted
Application number
CN202111087692.4A
Other languages
English (en)
Other versions
CN113839695B (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 CN202111087692.4A priority Critical patent/CN113839695B/zh
Publication of CN113839695A publication Critical patent/CN113839695A/zh
Application granted granted Critical
Publication of CN113839695B publication Critical patent/CN113839695B/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
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0204Channel estimation of multiple channels

Landscapes

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

Abstract

本发明公开了一种FDD大规模MIMO和速率最优统计预编码方法及设备。统计预编码较传统预编码方法可以有效地消除信道状态信息获取开销。基于波束域统计信道模型,将FDD预编码设计转化为总功率约束下的最大化遍历和速率优化问题。利用MM算法可以得到该优化问题的最优解。为了简化最优解中带有关于信道的期望项,引入确定性等同方法得到遍历和速率的近似,进而推导得到最优解的近似表达。相比于其它方法,本发明只需要统计信道信息就可以进行预编码设计,同时基于和速率最优的准则,可以有效地抑制用户间的干扰,从而提升整个通信系统的性能,具有较大的应用价值。

Description

FDD大规模MIMO和速率最优统计预编码方法及设备
技术领域
本发明属于通信技术领域,尤其涉及FDD大规模MIMO和速率最优统计预编码设计方法及设备。
背景技术
大规模多输入多输出(MIMO,Multiple-Input Multiple-Ouput)技术通过在基站侧配置大规模天线阵列,能显著地提升系统的吞吐量以及频谱效率,从而成为5G的关键技术之一。大规模MIMO在其它很多应用和服务上也已经有所应用,其中包括物联网以及无人机通信。基站利用信道信息(CSI)进行预编码设计,从而抑制用户间的干扰。目前大多数预编码相关研究把重点放在时分双工(TDD)系统上,利用上下行信道的互异性,基站可以直接从上行信道估计得到下行的信道信息,进而直接进行预编码设计。
频分双工(FDD)系统也在无线网络中起到重要作用,尤其在对时延敏感的系统中,FDD比TDD更具优势。同时,有较多可用的频率带宽可供FDD系统使用。由于FDD系统中不存在信道互异性,大多数研究都是基于上行反馈估计值的方式来获取信道信息然后进行预编码设计。由于反馈需要占用时频资源,这种方法会带来较大的开销。基于统计信道信息的预编码设计较好地解决了这个问题,利用上下行信道之间的统计互异性,不再需要上行反馈就可以进行预编码设计。目前有关大规模MIMO统计预编码的研究中,多数通过优化遍历和速率的上下界或其它准则如信号干扰加噪声比(SLNR)来设计最优预编码矩阵。这些准则在一些信噪比区域不能很好地逼近遍历和速率,所以会导致和速率性能的损失。目前为止,尚无直接利用遍历和速率以及统计信道信息进行预编码设计的研究。
发明内容
本发明目的在于提供一种FDD大规模MIMO和速率最优统计预编码设计方法及设备,以解决为FDD大规模MIMO预编码减轻信道信息获取开销的同时带来性能提升的技术问题。
为解决上述技术问题,本发明的具体技术方案如下:
一种FDD大规模MIMO和速率最优统计预编码设计方法,包括如下步骤:
步骤1、基站利用上行信道估计以及统计互异性获得各用户终端的下行统计信道信息;上行统计信道信息和下行统计信道信息是一致的;所述统计信道信息从波束域统计信道模型中提取得到,该波束域统计信道模型的导向矢量数量超过了天线数;上下行信道可以通过空间采样矩阵进行相互转换。
步骤2、获得遍历和速率的基于统计信道信息的确定性等同表达式;
步骤3、基站利用遍历和速率的确定性等同进行FDD大规模MIMO和速率最优统计预编码设计。
进一步的,所述步骤1中通过对基站和用户侧的导向矢量过采样得到物理信道的模型统计信道模型,该统计信道模型从波束域角度描述空间相关性。
进一步的,所述下行统计信道信息由基站通过对上行信道统计平均并利用统计互异性获得,不需要下行信道估计以及上行反馈过程。
进一步的,所述步骤1中基站侧的上行信道与下行信道空间采样矩阵进行相互转化,包括以下步骤:
步骤1.1、基站侧上下行信道的导向矢量的计算和中心频率有关,通过定义的上下行中心频率商因子进行相互转换;其中导向矢量中的采样方向余弦通过均匀采样得到;
步骤1.2、堆叠导向矢量分别得到水平和垂直方向的块矩阵,基站侧空间采样矩阵等于两个块矩阵的克罗内克积;
步骤1.3、利用基站侧上下行导向矢量的关系进而得到基站侧上下行空间采样矩阵的关系。
进一步的,所述步骤1.1中采样方向余弦的采样数大于对应的天线数,使得堆叠得到空间采样矩阵维度大于天线数。
进一步的,在所述的和速率最优统计预编码设计中,基站根据遍历和速率最大化的准则,进行各用户终端的线性预编码矩阵设计,遍历和速率表达式中的干扰加噪声协方差矩阵等效为高斯噪声的协方差矩阵,同时定义干扰加噪声协方差矩阵的瞬时表达。
进一步的,所述步骤2中将所述遍历和速率最大化准则替换为遍历和速率下界函数最大化准则;或者将所述遍历和速率最大化准则中和速率替换为其确定性等同。
进一步的,所述步骤3具体包括以下步骤:
通过MM算法将加权遍历和速率最大化预编码设计问题转化为迭代求解二次型优化问题得到最优预编码矩阵的迭代更新式;
通过确定性等同对最优迭代更新式中带有的关于随机矩阵的期望项作近似,得到最优预编码矩阵的闭式表达。
进一步的,所述方法结合MM算法以及确定性等同进行预编码设计,包括:
步骤3.1、初始化预编码矩阵并归一化使之满足功率条件;
步骤3.2、根据预编码矩阵以及统计信道信息计算遍历和速率确定性等同中的相关项;
步骤3.3、根据确定性等同的相关项计算预编码矩阵最优迭代更新式中带有的关于随机矩阵的期望项;
步骤3.4、更新预编码迭代更新的闭式表达,并归一化预编码矩阵;
重复步骤3.2到步骤3.4直到达到预设迭代次数或预编码收敛。
一种计算设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述计算机程序被加载至处理器时实现上述的FDD大规模MIMO和速率最优统计预编码设计方法。
本发明的FDD大规模MIMO和速率最优统计预编码设计方法及设备,具有以下优点:
本发明提出的FDD大规模MIMO和速率最优统计预编码设计方法及设备直接利用遍历和速率以及统计信道信息求解最优预编码矩阵,有效抑制用户间干扰的同时消除信道信息获取所带来的开销,为FDD系统下大规模MIMO预编码设计的实际应用提供支撑。
附图说明
图1为本发明的FDD大规模MIMO和速率最优统计预编码设计方法流程图;
图2(a)为用户1上行信道统计信息发布图;
图2(b)为用户1下行信道统计信息发布图;
图2(c)为用户2上行信道统计信息发布图;
图2(d)为用户2下行信道统计信息发布图;
图3为本发明的FDD大规模MIMO和速率最优统计预编码与BDMA预编码和速率性能对比图;
图4为本发明的FDD大规模MIMO和速率最优统计预编码与SLNR-SBF预编码以及ZF预编码和速率性能对比图。
具体实施方式
为了更好地了解本发明的目的、结构及功能,下面结合附图,对本发明一种FDD大规模MIMO和速率最优统计预编码设计方法及设备做进一步详细的描述。
如图1所示,本发明实施例公开的FDD大规模MIMO和速率最优统计预编码设计方法包括以下步骤:
步骤1、基站利用上行信道估计以及统计互异性获得各用户终端的下行统计信道信息;上行统计信道信息和下行统计信道信息是一致的;所述统计信道信息从波束域统计信道模型中提取得到,该波束域统计信道模型的导向矢量数量超过了天线数。
通过对基站和用户侧的导向矢量过采样得到物理信道的模型统计信道模型,该统计信道模型从波束域角度描述空间相关性。
所述下行统计信道信息由基站通过对上行信道统计平均并利用统计互异性获得,不需要下行信道估计以及上行反馈过程。
上下行信道可以通过空间采样矩阵进行相互转换,包括以下步骤:
步骤1.1、基站侧上下行信道的导向矢量的计算和中心频率有关,通过定义的上下行中心频率商因子进行相互转换;其中导向矢量中的采样方向余弦通过均匀采样得到。其中采样方向余弦的采样数大于对应的天线数,使得堆叠得到空间采样矩阵维度大于天线数。
步骤1.2、堆叠导向矢量分别得到水平和垂直方向的块矩阵,基站侧空间采样矩阵等于两个块矩阵的克罗内克积;
步骤1.3、利用基站侧上下行导向矢量的关系进而得到基站侧上下行空间采样矩阵的关系。
步骤2、获得遍历和速率的基于统计信道信息的确定性等同表达式;
在所述的和速率最优统计预编码设计中,基站根据遍历和速率最大化的准则,进行各用户终端的线性预编码矩阵设计,遍历和速率表达式中的干扰加噪声协方差矩阵等效为高斯噪声的协方差矩阵,同时定义干扰加噪声协方差矩阵的瞬时表达。
将所述遍历和速率最大化准则替换为遍历和速率下界函数最大化准则;或者将所述遍历和速率最大化准则中和速率替换为其确定性等同。
步骤3、基站利用遍历和速率的确定性等同进行FDD大规模MIMO和速率最优统计预编码设计,具体包括以下步骤:
通过MM算法将加权遍历和速率最大化预编码设计问题转化为迭代求解二次型优化问题得到最优预编码矩阵的迭代更新式;
通过确定性等同对最优迭代更新式中带有的关于随机矩阵的期望项作近似,得到最优预编码矩阵的闭式表达。
所述方法结合MM算法以及确定性等同进行预编码设计,包括:
步骤3.1、初始化预编码矩阵并归一化使之满足功率条件;
步骤3.2、根据预编码矩阵以及统计信道信息计算遍历和速率确定性等同中的相关项;
步骤3.3、根据确定性等同的相关项计算预编码矩阵最优迭代更新式中带有的关于随机矩阵的期望项;
步骤3.4、更新预编码迭代更新的闭式表达,并归一化预编码矩阵;
重复步骤3.2到步骤3.4直到达到预设迭代次数或预编码收敛。
一种计算设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述计算机程序被加载至处理器时实现上述的FDD大规模MIMO和速率最优统计预编码设计方法。
本发明方法主要适用于基站侧配备大规模天线阵列以同时服务多个用户的FDD大规模MIMO系统。下面结合具体的通信系统实例对本发明涉及和速率最优统计预编码设计方法的具体实现过程作详细说明,需要说明的是本发明方法不仅适用于下面示例所举的具体系统模型,也同样适用于其它配置的系统模型。
一、系统配置
考虑一个FDD大规模MIMO块平坦衰落系统模型,下行包括一个基站和K 个用户。基站配置大规模均匀面阵天线阵列(UPA),垂直方向的天线数量为 Mz,水平方向的天线数量为Mx,基站侧总天线数量为Mt=MzMx。简便起见,假设每个用户配置天线数为Mk的均匀线阵(ULA)。系统多普勒资源可以分为若干时隙,每一时隙包括Nb个多普勒块,而每个多普勒块则包含T个符号间隔。根据上下行信道中心频率的不同,可以将时隙结构分为两部分。对于上行部分,只在第一多普勒块传输上行导频信号,第2至Nb多普勒块则用于上行数据信号传输。对于下行部分,只在第一多普勒块传输下行导频信号,第2至Nb多普勒块则用于下行数据信号传输。
二、波束域统计信道模型
波束域统计信道模型在信道模型中引入比天线数量更多的空间导向矢量,从而更准确地描述信道统计特性。载波频率为f,光速为c。接收端线阵(ULA) 的天线间距为dr。为了表达方便,定义
Figure BDA0003266408670000071
ur=cosθr表示关于ULA的方向余弦,其中θr是信道路径在用户侧的角度。简便起见,定义(·)d表示下行信道对应的参数,(·)u上行信道对应的参数。用户侧导向矢量表示如下:
Figure BDA0003266408670000072
在基站侧,令方位角为φt,极角为θt。UPA放置在xz平面上。天线阵列的行之间的间隔为dx,天线阵列的列之间的间隔为dz。为了表达方便,定义
Figure BDA0003266408670000081
ut=cosθt表示关于UPA的水平天线阵列的方向余弦, vt=sinθtcosφt表示关于UPA的垂直天线阵列的方向余弦。基站侧导向矢量表示为如下克罗内克积形式
Figure BDA0003266408670000082
其中
Figure BDA0003266408670000083
Figure BDA0003266408670000084
定义上下行频率商因子为β=fd/fu,基站侧上下行导向矢量的关系可以表示为
Figure BDA0003266408670000085
Figure BDA0003266408670000086
其中pow(v,β)表示将v中每个元素进行β次幂运算,
Figure BDA0003266408670000087
Figure BDA0003266408670000088
为了得可用于实际大规模MIMO信道的信道模型,需要将ur,ut和vt离散化。令ur,i,ut,j和vt,l分别为ur,ut和vt的离散化采样方向余弦,Nk,Nz和Nx分别是其对应采样数。用户侧的空间采样矩阵可以表示为
Figure BDA0003266408670000089
定义变量Nt=NzNx,基站侧的下行空间采样矩阵表示为
Figure BDA00032664086700000810
Figure BDA00032664086700000811
其中
Figure BDA00032664086700000812
Figure BDA0003266408670000091
基站侧上行空间采样矩阵Vu可以通过转换式Vd=pow(Vu,β)得到基站侧下行空间采样矩阵Vd
基站到第k个用户的下行信道表示为Hkmn,其中m和n分别表示时隙和多普勒块的索引。简便起见,上标d被省略。信道矩阵可以作如下分解
Figure BDA0003266408670000092
其中
Figure BDA0003266408670000093
表示波束域信道,
Figure BDA0003266408670000094
是一个非负的确定矩阵,Wkmn为独立同分布复高斯随机矩阵。定义统计信道信息为波束域能量矩阵表示为
Ωk=Mk⊙Mk (13)
系数[Ωk]ij表示基站第j个特征向量与用户第i个特征向量的耦合平均能量。定义过采样因子
Figure BDA0003266408670000095
Figure BDA0003266408670000096
当Fk=Fz=Fx=1,矩阵U退化为酉矩阵,矩阵V可以表示为两个DFT矩阵的克罗内克积,信道模型退化为传统波束域信道。当Fk>1,Fz>1,Fx>1,由于导向矢量的数量超过了天线数,信道模型可以更精确地描述物理信道。下行的波束域能量矩阵可以通过利用上下行信道的统计互异性,从上行信道估计中得到。
下面给出预编码设计的传输模型。定义流数为dk
Figure BDA0003266408670000097
为维度是dk的单位矩阵。发送信号
Figure BDA0003266408670000098
满足
Figure BDA0003266408670000099
简便起见,忽略下标m和 n。第k个用户的接收信号可以表示为
Figure BDA00032664086700000910
其中
Figure BDA00032664086700000911
为第k个用户的归一化的预编码矩阵,基站到第k个用户的下行信道表示为Hk,zk是服从
Figure BDA00032664086700000912
的复高斯矢量,k表示用户的索引将干扰加噪声
Figure BDA00032664086700000913
看作高斯噪声,Rk表示干扰加噪声的协方差矩阵。瞬时的干扰加噪声的协方差矩阵表示如下
Figure BDA0003266408670000101
其中P\k=[P1 P2 … Pk-1 Pk+1 … PK]表示水平堆叠除第k个用户以外所有用户的预编码矩阵。(·)H表示矩阵的共轭转置。基于上面的定义与假设,用户速率可以表示为
Figure BDA0003266408670000102
其中P=[P1 P2 … PK]表示水平堆叠所有用户的预编码矩阵。
三、问题陈述
定义
Figure BDA0003266408670000103
为遍历和速率。将预编码设计问题转化为总功率约束下的最大化遍历和速率问题,表示如下
Figure BDA0003266408670000104
Figure BDA0003266408670000105
其中预编码矩阵的总功率约束为P。
在没有闭式解的情况下计算遍历速率比较复杂。从算子自由概率演化而来的著名分析方法被称为确定性等同,它能够通过推导收敛迭代公式来提供闭式的表达式。在下一节中,确定性等效方法将用于预编码设计。
四、统计预编码设计
4.1、基于MM算法的预编码设计
首先引入MM算法,在MM的当前迭代下,找到了一个与原始目标函数相切的下界函数。然后在下一次迭代过程中最大化下界函数。随着迭代次数的增加,原目标函数的值单调增加。算法的收敛性确保原优化问题可以取得稳定点。
当实值函数g(P1,P2,…,PK|P1[d],P2[d],…,PK[d])满足下面条件时,认为遍历和速率在第d次迭代时的下界函数
g(P1,P2,…,PK|P1[d],P2[d],…,PK[d])≤f(P1,P2,…,PK) (18)
g(P1[d],P2[d],…,PK[d]|P1[d],P2[d],…,PK[d])=f(P1[d],P2[d],…,PK[d]) (19)
Figure BDA0003266408670000111
其中
Figure BDA0003266408670000112
表示求导运算,(·)*表示矩阵的共轭。条件(20)可以使函数g在固定点与原目标函数相切,在得到下界函数之后,第d+1次迭代的预编码矩阵可以表示为
Figure BDA0003266408670000113
该最优解最终可以收敛到局部最优点。
定义总接收信号的协方差矩阵为
Figure BDA0003266408670000114
其中
Figure BDA0003266408670000115
为维度是Mk的单位矩阵下面的定理给出满足上面条件的一个下界函数。
定理1、定义函数g为
Figure BDA0003266408670000116
其中ck[d]为一个常量,
Figure BDA0003266408670000117
Figure BDA0003266408670000118
Figure BDA0003266408670000119
基于定理1中给出的下界函数,可以重新建模优化问题为
Figure BDA0003266408670000121
Figure BDA0003266408670000122
该问题是一个凹二次型优化问题,最优序列可以利用拉格朗日乘子法得到。拉格朗日方程可以表示为
Figure BDA0003266408670000123
进一步,可以得到优化问题的最优解为
Pk[d+1]=(D[d]+μ*I)-1Ak[d]Pk[d] (29)
观察上面的结构可以发现,Pk[d]先通过Ak[d]滤波之后再由(D[d]+ μ*I)-1滤波。Ak[d]包含了第k个用户的空间方向信息,(D[d]+μ*I)-1包含了干扰用户的空间方向信息。
4.2、基于确定性等同的预编码设计
在引入确定性等同之前,首先给出在下面推导会使用的定义。令
Figure BDA0003266408670000124
表示参数化信道协方差
Figure BDA0003266408670000125
进一步化简可以得到
Figure BDA0003266408670000126
令[·]ij表示矩阵中第i行第j列的元素,其中对角矩阵Λk计算如下
Figure BDA0003266408670000127
类似地,定义
Figure BDA0003266408670000128
Figure BDA0003266408670000129
其中对角矩阵
Figure BDA00032664086700001210
计算如下
Figure BDA00032664086700001211
根据上面定义的函数,定义
Figure BDA00032664086700001212
的确定性等同为
Figure BDA00032664086700001213
Figure BDA0003266408670000131
其中
Figure BDA0003266408670000132
Figure BDA0003266408670000133
计算如下
Figure BDA0003266408670000134
Figure BDA0003266408670000135
或者是
Figure BDA0003266408670000136
Figure BDA0003266408670000137
Figure BDA0003266408670000138
中的矩阵Γk
Figure BDA0003266408670000139
表示如下
Figure BDA00032664086700001310
Figure BDA00032664086700001311
以及Φk
Figure BDA00032664086700001312
表示如下
Figure BDA00032664086700001313
Figure BDA00032664086700001314
Figure BDA00032664086700001315
Figure BDA00032664086700001316
计算
Figure BDA00032664086700001317
的过程是类似的,将上面迭代式39到迭代式44中的P换成P\k得到Γ′k
Figure BDA00032664086700001318
Φ′k
Figure BDA00032664086700001319
Figure BDA00032664086700001326
以及
Figure BDA00032664086700001320
在下面的算法1中给出获取遍历和速率的确定性等同
Figure BDA00032664086700001321
的具体实现步骤。
步骤1、从算法2中获取预编码矩阵P1[d],P2[d],…,PK[d];
步骤2、堆叠所有用户的预编码矩阵得到P同时初始化Φk
Figure BDA00032664086700001322
步骤3、通过迭代式(39)计算
Figure BDA00032664086700001323
Figure BDA00032664086700001324
步骤4、根据表达式(37)计算
Figure BDA00032664086700001325
步骤5、堆叠除第k个用户以外所有用户的预编码矩阵得到P\k同时初始化Φ′k
Figure BDA0003266408670000141
Figure BDA00032664086700001421
以及
Figure BDA0003266408670000142
步骤6、通过如下迭代式计算
Figure BDA0003266408670000143
Figure BDA0003266408670000144
以及
Figure BDA0003266408670000145
步骤7、通过步骤6中的迭代式得到
Figure BDA0003266408670000146
Figure BDA0003266408670000147
步骤8、根据表达式(38)计算
Figure BDA0003266408670000148
步骤9、最后根据
Figure BDA0003266408670000149
得到遍历和速率的确定性等同。
观察遍历和速率的确定性等同,发现仅和预编码矩阵以及参数化协方差矩阵有关。基于已知的统计信道信息,可以直接计算得到
Figure BDA00032664086700001410
得到遍历和速率的确定性等同之后,利用
Figure BDA00032664086700001411
关于预编码矩阵的导数计算一些矩阵函数的确定性等同。观察发现Ak[d]和Bk[d]可以通过
Figure BDA00032664086700001412
Figure BDA00032664086700001413
关于Pk的导数联系。将
Figure BDA00032664086700001414
Figure BDA00032664086700001415
关于Pk求导得到Ak[d]和Bk[d]的确定性等同表达。在预编码迭代更新公式中代入Ak[d]和Bk[d]相应的确定性等同得到最终的闭式表达。
Figure BDA00032664086700001416
关于Pk求导可以得到
Figure BDA00032664086700001417
根据链式法则可以得到
Figure BDA00032664086700001418
其中Bl[d]表达如下
Figure BDA00032664086700001419
矩阵函数Ak[d]与遍历和速率的导数有如下关系
Figure BDA00032664086700001420
综合以上结论,得到下面定理。
定理2、Ak[d]和Bk[d]的确定性等同给出如下
Figure BDA0003266408670000151
Figure BDA0003266408670000152
矩阵函数D[d]的确定性等同定义为
Figure BDA0003266408670000153
将Ak[d]和D[d]替换为相应的确定性等同,得到最优预编码迭代式的闭式解为
Figure BDA0003266408670000154
在下面的算法2中给出FDD大规模MIMO线性统计预编码设计的具体步骤:
步骤1、初始化预编码矩阵并对它们归一化使之满足功率条件。设置初始迭代次数为d=0;
步骤2、根据算法1计算Γk和Γ′k
步骤3、根据步骤2中的结果计算
Figure BDA0003266408670000155
Figure BDA0003266408670000156
Figure BDA0003266408670000157
步骤4、更新预编码矩阵
Figure BDA0003266408670000158
并归一化使之满足能量条件;
步骤5、更新迭代次数d=d+1。
五、实施效果
为了使本技术领域的人员更好地理解本发明方案,下面给出一种具体系统配置下本实施例中采用的FDD大规模MIMO和速率最优统计预编码设计。
利用QuaDRiGa信道模型生成信道,考虑一个FDD大规模MIMO系统,其中基站天线配置Mx=16,Mz=8的UPA,用户侧配置ULA为Mk=4。采样数以及采样间隔分别设置为1000和2s。上下行中心频率分别设置为1.74GHz 和1.84GHz。假设K=20个用户随机均匀分布在小区,预编码的总功率限制为P=1。信噪比
Figure BDA0003266408670000161
取值范围设置为从-10dB到20dB。
图2(a),图2(b),图2(c),图2(d)给出FDD上下行信道统计信息的分布,观察发现统计信息的功率只在某些波束上有分布。相比于瞬时信道信息,统计信息的这种稀疏性在基站侧更容易被处理。对于同一个用户,上下行的统计信息分布几乎是相同的,其中微小的差别可以认为是信道模型引入的误差造成的。利用这个性质,FDD传统预编码方案中CSI的获取开销可以被消除。
图3给出FDD统计预编码和BDMA的性能对比曲线。两种算法的性能都随 SNR呈现线性增长。进一步,观察发现FDD统计预编码在所有SNR上都优于 BDMA而且两者的性能差异随着SNR的增加而线性增加。由于BDMA的优化函数使用的是遍历和速率的上界,该上界在高信噪比时不能很好的逼近遍历和速率,这使得FDD统计预编码算法在高信噪比时明显优于BDMA。具体来说,当SNR=20db时,FDD统计预编码相比于BDMA有22%左右的性能提升。与此同时,遍历和速率相对应的确定性等同也在图中给出,结果显示了所提算法的准确性。
图4给出FDD统计预编码和SLNR-SBF以及ZF预编码的性能对比曲线。相比于SLNR-SBF,FDD统计预编码有较显著的性能提升。具体来说,当过采样因子为2,SNR=20db时,FDD统计预编码相比于SLNR-SBF有63%左右的性能提升。该性能提升增长至153%当过采样因子设为1,即考虑了传统波束域信道。结果表明,直接对遍历和速率进行优化的方案会带来更好的性能增益。进一步,FDD统计预编码相比于ZF预编码有193%左右的性能提升。这些结果表明FDD统计预编码相比于SLNR-SBF以及ZF预编码,更有效地抑制了用户间的干扰噪声。
可以理解,本发明是通过一些实施例进行描述的,本领域技术人员知悉的,在不脱离本发明的精神和范围的情况下,可以对这些特征和实施例进行各种改变或等效替换。另外,在本发明的教导下,可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此,本发明不受此处所公开的具体实施例的限制,所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。

Claims (10)

1.一种FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,包括如下步骤:
步骤1、基站利用上行信道估计以及统计互异性获得各用户终端的下行统计信道信息;上行统计信道信息和下行统计信道信息是一致的;所述统计信道信息从波束域统计信道模型中提取得到,该波束域统计信道模型的导向矢量数量超过了天线数;上下行信道通过空间采样矩阵进行相互转换;
步骤2、获得遍历和速率的基于统计信道信息的确定性等同表达式;
步骤3、基站利用遍历和速率的确定性等同进行FDD大规模MIMO和速率最优统计预编码设计。
2.根据权利要求1所述的FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,所述步骤1中通过对基站和用户侧的导向矢量过采样得到统计信道模型,该模型从波束域角度描述空间相关性。
3.根据权利要求1所述的FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,所述下行统计信道信息由基站通过对上行信道统计平均并利用统计互异性获得,不需要下行信道估计以及上行反馈过程。
4.根据权利要求1所述的FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,所述步骤1中基站侧的上行信道与下行信道空间采样矩阵进行相互转化,包括以下步骤:
步骤1.1、基站侧上下行信道的导向矢量的计算和中心频率有关,通过定义的上下行中心频率商因子进行相互转换;其中导向矢量中的采样方向余弦通过均匀采样得到;
步骤1.2、堆叠导向矢量分别得到水平和垂直方向的块矩阵,基站侧空间采样矩阵等于两个块矩阵的克罗内克积;
步骤1.3、利用基站侧上下行导向矢量的关系进而得到基站侧上下行空间采样矩阵的关系。
5.根据权利要求4所述的FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,所述步骤1.1中采样方向余弦的采样数大于对应的天线数,使得堆叠得到空间采样矩阵维度大于天线数。
6.根据权利要求1所述的FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,在所述的和速率最优统计预编码设计中,基站根据遍历和速率最大化的准则,进行各用户终端的线性预编码矩阵设计,遍历和速率表达式中的干扰加噪声协方差矩阵等效为高斯噪声的协方差矩阵,同时定义干扰加噪声协方差矩阵的瞬时表达。
7.根据权利要求1所述的FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,所述步骤2中将所述遍历和速率最大化准则替换为遍历和速率下界函数最大化准则;或者将所述遍历和速率最大化准则中和速率替换为其确定性等同。
8.根据权利要求1所述的FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,所述步骤3具体包括以下步骤:
通过MM算法将加权遍历和速率最大化预编码设计问题转化为迭代求解二次型优化问题得到最优预编码矩阵的迭代更新式;
通过确定性等同对最优迭代更新式中带有的关于随机矩阵的期望项作近似,得到最优预编码矩阵的闭式表达。
9.根据权利要求8所述的FDD大规模MIMO和速率最优统计预编码设计方法,其特征在于,所述方法结合MM算法以及确定性等同进行预编码设计,包括:
步骤3.1、初始化预编码矩阵并归一化使之满足功率条件;
步骤3.2、根据预编码矩阵以及统计信道信息计算遍历和速率确定性等同中的相关项;
步骤3.3、根据确定性等同的相关项计算预编码矩阵最优迭代更新式中带有的关于随机矩阵的期望项;
步骤3.4、更新预编码迭代更新的闭式表达,并归一化预编码矩阵;
重复步骤3.2到步骤3.4直到达到预设迭代次数或预编码收敛。
10.一种计算设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述计算机程序被加载至处理器时实现根据权利要求1-9任一项所述的FDD大规模MIMO和速率最优统计预编码设计方法。
CN202111087692.4A 2021-09-16 2021-09-16 Fdd大规模mimo和速率最优统计预编码方法及设备 Active CN113839695B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111087692.4A CN113839695B (zh) 2021-09-16 2021-09-16 Fdd大规模mimo和速率最优统计预编码方法及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111087692.4A CN113839695B (zh) 2021-09-16 2021-09-16 Fdd大规模mimo和速率最优统计预编码方法及设备

Publications (2)

Publication Number Publication Date
CN113839695A true CN113839695A (zh) 2021-12-24
CN113839695B CN113839695B (zh) 2022-06-21

Family

ID=78959531

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111087692.4A Active CN113839695B (zh) 2021-09-16 2021-09-16 Fdd大规模mimo和速率最优统计预编码方法及设备

Country Status (1)

Country Link
CN (1) CN113839695B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117749255A (zh) * 2024-02-19 2024-03-22 成都本原星通科技有限公司 一种大规模mimo卫星通信的终端分组方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105071843A (zh) * 2015-07-29 2015-11-18 东南大学 大规模mimo系统低复杂度多项式展开矩阵求逆方法及应用
CN109981153A (zh) * 2019-04-11 2019-07-05 东南大学 一种人工噪声辅助的大规模mimo安全统计预编码方法
CN110635836A (zh) * 2019-09-12 2019-12-31 重庆大学 一种基于波束选择的毫米波大规模mimo系统mmse-pca信道估计方法
CN111181671A (zh) * 2019-12-27 2020-05-19 东南大学 一种基于深度学习的下行信道快速重建方法
US20200204226A1 (en) * 2017-08-31 2020-06-25 Southeast University Massive mimo robust precoding transmission method
CN111988069A (zh) * 2020-07-16 2020-11-24 东南大学 大规模mimo广义特征向量结构预编码求解方法、装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105071843A (zh) * 2015-07-29 2015-11-18 东南大学 大规模mimo系统低复杂度多项式展开矩阵求逆方法及应用
US20200204226A1 (en) * 2017-08-31 2020-06-25 Southeast University Massive mimo robust precoding transmission method
CN109981153A (zh) * 2019-04-11 2019-07-05 东南大学 一种人工噪声辅助的大规模mimo安全统计预编码方法
CN110635836A (zh) * 2019-09-12 2019-12-31 重庆大学 一种基于波束选择的毫米波大规模mimo系统mmse-pca信道估计方法
CN111181671A (zh) * 2019-12-27 2020-05-19 东南大学 一种基于深度学习的下行信道快速重建方法
CN111988069A (zh) * 2020-07-16 2020-11-24 东南大学 大规模mimo广义特征向量结构预编码求解方法、装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
AN-AN LU 等: "Free Deterministic Equivalents for the Analysis of MIMO Multiple Access Channel", 《IEEE TRANSACTIONS ON INFORMATION THEORY》, vol. 62, no. 8, 26 May 2016 (2016-05-26), pages 4604 - 4629, XP011616578, DOI: 10.1109/TIT.2016.2573309 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117749255A (zh) * 2024-02-19 2024-03-22 成都本原星通科技有限公司 一种大规模mimo卫星通信的终端分组方法及系统

Also Published As

Publication number Publication date
CN113839695B (zh) 2022-06-21

Similar Documents

Publication Publication Date Title
CN111181619A (zh) 基于深度强化学习的毫米波混合波束成形设计方法
WO2021109419A1 (zh) 大规模mimo波束域鲁棒预编码传输方法与系统
Huang et al. Beam-blocked channel estimation for FDD massive MIMO with compressed feedback
CN110138425B (zh) 低复杂度阵列天线多输入多输出系统混合预编码算法
CN112118033B (zh) 多用户大规模mimo系统的非线性混合预编码设计方法
Kreyndelin et al. Effective precoding and demodulation techniques for 5G communication systems
Xie et al. Dictionary learning for channel estimation in hybrid frequency-selective mmWave MIMO systems
JP7170863B2 (ja) アナログプリコーディング及びアナログ結合を可能にする方法
JP6472575B2 (ja) 正則化ゼロフォーシングプリコーダーを構成するためのパラメーターを決定する方法
CN113839695B (zh) Fdd大规模mimo和速率最优统计预编码方法及设备
Miretti et al. Team precoding towards scalable cell-free massive MIMO networks
Reddy et al. Literature review and research direction towards channel estimation and hybrid pre-coding in mmWave massive MIMO communication systems
Zarei et al. Robust MSE-balancing hierarchical linear/Tomlinson-Harashima precoding for downlink massive MU-MIMO systems
CN111988069B (zh) 大规模mimo广义特征向量结构预编码求解方法、装置
Ullah et al. Beyond traditional beamforming: Singular vector projection techniques for mu-mimo interference management
CN115065446B (zh) 一种基于流形学习的信道信息反馈方法
Dreifuerst et al. Machine Learning Codebook Design for Initial Access and CSI Type-II Feedback in Sub-6GHz 5G NR
CN107248878A (zh) 一种预编码的方法及装置
Zhou et al. Manifold learning inspired dynamic hybrid precoding with antenna partitioning algorithm for dual-hop hybrid fso-rf systems
JP2009268106A (ja) 信号対干渉電力と雑音比の決定方法およびその装置
CN109792271A (zh) 执行数据传输的方法和装置
Gautam et al. Hybrid precoding for millimeter wave MIMO: Trace optimization approach
Yu et al. Precoding for uplink distributed antenna systems with transmit correlation in Rician fading channels
CN112671435B (zh) 一种fdd大规模天线中基于字典迭代的信道估计方法
CN114793127B (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