CN114866117A - 天波大规模mimo波束结构预编码传输方法与系统 - Google Patents

天波大规模mimo波束结构预编码传输方法与系统 Download PDF

Info

Publication number
CN114866117A
CN114866117A CN202210317238.1A CN202210317238A CN114866117A CN 114866117 A CN114866117 A CN 114866117A CN 202210317238 A CN202210317238 A CN 202210317238A CN 114866117 A CN114866117 A CN 114866117A
Authority
CN
China
Prior art keywords
user
domain
precoder
beam domain
sky
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.)
Pending
Application number
CN202210317238.1A
Other languages
English (en)
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 CN202210317238.1A priority Critical patent/CN114866117A/zh
Priority to PCT/CN2022/086916 priority patent/WO2023184590A1/zh
Publication of CN114866117A publication Critical patent/CN114866117A/zh
Pending legal-status Critical Current

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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明公开了一种天波大规模MIMO波束结构预编码传输方法与系统。天波大规模MIMO通信基站利用波束结构预编码器生成发送信号，实现下行预编码传输。波束结构预编码器由各用户低维波束域预编码器、各用户波束映射模块、波束调制模块构成，各用户低维波束域预编码器为各用户波束集合上的预编码器，各用户波束映射将各用户低维波束域预编码信号映射成完整的波束域发送信号，波束调制为波束矩阵乘以波束域发送信号矢量，该波束域发送信号矢量为各用户波束域发送信号矢量之和。本发明可以解决天波大规模MIMO下行预编码传输的设计和实现复杂度问题，能够显著提升天波通信的频谱与功率效率、以及传输速率和传输距离。

Description

