一种基于稀疏张量分解的RIS-MIMO系统信道估计方法

技术领域

本发明属于无线通信技术领域，特别涉及一种基于稀疏张量分解的 RIS-MIMO系统信道估计方法。

背景技术

随着通信网络容量的数倍增长，海量的移动数据终端无处不在，各种新兴 技术不断被提出和实践，比如人工智能、物联网、车联网和云计算等高新发展 技术。5G致力于结合超密集异构网、毫米波(mmWave)和波束赋形等技术， 为用户提供更高速率和低时延的网络服务。但在此发展现实下，未来无线通信 所面临的挑战也随之而来，即复杂的网络结构、高成本的通信硬件部署及不断 增长的设备能耗等关键问题。因此，迫切需要研发高能效、低成本的绿色方案 来解决未来无线网络问题。在众多候选的新兴通信技术中，可重构智能表面 (RIS)因其极低功耗、低成本、可编程及易于部署等显著特点而被众多学者青 睐，成为脱颖而出的下一代移动通信研究热点。RIS平面往往是布置在大型的建 筑物表面，由大量的可调节无源反射元件构成，通过合理地数字系统设计，这 些元件可以改变自身的大小、形状和方向等，动态地改变入射信号的角度及幅 度等信息，从而实现人为干预信号传播环境，获得理想的传播信号，在改善通 信盲区信号质量差、无线传播距离有限、高速移动的数据交换、多径效应及海 量数据传输方面起着重要作用。但是，由于RIS中的单元大多数为无源的，即 只用来完成电磁波的反射，并不能实现相应的信号处理。同时，由于RIS往往 配备数量众多的无源元件，这导致相关联得信道往往具有很高的维度。因此， 对于RIS辅助的无线通信MIMO信道估计研究就具有一定的挑战性，如何设计 高效的信道估计算法也成为众多学者研究兴趣所在。

考虑到未来MIMO场景下的信道多数都会建模为mmWave信道，并在角度 域往往具有明显的角信道稀疏性。针对mmWave信道的稀疏特性，有研究人员 利用利用不同子载波上信道稀疏度，在减少导频开销的情况下，提出了一种基 于压缩感知的RIS辅助mmWave大规模MIMO系统的信道估计方案。有学者研 究了RIS辅助的大规模MIMO系统的半盲级联信道估计，利用近似消息传递原 理，开发了一种迭代算法来解决三线性推理问题，其中接收机同时估计信道系 数和部分未知发射信号。有学者研究了RIS辅助mmWave系统的信道估计和联合波束形成设计问题，利用级联信道固有的稀疏结构，提出了一种基于压缩感 知的信道估计方法。有学者利用RIS信道的单结构稀疏性，即不同用户的稀疏 角域通道仅在列维上具有相同的结构稀疏性，而在行维上没有相同的结构稀疏 性，提出了一种低开销降维信道估计方案。但上述信道估计算法依然具有较高 的复杂度和训练开销，因此，利用RIS辅助下MIMO系统中的mmWave信道的 本身特点，即高维度和显著的稀疏特性，设计相应的估计算法有可能会进一步 降低开销和复杂度，从而达到提高频谱效率和系统容量的目的。

发明内容

本发明为解决现有技术中的问题而提出了一种基于稀疏张量分解的 RIS-MIMO系统信道估计方法，改善了估计性能，降低了复杂度和训练开销， 且收敛速度快，稳定性好。为达到上述目的所采取的技术方案是：

一种基于稀疏张量分解的RIS-MIMO系统信道估计方法，包括如下步骤：

步骤A、对RIS辅助的多用户MIMO上行mmWave通信系统进行物理建模， 构造张量化观测模型；

步骤B、对张量化的系统模型进行平行因子分解，得到二维矩阵表达式；

步骤C、根据mmWave信道模型，将信道估计问题转化为结构化稀疏信号 恢复问题；

