CN112911610A - 一种面向无线通信网络优化的信道波束图样拓展方法 - Google Patents

一种面向无线通信网络优化的信道波束图样拓展方法 Download PDF

Info

Publication number
CN112911610A
CN112911610A CN202110102629.7A CN202110102629A CN112911610A CN 112911610 A CN112911610 A CN 112911610A CN 202110102629 A CN202110102629 A CN 202110102629A CN 112911610 A CN112911610 A CN 112911610A
Authority
CN
China
Prior art keywords
channel beam
beam pattern
gaussian process
neural network
equal
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
CN202110102629.7A
Other languages
English (en)
Other versions
CN112911610B (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 CN202110102629.7A priority Critical patent/CN112911610B/zh
Publication of CN112911610A publication Critical patent/CN112911610A/zh
Application granted granted Critical
Publication of CN112911610B publication Critical patent/CN112911610B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/24Cell structures
    • H04W16/28Cell structures using beam steering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/02Arrangements for optimising operational condition
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Molecular Biology (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种面向无线通信网络优化的信道波束图样拓展方法,首先,获取实际无线通信网络中的位置信息、工程参数以及通过波束训练得到的信道波束图样,对数据进行预处理,筛选合适的数据特征,并建立用于信道波束图样拓展的数据集。接着,建立由神经网络和多变量高斯过程回归模型串联而成的深度多变量高斯过程回归模型。然后,利用信道波束图样拓展数据集,根据负对数似然函数对深度多变量高斯过程回归模型进行训练,联合优化神经网络权重参数以及高斯过程回归超参数。最后,根据从无线通信网络中获取的新的位置信息和工程参数,利用得到的深度多变量高斯过程回归模型拓展得到对应的信道波束图样,可以有效实现信道波束图样拓展。

Description

一种面向无线通信网络优化的信道波束图样拓展方法
技术领域
本发明涉及一种面向无线通信网络优化数据增广技术,属于无线通信网络优化技术领域。
背景技术
为应对日益增长的数据高速传输需求和海量的终端连接,大规模多输入多输出技术广泛应用于第五代移动通信系统中。大规模天线阵列可以通过波束成形技术实现空分复用,从而大幅提高无线通信系统信道容量。而随着天线数目的增加,信道估计复杂度激增,因此一般通过波束训练获得波束域等效信道辅助信道估计和传输设计。在此基础上经过时间平均得到的信道波束图样在网络规划和优化的系统中长期性能提升方面具有重要作用。信道波束图样受到位置和工程参数的共同影响,在实际网络规划和优化的过程中,一般只能获得预设路径和工参组合下的信道波束图样。因此,信道波束图样拓展技术对网规网优有巨大价值。
信道波束图样拓展技术属于一种新型的信道信息拓展技术。传统信道信息拓展广泛采用论文(Di Taranto R,Muppirisetty S,Raulefs R,et al.Location-AwareCommunications for 5G Networks:How location information can improvescalability,latency,and robustness of 5G[J].IEEE Signal Processing Magazine,2014,31(6):102-112.参考文献1.)提出的基于高斯过程回归模型的信道插值方法。该方法利用信道大尺度衰落与位置之间的相关性,可以预测不同位置的路径损耗,但无法预测信道的空间特性,同时没有考虑基站工参配置对信道信息的影响。
近年来,为了克服基于模型的信道插值方法的局限性,深度学习方法被应用于信道信息拓展的研究中。论文(Ruichen Deng,Zhiyuan Jiang,Sheng Zhou*,et al.A Two-Step Learning and Interpolation Method for Location-Based Channel Database[C]//GLOBECOM 2018-2018IEEE Global Communications Conference.IEEE,2018.参考文献2.)提出了一种两阶段的基于位置的信道数据集构建方法。该方法先利用K近邻方法建立粗糙的基于位置的信道数据集,再利用卷积神经网络进行插值细化基于位置的信道数据集,即利用有限的信道数据拓展得到不同位置上的信道信息,但利用该方法构建的数据集中只包含信道的大尺度特征,且没有设计具体基站工参配置。论文(Alrabeiah M,Alkhateeb A .Deep Learning for TDD and FDD Massive MIMO:Mapping Channels inSpace and Frequency[C]//2019 53rd Asilomar Conference on Signals,Systems,andComputers.2019.参考文献3.)提出一种基于深度神经网络的信道映射方法。该方法可以根据当前频率下天线域和子载波域上的信道分布,利用深度神经网络学习当前环境特征,预测得到目标频率下天线域和子载波域上的信道分布。该方法考虑了多天线系统下信道的空间特性,但缺少波束域信道分布与位置信息以及工参配置之间的关系。
结合基于模型的信道信息拓展技术和利用深度学习基于数据的信道信息拓展技术,可以采用数据与模型协同驱动的深度多变量高斯过程回归模型来实现信道波束图样的拓展。论文(Wilson A G,Hu Z,Salakhutdinov R,et al.Deep Kernel Learning[J].Artificial Intelligence and Statistics,2016.参考文献4.)提出的深度核学习可用于高斯过程回归,通过神经网络对输入特征进行变换,之后经过高斯过程回归进行预测,从而减小模型失配带来的预测误差,同时可以联合优化神经网络权重参数和高斯回归超参数。论文(Chen Z,Wang B,Gorban A N.Multivariate Gaussian and Student-t ProcessRegression for Multi-output Prediction[J].Neural Computing and Applications,2019(3).参考文献5.)提出的多变量高斯过程回归引入了输出协方差函数来表征多输出预测之间的相关性,从而可用于解决多输出预测问题。因此,深度多变量高斯过程回归模型可用于信道波束图样拓展,从而提高无线通信网络优化的性能。
发明内容
技术问题:本发明旨在针对实际网络规划和优化的过程中只能获得预设路径和工参组合下的信道波束图样导致样本数量不足且分布不合理的问题,提出了一种基于深度多变量高斯过程回归模型的信道波束图样拓展方法,能够根据位置信息和工程参数拟合得到对应的信道波束图样,从而减小样本数量不足且分布不合理对无线通信网络优化算法的影响
技术方案:鉴于信道波束图样对网络规划和优化的重要性及信道波束图样拓展的实用价值,本发明提出一种面向无线通信网络优化的信道波束图样拓展方法,将深度多变量高斯过程回归模型用于信道波束图样拓展中,根据实际环境路测数据构建信道波束图样拓展数据集,并对所建立的深度多变量高斯过程回归模型进行训练,基于位置信息和工程参数利用该模型拟合对应的信道波束图样,拟合误差明显小于传统高斯过程回归和神经网络预测等方法。
本发明所提出的一种面向无线通信网络优化的信道波束图样拓展方法包括如下步骤:
步骤1:获取实际无线通信网络中的位置信息数据、工程参数数据以及通过波束训练得到的信道波束图样数据,对以上无线通信网络中获得的数据进行预处理并筛选合适的数据特征,建立用于信道波束图样拓展的数据集;
步骤2:建立由神经网络和多变量高斯过程回归模型串联而成的深度多变量高斯过程回归模型;
步骤3:利用信道波束图样拓展的数据集根据负对数似然函数对深度多变量高斯过程回归模型进行训练,联合优化神经网络权重参数和高斯过程回归超参数;
步骤4:根据从无线通信网络中获取的新的位置信息数据和工程参数数据利用已建立的并完成训练的深度多变量高斯过程回归模型拓展得到对应的信道波束图样。
其中,
步骤1中信道波束图样拓展数据集建立方法为:从无线通信网络中获取实测数据并进行数据预处理,选取信道波束图样相关的数据特征来建立信道波束图样拓展数据集;其中数据集的输入特征主要包括位置信息x=[xUE,yUE,hUE,xBS,yBS,hBS]和工程参数
Figure BDA0002914150880000021
xUE,yUE为用户位置,hUE为用户高度,xBS,yBS为基站位置,hBS为基站高度,
Figure BDA0002914150880000022
为基站天线阵列方位角、θ为基站天线阵列下倾角θ,输出标签为通过波束训练得到的信道波束图样
Figure BDA0002914150880000023
即基站与用户间信道在不同波束上投影值的时间平均值,受到位置信息x和工程参数Θ共同影响。
步骤2中深度多变量高斯过程回归模型建立方法如下:
深度多变量高斯过程回归模型由神经网络u=g(x,Θ;w)和多变量高斯过程回归模型
Figure BDA0002914150880000024
串联组成,其中x为位置信息,Θ为工程参数,而
Figure BDA0002914150880000025
为通过波束训练得到的信道波束图样,受到位置信息x和工程参数Θ共同影响;u为经过神经网络非线性变换后输出的中间特征,同时作为高斯回归过程的输入,w为神经网络权重参数,
Figure BDA0002914150880000026
为由均值函数μ,协方差核函数kγ以及输出相关矩阵Ω确定的多变量高斯过程,γ为高斯过程回归超参数,均值函数μ用以表征每一维信道波束图样均值与输入特征的关系,协方差核函数kγ用以表征不同输入特征对应的信道波束图样的相关性,输出相关函数Ω用以表征不同维信道波束图样的相关性。
步骤3中深度多变量高斯过程回归模型进行训练方法如下:
基于步骤1建立的信道波束图样拓展的数据集及步骤2建立的深度多变量高斯过程回归模型,以最小化负对数似然函数为目标优化神经网络权重参数w和高斯过程回归超参数γ,负对数似然函数计算方法如下:
Figure BDA0002914150880000027
其中N为数据集样本总数,D为信道波束图样的维数;
Figure BDA0002914150880000031
为数据集内样本对应信道波束图样标签组成的矩阵;M为均值矩阵,Md,n=μd(un),1≤d≤D,1≤n≤N,Σ为列协方差矩阵,Σi,j=kγ(ui,uj),1≤i≤N,1≤j≤N。kγ为协方差核函数,γ为高斯过程回归超参数;un=g(xnn;w),1≤n≤N,是位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征,w为神经网络权重参数,Ω为行协方差矩阵,
Figure BDA0002914150880000032
根据链式法则,关于神经网络权重参数w和高斯过程回归超参数γ的偏导数计算如下:
Figure BDA0002914150880000033
Figure BDA0002914150880000034
其中
Figure BDA0002914150880000035
由对数化似然函数计算公式获得,
Figure BDA0002914150880000036
Figure BDA0002914150880000037
由基本核函数具体定义获得,
Figure BDA0002914150880000038
通过神经网络的反向传播获得;根据以上偏导数,根据梯度下降的方法实现神经网络权重参数w和高斯过程回归超参数γ的联合优化,即完成深度高斯过程回归模型的联合训练。
步骤4中基于深度多变量高斯过程回归模型信道波束图样拓展方法如下:
根据步骤2的深度多变量高斯过程回归模型,信道波束图样满足多变量高斯过程:
Figure BDA0002914150880000039
其中μ为均值函数,kγ为协方差核函数,Ω为输出相关函数,w为神经网络权重参数,γ为高斯过程回归超参数;
根据多变量高斯过程的性质,步骤1建立的信道波束图样数据集中的所有样本对应信道波束图样满足多变量高斯分布:
Figure BDA00029141508800000310
其中
Figure BDA00029141508800000311
为数据集内样本对应信道波束图样标签组成的矩阵;M为均值矩阵,Md,n=μd(un),1≤d≤D,1≤n≤N,Σ为列协方差矩阵,Σi,j=kγ(ui,uj),1≤i≤N,1≤j≤N。kγ为协方差核函数。un=g(xnn;w),1≤n≤N,是位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征;Ω为行协方差矩阵,
Figure BDA00029141508800000312
N为数据集样本总数,D为信道波束图样的维数;
此时,重新获取无线通信网络中的路测数据并对数据进行预处理,获得新的位置信息和工程参数作为测试样本,其对应的信道波束图样未知,根据高斯过程模型的性质,测试样本对应的信道波束图样与数据集中样本对应的信道波束图样满足以下联合多变量高斯分布:
Figure BDA0002914150880000041
其中
Figure BDA0002914150880000042
为所有测试样本对应的信道波束图样组成的矩阵,N为数据集样本总数,M为测试样本总数;M为均值矩阵,Md,m=μd(uN+m),1≤d≤D,1≤m≤M,D为信道波束图样的维数;Σ_*,Σ*_,Σ**均为协方差矩阵,(Σ_*)i,j=kγ(ui,uN+j),1≤i≤N,1≤j≤M,(Σ*_)i,j=kγ(uN+i,uj),1≤i≤M,1≤j≤N,Σi,j=kγ(uN+i,uN+j),1≤i≤M,1≤j≤M。kγ为协方差核函数;un=g(xnn;w),1≤n≤N+M,位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征,Ω为行协方差矩阵,
Figure BDA0002914150880000043
根据联合高斯分布的性质,已知数据集中样本的位置信息、工程参数及其对应信道波束图样以及测试样本的位置信息和工程参数条件下,测试样本对应的信道波束图样的分布满足以下多变量高斯分布:
Figure BDA0002914150880000044
其中
Figure BDA0002914150880000045
N为数据集样本总数,M为测试样本总数,un=g(xnn;w),1≤n≤N+M,位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征,
Figure BDA0002914150880000046
Figure BDA0002914150880000047
分别表示均值矩阵、列协方差矩阵和行协方差矩阵计算方式如下:
Figure BDA0002914150880000048
Figure BDA0002914150880000049
Figure BDA00029141508800000410
根据高斯分布性质,拓展得到测试样本对应信道波束图样
Figure BDA00029141508800000411
的均值和协方差如下:
Figure BDA00029141508800000412
拓展得到的测试样本对应信道波束图样
Figure BDA00029141508800000413
的均值作为测试样本对应的信道波束图样的拓展结果,协方差则提供测试样本对应的信道波束图样的置信区间。
有益效果:与现有技术相比,本发明除位置信息外,考虑了工程参数对信道波束图样的影响,实现了信道波束图样的拓展,同时结合了基于模型的插值方法和深度学习方法,使得拟合误差明显小于传统高斯过程回归和神经网络预测等方法。本发明可用于无线通信网络优化中的数据增广,减小数据样本不足或分布不合理对无线通信网络优化算法的影响。
附图说明
下面结合附图和实施例对本发明做进一步说明。
图1是本发明所述面向无线通信网络优化的信道波束图样拓展方法算法示意图。
图2是示例性实施例一中使用本发明所述信道波束图样拓展方法实际通信环境示意图。
图3是示例性实施例一中使用本发明所述信道波束图样拓展方法平均绝对误差分布图。
图4是示例性实施例一中使用本发明所述信道波束图样拓展方法平均相对误差分布图。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求项要求所限定的范围。
以下为本发明的一种适合于面向无线通信网络优化的信道波束图样拓展方法的实施方式:
在本实施例中,如图2所示,图2为实际通信环境示意图。其中基站坐标为(242.42m,489.17m,6m),用户坐标为(x,y,z),x在242.42m≤x≤272.42m的范围内以2m为间隔选取,y在479.17m≤y≤499.17m范围内以2m为间隔选取,z=2m。基站采用8×4的均匀平面阵进行传输,基站天线阵列方位角
Figure BDA0002914150880000051
为天线阵列法线与xOz平面的夹角,以顺时针方向为正方向,在-5°,0°,5°范围内选取,基站天线阵列下倾角θ为天线阵列法线与xOy平面的夹角,以顺时针方向为正方向,在0°,5°,10°范围内选取。基站发射功率固定为1W。
本发明所提出的一种面向无线通信网络优化的信道波束图样拓展方法包括如下步骤:
步骤1:获取实际无线通信网络中的位置信息数据、工程参数数据以及通过波束训练得到的信道波束图样数据,对以上无线通信网络中获得的数据进行预处理并筛选合适的数据特征,建立用于信道波束图样拓展的数据集;
步骤2:建立由神经网络和多变量高斯过程多变量高斯过程回归模型串联而成的深度多变量高斯过程回归模型;
步骤3:利用信道波束图样拓展的数据集根据负对数似然函数对深度多变量高斯过程回归模型进行训练,联合优化神经网络权重参数和高斯过程回归超参数;
步骤4:根据从无线通信网络中获取的新的位置信息数据和工程参数数据利用已建立的并完成训练的深度多变量高斯过程回归模型拓展得到对应的信道波束图样。
在所述步骤1中,从无线通信网络中获取实测数据并进行数据预处理,选取信道波束图样相关的数据特征来建立信道波束图样拓展数据集。其中数据集的输入特征主要包括位置信息x=[xUE,yUE,hUE,xBS,yBS,hBS]和工程参数
Figure BDA0002914150880000052
xUE,yUE为用户位置,hUE为用户高度,xUE,yUE为用户位置,hUE为用户高度,
Figure BDA0002914150880000053
为基站天线阵列方位角、θ为基站天线阵列下倾角θ,输出标签为通过波束训练得到的信道波束图样
Figure BDA0002914150880000054
即基站与用户间信道在不同波束上投影值的时间平均值,受到位置信息x和工程参数Θ共同影响。特别地,在本实施例中,基站采用二维离散傅里叶矩阵(2D-DFT)作为波束码本进行波束训练并将所有波束上的参考信号接收强度(RSRP)作为信道波束图样
Figure BDA0002914150880000055
单位为dBm,其计算方式如下:
Figure BDA0002914150880000056
其中h为当前位置信息和工程参数下基站和用户间的信道矢量,D为2D-DFT矩阵,p为基站的发射功率。
Figure BDA0002914150880000057
Dk为D的第k列,即波束码本中的第k个波束,1≤k≤W×H。
从图2所示的实际无线通信网络中获取的位置信息和工程参数以及对应的信道波束图样,随机选取其中一定比例样本构成信道波束图样拓展数据集用于深度多变量高斯过程回归模型的训练,剩余样本作为测试样本,用于评估信道波束图样拓展的性能。
在所述步骤2中,基于步骤1建立的信道波束图样拓展数据集,建立由神经网络u=g(x,Θ;w)和多变量高斯过程回归模型
Figure BDA0002914150880000061
串联组成的深度多变量高斯过程回归模型,其中x为位置信息,Θ为工程参数,而
Figure BDA0002914150880000062
为通过波束训练得到的信道波束图样,受到位置信息x和工程参数Θ共同影响。u为经过神经网络非线性变换后输出的中间特征,同时作为高斯回归过程的输入。w为神经网络权重参数。
Figure BDA0002914150880000063
为由均值函数μ,协方差核函数kγ以及输出相关矩阵Ω确定的多变量高斯过程,γ为高斯过程回归超参数。均值函数μ用以表征每一维信道波束图样均值与输入特征的关系,协方差核函数kγ用以表征不同输入特征对应的信道波束图样的相关性,输出相关函数Ω用以表征不同维信道波束图样的相关性。本实施例中,神经网络采用全连接网络,多变量高斯回归过程均值函数采用常数函数,
Figure BDA0002914150880000064
N为数据集样本总数,输出相关函数则采用狄拉克函数,Ω=δi,j,当i=j时,δi,j=1,当i=j时,δi,j=0,协方差核函数则采用带随机噪声项的平方指数函数
Figure BDA0002914150880000065
σ2为随机噪声,用来表征信道波束图样自身的随机性。当i=j时,δi,j=1,当i=j时,δi,j=0。sf和l分别为平方指数核函数的幅度超参数和尺度超参数,γ=[sf,l,σ]。
为体现深度多变量高斯过程回归模型的优势,本实施例中,采用传统多变量高斯过程回归模型和神经网络模型作为对比算法,此处采用与深度多变量高斯过程回归模型相似的配置。
在所述步骤3中,基于步骤1建立的信道波束图样数据集及步骤2建立的深度多变量高斯过程回归模型,以最小化负对数似然函数为目标优化神经网络权重参数w和高斯过程回归超参数γ,负对数似然函数计算方法如下:
Figure BDA0002914150880000066
其中N为数据集样本总数,D为信道波束图样的维数。
Figure BDA0002914150880000067
为数据集内样本对应信道波束图样标签组成的矩阵。M为均值矩阵,Md,n=μd(un),1≤d≤D,1≤n≤N,Σ为列协方差矩阵,Σi,j=kγ(ui,uj),1≤i≤N,1≤j≤N。kγ为协方差核函数,γ为高斯过程回归超参数。un=g(xnn;w),1≤n≤N,是位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征,w为神经网络权重参数。Ω为行协方差矩阵,
Figure BDA0002914150880000068
根据链式法则,关于神经网络权重参数w和高斯过程回归超参数γ的偏导数计算如下:
Figure BDA0002914150880000071
Figure BDA0002914150880000072
其中
Figure BDA0002914150880000073
由对数化似然函数计算公式获得,
Figure BDA0002914150880000074
由基本核函数具体定义获得,
Figure BDA0002914150880000075
通过神经网络的反向传播获得。根据以上偏导数,根据梯度下降的方法实现神经网络权重参数w和高斯过程回归超参数γ的联合优化,即完成深度高斯过程回归模型的联合训练。
为体现深度多变量高斯过程回归模型的优势,本实施例中,采用多变量高斯过程回归模型和全连接网络模型作为对比算法,此处利用相似的方法对模型进行训练,全连接网络模型仅优化神经网络权重参数,多变量高斯过程回归模型仅优化高斯过程回归超参数。
在所述步骤4中,根据步骤2的深度多变量高斯过程回归模型,信道波束图样满足多变量高斯过程:
Figure BDA0002914150880000076
其中μ为均值函数,kγ为协方差核函数,Ω为输出相关函数,w为神经网络权重参数,γ为高斯过程回归超参数。
根据多变量高斯过程的性质,步骤1建立的信道波束图样数据集中的所有样本对应信道波束图样满足多变量高斯分布:
Figure BDA0002914150880000077
其中
Figure BDA0002914150880000078
为数据集内样本对应信道波束图样标签组成的矩阵。M为均值矩阵,Md,n=μd(un),1≤d≤D,1≤n≤N,Σ为列协方差矩阵,Σi,j=kγ(ui,uj),1≤i≤N,1≤j≤N。kγ为协方差核函数。un=g(xnn;w),1≤n≤N,是位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征。Ω为行协方差矩阵,
Figure BDA0002914150880000079
N为数据集样本总数,D为信道波束图样的维数。
此时,将测试样本的位置信息和工程参数作为模型输入,其对应的信道波束图样未知,根据高斯过程模型的性质,测试样本对应的信道波束图样与数据集中样本对应的信道波束图样满足以下联合多变量高斯分布:
Figure BDA00029141508800000710
其中
Figure BDA00029141508800000711
为所有测试样本对应的信道波束图样组成的矩阵,N为数据集样本总数,M为测试样本总数。M为均值矩阵,Md,m=μd(uN+m),1≤d≤D,1≤m≤M,D为信道波束图样的维数。Σ_*,Σ*_,Σ**均为协方差矩阵,(Σ_*)i,j=kγ(ui,uN+j),1≤i≤N,1≤j≤M,(Σ*_)i,j=kγ(uN+i,uj),1≤i≤M,1≤j≤N,Σi,j=kγ(uN+i,uN+j),1≤i≤M,1≤j≤M。kγ为协方差核函数。un=g(xnn;w),1≤n≤N+M,位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征。Ω为行协方差矩阵,
Figure BDA0002914150880000081
根据联合高斯分布的性质,已知数据集中样本的位置信息、工程参数及其对应信道波束图样以及测试样本的位置信息和工程参数条件下,测试样本对应的信道波束图样的分布满足以下多变量高斯分布:
Figure BDA0002914150880000082
其中
Figure BDA0002914150880000083
N为数据集样本总数,M为测试样本总数。un=g(xnn;w),1≤n≤N+M,位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征。
Figure BDA0002914150880000084
Figure BDA0002914150880000085
分别表示均值矩阵、列协方差矩阵和行协方差矩阵,计算方式如下:
Figure BDA0002914150880000086
Figure BDA0002914150880000087
Figure BDA0002914150880000088
根据高斯分布性质,拓展得到测试样本对应信道波束图样
Figure BDA0002914150880000089
的均值和协方差如下:
Figure BDA00029141508800000810
拓展得到的测试样本对应信道波束图样
Figure BDA00029141508800000811
的均值作为测试样本对应的信道波束图样的拓展结果,协方差则提供测试样本对应的信道波束图样的置信区间。
为体现深度多变量高斯过程回归模型的优势,本实施例中,采用多变量高斯过程回归模型和全连接网络模型作为对比算法,此处全连接网络模型将神经网络的输出作为拓展结果,而多变量高斯过程回归模型将模型输出均值作为拓展结果,模型输出的协方差则可用于置信度的评估。
信道波束图样拓展的性能可通过信道波束图样拓展结果相对于实测值在各波束上的平均绝对误差和平均相对误差进行评估,其计算方式如下:
Figure BDA00029141508800000812
Figure BDA00029141508800000813
其中
Figure BDA00029141508800000814
Figure BDA00029141508800000815
分别为第m个样本第n个波束上的信道波束图样的拓展值和实测值,1≤m≤M,1≤d≤D,M为测试样本总数,D为信道波束图样的维数。
以下是使用本发明所述信道波束图样拓展方法在各波束上平均绝对误差和平均相对误差的分布以及与对比算法的比较:
由图3和图4可以看出,使用本发明所述信道波束图样拓展方法得到的拓展结果相比于实测值能保持较低的平均绝对误差和平均相对误差,即基于有限的数据集可以实现信道波束图样拓展。拓展结果可用于无线通信网络优化中的数据增广,从而减小数据样本不足或分布不合理等对无线通信网络优化算法的影响而且信道波束图样拓展拟合误差明显小于传统高斯过程回归和神经网络预测。

Claims (5)

1.一种面向无线通信网络优化的信道波束图样拓展方法,其特征在于,该方法包括如下步骤:
步骤1:获取实际无线通信网络中的位置信息数据、工程参数数据以及通过波束训练得到的信道波束图样数据,对以上无线通信网络中获得的数据进行预处理并筛选合适的数据特征,建立用于信道波束图样拓展的数据集;
步骤2:建立由神经网络和多变量高斯过程回归模型串联而成的深度多变量高斯过程回归模型;
步骤3:利用信道波束图样拓展的数据集根据负对数似然函数对深度多变量高斯过程回归模型进行训练,联合优化神经网络权重参数和高斯过程回归超参数;
步骤4:根据从无线通信网络中获取的新的位置信息数据和工程参数数据利用已建立的并完成训练的深度多变量高斯过程回归模型拓展得到对应的信道波束图样。
2.根据权利要求1所述的一种面向无线通信网络优化的信道波束图样拓展方法,其特征在于,步骤1中信道波束图样拓展数据集建立方法为:从无线通信网络中获取实测数据并进行数据预处理,选取信道波束图样相关的数据特征来建立信道波束图样拓展数据集;其中数据集的输入特征主要包括位置信息x=[xUE,yUE,hUE,xBS,yBS,hBS]和工程参数
Figure FDA0002914150870000011
xUE,yUE为用户位置,hUE为用户高度,xBS,yBS为基站位置,hBS为基站高度,
Figure FDA0002914150870000012
为基站天线阵列方位角、θ为基站天线阵列下倾角θ,输出标签为通过波束训练得到的信道波束图样
Figure FDA0002914150870000013
即基站与用户间信道在不同波束上投影值的时间平均值,受到位置信息x和工程参数Θ共同影响。
3.根据权利要求1所述的一种面向无线通信网络优化的信道波束图样拓展方法,其特征在于,步骤2中深度多变量高斯过程回归模型建立方法如下:
深度多变量高斯过程回归模型由神经网络u=g(x,Θ;w)和多变量高斯过程回归模型
Figure FDA0002914150870000014
串联组成,其中x为位置信息,Θ为工程参数,而
Figure FDA0002914150870000015
为通过波束训练得到的信道波束图样,受到位置信息x和工程参数Θ共同影响;u为经过神经网络非线性变换后输出的中间特征,同时作为高斯回归过程的输入,w为神经网络权重参数,
Figure FDA0002914150870000016
为由均值函数μ,协方差核函数kγ以及输出相关矩阵Ω确定的多变量高斯过程,γ为高斯过程回归超参数,均值函数μ用以表征每一维信道波束图样均值与输入特征的关系,协方差核函数kγ用以表征不同输入特征对应的信道波束图样的相关性,输出相关函数Ω用以表征不同维信道波束图样的相关性。
4.根据权利要求1所述的一种面向无线通信网络优化的信道波束图样拓展方法,其特征在于,步骤3中深度多变量高斯过程回归模型进行训练方法如下:
基于步骤1建立的信道波束图样拓展的数据集及步骤2建立的深度多变量高斯过程回归模型,以最小化负对数似然函数为目标优化神经网络权重参数w和高斯过程回归超参数γ,负对数似然函数计算方法如下:
Figure FDA0002914150870000017
其中N为数据集样本总数,D为信道波束图样的维数;
Figure FDA0002914150870000021
为数据集内样本对应信道波束图样标签组成的矩阵;M为均值矩阵,Md,n=μd(un),1≤d≤D,1≤n≤N,Σ为列协方差矩阵,Σi,j=kγ(ui,uj),1≤i≤N,1≤j≤N。kγ为协方差核函数,γ为高斯过程回归超参数;un=g(xnn;w),1≤n≤N,是位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征,w为神经网络权重参数,Ω为行协方差矩阵,
Figure FDA0002914150870000022
根据链式法则,关于神经网络权重参数w和高斯过程回归超参数γ的偏导数计算如下:
Figure FDA0002914150870000023
Figure FDA0002914150870000024
其中
Figure FDA0002914150870000025
由对数化似然函数计算公式获得,
Figure FDA0002914150870000026
Figure FDA0002914150870000027
由基本核函数具体定义获得,
Figure FDA0002914150870000028
通过神经网络的反向传播获得;根据以上偏导数,根据梯度下降的方法实现神经网络权重参数w和高斯过程回归超参数γ的联合优化,即完成深度高斯过程回归模型的联合训练。
5.根据权利要求1所述的一种面向无线通信网络优化的信道波束图样拓展方法,其特征在于,步骤4中基于深度多变量高斯过程回归模型信道波束图样拓展方法如下:
根据步骤2的深度多变量高斯过程回归模型,信道波束图样满足多变量高斯过程:
Figure FDA0002914150870000029
其中μ为均值函数,kγ为协方差核函数,Ω为输出相关函数,w为神经网络权重参数,γ为高斯过程回归超参数;
根据多变量高斯过程的性质,步骤1建立的信道波束图样数据集中的所有样本对应信道波束图样满足多变量高斯分布:
Figure FDA00029141508700000210
其中
Figure FDA00029141508700000211
为数据集内样本对应信道波束图样标签组成的矩阵;M为均值矩阵,Md,n=μd(un),1≤d≤D,1≤n≤N,Σ为列协方差矩阵,Σi,j=kγ(ui,uj),1≤i≤N,1≤j≤N。kγ为协方差核函数。un=g(xnn;w),1≤n≤N,是位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征;Ω为行协方差矩阵,
Figure FDA00029141508700000212
N为数据集样本总数,D为信道波束图样的维数;
此时,重新获取无线通信网络中的路测数据并对数据进行预处理,获得新的位置信息和工程参数作为测试样本,其对应的信道波束图样未知,根据高斯过程模型的性质,测试样本对应的信道波束图样与数据集中样本对应的信道波束图样满足以下联合多变量高斯分布:
Figure FDA0002914150870000031
其中
Figure FDA0002914150870000032
为所有测试样本对应的信道波束图样组成的矩阵,N为数据集样本总数,M为测试样本总数;M为均值矩阵,Md,m=μd(uN+m),1≤d≤D,1≤m≤M,D为信道波束图样的维数;
Figure FDA00029141508700000312
Σ**均为协方差矩阵,
Figure FDA00029141508700000313
1≤i≤N,1≤j≤M,
Figure FDA00029141508700000314
1≤i≤M,1≤j≤N,Σi,j=kγ(uN+i,uN+j),1≤i≤M,1≤j≤M。kγ为协方差核函数;un=g(xnn;w),1≤n≤N+M,位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征,Ω为行协方差矩阵,
Figure FDA0002914150870000033
根据联合高斯分布的性质,已知数据集中样本的位置信息、工程参数及其对应信道波束图样以及测试样本的位置信息和工程参数条件下,测试样本对应的信道波束图样的分布满足以下多变量高斯分布:
Figure FDA0002914150870000034
其中
Figure FDA0002914150870000035
N为数据集样本总数,M为测试样本总数,un=g(xnn;w),1≤n≤N+M,位置信息xn和工程参数Θn经过神经网络非线性变换的中间特征,
Figure FDA0002914150870000036
Figure FDA0002914150870000037
分别表示均值矩阵、列协方差矩阵和行协方差矩阵计算方式如下:
Figure FDA0002914150870000038
根据高斯分布性质,拓展得到测试样本对应信道波束图样
Figure FDA0002914150870000039
的均值和协方差如下:
Figure FDA00029141508700000310
拓展得到的测试样本对应信道波束图样
Figure FDA00029141508700000311
的均值作为测试样本对应的信道波束图样的拓展结果,协方差则提供测试样本对应的信道波束图样的置信区间。
CN202110102629.7A 2021-01-25 2021-01-25 一种面向无线通信网络优化的信道波束图样拓展方法 Active CN112911610B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110102629.7A CN112911610B (zh) 2021-01-25 2021-01-25 一种面向无线通信网络优化的信道波束图样拓展方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110102629.7A CN112911610B (zh) 2021-01-25 2021-01-25 一种面向无线通信网络优化的信道波束图样拓展方法

Publications (2)

Publication Number Publication Date
CN112911610A true CN112911610A (zh) 2021-06-04
CN112911610B CN112911610B (zh) 2022-07-12

Family

ID=76119597

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110102629.7A Active CN112911610B (zh) 2021-01-25 2021-01-25 一种面向无线通信网络优化的信道波束图样拓展方法

Country Status (1)

Country Link
CN (1) CN112911610B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115882985A (zh) * 2022-11-28 2023-03-31 西安交通大学 一种基于高斯过程回归的低轨卫星信道预测方法及系统
WO2023097641A1 (zh) * 2021-12-03 2023-06-08 Oppo广东移动通信有限公司 信道数据的数据扩充方法、装置、设备及存储介质
WO2024051789A1 (zh) * 2022-09-08 2024-03-14 华为技术有限公司 一种波束管理方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109617584A (zh) * 2019-01-08 2019-04-12 南京邮电大学 一种基于深度学习的mimo系统波束成形矩阵设计方法
CN110598842A (zh) * 2019-07-17 2019-12-20 深圳大学 一种深度神经网络超参数优化方法、电子设备及存储介质
CN110995327A (zh) * 2019-12-17 2020-04-10 电子科技大学 一种多载波mimo系统的混合波束成形优化方法及系统
CN111817767A (zh) * 2020-06-30 2020-10-23 山西省信息产业技术研究院有限公司 一种基于动态优化策略的mvdr波束成形方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109617584A (zh) * 2019-01-08 2019-04-12 南京邮电大学 一种基于深度学习的mimo系统波束成形矩阵设计方法
CN110598842A (zh) * 2019-07-17 2019-12-20 深圳大学 一种深度神经网络超参数优化方法、电子设备及存储介质
CN110995327A (zh) * 2019-12-17 2020-04-10 电子科技大学 一种多载波mimo系统的混合波束成形优化方法及系统
CN111817767A (zh) * 2020-06-30 2020-10-23 山西省信息产业技术研究院有限公司 一种基于动态优化策略的mvdr波束成形方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
CHENG ZHANG等: "Channel Beam Pattern Extension for Massive MIMO via Deep Gaussian Process Regression", 《2021 IEEE/CIC INTERNATIONAL CONFERENCE ON COMMUNICATIONS IN CHINA (ICCC)》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023097641A1 (zh) * 2021-12-03 2023-06-08 Oppo广东移动通信有限公司 信道数据的数据扩充方法、装置、设备及存储介质
WO2024051789A1 (zh) * 2022-09-08 2024-03-14 华为技术有限公司 一种波束管理方法
CN115882985A (zh) * 2022-11-28 2023-03-31 西安交通大学 一种基于高斯过程回归的低轨卫星信道预测方法及系统

Also Published As

Publication number Publication date
CN112911610B (zh) 2022-07-12

Similar Documents

Publication Publication Date Title
CN112911610B (zh) 一种面向无线通信网络优化的信道波束图样拓展方法
Morais et al. Position-aided beam prediction in the real world: How useful GPS locations actually are?
Butt et al. RF fingerprinting and deep learning assisted UE positioning in 5G
CN113014295B (zh) 一种去蜂窝大规模mimo系统上行联合接收方法
Heng et al. Learning site-specific probing beams for fast mmWave beam alignment
CN113411110A (zh) 一种基于深度强化学习的毫米波通信波束训练方法
CN110708129B (zh) 一种无线信道状态信息获取方法
Jiang et al. Inferring remote channel state information: Cramér-Rae lower bound and deep learning implementation
CN113595666A (zh) 大规模mimo系统基于时频联合相关性的信道预测方法和系统
CN114826462A (zh) 面向大规模mimo统计端口选择的波束域信道增广方法
CN113438682A (zh) 一种基于波束成形的sage-bem 5g无线信道参数提取方法
Moosavi et al. Fingerprinting localization method based on clustering and Gaussian process regression in distributed massive MIMO systems
Kazemi et al. Beam SNR prediction using channel charting
Deng et al. Network-side localization via semi-supervised multi-point channel charting
US20210376886A1 (en) Method of Parameter Estimation for a Multi-input Multi-output System
CN107346985B (zh) 一种结合发射天线选择技术的干扰对齐方法
Prasad et al. A numerical approximation method for RSS-based user positioning in distributed massive MIMO
CN112752339B (zh) 一种基于大规模mimo单站系统的指纹库更新方法
CN113030853B (zh) 基于rss和aoa联合测量的多辐射源无源定位方法
CN112217750B (zh) 基于卷积神经网络的信道统计信息估计方法
Sankar et al. Learning to Precode for Integrated Sensing and Communication Systems
CN115310480A (zh) 一种基于深度神经网络的doa估计算法性能分析方法
Yang et al. Deep-learning-based bouncing-order prediction for propagation channel characterization using graph modeling
Long et al. Deep Learning-Based Realization of Millimeter-Wave Beam Alignment with Non-Uniform Distribution of Users
CN111277313A (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