用于双IRS辅助通信系统的信道估计方法及系统

用于双irs辅助通信系统的信道估计方法及系统
技术领域
1.本发明属于图像识别处理技术领域，本发明涉及一种用于双irs辅助通信系统的信道估计方法及系统。


背景技术：

2.智能反射面irs（intelligent reflecting surface）由大量准无源和低成本反射单元构成，是在6g无线通信系统中备受关注的备选技术。不同于传统的无线传输技术，irs可以通过适当调整反射单元的相移，重构无线传播环境。此外，由于具有重量轻且无射频链路的特点，irs可以密集部署在无线通信系统中，以较低能耗和硬件成本提高无线通信系统性能。相较于传统的单irs辅助通信系统，双irs辅助通信系统能实现更高的被动波束增益，更好地实现信道补盲。但是准无源irs并不具备信号处理能力，基站天线和irs反射单元的级联给信道估计带来巨大的压力。相较于单irs辅助通信系统而言，双irs辅助通信系统的信道估计技术面临更大的挑战。
3.申请号2021115230665公开了一种基于压缩感知的ris辅助毫米波系统的信道估计方法，基于毫米波的稀疏性将级联信道矩阵从空间域转换到虚拟角域，得到信道估计的稀疏表示；采用基于行-结构稀疏的改进压缩感知算法进行矩阵恢复，完成信道估计。但是该方法只能适用于单irs辅助通信系统，并不适用于双irs辅助通信系统的信道估计。


技术实现要素：