步骤D、使用LASSO方法构造相应的优化问题，设置初始化参数值，通过 交替迭代算法框架估计稀疏信道的相关参数，以实现对信道矩阵的重构

步骤E、重复步骤D，直到算法收敛。

优选的，步骤A具体包括：

根据单小区单用户MIMO上行链路通信系统，假设RIS平面有N个反射单 元，第m个用户到RIS和基站到RIS间的信道向量/矩阵分别表示为

Y

其中，

其中，

且

优选的，步骤B具体包括：

根据张量化观测模型，

其中，

模式-1：

模式-2：

模式-3：

优选的，步骤C具体包括：

根据mmWave信道模型特征，信道向量/矩阵

g

因此，M个用户和RIS间的信道矩阵

G＝F

其中，

Z

其中，

至此，关于信道矩阵H和G估计被转化为针对B和Q的稀疏信号恢复问题。

优选的，步骤D具体包括：

步骤D1、在最小二乘方法的基础上，使用l

步骤D2、设置合理的初始化参数值，通过交替迭代算法框架估计稀疏信道 的相关参数，以实现对信道矩阵的重构。

优选的，步骤D1具体包括：

向量化Z

其中，

s.t.||ξ

s.t.||ξ

优选的，步骤D2具体包括：

将(P1)与(P2)等价为：

通过交替求解目标函数来获得关于ξ

最终可获得信道矩阵G和H的估计值：

本发明所具有的有益效果为：本发明根据RIS辅助的多用户MIMO上行 mmWave通信系统进行物理建模，构造具有稀疏结构的张量化观测模型,对张量 进行平行因子分解将信道估计转化为稀疏信号恢复问题,在最小二乘方法的基础 上，利用LASSO方法提出一种基于稀疏张量分解的mmWave信道估计迭代算 法。数值仿真实验探究了稀疏信道中稀疏度及各个参数对估计性能的影响，验 证了所提算法的显著性能优势，复杂度和训练开销相对较少，且收敛速度快， 稳定性好。为进一步提升频谱效率提供了新思路。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面对实施例 或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中 的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出 创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明的流程图；

图2为相位配置状态数量K不同时，不同信道估计算法的NMSE随信噪比 SNR变化的示意图；

图3为RIS反射单元个数N不同时，不同信道估计算法的NMSE随信噪比 SNR变化的示意图；

图4为稀疏程度η不同时，不同信道估计算法的NMSE随信噪比SNR变化 的示意图；

图5为信道估计算法的NMSE随η变化示意图；

图6为本发明的系统结构图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明 实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然， 所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发 明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有 其他的实施例，都属于本发明保护的范围。

如图1所示，一种基于稀疏张量分解的RIS-MIMO系统信道估计方法，包 括如下步骤：

步骤A、对RIS辅助的多用户MIMO上行mmWave通信系统进行物理建模， 构造张量化观测模型；

步骤B、对张量化的系统模型进行平行因子分解，得到二维矩阵表达式；

步骤C、根据mmWave信道模型，将信道估计问题转化为结构化稀疏信号 恢复问题；

步骤D、使用LASSO方法构造相应的优化问题，设置初始化参数值，通过 交替迭代算法框架估计稀疏信道的相关参数，以实现对信道矩阵的重构

步骤E、重复步骤D，直到算法收敛。

步骤A具体包括：

考虑一个由RIS辅助的多用户上行毫米波MIMO无线通信系统。其中，为 了提高基站和用户之间的通信质量，部署了一个由N个无源反射元件组成的 RIS，基站(BS)配置L根天线，与M个单天线用户(User)间完成通信传输。 此外，假设RIS上的N个反射元件之间是相互独立的。由于传播环境中的不利 条件，BS和Users之间的直接链路信道可假设为被忽略，重点讨论如何有效地 估计BS-RIS和RIS-Users之间的信道。同时，假设所有的信道都为准静态平坦 衰落信道模型，即在T个时隙内信道状态信息保持不变。

