CN115395991A - 一种非线性多输入多输出信道估计方法和估计系统 - Google Patents

一种非线性多输入多输出信道估计方法和估计系统 Download PDF

Info

Publication number
CN115395991A
CN115395991A CN202210824124.6A CN202210824124A CN115395991A CN 115395991 A CN115395991 A CN 115395991A CN 202210824124 A CN202210824124 A CN 202210824124A CN 115395991 A CN115395991 A CN 115395991A
Authority
CN
China
Prior art keywords
baseband signal
neural network
channel estimation
representing
nonlinear
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
CN202210824124.6A
Other languages
English (en)
Other versions
CN115395991B (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.)
Beijing Information Science and Technology University
Original Assignee
Beijing Information Science and Technology 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 Beijing Information Science and Technology University filed Critical Beijing Information Science and Technology University
Priority to CN202210824124.6A priority Critical patent/CN115395991B/zh
Publication of CN115395991A publication Critical patent/CN115395991A/zh
Application granted granted Critical
Publication of CN115395991B publication Critical patent/CN115395991B/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
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/024—Channel estimation channel estimation algorithms
    • H04L25/0254—Channel estimation channel estimation algorithms using neural network algorithms
    • 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

本公开提出一种非线性多输入多输出信道估计方法和估计系统。该估计方法包括利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将复基带信号的幅度和预设相位作为观测值，进行复基带信号处理；将复基带信号处理后的信息输入到非线性多输入多输出信道中，并对非线性多输入多输出信道进行初始信道估计，得到初始信道估计值；利用GAMP算法，将初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对深度神经网络和卷积神经网络进行训练，得到最终信道估计值。该估计系统包括信号处理模块、初始信道估计模块和神经网络模块。本公开能够增强信道的均方误差性能和鲁棒性。

Description