4.本发明的目的在于提供一种用于双irs辅助通信系统的信道估计方法及系统，基于k-means聚类算法的自适应迭代终止阈值选取算法，进一步提出了一种基于k-means的自适应双结构正交匹配追踪（kads-omp）算法，用于估计和信道；在双irs辅助信道估计步骤中，提出了一种基于压缩感知的三重压缩信道估计（tc-ce）算法，用于估计信道，有效降低信道估计的导频开销。
5.实现本发明目的的技术解决方案为：一种用于双irs辅助通信系统的信道估计方法，双irs辅助通信系统包括至少一个基站、irs1反射单元和irs2反射单元，用户通过irs反射单元与基站通信，包括以下步骤：s01：若irs1和irs2中有一个关闭，进行单irs辅助信道估计；s02：若irs1和irs2同时打开，估计信道；所述信道由如下公式计算得到：
6.其中，为基站的字典酉矩阵，为信道的行支撑集，与信道的基站到达角支撑集相同，为第个基站到达角对应的信道，为
用户编号，为的元素个数，和分别为irs1和irs2反射单元的反射系数，，为导频开销，为用户到irs1的信道，为用户到irs2的信道，为irs1到基站的信道，为irs2到基站的信道，为irs1到irs2的信道；其中，由如下公式得到：其中，和为irs1和irs2的字典酉矩阵，为角域稀疏信道。
7.优选的技术方案中，所述s02中估计信道包括：s11：同时打开irs1和irs2，保持irs2反射单元的反射系数不变，在每个时隙改变irs1反射单元的反射系数，估计irs1的角度信息，得到每个用户的压缩矩阵；s12：保持irs1反射单元的反射系数不变，在每个时隙改变irs2反射单元的反射系数，估计irs2的角度信息，得到每个用户的压缩矩阵；s13：同时改变irs1和irs2反射单元的反射系数，估计所有用户每个基站到达角对应的信道；s14：得到信道的估计值。
8.优选的技术方案中，所述s11中得到每个用户的压缩矩阵的方法包括：s111：对所有用户的接收信号预处理得到；；
9.其中，和为噪声，右上角的是共轭转置；s112：对所有用户的预处理信号经omp算法得到所有用户irs1的角度支撑集；s113：由如下公式得到的压缩矩阵：
10.其中，为用户数。
11.优选的技术方案中，所述s13包括：s131：利用压缩矩阵、对公式(11)简化，由如下公式得到：
s132：经omp算法得到的估计值，进一步由公式(7)和公式(8)得到和的估计值和，是向量化。
12.优选的技术方案中，所述s01中若关闭irs2，打开irs1，通过kasd-omp算法估计信道，若打开irs2，关闭irs1通过kasd-omp算法估计信道，信道估计包括：s21：对所有用户的接收信号预处理得到；s22：经kasd-omp算法得到的估计值；s23：根据得到的，由公式(1)得到信道的估计值；
13.其中，为基站的字典酉矩阵，为角域稀疏信道，为用户编号，为irs1的字典酉矩阵；同理得到信道的估计值。
14.优选的技术方案中，所述s21中接收信号及预处理信号由如下公式得到：
[0015][0016]
其中，为感知矩阵，，，为irs1反射单元的反射系数，，为导频开销，和为噪声。
[0017]
优选的技术方案中，所述s22中kasd-omp算法包括如下步骤：s221：由行支撑估计算法得到所有用户的行支撑集和行稀疏度；s222：由列支撑估计算法得到所有用户的列支撑集；s223：根据ls算法得到所有用户的估计值。
[0018]
优选的技术方案中，所述s221中行支撑估计算法包括如下步骤：s2211：对所有的每一行做如下处理：其中，初始化为零向量；s2212：通过k-means算法对进行二分类，得到聚类中心和，通过如下式
orthogonal matching pursuit，kads-omp）算法，用于估计和信道；在双irs辅助信道估计步骤中，考虑信道的结构特性，本发明提出了一种基于压缩感知的三重压缩信道估计（triple compression channel estimation，tc-ce）算法，用于估计信道。与传统方案相比，本方案提出的kads-omp和tc-ce算法利用了毫米波信道的稀疏性和结构特性，有效降低信道估计的导频开销。
[0024]
实施例1：如图1所示，一种用于双irs辅助通信系统的信道估计方法，双irs辅助通信系统包括至少一个基站、irs1反射单元和irs2反射单元，用户通过irs反射单元与基站通信，包括以下步骤：s01：若irs1和irs2中有一个关闭，进行单irs辅助信道估计；s02：若irs1和irs2同时打开，估计信道；所述信道由如下公式计算得到：其中，为基站的字典酉矩阵，为信道的行支撑集，与信道的基站到达角支撑集相同，为第个基站到达角对应的信道，为用户编号，为的元素个数，和分别为irs1和irs2反射单元的反射系数，，为导频开销，为用户到irs1的信道，为用户到irs2的信道，为irs1到基站的信道，为irs2到基站的信道，为irs1到irs2的信道；其中，由如下公式得到：
[0025]
其中，和为irs1和irs2的字典酉矩阵，为角域稀疏信道。
[0026]
一较佳的实施例中，步骤s02中估计信道包括：s11：同时打开irs1和irs2，保持irs2反射单元的反射系数不变，在每个时隙改变irs1反射单元的反射系数，估计irs1的角度信息，得到每个用户的压缩矩阵；s12：保持irs1反射单元的反射系数不变，在每个时隙改变irs2反射单元的反射系数，估计irs2的角度信息，得到每个用户的压缩矩阵；s13：同时改变irs1和irs2反射单元的反射系数，估计所有用户每个基站到达角对应的信道；s14：得到信道的估计值。
[0027]
一较佳的实施例中，步骤s11中得到每个用户的压缩矩阵的方法包括：s111：对所有用户的接收信号预处理得到；
[0028][0029]
其中，和为噪声，右上角的是共轭转置；s112：对所有用户的预处理信号经omp算法得到所有用户irs1的角度支撑集；s113：由如下公式得到的压缩矩阵：
[0030]
其中，为用户数。
[0031]
一较佳的实施例中，步骤s13包括：s131：利用压缩矩阵、对公式(11)简化，由如下公式得到：
[0032]
s132：经omp算法得到的估计值，进一步由公式(7)和公式(8)得到和的估计值和，是向量化。
[0033]
一较佳的实施例中，步骤s01中若关闭irs2，打开irs1，通过kasd-omp算法估计信道，若打开irs2，关闭irs1通过kasd-omp算法估计信道，信道估计包括：s21：对所有用户的接收信号预处理得到；s22：经kasd-omp算法得到的估计值；s23：根据得到的，由公式(1)得到信道的估计值；
[0034]
其中，为基站的字典酉矩阵，为角域稀疏信道，为用户编号，为irs1的字典酉矩阵；同理得到信道的估计值。
[0035]
一较佳的实施例中，步骤s21中接收信号及预处理信号由如下公式得到：
[0036][0037]
其中，为感知矩阵，，，为irs1反射单元的反射系数，，为导频开销，和为噪声。
[0038]
一较佳的实施例中，步骤s22中kasd-omp算法包括如下步骤：s221：由行支撑估计算法得到所有用户的行支撑集和行稀疏度；s222：由列支撑估计算法得到所有用户的列支撑集；s223：根据ls算法得到所有用户的估计值。
[0039]
一较佳的实施例中，步骤s221中行支撑估计算法包括如下步骤：s2211：对所有的每一行做如下处理：
[0040]
其中，初始化为零向量；s2212：通过k-means算法对进行二分类，得到聚类中心和，通过如下式处理得到筛选阈值：
[0041]
其中，和为加权比重；s2213：将中大于的值对应的索引号作为行支撑集的元素，行稀疏度为元素的个数。
[0042]
一较佳的实施例中，步骤s222中所述列支撑估计算法包括如下步骤：s2221：对所有用户及所有中行支撑对应的稀疏行经omp算法得到对应的列支撑集；s2222：将在中累计出现过次的元素作为公共列支撑集的元素，其中，为向上取整运算符，为用户数，为松弛变量，；s2223：将作为omp算法迭代的初始值，重复步骤s2221，得到所有用户所有稀疏行的列支撑集。
[0043]
另一实施例中，一种用于双irs辅助通信系统的信道估计系统，包括一控制器，所
述控制器内存储有上述的用于双irs辅助通信系统的信道估计方法。
[0044]
具体的，下面以一较佳的实施例为例对用于双irs辅助通信系统的信道估计系统的工作流程说明如下：如图2所示，假定基站和用户之间的视线传播路径被阻断，基站和用户之间仅存在、和三条路径，为用户到irs1的信道，为用户到irs2的信道，为irs1到基站的信道，为irs2到基站的信道，为irs1到irs2的信道。本发明通过开关irs和改变irs反射单元的反射系数，分别估计上述三条路径，包括以下关键步骤：步骤1：关闭irs2，打开irs1，通过kasd-omp算法估计信道；步骤2：打开irs2，关闭irs1，通过kasd-omp算法估计信道；步骤3：同时打开irs1和irs2，通过tc-ce算法估计信道。
[0045]
较佳的，信道可由如下公式得到：
[0046]
其中，为基站的字典酉矩阵，为角域稀疏信道，为用户编号，为irs1的字典酉矩阵。
[0047]
步骤1进一步包括：步骤1.1：对所有用户的接收信号预处理得到；步骤1.2：经kasd-omp算法得到的估计值；步骤1.3：根据步骤1.2得到的，由公式(1)得到信道的估计值。
[0048]
接收信号及预处理信号可由如下公式得到：
[0049][0050]
其中，为感知矩阵，，，为irs1反射单元的反射系数，，为导频开销，和为噪声。
[0051]
kasd-omp算法包括如下步骤：步骤1.2.1：由行支撑估计算法得到所有用户的行支撑集和行稀疏度；步骤1.2.2：由列支撑估计算法得到所有用户的列支撑集；步骤1.2.3：根据ls（least square）算法得到所有用户的估计值。
[0052]
行支撑估计算法包括如下步骤：步骤1.2.1.1：对所有的每一行做如下处理：
其中，初始化为零向量；步骤1.2.1.2：通过k-means算法对进行二分类，得到聚类中心和，通过如下式处理得到筛选阈值：
[0053]
其中，和为加权比重；步骤1.2.1.3：将中大于的值对应的索引号作为行支撑集的元素, 行稀疏度为元素的个数。
[0054]
列支撑估计算法包括如下步骤：步骤1.2.2.1：对所有用户及所有中行支撑对应的稀疏行经omp算法得到对应的列支撑集；步骤1.2.2.2：将在中累计出现过次的元素作为公共列支撑集的元素，其中，为向上取整运算符，为用户数，为松弛变量，；步骤1.2.2.3：将作为omp算法迭代的初始值，重复步骤1.2.2.1，得到所有用户所有稀疏行的列支撑集。
[0055]
ls算法如下式所示：
[0056]
其中，为伪逆运算符。
[0057]
同理可知，步骤2进一步包括的步骤和步骤1类似。这里不再赘述。
[0058]
较佳的，信道可由如下公式得到：
[0059]
其中，为信道的行支撑集，与信道的基站到达角支撑集相同，为第个基站到达角对应的信道，为的元素个数，和分别为irs1和irs2反射单元的反射系数，, 为导频开销。可由如下公式得到：
[0060]
其中，和为irs1和irs2的字典酉矩阵，为角域稀疏信道。
[0061]
步骤3进一步包括：步骤3.1：同时打开irs1和irs2，保持irs2反射单元的反射系数不变，在每个时隙改变irs1反射单元的反射系数，估计irs1的角度信息，得到每个用户的压缩矩阵；步骤3.2：保持irs1反射单元的反射系数不变，在每个时隙改变irs2反射单元的反射系数，估计irs2的角度信息，得到每个用户的压缩矩阵；步骤3.3：同时改变irs1和irs2反射单元的反射系数，估计所有用户每个基站到达角对应的信道；步骤3.4：由公式(7)和公式(8)得到信道的估计值。
[0062]
步骤3.1进一步包括：步骤3.1.1：对所有用户的接收信号预处理得到；步骤3.1.2：对所有用户的预处理信号经omp算法得到所有用户irs1的角度支撑集；步骤3.1.3：由如下公式得到的压缩矩阵：
[0063]
较佳的，接收信号及预处理信号可由如下公式得到：可由如下公式得到：其中，和为噪声。
[0064]
同理可知，步骤3.2进一步包括的步骤和步骤3.1类似。这里不再赘述。
[0065]
步骤3.3进一步包括：步骤3.3.1：利用压缩矩阵、对公式(11)简化，可由如下公式得到：步骤3.3.2：经omp算法可得到的估计值，进一步由公式(7)和公式(8)得到和的估计值和。
[0066]
本发明的有效性可通过以下仿真实验进一步说明：在以下仿真中，基站天线为64根天线组成的均匀线性阵列，irs反射单元为16
×
16的均匀平面阵列，用户天线为1根天线。用户和irs1之间的距离为10米，用户和irs2之间的距离为100米，基站和irs1之间的距离为100米，基站和irs2之间的距离为10米，irs1和irs2之间的距离为100米。
[0067]
本发明在估计和信道上的性能如图3、图4所示。图3为导频开销为32时，本发明方案的kasd-omp算法在不同信噪比下与传统压缩感知算法和均值阈值正交匹配追踪算法的性能仿真图。图4为信噪比为0分贝时，本发明方案的kasd-omp算法在不同导频开销下与传统压缩感知算法和均值阈值正交匹配追踪算法的性能仿真图。由图可知，本发明方案在不同信噪比和导频开销下均优于两种对比方案。
[0068]
本发明在估计信道上的性能如图5、图6所示。图5为导频开销为64时，本发明方案的tc-ce算法在不同信噪比下的性能仿真图。图6为信噪比为0分贝时，本发明方案的tc-ce算法在不同导频开销下的性能仿真图。
[0069]
上述实施例为本发明优选地实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

