CN114338302B - 一种基于毫米波联合结构的智能反射面两级信道估计方法 - Google Patents

Abstract

本发明公开了一种基于毫米波联合结构的智能反射面两级信道估计方法，该所方法包括：构建智能反射面辅助的上行毫米波MIMO系统，采集不同时刻的接收信号；推导毫米波信道具有角扩展的波束空间模型，证明其联合稀疏与低秩结构；将时域接收信号模型重构为低秩矩阵采样形式，提出改进的秩算子，利用低秩矩阵近似算法完成接收信号的去噪处理；定义级联信道虚波束空间矩阵，引入vec算子及Kronecker积，将信道模型转换为稀疏信号恢复问题；联合恢复的虚波束空间矩阵与毫米波波束空间模型，完成直接信道与级联信道的最终估计。本发明在IRS辅助的毫米波MIMO系统中提高信道估计精度、降低信道估计开销，也同时实现了对IRS反射系数矩阵的优化。

Description

一种基于毫米波联合结构的智能反射面两级信道估计方法

技术领域

本发明涉及无线通信技术领域，具体涉及一种基于毫米波联合结构的智能反射面两级信道估计方法。

背景技术

随着6G移动通信网络对通信质量的更高需求，智能反射面作为增强通信的一种有效技术受到了广泛的关注。IRS是一种由大量微单元组成的无源反射面板，可以通过控制传输信号的相位和幅度来灵活配置无线传输环境，无需信号解码/放大处理地实现全双工无源波束赋形。因此，在理论上，IRS不仅可以有效地提高网络传输速率、通信覆盖范围、频谱效率等性能，还可以降低硬件成本和能耗，具有杰出的性能和广阔的前景。

现有研究分析了IRS系统中通信资源分配框架和无线携能通信系统等设计，验证了IRS在提高能效方面的重要理论作用。而广泛研究的基于IRS的波束赋形技术被证明其可有效提高通信系统的可靠性。IRS还应用于物理层服务集成系统中，辅助业务消息集成与接入点功率分配优化。值得注意的是，上述对IRS的研究都需要信道状态信息的支持。因此，有效的信道估计方法是IRS发挥其优异性能的重要条件之一。

然而，IRS系统的信道估计问题面临巨大的挑战。一方面，由于IRS反射单元是无源的，没有发送或接收导频信号的能力，传统有源设备广泛研究的传输训练序列方法不再适用。此外，IRS上丰富的反射微单元会带来更多的训练开销和更高的估计复杂度，需要被进一步降低。由于毫米波具有联合稀疏与低秩结构，可有效降低信道估计的导频开销，因此，验证这种联合结构并引入IRS系统，继而分别基于低秩性与稀疏性完成信道估计具有广阔的前景。

发明内容

针对现有技术的不足，本发明提供一种基于毫米波联合结构的智能反射面两级信道估计方法，构建了具有联合稀疏与低秩结构的毫米波角扩展信道模型，并应用于IRS系统进行两级信道估计，第一阶段基于低秩矩阵近似算法，提出改进的秩算子对接收信号进行去噪处理，第二阶段基于稀疏信号恢复算法，通过合并级联信道的虚波束空间对直接信道与级联信道进行估计，从而最终在IRS辅助的毫米波MIMO系统中提高信道估计精度、降低信道估计开销，也同时实现了对IRS反射系数矩阵的优化。

本发明的第二目的在于提供一种基于毫米波联合结构的智能反射面两级信道估计系统。

本发明的第三目的在于提供一种存储介质。

本发明的第四目的在于提供一种计算设备。

本发明的目的可以通过采取如下技术方案达到：

一种基于毫米波联合结构的智能反射面两级信道估计方法，包括下述步骤：

构建智能反射面辅助的上行毫米波MIMO系统，采集不同时刻的接收信号：

构建毫米波信道具有角扩展的波束空间模型，表现出联合稀疏与低秩结构；