一种非线性多输入多输出信道估计方法和估计系统
技术领域
本公开涉及无线通信技术领域，尤其涉及一种非线性多输入多输出信道估计方法和估计系统。
背景技术
非线性多输入多输出(Non-Linear Multiple Input Multiple Output， Non-Linear MIMO)通信技术能够将射频信号直接下变频到基带，而传统MIMO的接收机一般是将信号由高频变频到中频，再变频到低频后到基带，二者相比，非线性MIMO简化了信号处理流程，在功耗方面具有明显优势。
目前，对非线性MIMO系统的信道估计算法的研究，主要集中在对信息的传递算法上，但由于接收端相位观测值的均值、方差都没有封闭表达式，在进行采样技术迭代求解时，会导致收敛率变得不稳定，甚至在一些情况下该算法发散问题严重。此外，在实际应用中，信道特性会随着传输环境、工作频率等因素的变化而变化，这些因素通常又是未知的，无法用具体的分析方法来表示。因此，有必要改善上述相关技术方案中存在的一个或者多个问题。
需要说明的是，在上述背景技术部分公开的信息仅用于加强对本公开的背景的理解，因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。
发明内容
本公开实施例的目的在于提供一种，以提高信道的均方误差性能和鲁棒性。
根据本公开实施例的第一方面，提供一种非线性多输入多输出 (Non-LinearMultiple Input Multiple Output，Non-Linear MIMO)信道估计方法，包括：
利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行基带信号处理；
将基带信号处理后的信息输入到非线性多输入多输出信道中，并对所述非线性多输入多输出信道进行初始信道估计，得到初始信道估计值；
利用GAMP算法，将所述初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对所述深度神经网络和所述卷积神经网络进行训练，得到最终信道估计值。
本公开的一示例性实施例中，所述利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行基带信号处理的步骤中，所述调制解调器包括同相正交调制解调器；利用包络检波器获取所述复基带信号的幅度，利用相位检测器获取所述复基带信号的预设相位，所述预设相位为π相位或2π相位。
本公开的一示例性实施例中，所述利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行基带信号处理的步骤包括：
构建广义线性混合模型，所述广义线性混合模型包括：
Figure RE-GDA0003916870340000031
其中，z＝As表示噪声影响前的结果，
Figure RE-GDA0003916870340000032
表示待恢复的 M×N维度的观测矩阵，
Figure RE-GDA0003916870340000033
表示作为复基带信号的幅度和π相位的观测值，w表示加性复高斯噪声
Figure RE-GDA0003916870340000034
Figure RE-GDA0003916870340000035
表示服从均值为0，方差为
Figure RE-GDA0003916870340000036
的复高斯分布；
Figure RE-GDA0003916870340000037
表示分量算子，第m个元素gm的定义包括：
Figure RE-GDA0003916870340000038
其中，1(·)表示对复基带信号不做处理，∠π(·)表示复基带信号进行了π相位移动处理，T为矩阵转置。
本公开的一示例性实施例中，将基带信号处理后的信息输入到非线性多输入多输出信道中，并对所述非线性多输入多输出信道进行初步信道估计，得到初步信道估计值的步骤包括：
将复基带信号处理后的信息初始化，并利用先验模型和似然函数，得到先验概率，根据所述先验概率得到后验概率；
对信息输出通道进行线性处理和非线性处理，同时对信息输入通道进行线性处理和非线性处理；
重复对输出通道和输入通道进行线性和非线性处理，通过不断迭代过程，得到初始信道估计值。
本公开的一示例性实施例中，所述将经基带信号处理后的信息初始化，并利用似然函数，得到先验概率，根据所述先验概率得到后验概率的步骤包括：
所述先验概率包括：
Figure RE-GDA0003916870340000041
根据所述先验概率得到所述后验概率包括：
Figure RE-GDA0003916870340000042
其中，zk表示z＝As的第k个值，z＝As表示噪声影响前的结果，
Figure RE-GDA0003916870340000043
RF链为射频链，
Figure RE-GDA0003916870340000044
表示以
Figure RE-GDA0003916870340000045
为变量，均值为zk，方差为
Figure RE-GDA0003916870340000046
的复高斯分布。
本公开的一示例性实施例中，所述对信息输出通道进行线性处理和非线性处理，同时对信息输入通道进行线性处理和非线性处理得过程包括：
进行输出通道的线性处理，计算zk的输出通道的方差公式包括：
Figure RE-GDA0003916870340000047
计算zk的输出通道的估计值公式包括：
Figure RE-GDA0003916870340000048
其中，Akn表示信道传输中的矩阵；
Figure RE-GDA0003916870340000049
表示第t次迭代s向量中第n个元素的方差列的值；
Figure RE-GDA00039168703400000410
表示第t次迭代P向量中第k个元素的方差；uk(t-1)表示中间变量；
Figure RE-GDA00039168703400000411
表示t时刻sn的均值；
Figure RE-GDA00039168703400000412
表示t时刻pk的估计值；
进行输出通道的非线性处理包括：
Figure RE-GDA0003916870340000051
Figure RE-GDA0003916870340000052
其中，uk(t)是通过后验概率
Figure RE-GDA0003916870340000053
估计的期望，其中
Figure RE-GDA0003916870340000054
是概率归一化常数，
Figure RE-GDA0003916870340000055
是以zk为变量，均值为
Figure RE-GDA0003916870340000056
方差为
Figure RE-GDA0003916870340000057
的复高斯分布；
Figure RE-GDA0003916870340000058
是通过后验概率
Figure RE-GDA0003916870340000059
估计的方差，
Figure RE-GDA00039168703400000510
Figure RE-GDA00039168703400000511
Figure RE-GDA00039168703400000512
表示来自功率分路器和传统射频链的方差；
进行输入通道的线性处理，计算zk的输入通道的方差公式包括：
Figure RE-GDA00039168703400000513
计算zk的输入通道的均值公式包括
Figure RE-GDA00039168703400000514
其中，
Figure RE-GDA00039168703400000515
表示第t次迭代α向量中第k个元素的方差列的值；
Figure RE-GDA00039168703400000516
表示第t次迭代r向量中第n个元素的方差；uk(t)表示中间变量；
Figure RE-GDA00039168703400000517
表示t时刻sn的均值；
Figure RE-GDA00039168703400000518
表示t时刻rn的估计值；
进行输入通道的非线性处理包括：
Figure RE-GDA00039168703400000519
Figure RE-GDA00039168703400000520
其中，
Figure RE-GDA00039168703400000521
表示t+1个时刻的sn的均值，
Figure RE-GDA00039168703400000522
表示在
Figure RE-GDA0003916870340000061
和
Figure RE-GDA0003916870340000062
已知时sn的期望；
Figure RE-GDA0003916870340000063
表示t+1个时刻的sn的方差，
Figure RE-GDA0003916870340000064
表示在
Figure RE-GDA0003916870340000065
和
Figure RE-GDA0003916870340000066
已知时sn的方差。
本公开的一示例性实施例中，所述重复对输出通道和输入通道进行线性和非线性处理，通过不断迭代过程，得到初始信道估计值。
本公开的一示例性实施例中，所述利用GAMP算法，将所述初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对所述深度神经网络和所述卷积神经网络进行训练，得到最终信道估计值的步骤包括：
所述深度神经网络和所述卷积神经网络分别均包括输入层、隐含层和输出层，其中，所述深度神经网络中的隐含层和输出层均为全连接层；卷积神经网络中的输出层为全连接层；
所述全连接层的神经元的计算公式包括：
Figure RE-GDA0003916870340000067
其中，yl表示第l层神经元的值；wkl表示权重；bl表示偏差； xk表示前一层神经元的输出值；f表示激活函数，f＝max(0,wTx+b)；
所述输出层的损失函数的计算公式包括：
Figure RE-GDA0003916870340000068
其中，MSE表示离线训练过程中的均方误差；yd表示所述深度神经网络和所述卷积神经网络的输出；xd表示标签，n表示样本总数。
本公开的一示例性实施例中，还包括：利用RMSProp算法对深度神经网络的参数进行更新，利用Adam算法对卷积神经网络的参数进行更新；
所述RMSProp算法包括利用加权和偏差梯度的差分平方加权平均值，得到前一次迭代过程中损失函数累计的梯度动量，所述梯度动量包括：
Figure RE-GDA0003916870340000071
其中，σ表示梯度下降的衰减因子；dw和db表示梯度；w和b 表示更新参数，
Figure RE-GDA0003916870340000072
η表示学习率；μ表示累积变量；
所述Adam算法包括计算梯度的指数移动均值，包括：
Figure RE-GDA0003916870340000073
其中，Vdw表示梯度的一阶矩；Sdw表示梯度的二阶矩；
Figure RE-GDA0003916870340000074
表示Vdw的偏置矫正；
Figure RE-GDA0003916870340000075
表示Sdw的偏置矫正；
对所述Adam算法进行偏差修正得到包括：
Figure RE-GDA0003916870340000076
更新权重得到包括：
Figure RE-GDA0003916870340000081
其中，α表示学习率或步长因子；β1表示一阶矩估计的指数衰减率、β2表示二阶矩估计的指数衰减率；ε表示为数值很小的超参数。
本公开实施例的第二方面，提供了一种非线性多输入多输出信道估计系统。包括：
信号处理模块，所述信号处理模块包括利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行复基带信号处理；
初始信道估计模块，所述初始信道估计模块将经基带信号处理后的信息输入到非线性MIMO信道中，并对所述非线性MIMO信道进行初始信道估计，得到初始信道估计值；
神经网络模块，所述神经网络模块包括利用GAMP算法，将所述初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对所述深度神经网络和所述卷积神经网络进行训练，得到最终信道估计值。
本公开提供的技术方案可以包括以下有益效果：
本公开实施例提出一种非线性多输入多输出信道估计方法和估计系统，利用GAMP算法，将产生的初步信道估计值作为深度学习的输入，通过深度神经网络和卷积神经网络进一步细化估计值，以获得更为精确的信道估计结果。本公开的信道估计方法和估计系统具有更为优异的均方误差性能和更强的鲁棒性。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本公开。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理。显而易见的，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出本公开示例性实施例非线性MIMO信道估计方法的步骤示意图；
图2示出本公开示例性实施例非线性MIMO信道估计系统的结构示意图；
图3示出本公开示例性实施例非线性MIMO信道估计方法的流程图；
图4示出仿真实验中在不同信道估计算法下的均方误差曲线；
图5示出仿真实验中本公开的非线性MIMO信道估计方法在不同天线距离下的均方误差曲线；
图6示出仿真实验中本公开的非线性MIMO信道估计方法在不同天线数量时的均方误差曲线。
具体实施方式
现在将参考附图更全面地描述示例实施方式。然而，示例实施方式能够以多种形式实施，且不应被理解为限于在此阐述的范例；相反，提供这些实施方式使得本公开将更加全面和完整，并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。
本示例实施方式的第一方面提供了一种非线性MIMO信道估计方法，该信道估计方法包括以下步骤：
步骤S101：利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行基带信号处理；
步骤S102：将基带信号处理后的信息输入到非线性MIMO信道中，并对所述非线性MIMO信道进行初始信道估计，得到初始信道估计值；
步骤S103：利用GAMP算法，将所述初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对所述深度神经网络和所述卷积神经网络进行训练，得到最终信道估计值。
下面，将对本示例实施方式中的上述方法的各个步骤进行更详细的说明。
在步骤S101中，这里引入了低精度同相正交(Inphase Quadradure，IQ)调制解调器，将信道接收到的射频信号直接下变频到基带，根据基带信号获取复基带信号，采用包络检波器获取复基带信号的幅度，采用相位检测器获取复基带信号的预设相位，本示例实施方式中的预设相位为π相位，也可以是2π相位。然后将获取到的复基带信号的幅度和预设相位作为观测值，进行基带信号处理。
这里用到了广义近似消息传递算法(GAMP算法)，将非线性 MIMO的信道估计问题归类为广义线性混合问题，即由非线性 MIMO构建出广义线性混合模型，该模型的计算公式包括：
Figure RE-GDA0003916870340000111
其中，z＝As表示噪声影响前的结果，
Figure RE-GDA0003916870340000112
表示待恢复的 M×N维度的观测矩阵，
Figure RE-GDA0003916870340000113
表示作为复基带信号的幅度和π相位的观测值，w表示加性复高斯噪声
Figure RE-GDA0003916870340000114
Figure RE-GDA0003916870340000115
表示服从均值为0，方差为
Figure RE-GDA0003916870340000116
的复高斯分布；
Figure RE-GDA0003916870340000117
表示分量算子，第m个元素gm的定义包括：
Figure RE-GDA0003916870340000118
其中，1(·)表示对复基带信号不做处理，∠π(·)表示复基带信号进行了π相位移动处理，T为矩阵转置。
在步骤S102中，具体包括以下步骤：
步骤S1021：将基带信号处理后的信息进行初始化处理，并利用似然函数，得到先验概率，再根据先验概率得到后验概率；
这里，先验概率包括：
Figure RE-GDA0003916870340000119
根据所述先验概率得到所述后验概率包括：
Figure RE-GDA00039168703400001110
其中，其中，zk表示z＝As的第k个值，z＝As表示噪声影响前的结果，
Figure RE-GDA0003916870340000121
RF链为射频链，
Figure RE-GDA0003916870340000122
表示以
Figure RE-GDA0003916870340000123
为变量，均值为zk，方差为
Figure RE-GDA0003916870340000124
的复高斯分布。
步骤S1022：所述对信息输出通道进行线性处理和非线性处理，同时对信息输入通道进行线性处理和非线性处理得过程包括：
进行输出通道的线性处理，计算zk的输出通道的方差公式包括：
Figure RE-GDA0003916870340000125
计算zk的输出通道的估计值公式包括：
Figure RE-GDA0003916870340000126
其中，Akn表示信道传输中的矩阵；
Figure RE-GDA0003916870340000127
表示第t次迭代s向量中第n个元素的方差列的值；
Figure RE-GDA0003916870340000128
表示第t次迭代P向量中第k个元素的方差；uk(t-1)表示中间变量；
Figure RE-GDA0003916870340000129
表示t时刻sn的均值；
Figure RE-GDA00039168703400001210
表示t时刻pk的估计值；
进行输出通道的非线性处理包括：
Figure RE-GDA00039168703400001211
Figure RE-GDA00039168703400001212
其中，uk(t)是通过后验概率
Figure RE-GDA00039168703400001213
估计的期望，其中
Figure RE-GDA00039168703400001214
是概率归一化常数，
Figure RE-GDA00039168703400001215
是以zk为变量，均值为
Figure RE-GDA00039168703400001216
方差为
Figure RE-GDA00039168703400001217
的复高斯分布；
Figure RE-GDA00039168703400001218
是通过后验概率
Figure RE-GDA00039168703400001219
估计的方差，
Figure RE-GDA00039168703400001220
Figure RE-GDA0003916870340000131
σ2 conv表示来自功率分路器和传统射频链的方差；
进行输入通道的线性处理，计算zk的输入通道的方差公式包括：
Figure RE-GDA0003916870340000133
计算zk的输入通道的均值公式包括
Figure RE-GDA0003916870340000134
其中，
Figure RE-GDA0003916870340000135
表示第t次迭代α向量中第k个元素的方差列的值；
Figure RE-GDA0003916870340000136
表示第t次迭代r向量中第n个元素的方差；uk(t)表示中间变量；
Figure RE-GDA0003916870340000137
表示t时刻sn的均值；
Figure RE-GDA0003916870340000138
表示t时刻rn的估计值；
进行输入通道的非线性处理包括：
Figure RE-GDA0003916870340000139
Figure RE-GDA00039168703400001310
其中，
Figure RE-GDA00039168703400001311
表示t+1个时刻的sn的均值，
Figure RE-GDA00039168703400001312
表示在
Figure RE-GDA00039168703400001313
和
Figure RE-GDA00039168703400001314
已知时sn的期望；
Figure RE-GDA00039168703400001315
表示t+1个时刻的sn的方差，
Figure RE-GDA00039168703400001316
表示在
Figure RE-GDA00039168703400001317
和
Figure RE-GDA00039168703400001318
已知时sn的方差；
对后验概率密度
Figure RE-GDA00039168703400001319
求得的期望，其中
Figure RE-GDA00039168703400001320
是概率归一化常数，ps(sn)表示sn的先验分布，
Figure RE-GDA00039168703400001321
表示以sn为变量的均值为
Figure RE-GDA00039168703400001322
方差为
Figure RE-GDA00039168703400001323
的复高斯分布。
步骤S1023：重复对输出通道和输入通道进行线性和非线性处理，通过不断迭代过程，得到初始信道估计值。
在步骤S103中，利用广义近似消息传递算法(GAMP算法)，将步骤S102中得到的初始信道估计值作为深度神经网络(Deep Neural Network，DNN)和卷积神经网络(Convolutional Neural Network,CNN)的输入，将信道的实际响应作为标签，对DNN和 CNN进行离线训练，从而得到最终信道估计道。
参照图3，本示例实施方式中，该DNN的结构包含六层，包括一个输入层、四个隐含层和一个输出层；其中，DNN的隐含层和输出层均为全连接层，用于提取最后一层的所有特征。该CNN的结构也包含六层，包括一个输入层、四个隐含层和一个输出层；其中， CNN的输出层为全连接层。全连接层的每个神经元都可以表示为前一层中所有神经元输出值的加权和，公式包括：
Figure RE-GDA0003916870340000141
其中，yl表示第l层神经元的值；wkl表示权重；bl表示偏差； xk表示前一层神经元的输出值；f表示激活函数。这里在每个隐含层之后使用线性整流函数(Linearrectificationfunction,ReLU)作为激活函数，用于对前一层的加权输出值进行非线性决策，可以表示为：f＝max(0,wTx+b)。
在DNN和CNN的结构中使用回归层作为输出层。初始信道估计值由DAMP算法经过信道估计产生，将初始信道估计值H作为 DNN和CNN的输入，可以有多种方法来描述DNN和CNN的输出与发送数据之间的误差，在离线训练过程中，选择采用均方误差 (Mean SquareError，MSE)作为损失函数。其计算公式为：
Figure RE-GDA0003916870340000142
其中，L表示损失函数；MSE表示离线训练过程中的均方误差； yd表示所述深度神经网络和所述卷积神经网络的输出；xd表示标签， n表示样本总数。
该损失函数L就是最终信道估计值。
为了得到更精确的信道估计，需要不断重复上面论述的步骤 S101和步骤S102，这样可以得到不断更新的初始信道估计值H，将H不断作为输入从而更新为更为精确的最终信道估计值。
具体的，使用均方根传播(Root Mean Square Prop，RMSprop) 来更新DNN的参数。RMSprop利用加权和偏差梯度的差分平方加权平均值，在前一次迭代过程中损失函数累计的梯度动量，该梯度动量可以包括：
Figure RE-GDA0003916870340000151
其中，σ表示梯度动量下降的衰减因子；dw和db表示梯度；
RMSprop
更新参数w和b可以表示为：
Figure RE-GDA0003916870340000152
η表示学习率；μ当的值很小时，可以防止分母为零。
由于RMSprop是一种通过调整更新参数的步长来加速梯度下降的算法。在剧烈变化的方向上的速度较小，而在平缓变化的方向上速度较大，可以优化更新过度摆动问题的损失函数。因此，可以设置较大的学习率，来加快算法的收敛速度。
使用自适应动量算法(Adaptive momentum,Adam)来更新CNN 的参数。一般的随机梯度下降算法是在随机梯度下降时，保持单一的学习率更新所有的权重，学习率在训练过程中并不会改变；而 Adam算法通过计算梯度的一阶矩估计和二阶矩估计能为不同的参数涉及独立的自适应学习率。Adam算法可以看成是两种随机梯度下降扩展式的优点集合，即:(1)适应性梯度算法(AdaGrad)为每一个参数保留一个学习率以提升在稀疏梯度(即自然语言和计算机视觉问题)上的性能。(2)均方根传播(RMSprop)基于权重梯度最近量级的均值为每一个参数适应性的保留学习率。这意味着算法在非稳态和在线问题上有很多优秀的性能。
Adam算法同时获得了AdaGrad和RMSprop算法的优点。Adam 不仅如RMSprop算法那样基于一阶矩均值计算适应性参数学习率，它同时还充分利用了梯度的二阶矩均值。具体来说，该算法计算了梯度的指数移动均值，超参数β1和β2控制了这些移动均值的衰减率。该梯度的指数移动均值包括：
Figure RE-GDA0003916870340000161
其中，Vdw表示梯度的一阶矩；Sdw表示梯度的二阶矩；
Figure RE-GDA0003916870340000162
表示Vdw的偏置矫正；
Figure RE-GDA0003916870340000163
表示Sdw的偏置矫正；
对所述Adam算法进行偏差修正得到包括：
Vdw＝β1Vdw+(1-β1)dw (16)
Sdw＝β2Sdw+(1-β2)dw2
更新权重得到包括：
Figure RE-GDA0003916870340000171
其中，α为学习率或步长因子，它控制了权重的更新比率(如 0.001)，较大的值(如0.3)在学习率更新前会有更快的初始学习能力，而较小的值(如1.0E-5)会令训练收敛到更好的性能；β1为一阶矩估计的指数衰减率、β2为二阶矩估计的指数衰减率；ε表示为数值很小的超参数，目的是为了防止分母为0；
本示例实施方式的第二方面提供了一种非线性MIMO信道估计系统，该信道估计系统包括：
信号处理模块，所述信号处理模块包括利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行基带信号处理；
初始信道估计模块，所述初始信道估计模块将经基带信号处理后的信息输入到非线性MIMO信道中，并对所述非线性MIMO信道进行初始信道估计，得到初始信道估计值；
神经网络模块，所述神经网络模块包括利用GAMP算法，将所述初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对所述深度神经网络和所述卷积神经网络进行训练，得到最终信道估计值。
仿真实验
参照图4，在用户密度为(100，5)、基站天线数为8、用户数为4、天线距离为1.5λ的场景下，将本公开提出的非线性MIMO信道估计方法(记为A方法DNN-HPO)与B方法和C方法进行均方误差性能对比。PO-MIMO为纯相位MIMO。B方法是指基于GAMP 算法的非线性MIMO信道估计方法(HPO-MIMO)；C方法为传统 MIMO信道估计方法(Conv-MIMO)。由图4可以看到，对比结果显示了本公开提出的A方法的均方误差性能明显优于其他两种方法。参照图5，对本公开提出的非线性MIMO信道估计方法在不同天线距离情况下的均方误差性能，可以看到，当用户密度为(100，5)、基站天线数为32、用户数为4的场景下对比了5、15、30三种天线距离的均方误差曲线，由图5可以看到，但个信噪比增加时，本公开提出的非线性MIMO信道估计方法的均方误差性能随着天线距离的增加而增加。
参照图6，对本公开提出的非线性MIMO信道估计方法在不同天线数量情况下的均方误差性能。可以看到，应用的GAMP算法在 32x8场景时的均方误差性能要优于64x8，而且优于B方法和C方法的均方误差性能。
综上所述，本公开的实施例中提出了一种非线性MIMO信道估计方法和估计系统。该估计方法和估计系统利用GAMP算法，将产生的初步信道估计值作为深入学习的输入，应用深度神经网络和卷积神经网络做进一步的信道估计，从而获得更精确的最终信道估计值作为结果。并且作为一种模型驱动方案，所提出的DNN模型和CNN模型的输入可以直接从初步信道估计值中得到。本公开相较于现有算法具有更好的均方误差性能和更强的鲁棒性。
关于上述实施例中的系统，其中各个单元执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
应当注意，尽管在上文详细描述中提及了用于动作执行的系统的若干单元，但是这种划分并非强制性的。实际上，根据本公开的实施方式，上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。反之，上文描述的一个单元的特征和功能可以进一步划分为由多个单元来具体化。可以根据实际的需要选择其中的部分或者全部单元来实现本公开方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本公开的其它实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本公开的真正范围和精神由所附的权利要求指出。