天波大规模MIMO波束结构预编码传输方法与系统
技术领域
本发明涉及天波通信领域，特别是一种涉及天波大规模MIMO波束结构预编码传输方法 与系统。
背景技术
天波通信通常利用1.6～30MHz的短波波段，通过电离层反射的方式实现长达几千公里的 远距离信号传输。由于电离层信道传播特性异常复杂，传统的单入单出天波通信系统往往仅 具有很低的系统数据率。大规模MIMO技术通过在基站侧配置大量天线，在同一时频资源上 服务大量用户，可以极大提升系统容量和可靠性。
现有地面蜂窝通信系统一般采用在空间域设计用户预编码器，其维度等于基站天线数。 在大规模MIMO系统中，这个维度是相当大的，因此求解高维的空间域预编码设计问题的计 算量非常可观。另一方面，天波大规模MIMO信道的角度扩展典型的非常小，也就是说波束 域信道是稀疏的，发射端发送的无线信号仅通过有限的方向到达接收端，所以可以选择在部 分波束方向上发送信号，并可取得近乎最优的和速率性能。
发明内容
鉴于此，本发明的目的在于提供一种天波大规模MIMO波束结构预编码传输方法与系统， 该方法利用具有波束结构的空间域预编码器，通过灵活选择波束集合，进行相应的低维波束 域预编码器设计，可以大幅提升天波通信的频谱效率，进而大幅提高传输速率和传输距离， 并显著降低空间域预编码器设计的复杂度。
为了实现上述目的，本发明采用如下技术方案：
天波大规模MIMO波束结构预编码传输方法，所述波束结构预编码传输方法包括：天波 大规模MIMO通信基站利用具有波束结构的预编码器生成发送信号，实现与一组用户的下行 预编码传输；波束结构预编码器由各用户低维波束域预编码器、各用户波束映射模块、波束 调制模块构成，各用户低维波束域预编码器为各用户波束集合上的预编码器，各用户波束映 射将各用户低维波束域预编码信号映射成完整的波束域发送信号，波束调制为波束矩阵乘以 波束域发送信号矢量，该波束域发送信号矢量为各用户波束域发送信号矢量之和；基站依据 各用户波束基信道表示及波束域信道信息，设计各用户低维波束域预编码器。
作为本发明的一种改进，所述波束矩阵为选定的一组空间角度采样格点所对应的阵列方 向矢量构成的矩阵，每个阵列方向矢量称为一个波束。
作为本发明的一种改进，所述各用户波束集合为各用户波束基信道表示中波束域信道非 零元素所对应波束的集合或包含该波束集合的选定集合。
作为本发明的一种改进，所述波束基信道表示为波束矩阵乘以波束域信道矢量；所述波 束域信道信息包括波束域信道矢量的估计值及估计误差的方差。
作为本发明的一种改进，所述波束域预编码器的设计包括：优化目标为最大化系统和速 率的设计、最大化系统遍历和速率的设计及最大化系统遍历和速率上界的设计，其中，
(1)所述优化目标为最大化系统和速率的设计，利用波束矩阵、各用户波束映射矩阵、 各用户波束域预编码器和各用户信道估计值更新系统和速率表达式，将空间域预编码器设计 问题转化为波束域预编码器设计问题，且波束域预编码器的迭代设计包括如下步骤：
①初始化各用户波束域预编码器，使其满足功率约束；
②利用MM算法框架得到系统和速率在当前迭代下的凸替代函数；
③使用拉格朗日乘子法求解当前迭代的凸问题；
重复步骤②-③直到达到预设迭代次数或预编码收敛，得到各用户最优波束域预编码器。
(2)所述优化目标为最大化遍历和速率的设计，利用波束矩阵、各用户波束映射矩阵、 各用户波束域预编码器、各用户波束域信道和各用户波束域统计信道信息更新系统遍历和速 率表达式，将空间域预编码器设计问题转化为波束域预编码器设计问题，且波束域预编码器 的迭代设计包括如下步骤：
①初始化各用户波束域预编码器，使其满足功率约束；
②利用MM算法框架得到系统遍历和速率在当前迭代下的凸替代函数；
③使用拉格朗日乘子法求解当前迭代的凸问题；
重复步骤②-③直到达到预设迭代次数或预编码收敛，得到各用户最优波束域预编码器。
(3)所述优化目标为最大化系统遍历和速率上界的设计，系统遍历和速率上界为系统遍 历和速率利用詹森不等式得到，表达式包括波束矩阵、各用户波束映射矩阵、波束域预编码 器和各用户波束域统计信道信息，将空间域预编码器设计问题转化为波束域预编码器设计问 题，且波束域预编码器的迭代设计包括如下步骤：
①初始化各用户波束域预编码器，使其满足功率约束；
②利用MM算法框架得到系统遍历和速率上界在当前迭代下的凸替代函数；
③使用拉格朗日乘子法求解当前迭代的凸问题；
重复步骤②-③直到达到预设迭代次数或预编码收敛，得到各用户最优波束域预编码器。
作为本发明的一种改进，所述根据设计生成的波束域预编码器实施与用户的下行信号传 输包括如下步骤：
(1)利用用户波束域预编码器与其发送数据符号相乘生成低维波束域发送信号；
(2)利用波束映射矩阵和低维波束域发送信号相乘得到用户波束域发送信号；
(3)将各用户波束域发送信号进行叠加，得到包含所有用户的波束域发送信号；
(4)利用波束矩阵和所有用户的波束域发送信号相乘生成空间域发送信号。
其中，所述步骤(4)利用Chirp-z变换进行有效实现。
作为本发明的一种改进，所述天波大规模MIMO通信基站包含大规模天线阵列，工作载 频为1.6～30MHz的短波波段，基站通过电离层反射的方式与用户进行信号传输。
天波大规模MIMO波束结构预编码传输系统，包括基站和多个用户，其特征在于，所述 基站实现根据权利要求1-7任一项所述的天波大规模MIMO波束结构预编码传输方法。
本发明的有益效果是：
本发明能够大幅提升天波大规模MIMO通信在各种典型通信场景下的频谱效率，进而大 幅提高传输速率和传输距离；基站侧预编码器的设计仅涉及各用户低维波束域预编码器设计， 设计复杂度可以大幅度降低，而所形成的波束结构预编码则能够显著降低基站侧预编码实现 复杂度。本发明充分利用天波通信信道波束域稀疏特性，进行低维的波束域预编码器的最优 设计，显著降低空间域预编码器设计的复杂度，且具有近乎最优的(遍历)和速率性能。该 方案可以根据需要灵活调整波束映射矩阵来设计生成用户预编码器。
附图说明
图1为波束结构预编码系统框图；
图2为实施例1中提供的一种天波大规模MIMO波束结构预编码传输方法的流程示意图；
图3为实施例1中提供的一种天波大规模MIMO波束结构预编码传输方法与基于MMSE (minimum mean-squared error)预编码器传输方法的遍历和速率结果比较图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附 图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明 一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有 做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1：
参见图1-图3，本实施例提供一种天波大规模MIMO波束结构预编码传输方法，该方法 主要适用于基站配备大规模天线阵列以同时服务大量单天线用户的天波大规模MIMO通信系 统。该方法具体包括：
天波大规模MIMO通信基站利用具有波束结构的预编码器生成发送信号，实现与一组用 户的下行预编码传输。
波束结构预编码器由各用户低维波束域预编码器、各用户波束映射模块、波束调制模块 构成，各用户低维波束域预编码器为各用户波束集合上的预编码器，各用户波束映射将各用 户低维波束域预编码信号映射成完整的波束域发送信号，波束调制为波束矩阵乘以波束域发 送信号矢量，该波束域发送信号矢量为各用户波束域发送信号矢量之和。
基站依据各用户波束基信道表示及波束域信道信息，设计各用户低维波束域预编码器。
在本实施例中通过具体的通信系统实例对该一种天波大规模MIMO波束结构预编码传输 方法的具体实现过程作详细说明，需要说明的是该方法不仅适用于本实施例中所举的具体系 统模型，也同样适用于其它配置的系统模型。
1、系统模型
1.1、系统设置和信号模型
考虑一个天波大规模MIMO通信系统，考虑OFDM调制且工作在时分双工TDD模式。基站配备大规模天线阵列，为M天线、天线间距为d的均匀线阵，服务U个单天线用户。系 统载频为fc且位于1.6～30MHz的短波波段。基站通过电离层反射的方式与用户进行信号传输。 天波大规模MIMO-OFDM通信系统每个信号帧包含L个OFDM符号且分为
Figure BDA0003570193230000041
个上行数据符号、1个上行训练符号及
Figure BDA0003570193230000042
个下行数据符号。
记
Figure BDA0003570193230000043
表示发给用户u的模拟基带信号。用户u接收到的模拟基带信号表示为
Figure BDA0003570193230000044
其中
Figure BDA0003570193230000045
表示用户u的下行信道冲激响应，zu(t)表示复白高斯噪声过程。
记Nc、Nc和Ts分别表示子载波数目、循环前缀长度和系统采样间隔。用
Figure BDA0003570193230000046
表示 发给用户u的子载波k符号
Figure BDA0003570193230000047
上的信号。用户u处解调的子载波k符号
Figure BDA0003570193230000048
上的信号表示为
Figure BDA0003570193230000049
其中，
Figure BDA00035701932300000410
表示用户u的下行信道频率响应，
Figure BDA00035701932300000411
其中，gu(t,f)是hu(t,τ)的傅里叶变换，Δf＝1/Tc为子载波间隔，
Figure BDA00035701932300000412
是循环对称复高斯噪 声且均值为零、方差为
Figure BDA0003570193230000051
1.2、波束基信道模型
假设基站和用户u之间有Pu条可分辨的传播路径。用τu,p表示基站到用户u的路径p的传 播时延。记Δτ＝d/c，其中c表示光速。基站第m根天线到用户u的信道冲激响应为
Figure BDA0003570193230000052
其中，αu,p(t)和Ωu,p分别表示路径p的复值增益和方向余弦，τu,p为基站第1根天线到用户u 的路径p的时延。
假设路径p包含Qp条子径，且αu,p(t)可以被建模为
Figure BDA0003570193230000053
其中，βu,p,q、φu,p,q和υu,p,q分别表示子径q的增益、初始相位和多普勒频移。假设φu,p,q在区间 [0,2π)上均匀分布，当Qp趋于无穷大时，αu,p(t)是一个零均值复高斯随机过程。
根据式(3)，子载波k符号
Figure BDA0003570193230000054
上的用户u的信道频率响应可以表示为
Figure BDA0003570193230000055
其中，
Figure BDA0003570193230000056
是一个零均值复高斯随机变量且
Figure BDA0003570193230000057
表示子载波k上对应方向余弦Ω的阵列方向矢量。上标T表示矩阵或矢量的转置。
对方向余弦Ω进行均匀采样。记
Figure RE-GDA0003731711050000058
1≤n≤N，其中N≥M表 示采样数目。集合
Figure RE-GDA0003731711050000059
中的路径方向余弦被近似为ξn＝(2n-1)/N-1。记集合
Figure RE-GDA00037317110500000510
可以得到波束基信道模型为
Figure BDA0003570193230000059
其中，v(ξn,k)是采样的方向余弦，其对应一个物理的空间波束，且有
Figure BDA0003570193230000061
是独立的零均值复高斯分布波束域信道元素。∩表示集合的交集。
记
Figure BDA0003570193230000062
为尺度化的用户u在子载波k符号
Figure BDA0003570193230000063
上的波束域 信道，定义子载波k上的波束矩阵为
Figure BDA0003570193230000064
可以看出，式(10)给出的波束矩阵为选定的一组空间角度采样格点所对应的阵列方向矢量 构成的矩阵，每个阵列方向矢量称为一个波束。可以改写式(8)为
Figure BDA0003570193230000065
其可以被认为是信道估计前的先验信道模型。
定义用户u的波束域信道功率矢量为
Figure BDA0003570193230000066
其中
Figure BDA0003570193230000067
其独立于符号与子载波。可以简记
Figure BDA0003570193230000068
上标*表示矩阵或矢量的共轭，E{·}表示求 数学期望，
Figure BDA0003570193230000069
表示哈达玛积，|·|表示求模值运算。
1.3、波束基后验信道模型
记
Figure BDA00035701932300000610
表示用户u在子载波k训练符号
Figure BDA00035701932300000611
上的估计信道，其中
Figure BDA00035701932300000612
表示估计的波 束域信道。由于信道老化的影响，通过用一阶高斯马尔可夫过程描述信道时变特性，可以建 模子载波k符号
Figure BDA00035701932300000613
上的波束基信道为
Figure BDA00035701932300000614
其可以被认为是后验信道模型并可描述不同移动场景下的不完美信道状态信息。
Figure BDA00035701932300000615
是训练符 号
Figure BDA00035701932300000616
和符号
Figure BDA00035701932300000617
之间的时间相关系数。一般来说，
Figure BDA00035701932300000618
和信道多普勒扩展有关，可以通过选择 合适的
Figure BDA00035701932300000619
来描述相应的信道不确定性。特别地，当
Figure BDA00035701932300000620
非常接近于1时，信道可以被认为是 准静态的，当
Figure BDA0003570193230000071
接近0时，信道变化非常剧烈。式(14)给出的波束基信道表示为波束矩阵 乘以波束域信道矢量。波束域信道信息包括波束域信道矢量的估计值及估计误差的方差。
天波大规模MIMO通信中波束域信道一般是空间稀疏的，这表明了波束域信道的大部分 元素接近0。记波束域信道非零元素所对应波束的集合为
Figure BDA0003570193230000072
且有Nu＝|Nu|。 定义集合N＝N1∪…∪Nu且有
Figure BDA0003570193230000073
其中∪表示集合的并集。用
Figure BDA0003570193230000074
表示 降维波束矩阵。可以改写式(14)为
Figure BDA0003570193230000075
其中
Figure BDA0003570193230000076
且
Figure BDA0003570193230000077
个波束就可以得到信道频率响应的完整表示。
2、鲁棒预编码
2.1、问题形成
考虑子载波k符号
Figure BDA0003570193230000078
上的下行传输。为了表示简洁，后续省略符号下标k和
Figure BDA0003570193230000079
考虑线性 预编码，式(2)给出的信号模型可以改写为
Figure BDA00035701932300000710
其中，pu是用户u的预编码器，su是发送给用户u的数据符号且具有零均值和单位方差。
用户u处的干扰+噪声
Figure RE-GDA00037317110500000713
被当作高斯噪声，其协方差为
Figure RE-GDA00037317110500000714
当用户u可以获得该方差时，其遍历速率可以表示为
Figure BDA00035701932300000713
最大化遍历和速率的鲁棒预编码设计问题表述为
Figure BDA00035701932300000714
Figure BDA0003570193230000081
其中P为总的发送功率约束。由于大规模MIMO通信的基站天线数是非常多的，所以上述优 化问题是大维的，其最优解的求取需要复杂的计算。
2.2、波束结构鲁棒预编码
定义矢量
Figure BDA0003570193230000082
u＝1,…,U。可以改写遍历速率ru为
Figure BDA0003570193230000083
其中包含波束域信道
Figure BDA0003570193230000084
而不是空间域信道h1,…,hU。因此，q1,…,qU可以被认为是波束 域信道下的预编码器。另一方面，q1,…,qU可以由空间域预编码器p1,…,pU通过波束矩阵
Figure BDA0003570193230000085
转 换得到。因此，称q1,…,qU为波束域预编码器。
简记
Figure BDA0003570193230000086
和
Figure BDA0003570193230000087
其中上标
Figure BDA0003570193230000088
表示伪逆运算符。接下来，考虑优化矢 量q1,…,qU来最大化遍历和速率。可以证明对于任意的矢量
Figure BDA0003570193230000089
使得
Figure BDA00035701932300000810
成立且
Figure BDA00035701932300000811
的充分必要条件为
Aa＝0 (21)
由此可得，如果qu满足Aqu＝0，空间域预编码器pu根据关系式
Figure BDA00035701932300000812
总是存在的。可以 得到优化问题(19)的最优解为
Figure BDA00035701932300000813
当u＝1,…,U (22)
其中
Figure BDA00035701932300000814
Figure BDA00035701932300000815
当u＝1,…,U
式(22)给出了最优的空间域预编码结构。当求解优化问题(23)得到最优解后，最优 的空间域预编码器在式(22)中给出。可是问题(23)的维度可能依然相当大，其求解过程是计算复杂的。进一步假设用户数是有限的，且用户的方向余弦是离散且有限的。当M趋于无穷大时，有
Figure BDA0003570193230000091
进一步有
Figure BDA0003570193230000092
和
Figure BDA0003570193230000093
此时最优解(22)变为
Figure BDA0003570193230000094
当u＝1,…,U (24)
其中
Figure BDA0003570193230000095
Figure BDA0003570193230000096
得到波束域优化问题(25)的最优解后，渐近最优空间域预编码器表达为式(24)。特 别地，观察到(22)式的空间域预编码器结构变为一个如式(24)所示的简单波束结构预编 码器。换句话说，当基站侧具有足够多的天线数时，波束结构预编码器设计是渐近最优的。 进一步得到问题(25)的最优解满足
Figure BDA0003570193230000097
当
Figure BDA0003570193230000098
上式揭示了当M足够大时，波束域预编码器qu的非零波束集合外的元素都为0。基于此 结论，可以仅关注q1,…,qU的非零元素，其维度是相当小的。
定义矩阵
Figure BDA0003570193230000099
为
Figure BDA00035701932300000910
此时波束结构预编码器(24)变为
pu＝VΦuwu (28)
其中wu是对应qu的非零元素的低维波束域预编码器，Φu是波束映射矩阵用来映射用户u的 非零波束。
根据式(28)，发送信号矢量x可以被改写为
Figure BDA00035701932300000911
根据式(29)，利用设计生成的波束域预编码器实施与用户的下行信号传输包括如下步骤：
步骤1：wusu为利用用户u的波束域预编码器与其发送数据符号相乘生成低维波束域发 送信号；
步骤2：Φuwusu表示利用波束映射矩阵Φu和低维波束域发送信号相乘得到用户u的波束 域发送信号；
步骤3：将各用户波束域发送信号进行叠加，得到包含所有用户的波束域发送信号
Figure BDA0003570193230000101
步骤4：利用波束矩阵V和所有用户的波束域发送信号
Figure BDA0003570193230000102
相乘生成空间域发送信 号x，即为波束调制过程。
特别地，步骤4可以利用Chirp-z变换进行有效实现。
简记
Figure BDA0003570193230000103
和
Figure BDA0003570193230000104
此时有
Figure BDA0003570193230000105
定义
Figure BDA0003570193230000106
为
Figure BDA0003570193230000107
式(18)给出的用户u的遍历速率表达式变为
Figure BDA0003570193230000108
遍历速率表达式(32)包含了波束矩阵、用户波束映射矩阵、用户波束域预编码器和用户信 道估计值。
接下来考虑优化波束域预编码器w1,…,wU来最大化遍历和速率，问题表述为
Figure BDA0003570193230000109
Figure BDA00035701932300001010
这样空间域预编码器设计问题(19)转化为波束域预编码器设计问题(33)。和空间域预编 码器设计问题(19)式相比，波束域预编码器设计问题(33)的维度是非常小的。特别地， 当M足够大且用户的非零波束变得正交时，问题(33)和问题(25)是相同的。
2.3、基于Chirp-z变换的波束调制
式(29)给出的波束调制可以改写为
x＝Vs (34)
其中
Figure BDA0003570193230000111
在式(34)中，注意到波束矩阵V可以被认为是CZT矩阵，表示为三个矩阵的乘积：分别为对角阵
Figure BDA0003570193230000112
托普利兹矩阵
Figure BDA0003570193230000113
和另一个对角阵
Figure BDA0003570193230000114
即
V＝ΓTΛ (36)
其中
Figure BDA0003570193230000115
Figure BDA0003570193230000116
Figure BDA0003570193230000117
其中
Figure BDA0003570193230000118
进一步表示托普利兹矩阵T为
Figure BDA0003570193230000119
其中S≥M+N-1，Π是一个对角阵，FS×M表示一个矩阵包含S点DFT矩阵的前M列。将式(40)和(36)代入式(34)，可以得到
Figure BDA00035701932300001110
其可通过快速傅里叶变换(FFT)和逆FFT来有效实现。
3、波束域鲁棒预编码器设计
3.1、基于MM算法框架的波束域预编码器设计
波束域预编码器设计问题(33)是非凸的，其全局最优解很难获得。基于MM算法框架， 将推导一个迭代的局部最优解。在第d次迭代中，用
Figure BDA0003570193230000121
表示wu，定义函数
Figure BDA0003570193230000122
来最小化
Figure BDA0003570193230000123
当
Figure BDA0003570193230000124
Figure BDA0003570193230000125
进一步地，式(42)和(43)表明
Figure BDA0003570193230000126
当u＝1,…,U (44)
接下来寻找替代函数f其在任意点上最小化
Figure BDA0003570193230000127
然后最大化f来得到原问题的迭代解。 特别地，记
Figure BDA0003570193230000128
为最大化f的解，根据式(42)和(43)可以得到
Figure BDA0003570193230000129
条件(44)和(45)可以确保生成序列可以收敛到局部最优解。
简记ρu为
Figure BDA00035701932300001210
可以得到下面的关于
Figure BDA00035701932300001211
的一个替代函数为
Figure BDA00035701932300001212
其中
Figure BDA00035701932300001213
为常数且
Figure BDA00035701932300001214
Figure BDA00035701932300001215
Figure BDA00035701932300001216
Figure BDA0003570193230000131
基于替代函数f，第(d+1)次迭代的最优解可以通过求解下面的优化问题得到
Figure BDA0003570193230000132
Figure BDA0003570193230000133
该问题是一个凹二次型优化问题，其最优解可通过拉格朗日乘子法求解。
优化问题(52)的拉格朗日函数可以表示为
Figure RE-GDA0003731711050000134
其中μ≥0表示拉格朗日乘子。根据一阶最优条件，问题(52)的最优解为
Figure BDA0003570193230000135
其中μop是最优拉格朗日乘子。注意到
Figure BDA0003570193230000136
是μ的单调递减函数。因此，如 果μop＝0且
Figure BDA0003570193230000137
最优解变为
Figure BDA0003570193230000138
否则的话，可以 利用二分法来求取最优μop。下面给出了波束域鲁棒预编码设计的具体过程：
步骤1：初始化各用户波束域预编码器
Figure BDA0003570193230000139
使其满足功率约束
Figure BDA00035701932300001310
设置d＝0；
步骤2：利用MM算法框架得到系统遍历和速率在当前迭代下的凸替代函数，即分别根 据式(48)、(49)、(50)和(51)计算
Figure BDA00035701932300001311
和
Figure BDA00035701932300001312
步骤3：使用拉格朗日乘子法求解当前迭代的凸问题，即根据式(54)更新
Figure BDA00035701932300001313
并设置 d＝d+1。重复步骤2和步骤3直到达到预设迭代次数或预编码收敛；
得到各用户最优波束域预编码器
Figure BDA00035701932300001314
当u＝1,…,U。
3.2、低复杂度波束域预编码器设计
一般来说，式(32)给出的遍历速率没有闭式表达式，往往需要复杂的蒙特卡洛平均来 计算相应的期望值。进一步研究了基于遍历和速率上界的低复杂度波束域鲁棒预编码器设计。 系统遍历速率上界为系统遍历速率ru利用詹森不等式得到
Figure BDA0003570193230000141
该速率上界对单天线用户是非常紧致的且表达式包括波束矩阵、用户波束映射矩阵、波束域 预编码器和用户波束域统计信道信息。
考虑一个新的波束域预编码器设计问题来最大化遍历和速率上界，表述为
Figure BDA0003570193230000142
Figure BDA0003570193230000143
其目标函数不包含期望运算。但是，问题(56)仍然是非凸的，其全局最优解很难获得。
观察到问题(56)和(33)的不同在于目标函数。在(56)中存在
Figure BDA0003570193230000144
而不是(33)中的
Figure BDA0003570193230000145
且不包含任何期望运算。基于MM算法框架，类似地可以推导得到问题(56)的迭代 局部最优解。在第(d+1)次迭代时，最优解可以表示为
Figure BDA0003570193230000146
其中μop是最优拉格朗日乘子，
Figure BDA0003570193230000147
Figure BDA0003570193230000148
Figure BDA0003570193230000149
通过替代定理3中
Figure BDA00035701932300001410
和
Figure BDA00035701932300001411
中的
Figure BDA00035701932300001412
为
Figure BDA00035701932300001413
且略去相应的期望运算操作。注意到
Figure BDA00035701932300001414
进一步得到
Figure BDA0003570193230000151
下面给出低复杂度波束域鲁棒预编码器设计过程：
步骤1：初始化各用户波束域预编码器
Figure BDA0003570193230000152
使其满足功率约束
Figure BDA0003570193230000153
设置d＝0；
步骤2：利用MM算法框架得到系统遍历和速率上界在当前迭代下的凸替代函数，即分 别根据式(58)和(60)计算
Figure BDA0003570193230000154
和
Figure BDA0003570193230000155
步骤3：使用拉格朗日乘子法求解当前迭代的凸问题，即根据式(61)更新
Figure BDA0003570193230000156
并设置 d＝d+1。重复步骤2和步骤3直至收敛；
得到各用户最优波束域预编码器
Figure BDA0003570193230000157
当u＝1,…,U。
为了验证本实施例提供的一种利用不完美信道状态信息的天波大规模MIMO波束结构预 编码传输方法的先进性和优越性，本实施例对该方法以及基于瞬时信道状态信息的MMSE预 编码下行传输方法进行了仿真对比试验。
具体的说，考虑天波大规模MIMO-OFDM通信系统，系统参数配置如下：载频fc＝25MHz， 天波通信基站天线阵列间距d＝5.8m，系统带宽B＝192kHz，系统采样间隔Ts＝3.9μs，子 载波间隔Δf＝125Hz，子载波个数Nc＝2048，CP点数Ng＝512。设置天波大规模MIMO通信 基站天线数M＝256，采样波束个数N＝512，用户数U＝64。定义总发送功率为64个用户在系统带宽B＝192kHz上的发送功率之和，遍历和速率为所有有效子载波上的遍历和速率的平 均。
图3给出了本实施例方法在处于不同总发送功率下的最大化遍历和速率的空间域预编码 器设计、最大化遍历和速率的波束域预编码器设计及最大化遍历和速率上界的波束域预编码 器(低复杂度波束域预编码器)设计的遍历和速率结果比较。从图3中可以看出，系统遍历和 速率结果随总发送功率的增加而增大。与空间域预编码器传输方法相比，本实施例中天波大 规模MIMO通信波束结构预编码传输方法能够取得近乎最优的系统遍历和速率性能且具有相 当低的复杂度。
本发明未详述之处，均为本领域技术人员的公知技术。
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造 性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本 发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案， 皆应在由权利要求书所确定的保护范围内。