第一阶段：将时域接收信号模型重构为低秩矩阵采样形式，利用低秩矩阵近似算法完成接收信号的去噪处理：

第二阶段：定义级联信道的虚波束空间矩阵，引入vec算子及Kronecker积，将信道模型转换为稀疏信号恢复问题：

从联合恢复的虚波束空间矩阵与毫米波波束空间模型完成直接信道与级联信道的最终估计。

作为优选的技术方案，所述构建智能反射面辅助的上行毫米波MIMO系统，具体步骤包括：

在配备IRS的点对点上行mmWave MIMO系统中，基站端和用户端分别设置N_BS和N_US个天线，IRS有M个反射单元，系统中有US-BS、US-IRS和IRS-BS三条信道，在每条信道的发射端和接收端之间，采用仅需一个射频链的模拟发射波束赋形和接收组合结构；

在时刻t，发射端通过波束赋形矩阵发射元素全为1的符号向量经过直接信道H_d与级联信道GΘH_r后与接收组合向量/>结合成为最终时域接收信号：

y(t)＝z^H(t)(H_d+GΘH_r)F(t)s(t)+ω(t)

其中，为IRS的反射系数矩阵，ω(t)为均值0，方差σ²的加性高斯噪声。

作为优选的技术方案，所述构建毫米波信道具有角扩展的波束空间模型，具体步骤包括：

构建毫米波信道的角扩展模型：当AoA域上的角分布来自单个公共AoD射线时，信道模型为：

其中，a_A和a_D是AoA和AoD对应的阵列响应向量，θ是平均AoA，φ是该公共AoD，I是AoD引起的AoA域路径数目，α_i是第i条路径的增益，ν_i是第i条路径相对于平均AOA的偏移量；

信道表示为：

其中a_A和a_D是AoA与AoD对应的阵列响应向量，θ和φ是平均AoA与AoD，I和I′是AOA域的路径数目，α_i和α′分别是第i条路径的增益，ν_i是第i条路径相对于平均AOA的偏移量，两个AoD间距具有相似的功率角模式，即信道模型简化为：

毫米波的几何信道模型表示为：

其中a_A和a_D是AoA与AoD对应的阵列响应向量，L为毫米波系统内聚类数，在每个聚类l中，θ和φ是平均AoA与AoD，I和I′是AOA域的路径数目，α_i和α′分别是第i条路径的增益，ν_i和分别是第i和j条路径相对于平均AoA和AoD的偏移量；

a_A和a_D为：

其中λ是载波的波长，d是天线间距等于波长的一半；

在构建了毫米波具有角扩展的几何信道模型后，将几何信道模型转换为波束空间形式：

其中和/>分别表示AoA域和AoD域上的路径增益的集合，和/>为过完备矩阵，其每列分别对应一个由预离散AoA/AoD参数化的导向向量，虚波束空间信道H_v是L个稀疏矩阵的和。

作为优选的技术方案，所述将时域接收信号模型重构为低秩矩阵采样形式，具体步骤包括：

假设和/>为预先准备的波束赋形矢量和接收组合矢量的码本，分别有N_Z和N_F个元素，z(t)和f(t)从中随机不重复地选择，收集/>和/>中所有元素，分别用/>和表示，接收信号用低秩矩阵抽样形式重构为：

其中Y_ij是Y的第ij个元素，Ω表示观察到的集合；

定义一种改进的算子：P_δ(Γ)＝Γ(Γ^HΓ+δI)^-1Γ^H,δ≥0，其中为待优化的矩阵，tr()为矩阵的秩，当矩阵满秩时，δ＝0；当矩阵低秩时，具有性质：

即当系数δ无限趋近于0时，该改进的秩算子无限趋近于原矩阵的秩；

通过改进的秩算子完成低秩矩阵近似运算，完成接收信号的去噪处理：

其中，ε₁为精确阈值，Y_Ω为观测到的待去噪接收信号，为去噪后的纯净接收信号。

作为优选的技术方案，所述将信道模型转换为稀疏信号恢复问题，具体步骤包括：