Claims (10)

1.一种非线性多输入多输出信道估计方法，其特征在于，包括：
利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行基带信号处理；
将基带信号处理后的信息输入到非线性多输入多输出信道中，并对所述非线性多输入多输出信道进行初始信道估计，得到初始信道估计值；
利用GAMP算法，将所述初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对所述深度神经网络和所述卷积神经网络进行训练，得到最终信道估计值。
2.根据权利要求1所述非线性多输入多输出信道估计方法，其特征在于，所述利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行基带信号处理的步骤中，所述调制解调器包括同相正交调制解调器；利用包络检波器获取所述复基带信号的幅度，利用相位检测器获取所述复基带信号的预设相位，所述预设相位为π相位或2π相位。
3.根据权利要求2所述非线性多输入多输出信道估计方法，其特征在于，所述利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行基带信号处理的步骤包括：
构建广义线性混合模型，所述广义线性混合模型包括：
Figure RE-FDA0003916870330000021
其中，z＝As表示噪声影响前的结果，
Figure RE-FDA0003916870330000022
表示待恢复的M×N维度的观测矩阵，
Figure RE-FDA0003916870330000023
表示作为复基带信号的幅度和π相位的观测值，w表示加性复高斯噪声
Figure RE-FDA0003916870330000024
Figure RE-FDA0003916870330000025
表示服从均值为0，方差为
Figure RE-FDA0003916870330000026
的复高斯分布；
Figure RE-FDA0003916870330000027
表示分量算子，第m个元素gm的定义包括：
Figure RE-FDA0003916870330000028
其中，1(·)表示对复基带信号不做处理，∠π(·)表示复基带信号进行了π相位移动处理，T为矩阵转置。
4.根据权利要求1所述非线性多输入多输出信道估计方法，其特征在于，将基带信号处理后的信息输入到非线性多输入多输出信道中，并对所述非线性多输入多输出信道进行初步信道估计，得到初步信道估计值的步骤包括：
将复基带信号处理后的信息初始化，并利用先验模型和似然函数，得到先验概率，根据所述先验概率得到后验概率；
对信息输出通道进行线性处理和非线性处理，同时对信息输入通道进行线性处理和非线性处理；
重复对输出通道和输入通道进行线性和非线性处理，通过不断迭代过程，得到初始信道估计值。
5.根据权利要求4所述非线性多输入多输出信道估计方法，其特征在于，所述将经基带信号处理后的信息初始化，并利用似然函数，得到先验概率，根据所述先验概率得到后验概率的步骤包括：
所述先验概率包括：
Figure RE-FDA0003916870330000031
根据所述先验概率得到所述后验概率包括：
Figure RE-FDA0003916870330000032
其中，zk表示z＝As的第k个值，z＝As表示噪声影响前的结果，
Figure RE-FDA0003916870330000033
RF链为射频链，
Figure RE-FDA0003916870330000034
表示以
Figure RE-FDA0003916870330000035
为变量，均值为zk，方差为
Figure RE-FDA0003916870330000036
的复高斯分布。
6.根据权利要求5所述非线性多输入多输出信道估计方法，其特征在于，所述对信息输出通道进行线性处理和非线性处理，同时对信息输入通道进行线性处理和非线性处理得过程包括：
进行输出通道的线性处理，计算zk的输出通道的方差公式包括：
Figure RE-FDA0003916870330000037
计算zk的输出通道的估计值公式包括：
Figure RE-FDA0003916870330000038
其中，Akn表示信道传输中的矩阵；
Figure RE-FDA0003916870330000039
表示第t次迭代s向量中第n个元素的方差列的值；
Figure RE-FDA00039168703300000310
表示第t次迭代P向量中第k个元素的方差；uk(t-1)表示中间变量；
Figure RE-FDA00039168703300000311
表示t时刻sn的均值；
Figure RE-FDA00039168703300000312
表示t时刻pk的估计值；
进行输出通道的非线性处理包括：
Figure RE-FDA0003916870330000041
Figure RE-FDA0003916870330000042
其中，uk(t)是通过后验概率
Figure RE-FDA0003916870330000043
估计的期望，其中
Figure RE-FDA0003916870330000044
是概率归一化常数，
Figure RE-FDA0003916870330000045
是以zk为变量，均值为
Figure RE-FDA0003916870330000046
方差为
Figure RE-FDA0003916870330000047
的复高斯分布；
Figure RE-FDA0003916870330000048
是通过后验概率
Figure RE-FDA0003916870330000049
估计的方差，
Figure RE-FDA00039168703300000410
Figure RE-FDA00039168703300000411
σ2 conv表示来自功率分路器和传统射频链的方差；
进行输入通道的线性处理，计算zk的输入通道的方差公式包括：
Figure RE-FDA00039168703300000412
计算zk的输入通道的均值公式包括：
Figure RE-FDA00039168703300000413
其中，
Figure RE-FDA00039168703300000414
表示第t次迭代α向量中第k个元素的方差列的值；
Figure RE-FDA00039168703300000415
表示第t次迭代r向量中第n个元素的方差；uk(t)表示中间变量；
Figure RE-FDA00039168703300000416
表示t时刻sn的均值；
Figure RE-FDA00039168703300000417
表示t时刻rn的估计值；
进行输入通道的非线性处理包括：
Figure RE-FDA00039168703300000418
Figure RE-FDA00039168703300000419
其中，
Figure RE-FDA00039168703300000420
表示t+1个时刻的sn的均值，
Figure RE-FDA00039168703300000421
表示在
Figure RE-FDA00039168703300000422
和
Figure RE-FDA00039168703300000423
已知时sn的期望；
Figure RE-FDA00039168703300000424
表示t+1个时刻的sn的方差，
Figure RE-FDA0003916870330000051
表示在
Figure RE-FDA0003916870330000052
和
Figure RE-FDA0003916870330000053
已知时sn的方差。
7.根据权利要求6所述非线性多输入多输出信道估计方法，其特征在于，所述重复对输出通道和输入通道进行线性和非线性处理，通过不断迭代过程，得到初始信道估计值。
8.根据权利要求1所述非线性多输入多输出信道估计方法，其特征在于，所述利用DAMP算法，将所述初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对所述深度神经网络和所述卷积神经网络进行训练，得到最终信道估计值的步骤包括：
所述深度神经网络和所述卷积神经网络分别均包括输入层、隐含层和输出层，其中，所述深度神经网络中的隐含层和输出层均为全连接层；卷积神经网络中的输出层为全连接层；
所述全连接层的神经元的计算公式包括：
Figure RE-FDA0003916870330000054
其中，yl表示第l层神经元的值；wkl表示权重；bl表示偏差；xk表示前一层神经元的输出值；f表示激活函数，f＝max(0,wTx+b)；
所述输出层的损失函数的计算公式包括：
Figure RE-FDA0003916870330000055
其中，MSE表示离线训练过程中的均方误差；yd表示所述深度神经网络和所述卷积神经网络的输出；xd表示标签，n表示样本总数。
9.根据权利要求1所述非线性多输入多输出信道估计方法，其特征在于，还包括：利用RMSProp算法对深度神经网络的参数进行更新，利用Adam算法对卷积神经网络的参数进行更新；
所述RMSProp算法包括利用加权和偏差梯度的差分平方加权平均值，得到前一次迭代过程中损失函数累计的梯度动量，所述梯度动量包括：
Figure RE-FDA0003916870330000061
其中，σ表示梯度下降的衰减因子；dw和db表示梯度；w和b表示更新参数，
Figure RE-FDA0003916870330000062
η表示学习率；μ表示累积变量；
所述Adam算法包括计算梯度的指数移动均值，包括：
Figure RE-FDA0003916870330000063
其中，Vdw表示梯度的一阶矩；Sdw表示梯度的二阶矩；
Figure RE-FDA0003916870330000064
表示Vdw的偏置矫正；
Figure RE-FDA0003916870330000065
表示Sdw的偏置矫正；
对所述Adam算法进行偏差修正得到包括：
Figure RE-FDA0003916870330000066
更新权重得到包括：
Figure RE-FDA0003916870330000071
其中，α表示学习率或步长因子；β1表示一阶矩估计的指数衰减率、β2表示二阶矩估计的指数衰减率；ε表示为数值很小的超参数。
10.一种非线性多输入多输出信道估计系统，其特征在于，包括：
信号处理模块，所述信号处理模块包括利用调制解调器，将接收到的射频信号下变频到基带，根据基带信号获取复基带信号，并将所述复基带信号的幅度和预设相位作为观测值，进行复基带信号处理；
初始信道估计模块，所述初始信道估计模块将经复基带信号处理后的信息输入到非线性多输出多输出信道中，并对所述非线性多输入多输出信道进行初始信道估计，得到初始信道估计值；
神经网络模块，所述神经网络模块包括利用GAMP算法，将所述初始信道估计值作为深度神经网络和卷积神经网络的输入，将信道的实际响应作为标签，对所述深度神经网络和所述卷积神经网络进行训练，得到最终信道估计值。
CN202210824124.6A 2022-07-13 2022-07-13 一种非线性多输入多输出信道估计方法和估计系统 Active CN115395991B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210824124.6A CN115395991B (zh) 2022-07-13 2022-07-13 一种非线性多输入多输出信道估计方法和估计系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210824124.6A CN115395991B (zh) 2022-07-13 2022-07-13 一种非线性多输入多输出信道估计方法和估计系统