Claims (10)

1.天波大规模MIMO波束结构预编码传输方法，其特征在于：所述方法包括以下步骤：
天波大规模MIMO通信基站利用具有波束结构的预编码器生成发送信号，实现与一组用户的下行预编码传输；
波束结构预编码器由各用户低维波束域预编码器、各用户波束映射模块、波束调制模块构成，各用户低维波束域预编码器为各用户波束集合上的预编码器，各用户波束映射将各用户低维波束域预编码信号映射成完整的波束域发送信号，波束调制为波束矩阵乘以波束域发送信号矢量，该波束域发送信号矢量为各用户波束域发送信号矢量之和；
基站依据各用户波束基(Beambased)信道表示及波束域信道信息，设计各用户低维波束域预编码器。
2.根据权利要求1所述的天波大规模MIMO波束结构预编码传输方法，其特征在于，所述波束矩阵为选定的一组空间角度采样格点所对应的阵列方向矢量构成的矩阵，每个阵列方向矢量称为一个波束。
3.根据权利要求1所述的天波大规模MIMO波束结构预编码传输方法，其特征在于，所述各用户波束集合为各用户波束基信道表示中波束域信道非零元素所对应波束的集合或包含该波束集合的选定集合。
4.根据权利要求3所述的天波大规模MIMO波束结构预编码传输方法，其特征在于，所述波束基信道表示为波束矩阵乘以波束域信道矢量；所述波束域信道信息包括波束域信道矢量的估计值及估计误差的方差。
5.根据权利要求1所述的天波大规模MIMO波束结构预编码传输方法，其特征在于，所述波束域预编码器的设计包括：优化目标为最大化系统和速率的设计、最大化系统遍历和速率的设计及最大化系统遍历和速率上界的设计。
6.根据权利要求1所述的天波大规模MIMO波束结构预编码传输方法，其特征在于，
其中，
(1)所述优化目标为最大化系统和速率的设计，利用波束矩阵、各用户波束映射矩阵、各用户波束域预编码器和各用户信道估计值更新系统和速率表达式，将空间域预编码器设计问题转化为波束域预编码器设计问题，且波束域预编码器的迭代设计包括如下步骤：
①初始化各用户波束域预编码器，使其满足功率约束；
②利用MM算法框架得到系统和速率在当前迭代下的凸替代函数；
③使用拉格朗日乘子法求解当前迭代的凸问题；
重复步骤②-③直到达到预设迭代次数或预编码收敛，得到各用户最优波束域预编码器。
(2)所述优化目标为最大化遍历和速率的设计，利用波束矩阵、各用户波束映射矩阵、各用户波束域预编码器、各用户波束域信道和各用户波束域统计信道信息更新系统遍历和速率表达式，将空间域预编码器设计问题转化为波束域预编码器设计问题，且波束域预编码器的迭代设计包括如下步骤：
①初始化各用户波束域预编码器，使其满足功率约束；
②利用MM算法框架得到系统遍历和速率在当前迭代下的凸替代函数；
③使用拉格朗日乘子法求解当前迭代的凸问题；
重复步骤②-③直到达到预设迭代次数或预编码收敛，得到各用户最优波束域预编码器。
(3)所述优化目标为最大化系统遍历和速率上界的设计，系统遍历和速率上界为系统遍历和速率利用詹森不等式得到，表达式包括波束矩阵、各用户波束映射矩阵、波束域预编码器和各用户波束域统计信道信息，将空间域预编码器设计问题转化为波束域预编码器设计问题，且波束域预编码器的迭代设计包括如下步骤：
①初始化各用户波束域预编码器，使其满足功率约束；
②利用MM算法框架得到系统遍历和速率上界在当前迭代下的凸替代函数；
③使用拉格朗日乘子法求解当前迭代的凸问题；
重复步骤②-③直到达到预设迭代次数或预编码收敛，得到各用户最优波束域预编码器。
7.根据权利要求1所述的天波大规模MIMO波束结构预编码传输方法，其特征在于，所述根据设计生成的波束域预编码器实施与用户的下行信号传输包括如下步骤：
(1)利用用户波束域预编码器与其发送数据符号相乘生成低维波束域发送信号；
(2)利用波束映射矩阵和低维波束域发送信号相乘得到用户波束域发送信号；
(3)将各用户波束域发送信号进行叠加，得到包含所有用户的波束域发送信号；
(4)利用波束矩阵和所有用户的波束域发送信号相乘生成空间域发送信号。
8.根据权利要求6所述的天波大规模MIMO波束结构预编码传输方法，其特征在于，所述步骤(4)利用Chirp-z变换进行有效实现。
9.天波大规模MIMO波束结构预编码传输系统，包括天波大规模MIMO通信基站和多个用户，其特征在于，所述天波大规模MIMO通信基站包含大规模天线阵列，工作载频为1.6～30MHz的短波波段，基站通过电离层反射的方式与用户进行信号传输。
10.天波大规模MIMO波束结构预编码传输系统，包括基站和多个用户，其特征在于，所述基站实现根据权利要求1-7任一项所述的天波大规模MIMO波束结构预编码传输方法。
CN202210317238.1A 2022-03-29 2022-03-29 天波大规模mimo波束结构预编码传输方法与系统 Pending CN114866117A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210317238.1A CN114866117A (zh) 2022-03-29 2022-03-29 天波大规模mimo波束结构预编码传输方法与系统
PCT/CN2022/086916 WO2023184590A1 (zh) 2022-03-29 2022-04-14 天波大规模mimo波束结构预编码传输方法与系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210317238.1A CN114866117A (zh) 2022-03-29 2022-03-29 天波大规模mimo波束结构预编码传输方法与系统