将毫米波的角扩展波束空间模型引入IRS系统：

定义级联信道的虚波束空间模型：

即合并信道H_d、信道G和IRS的反射系数矩阵Θ；

联合毫米波的角扩展波束空间模型与级联信道的虚波束空间模型，将低秩矩阵抽样模型表示为：

其中，为级联信道虚波束空间矩阵，H_dv为直接信道的虚波束空间矩阵，定义系数矩阵/>从而简化模型；

使用vec算子及Kronecker积对该模型进行处理：

其中，表示Kronecker积，vec算子将直接信道H_dv和级联信道H_cv的矩阵转化为稀疏向量形式：

y_vec＝ψh

其中并且/>为待恢复的稀疏向量，用稀疏信号恢复算法直接求解：

作为优选的技术方案，从联合恢复的虚波束空间矩阵与毫米波波束空间模型完成直接信道与级联信道的最终估计，具体步骤包括：

从h中恢复出H_dv及H_cv，从系数矩阵中恢复A_dA，并带入波束空间形式：

最终完成两阶段的直接信道与级联信道的估计。

为了达到上述第二目的，本发明采用以下技术方案：

一种基于毫米波联合结构的智能反射面两级信道估计系统，包括：上行毫米波MIMO系统构建模块、波束空间模型构建模块、重构模块、虚波束空间矩阵定义模块、信道模型转换模块、最终估计模块；

所述上行毫米波MIMO系统构建模块用于构建智能反射面辅助的上行毫米波MIMO系统，采集不同时刻的接收信号：

所述波束空间模型构建模块用于构建毫米波信道具有角扩展的波束空间模型，表现出联合稀疏与低秩结构；

所述重构模块用于将时域接收信号模型重构为低秩矩阵采样形式，利用低秩矩阵近似算法完成接收信号的去噪处理：

所述虚波束空间矩阵定义模块用于定义级联信道的虚波束空间矩阵；

所述信道模型转换模块用于引入vec算子及Kronecker积，将信道模型转换为稀疏信号恢复问题：

所述最终估计模块用于从联合恢复的虚波束空间矩阵与毫米波波束空间模型完成直接信道与级联信道的最终估计。

为了达到上述第三目的，本发明采用以下技术方案：

一种计算机可读存储介质，存储有程序，所述程序被处理器执行时实现上述基于毫米波联合结构的智能反射面两级信道估计方法。

为了达到上述第四目的，本发明采用以下技术方案：

一种计算设备，包括处理器和用于存储处理器可执行程序的存储器，所述处理器执行存储器存储的程序时，实现如上述基于毫米波联合结构的智能反射面两级信道估计方法。

本发明与现有技术相比，具有如下优点和有益效果：

1)本发明构建了具有联合稀疏与低秩结构的毫米波角扩展信道模型，并引入IRS系统辅助两级信道估计。

2)本发明提出了改进的秩算子，得到低秩矩阵抽样模型后，基于低秩矩阵近似算法对接收信号进行去噪处理。

3)本发明方法利用vec算子和Kronecker积的性质，找到了级联信道的稀疏表示，将级联信道估计问题转化为稀疏信号恢复问题，可以直接使用现有的压缩感知方法。

4)本发明方法定义了的级联信道虚波束空间模型，在完成准确、有效信道估计的同时优化了IRS的反射系数矩阵。

附图说明

图1为本发明基于毫米波联合结构的智能反射面两级信道估计方法的流程图；

图2为本发明上行毫米波MIMO系统模型示意图；

图3为本发明的信道估计最小均方误差效果图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

实施例1

如图1所示，本发明公开了一种基于毫米波联合结构的智能反射面两级信道估计的方法，步骤如下：

S1、构建智能反射面辅助的上行毫米波MIMO系统，采集不同时刻的接收信号；

本实施例中，步骤S1过程如下：

