可重构智能表面辅助的多用户通信宽带信道估计.pdf

可重构智能表面辅助的多用户通信宽带信道估计李贵勇杜一舟*王丹(重庆邮电大学通信与信息工程学院重庆400065)摘要：针对太赫兹(THz)链路的严重传输衰减和宽带系统中波束斜视导致传统信道估计方案性能下降的问题，该文构建了可重构智能表面(RIS)辅助多用户THz通信模型，并提出一种低复杂度的两阶段级联信道估计方案。在第1阶段，利用THz的稀疏性和对数和函数，将信道估计问题转化为目标优化问题，通过梯度下降法优化目标函数，使待估信道参数迭代逼近最优解，从而估计出典型用户级联信道；在第2阶段，利用其他用户的级联信道与典型用户信道的强相关性，以较低的导频开销来估计其他用户的级联信道。仿真结果表明，所提方案相较于其他方案具有更好的性能。关键词：可重构智能反射表面；太赫兹通信；信道估计；波束斜视中图分类号：TN929.5文献标识码：A文章编号：1009-5896(2023)07-2443-08DOI:10.11999/JEIT220775Wideband Channel Estimation for Multiuser Communication Based onReconfigurable Intelligent Surface AssistedLIGuiyongDUYizhouWANGDan(School of Communication and Information Engineering,Chongqing University of Posts andTelecommunications,Chongqing 400065,China)Abstract:TosolvetheproblemofseveretransmissionattenuationofTeraHertz(THz)links,andtheperformancedegradationoftraditionalchannelestimationschemescausedbybeamsquintinwidebandsystems,amulti-userTHzcommunicationmodelassistedbyReconfigurableIntelligentSurfaces(RIS)isconstructedinthispaper,andalowcomplexitytwo-stagecascadedchannelestimationschemeisproposed.Inthefirststage,thechannelestimationproblemistransformedintoanobjectiveoptimizationproblembyusingthesparsityoftheTHzandlog-sumfunction,andtheobjectivefunctionisoptimizedbygradientdescentmethod,sothatthechannelparameterstobeestimatedareiterativelyclosetotheoptimalsolution,thusthetypicalusercascadechannelisestimated.Inthesecondstage,thecascadechannelsofotherusersareestimatedwithlowerpilotoverheadbyusingthestrongcorrelationbetweenthecascadechannelsofotherusersandthetypicaluserchannel.Thesimulationresultsshowthattheproposedschemehasbetterperformancethanotherschemes.Key words:ReconfigurableIntelligentSurface(RIS);TeraHertz(THz)communication;Channelestimation;Beamsquint1 引言太赫兹(TeraHertz,THz)频谱的大带宽和高数据速率，有望满足无线通信系统中爆炸式增长的数据流量需求；然而，THz的极高传播衰减和分子吸收通常会限制信号的传输距离和覆盖范围1。幸运的是，可重构智能表面(ReconfigurableIntelligentSurface,RIS)能以较低的硬件成本和能耗提高无线通信系统的容量和覆盖率2,3。然而，RIS提供的性能增益很大程度上依赖于信道状态信息(ChannelStateInformation,CSI)的准确性，且RIS的无源性和高维级联信道使得CSI的准确获取具有挑战性3。针对RIS辅助单用户通信的信道估计方案已取得很多有价值的研究成果46。然而，文献46这类针对RIS辅助单用户通信系统的信道估计方案，直接推广到RIS辅助的多用户通信系统中是低效的7。在已有的RIS辅助多用户通信的宽带信道估计方案收稿日期：2022-06-13；改回日期：2022-10-20；网络出版：2022-10-25*通信作者：杜一舟S基金项目：重庆市自然科学基金(cstc2021jcyj-msxmX0454)FoundationItem:TheNaturalScienceFoundationofChongqing(cstc2021jcyj-msxmX0454)第45卷第7期电子与信息学报Vol.45No.72023年7月JournalofElectronics&InformationTechnologyJul.2023中，文献8利用基站(BaseStation,BS)到RIS的信道高维准静态特性和RIS到用户设备(UserEquip-ment,UE)的信道移动低维特性，提出一种低导频开销的双时间尺度信道估计方案，利用坐标下降算法估计BS-RIS信道，利用最小二乘(LeastSquare,LS)算法估计RIS-UE和BS-UE信道。文献9研究载波频偏(CarrierFrequencyOffset,CFO)对RIS辅助的正交频分复用(OrthogonalFrequencyDivi-sionMultiplexing,OFDM)无线通信系统中信道估计的影响，并提出一种联合CFO和信道冲激响应的信道估计方法。文献10提出一种基于深度去噪神经网络的压缩信道估计方法，以减少训练开销。文献11将训练信号建模为张量模型，并利用代数算法解决张量补全问题，从而恢复出信道的多径参数。上述宽带信道估计方案存在以下局限性：首先，文献811的宽带信道估计方案考虑了宽带信道的频率选择性，但忽略了波束斜视的影响12。在OFDM宽带系统中，波束斜视使阵列响应变得频率相关，导致传统的信道估计方案失效13。其次，文献811的宽带信道估计方案未利用信道间的强相关性，虽然文献14利用了信道间的强相关性，但所提方案是针对窄带系统的。最后，针对RIS辅助的宽带通信系统且考虑波束斜视的研究仅有文献15，但文献15是针对单用户系统的。因此，针对波束斜视影响下RIS辅助多用户通信的宽带信道估计方案的研究具有非常重要的理论意义与现实价值。本文研究了波束斜视影响下RIS辅助多用户THz通信的上行信道估计问题，并提出一种低复杂度的稳健两阶段信道估计方案。在第1阶段，利用THz信道稀疏性，将信道估计问题转化为一个非凸目标优化问题，提出一种基于梯度下降法的迭代加权的对数和约束信道估计算法，本算法以对数和函数为约束条件，用优化最小化(Majorization-Mini-mization,MM)方法得到一个代理目标函数，然后通过梯度下降法优化目标函数，从而准确估计出信道增益、角度信息和时延参数。在第2阶段，利用所有用户在各自的级联信道(UE-RIS-BS)中共享公共信道RIS-BS信道这一特性，将其他用户的级联信道表示为典型用户CSI的缩放版本，利用LS算法进行信道估计。仿真结果表明，本文算法具有较好的性能。diag()blkdiag()Fvec()aarank()()符号说明：为克罗内克积，为哈达玛积，为对角矩阵，为块对角矩阵，为弗罗贝尼乌斯范数，为列序向量化，为大于或等于 的最小整数，为矩阵的秩，为伪逆。2 系统模型UM=Mr MoNRFNd=/2QTallTupP本文研究RIS辅助多用户大规模多输入单输出THz无线通信系统中的信道估计问题。在本系统中，BS在RIS的辅助下同时为 个单天线UE提供服务。RIS为配备个无源反射元件的均匀平面阵列，BS为配备个射频链(RadioFre-quency,RF)和根天线的均匀线性阵列，RIS上反射元件之间的间距或天线阵列之间的间距为半波长，即，为载波波长。为对抗多径效应，采用具有 个子载波的OFDM技术。BS和RIS安装位置固定，且考虑低移动性UE。考虑准静态块衰落信道，每个信道在个时隙的信道相干块中保持近似恒定16。假设每个相干块中个时隙用于上行信道估计，每个时隙包含 个导频符号。假设UE-BS之间的直接路径被大的障碍物阻塞6。qb在上行信道估计期间，在第 个子载波的第 个时隙BS的接收信号为yb,q=wHbHBR,qdiag(b)Uu=1HRU,u,qsu,b,q+nb,q(1)HRU,u,quqHBR,q=1ej1,2ej2,.,MejMTm 0,1m 0,2)mwb CNNRFbsu,b,q=su,b,q,1,su,b,q,2,.,su,b,q,Puqbnb,qsu,b,qsHu,b,q=RPsu1,b,qsHu2,b,q=0,u1=u2 Rqb其中，为第 个UE的第 个子载波的UE到RIS的信道矩阵，为RIS到BS的信道矩阵，为RIS的相移系数，和为第个反射元件的幅度和相移,为BS第 个时隙处的组合矩阵；为第 个UE的第个子载波在第 个时隙处发送的导频序列,为加性高斯白噪声矩阵。假设采用正交导频传输策略，即,为发射功率15。因此，在第 个子载波的时隙 处BS接收到的信号为e yu,b,q=wHbHBR,qdiag(b)HRU,u,q+e nu,b,q(2)e yu,b,q1RPyu,b,qsHu,b,qe nu,b,q1RPnu,b,qsHu,b,q其 中，,。3 宽带信道模型Bfcqfq=fc+B(q 1 (Q 1)/2)/Q,q QLRUu,l,mumlBRu,m,nmnum假设该系统的传输带宽为，载波频率为，第个子载波的频率为。本文只考虑BS和RIS之间的视距传输3,17。假设UE和RIS之间存在 条路径，表示第 个UE到RIS第个反射元件在第 条路径上的路径时延，表示RIS的第个反射元件到BS的第 根天线的路径时延。因此，第 个UE到RIS的第个反射元件的上行链路信道的时域响应为132444电子与信息学报第45卷hRU,u,m(t)=Ll=1RU,lej2fcRUu,l,m(t RUu,l,m)(3)RU,ll其中，为UE-RIS信道的第 条路径的信道增益。uq对式(3)进行傅里叶变换，则第 个UE的第 个子载波在UE-RIS信道整体频域信道响应为HRU,u,q=Ll=1RU,laRIS(l,l,fq)ej2fqRUu,l(4)RU,l RU,lej2fcRUu,llRUu,l RUu,l,0lllaRIS(l,l,fq)aRIS(,fq)其中，为第 条路径等效路径增益，且，和 分别为RIS第 路径的到达角(Angle-Of-Arrival,AOA)的仰角和方位角，为 R I S 的 阵 列 导 向 矢 量 矩 阵。表示为13aRIS(,fq)=1M1,ej2(1+fqfc)d(sin cos+sin sin),.,ej2(1+fqfc)d(Mr1)sin cos+(Mo1)sin sinT(5)uq同理，第 个UE的第 个子载波在RIS-BS信道整体频域信道响应为HBR,u,q=BRaBS(,fq)aHRIS(,fq)ej2fqBRu(6)BR BRej2fcBRuBRBRu BRu,0,0aBS(,fq)aBS(,fq)其中，为等效路径增益，为RIS-BS信道的路径增益，和 分别为RIS的离开角(Angle-Of-Departure,AOD)的仰角和方位角，为BS的AOD，为BS的阵列导向矢量矩阵。表示为13aBS(,fq)=1N1,ej2(1+fqfc)dsin,.,ej2(1+fqfc)(Q1)dsinT(7)定义HURB,u,q=HBR,qdiag(HRU,u,q)(8)HURB,u,quq其中，为第 个UE在第 个子载波处的级联信道。4 两阶段级联信道估计方案T1Tup T1针对RIS辅助的宽带多用户系统的两阶段信道估计方案的主要思想如下：在第1阶段中，选择一个典型用户，其传输了个时隙的正交导频序列，通过一种基于梯度下降法的迭代优化算法来估计AOA、信道增益和时延参数，得到典型用户的级联信道。在第2阶段中，其他用户将发送个时隙的正交导频序列，该阶段利用RIS信道的空间相关性和用户信道之间的相关性减少信道估计的训练开销。4.1 第1阶段：典型用户级联信道估计4.1.1 问题描述u=1T1q在这一阶段，选择最接近BS的UE作为典型用户，记为UE1()。该阶段只有典型用户在个时隙中向BS发送正交导频信号，则BS在第 个子载波处的接收信号为y1,qe yT1,q,1,e yT1,q,2,.,e yT1,q,T1T=WHLl=1glA1,q(,fq)aRIS(l,l,fq)aBS(,fq)ej2fq1,l+n1,q(9)gl RU,lBRl1,l BR1+RU1,llW blkdiag(w1,w2,.,wT1)A1,q(,fq)aHRIS(,fq)diag(1)T,aHRIS(,fq)diag(2)T,.aHRIS(,fq)diag(T1)TTn1,qe nT1,q,1e nT1,q,2,.,e nT1,q,T1T其中，为第 条路径的级联信道增益，为UE1的第 条路径的级联信道时延，,。QBS处 个子载波的接收信号为Y1yT1,1,yT1,2,.,yT1,QT=WHLl=1gll+N1(10)W blkdiag(W,W,.,W)CQT1NQT1NN1nT1,1,nT1,2,.,nT1,QTlA1,1(,f1)aRIS(l,l,f1)aBS(,f1)ej2f11,lTA1,2(,f2)aRIS(l,l,f2)aBS(,f2)ej2f21,lT,.A1,Q(,fQ)aRIS(l,l,fQ)aBS(,fQ)ej2fQ1,lTT其中，,。gl,l,l,1,lLl=1由式(10)可知，级联信道的CSI由3种参数决定，即路径增益，角度信息和时延。此外，BS和RIS的位置固定，即,和 完全可以由其几何位置确定17。因此，可将级联信道估计问题转化为参数恢复问题，即估计参数集。将式(10)改写为Y1=WHKg+N1(11)g g1,g2,.,gLTK 1,2,.,L其中，为信道增益矩阵，为角度参数和时延参数组成的矩阵。gl,l,l,1,lLl=1Y1由于信道参数的数量远小于的维数，可将式(11)的信道估计问题表示为minb g,b,b,b b g0,s.t.?Y1 WHKg?F(12)b g g1,g1,.,gbLT其中，为估计的级联信道增益，第7期李贵勇等：可重构智能表面辅助的多用户通信宽带信道估计2445b 1,2,.,bLTb 1,2,.,bLTb v1,1,v1,2,.,v1,bLTbL 0,和为估计的角度参数矩阵，为估计的路径数，为误差容错参数。4.1.2 基于梯度下降的迭代加权算法求解式(12)的主要困难在于优化是一个NP-hard问题。文献18提出对数和函数对于稀疏信号恢复问题具有更强的稀疏性。利用对数和函数替换式(12)中的0范数，得到对数和最小化目标函数为ming,Ll=1lg(|gl|2F+),s.t.?Y1 WHKg?F(13)0其中，是一个较小的正则化常数。显然，式(13)中的优化问题是非凸的18。通过添加数据拟合项将问题转化为如式(14)的无约束优化问题19ming,(g,)=Ll=1lg(|gl|2F+)+?Y1 WHKg?2F(14)0其中，为控制稀疏性和数据拟合误差之间平衡的正则化参数。(g,)为得到不增的目标函数值，利用优化最小化(MajorizationMinimization,MM)方法求解式(14)，该方法通过迭代最小化一个简单的凸代理函数以优化给定的目标函数20。借鉴文献18的思想，将的代理函数建立为J(g|b g(i)=Ll=1lg(?g(i)l?2F+)+|gl|2F+?g(i)l?2F+1Ll=1lg(|gl|2F+)(15)b g(i)gi其中，为 第 次迭代时的估计值。J(g|b g(i)C=|gl|2F+0 x=?g(i)l?2F+0k(x)=Ll=1lgx+C/x 1k(x)=(x C)/x3k(x)=(2C x)/x3k(x)0J(g|b g(i)证明为证明式(15)中为凸函数，令常 数 项，变 量，则有,，易证，故为凸函数。证毕(g,)将式(15)代入式(14)中，的最小化等价于代理函数的最小化，即ming,G(i)(g,)=1gHD(i)g+?Y1WHKg?2F(16)G(i)(g,)i其中，为第 次迭代时的迭代代理函D(i)=diag(1/(?b g(i)1,1?2+),1/(?b g(i)1,2?2+).,1/(?b g(i)1,L?2+)数，。gb gG(i)与文献19相似，可以优化式(16)中的，以找到最佳点 和对应的的最佳值，则有g(i)opt(,)=(1D(i)+(WHK)HWHK)1(WHK)HY1(17)将式(17)代入式(16)中，最优目标函数可表示为G(i)opt(,)=YH1WHK(1D(i)+(WHK)HWHK)1(WHK)HY1+Y1HY1(18)ib(i+1)b(i+1)b(i+1)G(i)opt(,)最终，得到以最小化式(18)中的目标函数为目标的迭代算法。在所提参数估计算法中，在第 次迭代时，利用梯度下降法找到新的参数,和使得目标函数变得更小。在迭代过程中，动态正则化参数 随着迭代的进行不断更新，更新方式为=max(i),min)(19)(i)(i)=?Y1 WHK(i),(i),(i)g(i)?2F其中，为上一步的平方残差，为常数缩放因子，。LL?g(i+1)l?2F gth,l 1,2,.,bL(i)gthb g(i+1)g(i+1)blb(i+1)b(i+1)b(i+1)b(i+1)blb(i+1)blb v(i+1)1,blbL(i+1)b(i+1)考虑到实际场景中，路径数是未知的。故所提信道估计算法中，将路径数 初始化为一个较大的值，且与真实路径数保持相同的数量级。在迭代过程中，将路径增益过小的路径视为噪声。基于此的 一 种 剪 枝 策 略 为：如 果,为修剪阈值，则从中删除，从,和中删除,和，并更新使其与的长度相等。基于梯度下降的迭代加权算法的流程如下：thgthb g(0)=g(0)opt(b(0),b(0),b(0)步骤1设置终止阈值，修剪阈值，根据式(17)初始化。步骤2根据式(19)更新。G(i)opt(,)步骤3根据式(18)构造函数。b(i+1)b(i+1)b(i+1)步骤4利用梯度下降算法更新参数,和。b g(i+1)步骤5根据式(17)估计路径增益。?g(i+1)l?2F gthbL(i+1)步骤6如果，剪枝并更新路径数。?b g(i+1)b g(i)?2 th步骤7如果，则返回步2446电子与信息学报第45卷b gbbb cHURB,1,q,q 1,2,.,Q骤2；否则，得到,，由式(8)得到。4.2 第2阶段：其他用户级联信道估计HRU,u,q,mumqu,mm由于每个RIS反射元件通过相同的RIS-BS信道将来自不同用户的信号反射到BS，所以可将任意用户的反射信道视为其他用户的反射信道的缩放版本21。因此，在信道估计期间，只需要估计缩放因子(标量)，而不需要估计整个信道矢量。在该阶段，令表示第 个UE到RIS的第个反射元件在第 个子载波处的频域信道响应，表示与RIS的第个反射元件相关的缩放因子，利用不同UE的UE-RIS信道之间的强相关性14，可得u,m,q=HRU,u,q,mHRU,1,q,m,u=2,3,.,U,m=1,2,.,M(20)HRU,u,q对整个信道，则有HRU,u,q=HRU,1,q u,q(21)u,q u,1,q,u,2,q,.,u,M,qTHRU,u,q其中，为与相关的缩放矢量。由式(8)和式(21)可得HBR,qdiag(b)HRU,u,q=HURB,1,b,qu,q(22)HURB,1,b,q HURB,1,qdiag(b)其中，。HURB,u,qUu=2HURB,1,qu,qUu=2HURB,1,q,u,qUu=2HURB,u,qUu=1U 1T2=Tup T1q利用式(21)和式(22)的关系，可以由和的信息估计得到。因此，通过估计而不是估计，从而显著减少待估计的未知信道系数的数量。将式(22)代入式(1)并向量化，在个用户发送了时隙的正交导频序列后，BS接收到第 个子载波的信号为yqvec(y1,q),vec(y2,q),.,vec(yT2,q)T=qq+zq(23)qB2,1,q,B2,2,q,.,B2,T2,qT,.,BU,1,q,BU,2,q,.,BU,T2,qTBu,b,qsTu,b,q(wHbHURB,1,b,q)q 2,q,3,q,.,U,qTzq vec(n1,q),vec(n2,q),.,vec(nT2,q)T其中，,。QBS处 个子载波的接收信号为Y y1,y2,.,yQT=+Z(24)blkdiag1,2,.,QCT2PQNRF(U1)QM 1,2,.,PT C(U1)QM1Z z1,z2,.,zQT其中，,。T2PQNRF(U 1)QMrank()=文献22表明，使得式(24)存在唯一解的条件为：(1)；(2)(U 1)QMT2P (U 1)M/NRF。因此，该阶段传输的导频数应满足：。由LS算法可得b=()Y(25)b其中，为 的估计值。5 最小导频开销与计算复杂度分析5.1 最小导频开销分析Lb g3Lbbb T1 4L(U 1)M/NRF4LP+(U 1)M/NRF在第1阶段中，基于梯度下降的迭代加权算法估计了个级联信道增益，及个的角度参数,和，因此，在该阶段的导频开销为；在第2阶段中，因为需要确保式(24)存在唯一解，所需最小导频数目为。因此，提出的两阶段级联信道估计方案最小导频开销为。ULlgMULlg(L2EUEERIS)EUEERISUTupP N/NRFM/MRFMRFM+max(U 1,(U 1)M/NRF)2(M+1)/+U M/N+U 1对比方案中，文献4最小导频开销为，文献5最小导频开销为，其中和为角分辨率，文献6最小导频开销为，其中为RIS处RF个数，文 献 8 最 小 导 频 开 销 为，文 献 1 4 最 小 导 频 开 销 为，其中，表示一个常量系数。相比之下，所提方案的最小导频开销小于文献4和文献6。5.2 计算复杂度分析O(QT1NRFL2(QT1N+L)O(Q3U3M3+2Q3U2M2T2PNRF+Q2T2PNRFUM)O(QT1NRFL2(QT1N+L)+O(Q3U3M3+2Q3U2M2T2PNRF+Q2T2PNRFUM)在第1阶段的基于梯度下降的迭代加权算法中，复杂度主要来自步骤4中梯度计算，计算复杂度为；第2阶段利用了LS算法，复杂度主要来自式(25)中的伪逆计算，计 算 复 杂 度 为。因此，所提两阶段级联信道估计 方 案 总 复 杂 度 为。O(t3U+2P3t2UT2PNRF+Q2UT2PNRFt)O(Q2UT2PNRFt)+O(Q2UT2PNRFI)tIO(2t3UQTP+4Q2t2T2P2UN+2Q2T2P2UNt)O(T2P2QUNEUEERISEBS)O(T2P2QUEBS)EBSO(UQTP(M+N)3.5)O(Q3N3M3+2Q3N2M2T1PNRF+P2NMT1PNRF)O(Q3U3M3+2Q3U2M2T2PNRF+Q2T2PNRFUM)O(Q3N3M3U+2Q3N2M2UT2PNRF+Q2UT2PNRFNM)O(M3+M2NL+MNLt)对比方案中，文献4主要计算复杂度为+，其中 为迭代次数，为量化网格数，文献5主要计算复杂度为+，其中为角分辨率，文献6主要计算复杂度为，文献8主要计算复杂度为+，文 献 1 4 主 要 计 算 复 杂 度 为+。虽然本文提出的两阶段级联信第7期李贵勇等：可重构智能表面辅助的多用户通信宽带信道估计2447LO(PT1NRFL2(QT1N+L)道估计方案在第2阶段也涉及了伪逆计算，但在第1 阶 段 中，因 为 路 径 数较 小，所 以 复 杂 度较小。因此，相比于对比方案，所提方案的总的复杂度更低。6 仿真分析N=256M=256U=8L=3B=20 GHzfc=150 GHzQ=256本节探究所提两阶段级联信道估计方案的可行性和优越性。参数设置如下：BS处的天线数，R I S 处 的 天 线 数，U E 数 量，UE和RIS之间的路径数，传输带宽，载波频率，子载波数，信噪比(SignaltoNoiseRatio,SNR)为5dB。与本文所提两阶段信道估计方案的对比算法中，选择未考虑波束斜视影响的RIS辅助通信的窄带和宽带信道估计方案：基于压缩感知的分布式正交匹配追踪算法的文献4、基于压缩感知的3阶段信道估计方案的文献5、基于混合RIS结构的2阶段信道估计方案的文献6、多用户信道估计方案的文献8和基于坐标下降信道估计的文献14。图1给出了不同SNR下不同信道估计方案的NMSE性能对比。由于文献4和文献5中的信道估计方案从离散网格中选择估计角度，网格失配使得其性能不如其他算法，且网格字典大小的限制使得随着SNR的增加估计精度逐渐达到饱和。相比之下，所提2阶段级联信道估计方案中基于梯度下降的迭代加权算法可以避免量化误差，NMSE性能明显好于对比的5种信道估计方案。BBB图2给出了不同信道估计方案的NMSE性能与传输带宽 的关系。由图2可知，随着传输带宽 的增加，波束斜视的影响逐渐增大，文献4和文献5中的信道估计方案的性能越来越差。这是因为，对于较大的传输带宽，不同子载波处的宽带波束空间信道的支撑集差异性更大，导致公共稀疏信道支BB撑集假设不再成立。文献6、文献8和文献14中的信道估计方案的性能不随传输带宽 的增加而降低，因为这3种信道估计方案是逐个估计每个子载波对应的信道。相反，所提2阶段级联信道估计方案对传输带宽 是鲁棒的。U图3给出了不同UE数量下各种信道估计方案的信道估计性能关系。由图3可知，随着UE数量的逐渐增加，所提信道估计方案的性能有所下降，当UE数量大于20时，所提方案的性能不如文献6和文献14的信道估计方案。这是因为，所提2阶段级联信道估计方案的性能与典型用户的估计准确度密切相关，随着用户数量的增加，其他用户利用信道间的强相关性进行信道估计的信道误差将会逐渐增加。PTup图4给出了不同信道估计方案的导频长度与NMSE性能的关系。由图4可知，随着导频符号长度的增加，所有信道估计方案的准确度都会提高，因为通过大量的测量值可以获得更好的恢复和图1不同信噪比下的归一化均方误差图3不同UE数量下的归一化均方误差图2不同传输带宽下的归一化均方误差2448电子与信息学报第45卷PTup估计精度。其次，对于相同导频长度，所提两阶段级联信道估计方案始终比其他信道估计方案具有更好的NMSE性能。7 结束语本文研究了波束斜视影响下RIS辅助宽带多用户通信系统中的上行信道估计问题，并利用THz信道的稀疏性，提出一种两阶段级联信道估计方案。首先，将信道估计问题转化为无约束优化问题，并利用基于梯度下降法的迭代优化算法来估计AOA、信道增益和时延参数；然后，基于信道间的强相关性估计其他用户的级联反射信道。仿真结果表明，与传统宽带算法相比，本算法具有更好、更稳定的信道估计性能。针对信道估计中用户数增多导致性能下降和矩阵伪逆的计算导致较高复杂度的问题，后续工作将展开对用户进行分组和降低复杂度的分析与研究。参 考 文 献CHENHui,SARIEDDEENH,BALLALT,et al.Atutorialonterahertz-bandlocalizationfor6GcommunicationsystemsJ.IEEE Communications Surveys&Tutorials,2022,24(3):17801815.doi:10.1109/COMST.2022.3178209.1李斌,刘文帅,谢万城,等.智能超表面赋能移动边缘计算部分任务卸载策略J.电子与信息学报,2022,44(7):23092316.doi:10.11999/JEIT211595.LIBin,LIUWenshuai,XIEWancheng,et al.Partialcomputationoffloadingfordouble-RISassistedmulti-usermobileedgecomputingnetworksJ.Journal of Electronics&Information Technology,2022,44(7):23092316.doi:10.11999/JEIT211595.2WANZiwei,GAOZhen,GAOFeifei,et al.TerahertzmassiveMIMOwithholographicreconfigurableintelligentsurfacesJ.IEEE Transactions on Communications,2021,369(7):47324750.doi:10.1109/TCOMM.2021.3064949.WANZiwei,GAOZhen,andALOUINIMS.BroadbandchannelestimationforintelligentreflectingsurfaceaidedmmWavemassiveMIMOsystemsC.ICC2020-2020IEEEInternationalConferenceonCommunications(ICC),Dublin,Ireland,2020:16.doi:10.1109/ICC40277.2020.9149146.4LIN Tian,YU Xianghao,ZHU Yu,et al.ChannelestimationforIRS-assistedmillimeter-waveMIMOsystems:Sparsity-inspiredapproachesJ.IEEE Transactions onCommunications,2022,70(6):40784092.doi:10.1109/TCOMM.2022.3168876.5SCHROEDERR,HEJiguang,BRANTEG,et al.Two-stagechannelestimationforhybridRISassistedMIMOsystemsJ.IEEE Transactions on Communications,2022,70(7):47934806.doi:10.1109/TCOMM.2022.3176654.6ZHENGBeixiong,YOUChangsheng,MEIWeidong,et al.AsurveyonchannelestimationandpracticalpassivebeamformingdesignforintelligentreflectingsurfaceaidedwirelesscommunicationsJ.IEEE Communications Surveys&Tutorials,2022,24(2):10351071.doi:10.1109/COMST.2022.3155305.7HUChen,DAILinglong,HANShuangfeng,et al.Two-timescalechannelestimationforreconfigurableintelligentsurface aided wireless communicationsJ.IEEETransactions on Communications,2021,69(11):77367747.doi:10.1109/TCOMM.2021.3072729.8JEONGS,FARHANGA,PEROVICNS,et al.Low-complexityjointCFOandchannelestimationforRIS-aidedOFDMsystemsJ.IEEE Wireless Communications Letters,2022,11(1):203207.doi:10.1109/LWC.2021.3124049.9LIU Shicong,GAO Zhen,ZHANG Jun,et al.DeepdenoisingneuralnetworkassistedcompressivechannelestimationformmWaveintelligentreflectingsurfacesJ.IEEE Transactions on Vehicular Technology,2020,69(8):92239228.doi:10.1109/TVT.2020.3005402.10LINYuxing,JINShi,MATTHAIOUM,et al.Tensor-basedalgebraicchannelestimationforhybridIRS-assistedMIMO-OFDMJ.IEEE Transactions on WirelessCommunications,2021,20(6):37703784.doi:10.1109/TWC.2021.3053580.11SONG Yuhui,GONG Zijun,CHEN Yuanzhu,et al.Efficientchannelestimationforwidebandmillimeterwavemassive MIMO systems with beam squintJ.IEEETransactions on Communications,2022,70(5):34213435.doi:10.1109/TCOMM.2022.3156065.12WANGBolei,GAOFeifei,JINShi,et al.Spatial-andfrequency-widebandeffectsinmillimeter-wavemassiveMIMOsystemsJ.IEEE Transactions on Signal Processing,13图4不同导频数量下的归一化均方误差第7期李贵勇等：可重构智能表面辅助的多用户通信宽带信道估计24492018,66(13):33933406.doi:10.1109/TSP.2018.2831628.WANGZhaorui,LIULiang,andCUIShuguang.Channelestimation for intelligent reflecting surface assistedmultiusercommunications:Framework,algorithms,andanalysisJ.IEEE Transactions on WirelessCommunications,2020,19(10):66076620.doi:10.1109/TWC.2020.3004330.14MASiqi,SHENWenqian,ANJianping,et al.WidebandchannelestimationforIRS-aidedsystemsinthefaceofbeam squintJ.IEEE Transactions on WirelessCommunications,2021,20(10):62406253.doi:10.1109/TWC.2021