技术特征：
1.一种用于双irs辅助通信系统的信道估计方法，双irs辅助通信系统包括至少一个基站、irs1反射单元和irs2反射单元，用户通过irs反射单元与基站通信，其特征在于，包括以下步骤：s01：若irs1和irs2中有一个关闭，进行单irs辅助信道估计；s02：若irs1和irs2同时打开，估计信道；所述信道由如下公式计算得到：，其中，为基站的字典酉矩阵，为信道的行支撑集，与信道的基站到达角支撑集相同，为第个基站到达角对应的信道，为用户编号，为的元素个数，和分别为irs1和irs2反射单元的反射系数，，为导频开销，为用户到irs1的信道，为用户到irs2的信道，为irs1到基站的信道，为irs2到基站的信道，为irs1到irs2的信道；其中，由如下公式得到：，其中，和为irs1和irs2的字典酉矩阵，为角域稀疏信道。2.根据权利要求1所述的用于双irs辅助通信系统的信道估计方法，其特征在于，所述s02中估计信道包括：s11：同时打开irs1和irs2，保持irs2反射单元的反射系数不变，在每个时隙改变irs1反射单元的反射系数，估计irs1的角度信息，得到每个用户的压缩矩阵；s12：保持irs1反射单元的反射系数不变，在每个时隙改变irs2反射单元的反射系数，估计irs2的角度信息，得到每个用户的压缩矩阵；s13：同时改变irs1和irs2反射单元的反射系数，估计所有用户每个基站到达角对应的信道；s14：得到信道的估计值。3.根据权利要求2所述的用于双irs辅助通信系统的信道估计方法，其特征在于，所述s11中得到每个用户的压缩矩阵的方法包括：s111：对所有用户的接收信号预处理得到；，，其中，和为噪声，右上角的是共轭转置；s112：对所有用户的预处理信号经omp算法得到所有用户irs1的角度支撑集；s113：由如下公式得到的压缩矩阵：
，其中，为用户数。4.根据权利要求2所述的用于双irs辅助通信系统的信道估计方法，其特征在于，所述s13包括：s131：利用压缩矩阵、对公式(11)简化，由如下公式得到：，s132：经omp算法得到的估计值，进一步由公式(7)和公式(8)得到和的估计值和，是向量化。5.根据权利要求1所述的用于双irs辅助通信系统的信道估计方法，其特征在于，所述s01中若关闭irs2，打开irs1，通过kasd-omp算法估计信道，若打开irs2，关闭irs1通过kasd-omp算法估计信道，所述信道估计包括：s21：对所有用户的接收信号预处理得到；s22：经kasd-omp算法得到的估计值；s23：根据得到的，由公式(1)得到信道的估计值；，其中，为基站的字典酉矩阵，为角域稀疏信道，为用户编号，为irs1的字典酉矩阵；同理得到信道的估计值。6.根据权利要求5所述的用于双irs辅助通信系统的信道估计方法，其特征在于，所述s21中接收信号及预处理信号由如如下公式得到：，，其中，为感知矩阵，，，为irs1反射单元的反射系数，，为导频开销，和为噪声。7.根据权利要求5所述的用于双irs辅助通信系统的信道估计方法，其特征在于，所述s22中kasd-omp算法包括如下步骤：s221：由行支撑估计算法得到所有用户的行支撑集和行稀疏度；s222：由列支撑估计算法得到所有用户的列支撑集；s223：根据ls算法得到所有用户的估计值。8.根据权利要求7所述的用于双irs辅助通信系统的信道估计方法，其特征在于，所述s221中行支撑估计算法包括如下步骤：
s2211：对所有的每一行做如下处理：，其中，初始化为零向量；s2212：通过k-means算法对进行二分类，得到聚类中心和，通过如下式处理得到筛选阈值：，其中，和为加权比重；s2213：将中大于的值对应的索引号作为行支撑集的元素，行稀疏度为元素的个数。9.根据权利要求7所述的用于双irs辅助通信系统的信道估计方法，其特征在于，所述s222中所述列支撑估计算法包括如下步骤：s2221：对所有用户及所有中行支撑对应的稀疏行经omp算法得到对应的列支撑集；s2222：将在中累计出现过次的元素作为公共列支撑集的元素，其中，为向上取整运算符，为用户数，为松弛变量，；s2223：将作为omp算法迭代的初始值，重复步骤s2221，得到所有用户所有稀疏行的列支撑集。10.一种用于双irs辅助通信系统的信道估计系统，其特征在于，包括一控制器，所述控制器内存储有权利要求1-9任一项所述的用于双irs辅助通信系统的信道估计方法。

技术总结
本发明公开了一种用于双IRS辅助通信系统的信道估计方法及系统，具体步骤可分为单IRS辅助信道估计和双IRS辅助信道估计两个步骤。在单IRS辅助信道估计步骤中，提出了一种基于k-means聚类算法的自适应迭代终止阈值选取算法，进一步提出了一种基于k-means的自适应双结构正交匹配追踪（KADS-OMP）算法，用于估计h1—G1和h2—G2信道；在双IRS辅助信道估计步骤中，提出了一种基于压缩感知的三重压缩信道估计（TC-CE）算法，用于估计h1—D—G2信道。与传统方案相比，本方案提出的KADS-OMP和TC-CE算法利用了毫米波信道的稀疏性和结构特性，有效降低信道估计的导频开销。效降低信道估计的导频开销。效降低信道估计的导频开销。


技术研发人员：梁彦 涂威洋 李飞
受保护的技术使用者：南京邮电大学
技术研发日：2023.06.29
技术公布日：2023/8/1