如图2所示，在配备IRS的点对点上行mmWave MIMO系统中，基站端和用户端分别设置N_BS和N_US个天线，IRS有M个反射单元。系统中有US-BSUS-IRS和IRS-BS/>三条信道，在每条信道的发射端和接收端之间，为降低硬件成本，采用仅需一个射频链的模拟发射波束赋形和接收组合结构，且该结构内射频链数目在必要时可灵活扩展；

在时刻t，发射端通过波束赋形矩阵发射元素全为1的符号向量经过直接信道(H_d)与级联信道(GΘH_r)后与接收组合向量/>结合成为最终时域接收信号：

y(t)＝z^H(t)(H_d+GΘH_r)F(t)s(t)+ω(t) (1)

其中为IRS的反射系数矩阵，ω(t)为均值0，方差σ²的加性高斯噪声。

S2：推导毫米波信道具有角扩展的波束空间模型，证明其联合稀疏与低秩结构；

本实施例中，步骤S2过程如下：

研究表明，毫米波信道在AoA、AoD和高程域上采用角扩展的形式，而存在角度扩展时，信道可能表现出联合稀疏与低秩结构；

首先，从最简单情况着手，推导毫米波信道的角扩展模型。当AoA域上的角分布来自单个公共AoD射线时，信道模型为：

其中a_A和a_D是AoA和AoD对应的阵列响应向量，θ是平均AoA，φ是该公共AoD，I是AoD引起的AoA域路径数目，α_i是第i条路径的增益，ν_i是第i条路径相对于平均AOA的偏移量；

接着，扩展为AoA域上的角分布由两个间距很近的AoD引起，信道可表示为：

其中a_A和a_D是AoA与AoD对应的阵列响应向量，θ和φ是平均AoA与AoD，I和I′是AOA域的路径数目，α_i和α′分别是第i条路径的增益，ν_i是第i条路径相对于平均AOA的偏移量。由于两个AoD间距微小，可视为具有相似的功率角模式，即信道模型可简化为：

最后，对于毫米波信道，可视为AoA域上的角分布由J个间距很近的AoD引起，并且该毫米波系统内有L个此种小聚类。因此，毫米波的几何信道模型可表示为：

其中a_A和a_D是AoA与AoD对应的阵列响应向量，L为毫米波系统内聚类数，在每个聚类l中，θ和φ是平均AoA与AoD，I和I′是AOA域的路径数目，α_i和α′分别是第i条路径的增益，ν_i和分别是第i和j条路径相对于平均AoA和AoD的偏移量。具体地，a_A和a_D为：

其中λ是载波的波长，d是天线间距等于波长的一半；

在构建了毫米波具有角扩展的几何信道模型后，为了将信道估计表示为稀疏信号恢复问题，需要将几何信道模型进一步转换为更紧凑的波束空间形式：

其中和/>分别表示AoA域和AoD域上的路径增益的集合，和/>为过完备矩阵，其每列分别对应一个由预离散AoA/AoD参数化的导向向量。此外，由于角扩展只占整个角域的一小部分，所以α_l和β_l都是稀疏向量，只有少数非零项集中在与第l簇相关的平均AoA和AoD周围。因此虚波束空间信道H_v是L个稀疏矩阵的和。一方面，假设在每对{α_l,β_l}_l中，任意稀疏向量最多包含p个非零元素，则H_v最多有pL非零列和最多pL非零行。由于有限的散射性质和小的角扩散，通常有pL＜＜min{Ν₁,Ν₂}_l，所以H_v是稀疏的。另一方面，对于虚波束空间信道H_v，具有rank(H_v)＝L的低秩结构，所以H_v是低秩的。

S3：第一阶段：将时域接收信号模型重构为低秩矩阵采样形式，提出改进的秩算子，利用低秩矩阵近似算法完成接收信号的去噪处理；

本实施例中，步骤S3过程如下：

假设和/>为预先准备的波束赋形矢量和接收组合矢量的码本，分别有N_Z和N_F个元素，z(t)和f(t)从中随机不重复地选择。收集/>和/>中所有元素，分别用/>和表示。为了分析方便，忽略了噪声w的影响。因此，接收信号可以用低秩矩阵抽样形式重构为：