第m个用户到RIS和基站到RIS间的信道向量/矩阵分别表示为

其中，t∈[1,T]，x

其中，s

y

Y

其中，

为了更简洁地描述系统模型，假设导频信号矩阵X是一个正交的导频信号 矩阵，且满足XX

其中，

且

步骤B具体包括：

根据张量化观测模型，

利用n-模式积运算，

其中，

模式-1：

模式-2：

模式-3：

步骤C具体包括：

根据mmWave信道模型特征，BS-RIS间的信号矩阵建模为：

其中，P为BS-RIS间总的传输路径数，ρ

那么，与BS的天线阵列几何排列相关的导向矢量可表示为:

同样，由于RIS配置均匀平面天线阵列，所以相应的传输天线阵列响应为：

其中，

信道矩阵

其中，

其中，

其中，P′和ε

g

同时，

G＝F

其中，

Z

其中，

至此，关于信道矩阵H和G估计被转化为针对B和Q的稀疏信号恢复问题。

步骤D具体包括：

步骤D1、在最小二乘方法的基础上，使用l

步骤D2、设置合理的初始化参数值，通过交替迭代算法框架估计稀疏信道 的相关参数，以实现对信道矩阵的重构。

步骤D1具体包括：

向量化Z

其中，

s.t.||ξ

s.t.||ξ

步骤D2具体包括：

将(P1)与(P2)等价为：

其中，τ和β均为大于等于0的正则化参数，且τ和β可以调节稀疏程度，即正 则化参数的值越大，越多的参数将会被收缩到0，因此，正则化参数的取值将会 影响估计算法的估计性能。通过交替求解目标函数来获得关于ξ

最终可获得信道矩阵G和H的估计值：

下面将给出本发明的传输方案与现有的其它传输方案的比较，以使本发明 的优势及特征更加明显。

仿真参数的设定：最大迭代次数设置为20次，阙值ξ＝10

图2为相位配置状态数量K不同时，不同信道估计算法的NMSE随信噪比 SNR变化的示意图，其中，η＝1.5，L′＝24，M＝L＝K＝T＝16，N

图3为RIS反射单元个数N不同时，不同信道估计算法的NMSE随信噪比SNR 变化的示意图。其中，图3中的左边图，即图(a)，为

图4是η不同时，信道估计算法的NMSE随SNR变化的曲线示意图，η分别取 {1.25,1.5,1.75}，M＝L＝T＝K＝N＝16。从图中可以看出，对于估计信道矩阵

图5表示的是当SNR＝25dB，η分别取{1,1.25,1.5,1.75,2}，M＝L＝T＝K＝N＝16时，STF-CE算法和基于ALS的算法的NMSE曲线对比图，从图中可以看出，在 不同的η下，所提STF-CE算法的估计性能也是要显著优于其他估计方法，这与 图4形成呼应，进一步充分说明了不同的η下STF-CE算法的估计性能的优势。

图6为本发明的系统结构图，是一种基于UAMP的噪声精度和信道估计的系 统结构图，包括：初始化模块1、信道估计模块2以及噪声精度估计模块3。具体 来说：

初始化模块1，为了上述迭代系统能顺利启动，对相关参数的初始化设置；

参数估计模块2，基于预设的初始值，通过交替迭代算法框架估计稀疏信道 的相关参数；

信道重构模块3，基于上一模块得到的估计结果，对信道矩阵完成重构，并 更新信道估计值以便下一次迭代。

具体的基于稀疏张量分解的RIS-MIMO系统中信道估计方法中的参数估计 和信道信道重构过程可参见上述实施例，本发明实施例在此不再赘述。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限 制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员 应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其 中部分技术特征进行等同替换，但这些修改或者替换，并不使相应技术方案的 本质脱离本发明各实施例技术方案的精神和范围。