Publications (2)

Publication Number Publication Date
CN115395991A true CN115395991A (zh) 2022-11-25
CN115395991B CN115395991B (zh) 2023-08-25

Family

ID=84115985

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210824124.6A Active CN115395991B (zh) 2022-07-13 2022-07-13 一种非线性多输入多输出信道估计方法和估计系统

Country Status (1)

Country Link
CN (1) CN115395991B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104022978A (zh) * 2014-06-18 2014-09-03 中国联合网络通信集团有限公司 一种半盲信道估计方法和系统
DE102013110287B3 (de) * 2013-09-18 2014-10-16 Intel IP Corporation Verfahren und Mobilvorrichtung für die Kanalqualitätsschätzung für nichtlineare MIMO-Empfänger
CN109714086A (zh) * 2019-01-23 2019-05-03 上海大学 基于深度学习的优化mimo检测方法
WO2020253690A1 (zh) * 2019-06-17 2020-12-24 浙江大学 一种基于近似消息传递算法的深度学习波束域信道估计方法
CN112242969A (zh) * 2020-10-19 2021-01-19 南京爱而赢科技有限公司 一种基于模型驱动的深度学习的新型单比特ofdm接收机
CN112637094A (zh) * 2020-12-17 2021-04-09 南京爱而赢科技有限公司 一种基于模型驱动深度学习的多用户mimo接收方法
CN112737985A (zh) * 2020-12-25 2021-04-30 东南大学 基于深度学习的大规模mimo信道联合估计和反馈方法
CN113472706A (zh) * 2021-07-12 2021-10-01 南京大学 一种基于深度神经网络的mimo-ofdm系统信道估计方法
WO2021228743A1 (en) * 2020-05-14 2021-11-18 Nokia Technologies Oy Channel estimation for an antenna array
EP3994811A1 (en) * 2019-07-03 2022-05-11 Nokia Technologies Oy Transmission system with channel estimation based on a neural network

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102013110287B3 (de) * 2013-09-18 2014-10-16 Intel IP Corporation Verfahren und Mobilvorrichtung für die Kanalqualitätsschätzung für nichtlineare MIMO-Empfänger
CN104022978A (zh) * 2014-06-18 2014-09-03 中国联合网络通信集团有限公司 一种半盲信道估计方法和系统
CN109714086A (zh) * 2019-01-23 2019-05-03 上海大学 基于深度学习的优化mimo检测方法
WO2020253690A1 (zh) * 2019-06-17 2020-12-24 浙江大学 一种基于近似消息传递算法的深度学习波束域信道估计方法
EP3994811A1 (en) * 2019-07-03 2022-05-11 Nokia Technologies Oy Transmission system with channel estimation based on a neural network
WO2021228743A1 (en) * 2020-05-14 2021-11-18 Nokia Technologies Oy Channel estimation for an antenna array
CN112242969A (zh) * 2020-10-19 2021-01-19 南京爱而赢科技有限公司 一种基于模型驱动的深度学习的新型单比特ofdm接收机
CN112637094A (zh) * 2020-12-17 2021-04-09 南京爱而赢科技有限公司 一种基于模型驱动深度学习的多用户mimo接收方法
CN112737985A (zh) * 2020-12-25 2021-04-30 东南大学 基于深度学习的大规模mimo信道联合估计和反馈方法
CN113472706A (zh) * 2021-07-12 2021-10-01 南京大学 一种基于深度神经网络的mimo-ofdm系统信道估计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YI GONG ET AL.: "Nonlinear MIMO for Industrial Internet of Things in Cyber-Physical Systems", 《IEEE TRANSCTION ON INDUSTRIAL INFORMATICS》 *
巩译: "面向车联网的非线性MIMO高能效理论和关键技术研究", 《中国博士学位论文全文数据库（电子期刊）》, pages 5 *
张体操: "低成本大规模MIMO通信系统中的多用户检测技术研究", 《硕士论文》 *