其中Y_ij是Y的第ij个元素，Ω表示观察到的集合；

定义一种改进的算子：P_δ(Γ)＝Γ(Γ^HΓ+δI)^-1Γ^H,δ≥0，其中为待优化的矩阵，tr()为矩阵的秩。当矩阵满秩时，δ＝0；当矩阵低秩时，具有性质：

即当系数δ无限趋近于0时，该改进的秩算子无限趋近于原矩阵的秩；

因此，可通过改进的秩算子完成低秩矩阵近似运算，完成接收信号的去噪处理：

其中ε₁为精确阈值，Y_Ω为观测到的待去噪接收信号，为去噪后的纯净接收信号。

S4：第二阶段：定义级联信道的虚波束空间矩阵，引入vec算子及Kronecker积，将信道模型转换为稀疏信号恢复问题；

本实施例中，步骤S4过程如下：

将毫米波的角扩展波束空间模型引入IRS系统：

定义级联信道的虚波束空间模型：

即合并信道H_d、信道G和IRS的反射系数矩阵Θ，级联信道的合并估计既可降低信道估计复杂度，也可达到在估计时优化反射系数矩阵Θ的目的；

联合毫米波的角扩展波束空间模型与级联信道的虚波束空间模型，可将低秩矩阵抽样模型表示为：

其中为级联信道虚波束空间矩阵，H_dv为直接信道的虚波束空间矩阵，定义系数矩阵/>从而简化模型；

使用vec算子及Kronecker积对该模型进行处理：

其中表示Kronecker积，vec算子将直接信道H_dv和级联信道H_cv的矩阵转化为稀疏向量形式：

y_vec＝ψh (18)

其中y_vec＝vec(Y)，并且/>为待恢复的稀疏向量，可用稀疏信号恢复算法直接求解：

S5：从联合恢复的虚波束空间矩阵与毫米波波束空间模型，完成直接信道与级联信道的最终估计；

本实施例中，步骤S5过程如下：

从h中恢复出H_dv及H_cv，从系数矩阵中恢复并带入波束空间形式：

至此，完成了两阶段的直接信道与级联信道的估计。

本实施例公开的基于毫米波联合结构的智能反射面两级信道估计方法的效果通过以下实验进一步说明：

如图3所示，从所显示的信道估计的均方误差可以看出，本实施例公开的基于毫米波联合结构的智能反射面两级信道估计方法得到的估计结果误差小、精度高。

综上所述，本实施例公开了一种基于毫米波联合结构的智能反射面两级信道估计方法，通过毫米波联合结构提出两级信道估计，并定义了的级联信道虚波束空间模型，在完成准确、有效信道估计的同时优化了IRS的反射系数矩阵。

实施例2

一种基于毫米波联合结构的智能反射面两级信道估计系统，包括：上行毫米波MIMO系统构建模块、波束空间模型构建模块、重构模块、虚波束空间矩阵定义模块、信道模型转换模块、最终估计模块；

在本实施例中，上行毫米波MIMO系统构建模块用于构建智能反射面辅助的上行毫米波MIMO系统，采集不同时刻的接收信号：

在本实施例中，波束空间模型构建模块用于构建毫米波信道具有角扩展的波束空间模型，表现出联合稀疏与低秩结构；

在本实施例中，重构模块用于将时域接收信号模型重构为低秩矩阵采样形式，利用低秩矩阵近似算法完成接收信号的去噪处理：

在本实施例中，虚波束空间矩阵定义模块用于定义级联信道的虚波束空间矩阵；

在本实施例中，信道模型转换模块用于引入vec算子及Kronecker积，将信道模型转换为稀疏信号恢复问题：

在本实施例中，最终估计模块用于从联合恢复的虚波束空间矩阵与毫米波波束空间模型完成直接信道与级联信道的最终估计。

实施例3