Publications (1)

Publication Number Publication Date
CN114866117A true CN114866117A (zh) 2022-08-05

Family

ID=82628968

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210317238.1A Pending CN114866117A (zh) 2022-03-29 2022-03-29 天波大规模mimo波束结构预编码传输方法与系统

Country Status (2)

Country Link
CN (1) CN114866117A (zh)
WO (1) WO2023184590A1 (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114172549A (zh) * 2021-12-07 2022-03-11 东南大学 天波大规模mimo通信下行传输方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11394436B2 (en) * 2004-04-02 2022-07-19 Rearden, Llc System and method for distributed antenna wireless communications
CN107888249A (zh) * 2016-09-30 2018-04-06 北京信威通信技术股份有限公司 一种波束域参考信号发送的方法及装置
CN211554155U (zh) * 2019-12-04 2020-09-22 深圳有力科技有限公司 一种短波天线增益测量系统
CN112511201B (zh) * 2020-11-19 2021-10-26 东南大学 天波大规模mimo通信方法及模型和系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114172549A (zh) * 2021-12-07 2022-03-11 东南大学 天波大规模mimo通信下行传输方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XIANGLONG YU等: "Massive MIMO Communication Over HF Skywave Channels", 《2021 IEEE GLOBAL COMMUNICATIONS CONFERENCE (GLOBECOM)》, 2 February 2022 (2022-02-02), pages 1 - 6 *

Also Published As

Publication number Publication date
WO2023184590A1 (zh) 2023-10-05

Similar Documents

Publication Publication Date Title
CN107483088B (zh) 大规模mimo鲁棒预编码传输方法
US11258495B2 (en) Hybrid beamforming method for wireless multi-antenna and frequency-division duplex systems
US8594221B2 (en) Model-based channel estimator for correlated fading channels and channel estimation method thereof
US7385617B2 (en) Methods for multi-user broadband wireless channel estimation
Han et al. FDD massive MIMO based on efficient downlink channel reconstruction
CN107465636B (zh) 一种毫米波大规模阵列空频双宽带系统的信道估计方法
Zhang et al. Channel estimation and training design for hybrid multi-carrier mmwave massive MIMO systems: The beamspace ESPRIT approach
Zhou et al. Active terminal identification, channel estimation, and signal detection for grant-free NOMA-OTFS in LEO satellite Internet-of-Things
US11063724B1 (en) Reduced channel-sounding in MU-MIMO WLANS
Attiah et al. Deep learning for channel sensing and hybrid precoding in TDD massive MIMO OFDM systems
CN110034916B (zh) 基于终端反馈的天线相位同步与信道互易性校准的方法
Kim et al. Two-step approach to time-domain channel estimation for wideband millimeter wave systems with hybrid architecture
CN114172549B (zh) 天波大规模mimo通信下行传输方法
Uwaechia et al. Compressed channel estimation for massive MIMO-OFDM systems over doubly selective channels
CN115065432A (zh) 天波大规模mimo三重波束基信道建模及信道信息获取
CN113193895B (zh) 大规模mimo信道状态信息获取方法、系统及计算机存储介质
Srivastava et al. Sparse Bayesian learning (SBL)-based frequency-selective channel estimation for millimeter wave hybrid MIMO systems
Noh et al. Joint direct and indirect channel estimation for RIS-assisted millimeter-wave systems based on array signal processing
Mollén High-end performance with low-end hardware: Analysis of massive MIMO base station transceivers
CN109787672B (zh) 基于参数学习的大规模mimo格点偏移信道估计方法
Cuji et al. Multi-user Communications for Acoustic OFDM: A Broadband Beamforming Approach
CN114553640B (zh) 多频段大规模mimo系统中的跨频段统计信道状态信息估计方法
CN114866117A (zh) 天波大规模mimo波束结构预编码传输方法与系统
Ji et al. Multipath extraction based UL/DL channel estimation for FDD massive MIMO-OFDM systems
CN114629533A (zh) 大规模mimo信道估计的信息几何方法及系统

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