Also Published As

Publication number Publication date
CN115395991B (zh) 2023-08-25

Similar Documents

Publication Publication Date Title
Liao et al. A model-driven deep learning method for massive MIMO detection
CN111181619B (zh) 基于深度强化学习的毫米波混合波束成形设计方法
CN114051701B (zh) 用于机器学习辅助预编码的设备和方法
CN112637094A (zh) 一种基于模型驱动深度学习的多用户mimo接收方法
CN113162665B (zh) 一种基于深度学习信道预测的预编码方法
CN113541747B (zh) 一种大规模mimo检测方法、装置和存储介质
Kumar et al. DLNet: Deep learning-aided massive MIMO decoder
Liao et al. Deep learning aided low complex sphere decoding for MIMO detection
Nguyen et al. Deep learning for estimation and pilot signal design in few-bit massive MIMO systems
CN114499601A (zh) 一种基于深度学习的大规模mimo信号检测方法
Yang et al. Dynamic neural network for MIMO detection
Arjoune et al. Experience-driven learning-based intelligent hybrid beamforming for massive MIMO mmWave communications
Long et al. Broad learning based hybrid beamforming for mm-wave MIMO in time-varying environments
Yang et al. MIMO detector selection with federated learning
CN114938232B (zh) 基于lstm的同时同频全双工数字域自干扰抑制方法
CN115395991A (zh) 一种非线性多输入多输出信道估计方法和估计系统
Yang et al. Federated dynamic neural network for deep MIMO detection
Wang et al. Deep Transfer Learning for Model-Driven Signal Detection in MIMO-NOMA Systems
Long et al. Few-shot learning based hybrid beamforming under birth-death process of scattering paths
CN110868244B (zh) 一种基于信道穿刺的低复杂度通信信号检测方法
CN117060952A (zh) 多输入多输出系统中的信号检测方法及装置
CN113037409B (zh) 基于深度学习的大规模mimo系统信号检测方法
Hsu et al. Deep Learning Based Adaptive Hybrid Beamforming for mmWave MIMO Systems
CN113630160B (zh) 一种大规模mimo检测方法、装置、设备和存储介质
CN113438190B (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
GR01 Patent grant
GR01 Patent grant