本实施例提供一种存储介质，存储介质可以是ROM、RAM、磁盘、光盘等储存介质，该存储介质存储有一个或多个程序，程序被处理器执行时，实现实施例1的基于毫米波联合结构的智能反射面两级信道估计方法。

实施例4

本实施例提供一种计算设备，该计算设备可以是台式电脑、笔记本电脑、智能手机、PDA手持终端、平板电脑或其他具有显示功能的终端设备，该计算设备包括处理器和存储器，存储器存储有一个或多个程序，处理器执行存储器存储的程序时，实现实施例1的基于毫米波联合结构的智能反射面两级信道估计方法。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (9)

1.一种基于毫米波联合结构的智能反射面两级信道估计方法，其特征在于，包括下述步骤：

构建智能反射面辅助的上行毫米波MIMO系统，采集不同时刻的接收信号：

构建毫米波信道具有角扩展的波束空间模型，表现出联合稀疏与低秩结构；

第一阶段：将时域接收信号模型重构为低秩矩阵采样形式，利用低秩矩阵近似算法完成接收信号的去噪处理：

第二阶段：定义级联信道的虚波束空间矩阵，引入vec算子及Kronecker积，将信道模型转换为稀疏信号恢复问题：

从联合恢复的虚波束空间矩阵与毫米波波束空间模型完成直接信道与级联信道的最终估计。

2.根据权利要求1所述的基于毫米波联合结构的智能反射面两级信道估计方法，其特征在于，所述构建智能反射面辅助的上行毫米波MIMO系统，具体步骤包括：

在配备IRS的点对点上行mmWave MIMO系统中，基站端和用户端分别设置N_BS和N_US个天线，IRS有M个反射单元，系统中有US-BS、US-IRS和IRS-BS三条信道，在每条信道的发射端和接收端之间，采用仅需一个射频链的模拟发射波束赋形和接收组合结构；

在时刻t，发射端通过波束赋形矩阵发射元素全为1的符号向量经过直接信道H_d与级联信道GΘH_r后与接收组合向量/>结合成为最终时域接收信号：

y(t)＝z^H(t)(H_d+GΘH_r)F(t)s(t)+ω(t)

其中，为IRS的反射系数矩阵，ω(t)为均值0，方差σ²的加性高斯噪声。

3.根据权利要求1所述的基于毫米波联合结构的智能反射面两级信道估计方法，其特征在于，所述构建毫米波信道具有角扩展的波束空间模型，具体步骤包括：

构建毫米波信道的角扩展模型：当AoA域上的角分布来自单个公共AoD射线时，信道模型为：

其中，a_A和a_D是AoA和AoD对应的阵列响应向量，θ是平均AoA，φ是该公共AoD，I是AoD引起的AoA域路径数目，α_i是第i条路径的增益，ν_i是第i条路径相对于平均AOA的偏移量；

信道表示为：

其中a_A和a_D是AoA与AoD对应的阵列响应向量，θ和φ是平均AoA与AoD，I和I′是AOA域的路径数目，α_i和α′分别是第i条路径的增益，ν_i是第i条路径相对于平均AOA的偏移量，两个AoD间距具有相似的功率角模式，即信道模型简化为：

毫米波的几何信道模型表示为：

其中a_A和a_D是AoA与AoD对应的阵列响应向量，L为毫米波系统内聚类数，在每个聚类l中，θ和φ是平均AoA与AoD，I和I′是AOA域的路径数目，α_i和α′分别是第i条路径的增益，ν_i和分别是第i和j条路径相对于平均AoA和AoD的偏移量；

a_A和a_D为：

其中λ是载波的波长，d是天线间距等于波长的一半；

在构建了毫米波具有角扩展的几何信道模型后，将几何信道模型转换为波束空间形式：

其中和/>分别表示AoA域和AoD域上的路径增益的集合，和/>为过完备矩阵，其每列分别对应一个由预离散AoA/AoD参数化的导向向量，虚波束空间信道H_v是L个稀疏矩阵的和。

4.根据权利要求1所述的基于毫米波联合结构的智能反射面两级信道估计方法，其特征在于，所述将时域接收信号模型重构为低秩矩阵采样形式，具体步骤包括：

假设和/>为预先准备的波束赋形矢量和接收组合矢量的码本，分别有N_Z和N_F个元素，z(t)和f(t)从中随机不重复地选择，收集/>和/>中所有元素，分别用/>和表示，接收信号用低秩矩阵抽样形式重构为：

其中Y_ij是Y的第ij个元素，Ω表示观察到的集合；

定义一种改进的算子：P_δ(Γ)＝Γ(Γ^HΓ+δI)^-1Γ^H,δ≥0，其中为待优化的矩阵，tr()为矩阵的秩，当矩阵满秩时，δ＝0；当矩阵低秩时，具有性质：

即当系数δ无限趋近于0时，该改进的秩算子无限趋近于原矩阵的秩；

通过改进的秩算子完成低秩矩阵近似运算，完成接收信号的去噪处理：

其中，ε₁为精确阈值，Y_Ω为观测到的待去噪接收信号，为去噪后的纯净接收信号。

5.根据权利要求1所述的基于毫米波联合结构的智能反射面两级信道估计方法，其特征在于，所述将信道模型转换为稀疏信号恢复问题，具体步骤包括：

将毫米波的角扩展波束空间模型引入IRS系统：

定义级联信道的虚波束空间模型：

即合并信道H_d、信道G和IRS的反射系数矩阵Θ；

联合毫米波的角扩展波束空间模型与级联信道的虚波束空间模型，将低秩矩阵抽样模型表示为：

其中，为级联信道虚波束空间矩阵，H_dv为直接信道的虚波束空间矩阵，定义系数矩阵/>从而简化模型；

使用vec算子及Kronecker积对该模型进行处理：

其中，表示Kronecker积，vec算子将直接信道H_dv和级联信道H_cv的矩阵转化为稀疏向量形式：

y_vec＝ψh

其中并且/>为待恢复的稀疏向量，用稀疏信号恢复算法直接求解：

6.根据权利要求1所述的基于毫米波联合结构的智能反射面两级信道估计方法，其特征在于，从联合恢复的虚波束空间矩阵与毫米波波束空间模型完成直接信道与级联信道的最终估计，具体步骤包括：

从h中恢复出H_dv及H_cv，从系数矩阵中恢复A_dA，A_GA，/>并带入波束空间形式：

最终完成两阶段的直接信道与级联信道的估计。

7.一种基于毫米波联合结构的智能反射面两级信道估计系统，其特征在于，包括：上行毫米波MIMO系统构建模块、波束空间模型构建模块、重构模块、虚波束空间矩阵定义模块、信道模型转换模块、最终估计模块；

所述上行毫米波MIMO系统构建模块用于构建智能反射面辅助的上行毫米波MIMO系统，采集不同时刻的接收信号：

所述波束空间模型构建模块用于构建毫米波信道具有角扩展的波束空间模型，表现出联合稀疏与低秩结构；

所述重构模块用于将时域接收信号模型重构为低秩矩阵采样形式，利用低秩矩阵近似算法完成接收信号的去噪处理：

所述虚波束空间矩阵定义模块用于定义级联信道的虚波束空间矩阵；

所述信道模型转换模块用于引入vec算子及Kronecker积，将信道模型转换为稀疏信号恢复问题：

所述最终估计模块用于从联合恢复的虚波束空间矩阵与毫米波波束空间模型完成直接信道与级联信道的最终估计。

8.一种计算机可读存储介质，存储有程序，其特征在于，所述程序被处理器执行时实现如权利要求1-6任一项所述基于毫米波联合结构的智能反射面两级信道估计方法。

9.一种计算设备，包括处理器和用于存储处理器可执行程序的存储器，其特征在于，所述处理器执行存储器存储的程序时，实现如权利要求1-6任一项所述基于毫米波联合结构的智能反射面两级信道估计方法。