CN110380997B - 基于自适应压缩感知的毫米波信道估计方法 - Google Patents
基于自适应压缩感知的毫米波信道估计方法 Download PDFInfo
- Publication number
- CN110380997B CN110380997B CN201910635118.4A CN201910635118A CN110380997B CN 110380997 B CN110380997 B CN 110380997B CN 201910635118 A CN201910635118 A CN 201910635118A CN 110380997 B CN110380997 B CN 110380997B
- Authority
- CN
- China
- Prior art keywords
- millimeter wave
- iteration
- channel estimation
- wave channel
- optimization problem
- 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Radio Transmission System (AREA)
Abstract
一种基于自适应压缩感知的毫米波信道估计方法,所述方法包括:构建毫米波信道估计的优化问题模型;采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解;基于所述优化问题模型的最优解,计算得到毫米波信道估计结果。上述的方案,可以在信道中有效路径数目未知的条件下对毫米波信道进行估计。
Description
技术领域
本发明属于通信技术领域,特别是涉及一种基于自适应压缩感知的毫米波信道估计方法。
背景技术
毫米波MIMO***被认为是下一代无线通信技术中前景最好的一种通信技术,其原因有以下两点。一、由于毫米波的波长较短,因此在毫米波MIMO***中可以将大量天线集成在一个天线阵列之中,以此降低毫米波所带来的路径损失,并提高***的波束赋形增益。二、由于毫米波所处的30-300GHz的频段中存在着富裕的频谱资源,因而可以满足日益增长的频谱需求。
信道估计是毫米波MIMO***进行有效通信和预编码的基础,由于***中基站端(接收端)的天线数目庞大,如果将传统的信道估计算法直接应用于毫米波MIMO***,将产生较高的复杂度。因而,为了降低信道估计的复杂度并充分利用毫米波信道的稀疏特性,压缩感知算法被广泛应用于毫米波***的信道估计中。
为了解决毫米波信道在频域中的块状频谱泄露问题,研究者使用基于块稀疏的正交匹配追踪算法(BOMP)对毫米波的信道信息进行估计。然而,在使用BOMP算法估计毫米波信道信息的过程中,***需要知晓信道中有效路径的数目,使得该算法在实际场景的应用中具有一定的局限性。
发明内容
本发明解决的技术问题是如何在信道中有效路径数目未知的条件下对毫米波信道进行估计。
为了达到上述目的,本发明实施例提供了一种基于自适应压缩感知的毫米波信道估计方法,所述方法包括:
构建毫米波信道估计的优化问题模型;
采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解;
基于所述优化问题模型的最优解,计算得到毫米波信道估计结果。
可选地,所述毫米波信道估计的优化问题模型为:
h=argmin||y-Ah||2;
其中,h表示为估计所得信道信息,y表示压缩感知中的观测向量,A表示压缩感知中的感知矩阵,h=vec(H)表示,H表示基站与用户端之间的信道信息。
可选地,所述采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解,包括:
设置毫米波信道中的有效路径数的初始值L=1,毫米波信道在频域中的稀疏块大小S=2,残差初始值r0等于观测向量y,测量矩阵Ρ=ΑΨ,块索引为R(2),初始的非零元素的索引集合迭代次数的初始值t=1,算法结束门限值T;
计算测量矩阵Ρ与第t-1次迭代得到的残差rt-1的乘积向量,对所述乘积向量中的每一个元素分别求其二范数,将求得的二范数中L个最大值对应的L个块索引组成对应的块索引集合,将所述块索引集合与第t-1次的非零元素的索引集合合并,得到第t次迭代的非零元素的索引集合;
基于第t次迭代的非零元素的索引集合,求解得到稀疏信道的估计值;
将从求得的稀疏信道的估计值中选取的L个最大值对应的索引,组成对应的第二索引集合,并在测量矩阵P中选取所述第二索引集合中的索引对应的L列,组成对应的L列矩阵;
基于所述稀疏信道的估计值和L列矩阵,求解得到第t次迭代的残差rt;
当确定第t次迭代的残差rt的二范数大于第t-1次迭代的残差rt-1的二范数之差大于所述结束门限值时,设置L=L+1,t=t+1,执行下一次迭代;
当确定第t次迭代的残差rt的二范数大于第t-1次迭代的残差rt-1的二范数之差小于所述结束门限值时,将求得的稀疏信道的估计值作为所述优化问题模型的最优解。
可选地,采用如下的公式基于所述优化问题模型的最优解,计算得到毫米波信道估计结果,包括:
与现有技术相比,本发明的有益效果为:
上述的方案,通过构建毫米波信道估计的优化问题模型,并采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解,再基于所述优化问题模型的最优解,计算得到毫米波信道估计结果,可以在有效信道路径数未知的情况下,实现对毫米波信道进行估计。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例的一种基于自适应压缩感知的毫米波信道估计方法的流程示意图;
图2是本发明实施例中采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解的流程示意图;
图3是本发明实施例中的基于自适应压缩感知的毫米波信道估计方法与现有的基于块稀疏的正交匹配追踪算法对***进行信道估计所得的波束赋形增益与信噪比的曲线比较图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。本发明实施例中有关方向性指示(诸如上、下、左、右、前、后等)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等,如果该特定姿态发生改变时,则该方向性指示也相应地随之改变。
如背景技术所述,现有技术中,在使用BOMP算法估计毫米波信道信息的过程中,***需要知晓信道中有效路径的数目,使得该算法在实际场景的应用中具有一定的局限性。
本发明的技术方案通过构建毫米波信道估计的优化问题模型,并采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解,再基于所述优化问题模型的最优解,计算得到毫米波信道估计结果,可以在有效信道路径数未知的情况下,实现对毫米波信道进行估计。
为使本发明的上述目的、特征和有益效果能够更为明显易懂,下面结合附图对本发明的具体实施例做详细的说明。
图1是本发明实施例的一种基于自适应压缩感知的毫米波信道估计方法的流程示意图。参见图1,一种基于自适应压缩感知的毫米波信道估计方法,具体可以包括如下的步骤:
步骤S101:构建毫米波信道估计的优化问题模型。
在具体实施中,本发明实施例的毫米波MIMO***,***包括全连接结构的基站发射端以及全连接结构的用户接收端。其中,所述的基站发射端配有Nb根天线,用户接收端配有Nu根天线和一根射频链。基站与用户之间采用的是单流通信,即基站在与用户的下行信道估计中只使用了一根射频链。
在第m个时隙内,用户端的接收信号可以表示为:
y(m)=w(m)HHf(m)s(m)+n(m) (1)
其中,y(m)表示第m个时隙内用户端的接收信号,表示第m个时隙内基站端所使用的预编码矩阵,s(m)表示第m个时隙内基站端所发送的导频序列,且s(m)=1,w(m)表示第m个时隙内用户端所使用的预编码矩阵。
基站与用户端之间的信道信息H可以表示为:
其中,αl为基站与用户间的第l个散射信道的信道增益,L为基站与用户间散射信道的数目,a(θkl)为基站端和用户端的导向矢量,k∈{b,u},且:
将上式(1)经过一次Kronecker积变换可得:
其中,h=vec(H)。
那么,在连续M个时隙内,用户端的接收信号可以表示为:
y=Ah+n (5)
且:
y=[y(1) y(2) ... y(M)]T (7)
信道信息H经过傅里叶变换可以表示为:
其中,Du和Db分别为Nu维和Nb维的傅里叶矩阵。
将上式(8)经过一次Kronecker变换之后可得:
h=Ψhω (9)
由此可得:
y=AΨhω+n (10)
综上,所述毫米波信道估计的优化问题模型为:
h=argmin||y-Ah||2s.t.h=Ψhw (11)
步骤S102:采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解。
参见图2,在本发明一实施例中,采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解时,可以具体包括如下的步骤:
步骤S201:参数初始化设置。设置毫米波信道中的有效路径数的初始值L=1,毫米波信道在频域中的稀疏块大小S=2,残差初始值r0等于观测向量y,测量矩阵Ρ=ΑΨ,块索引为R(2),初始的非零元素的索引集合迭代次数的初始值t=1,算法结束门限值T。
步骤S202:计算测量矩阵Ρ与第t-1次迭代得到的残差rt-1的乘积向量,对所述乘积向量中的每一个元素分别求其二范数,采用求得的二范数中L个最大值对应的L个块索引组成对应的块索引集合ΙL={R(2){li}:i=1,…L},将所述块索引集合ΙL与第t-1次的非零元素的索引集合合并,得到第t次迭代的非零元素的索引集合Λt,即Λt=Λt-1∪ΙL。
步骤S204:将从求得的第t次迭代对应的稀疏信道hω的估计值中选取的L个最大值对应的索引,组成对应的第二索引集合ΛtL,并在测量矩阵P中选取所述第二索引集合ΛtL中的索引对应的L列,组成对应的L列矩阵
步骤S206:判断第t次迭代的残差rt的二范数与第t-1次迭代的残差rt-1的二范数之差是否小于所述结束门限值;当判断结果为否时,可以执行步骤S207;反之,则可以执行步骤S208。
步骤S207:当确定第t次迭代的残差rt的二范数与第t-1次迭代的残差rt-1的二范数之差大于所述结束门限值,即||rt||2-||rt-1||2>T时,设置L=L+1,t=t+1,并从步骤S202开始执行第t+1次迭代。
步骤S208:当确定第t次迭代的残差rt的二范数与第t-1次迭代的残差rt-1的二范数之差小于所述结束门限值,即||rt||2-||rt-1||2<T时,将求得的第t次迭代对应的稀疏信道hω的估计值作为所述优化问题模型的最优解。
当采用上述的采用基于块稀疏的自适应正交匹配追踪算法求解得到所述优化问题模型的最优解时,接着执行步骤S103。
步骤S103:基于所述优化问题模型的最优解,计算得到毫米波信道估计结果。
下面结合图3的仿真实验对本发明方法的性能进行分析。
本发明的仿真实验采用均匀线性天线阵列,基站端配有64根天线,即Nb=64,用户接收端配有32根天线,即Nu=32,基站与用户间的有效路径数L=4,且用户接收端估计信道时,未知信道中的有效路径数目。
本仿真实验所研究的是波束赋形增益(beaming gain)与信噪比(SNR)的变化关系。其中,波束赋形增益定义为
γ=|woptHfopt|2 (13)
且:
H=UΣVH (14)
其中,γ为波束赋形增益,H为***估计所得的信道矩阵,U和V分别为H通过SVD分解后所得左奇异矩阵和右奇异矩阵,wopt和fopt为别表示左奇异矩阵U和右奇异矩阵V的第一列。
图3中的BOMP算法曲线为使用现有的基于块稀疏的正交匹配追踪算法对***进行信道估计所得的波束赋形增益与信噪比的曲线图,ABOMP算法为本发明的基于块稀疏的自适应正交匹配追踪算法,ABOMP算法曲线为使用本发明的基于块稀疏的自适应正交匹配追踪算法对***进行信道估计所得的波束赋形增益与信噪比的曲线图。
由图3可见,在相同信噪比的情况下,使用本发明的算法估计的信道信息所产生的波束赋形增益要大于使用传统的基于块稀疏的正交匹配追踪算法估计信道时所产生的波束赋形增益。
上述的方案,通过构建毫米波信道估计的优化问题模型,并采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解,再基于所述优化问题模型的最优解,计算得到毫米波信道估计结果,可以在有效信道路径数未知的情况下,实现对毫米波信道进行估计。
采用本发明实施例中的上述方案,通过构建毫米波信道估计的优化问题模型,并采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解,再基于所述优化问题模型的最优解,计算得到毫米波信道估计结果,可以在有效信道路径数未知的情况下,实现对毫米波信道进行估计。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,本发明要求保护范围由所附的权利要求书、说明书及其等效物界定。
Claims (3)
1.一种基于自适应压缩感知的毫米波信道估计方法,其特征在于,包括:
构建毫米波信道估计的优化问题模型;
采用基于块稀疏的自适应正交匹配追踪算法求解所述优化问题模型的最优解,具体包括:设置毫米波信道中的有效路径数的初始值L=1,毫米波信道在频域中的稀疏块大小S=2,残差初始值r0等于观测向量y,测量矩阵Ρ=ΑΨ,块索引为R(2),初始的非零元素的索引集合迭代次数的初始值t=1,算法结束门限值T;计算测量矩阵Ρ与第t-1次迭代得到的残差rt-1的乘积向量,对所述乘积向量中的每一个元素分别求其二范数,将求得的二范数中L个最大值对应的L个块索引组成对应的块索引集合,将所述块索引集合与第t-1次的非零元素的索引集合合并,得到第t次迭代的非零元素的索引集合;基于第t次迭代的非零元素的索引集合,求解得到稀疏信道的估计值;将从求得的稀疏信道的估计值中选取的L个最大值对应的索引,组成对应的第二索引集合,并在测量矩阵P中选取所述第二索引集合中的索引对应的L列,组成对应的L列矩阵;基于所述稀疏信道的估计值和L列矩阵,求解得到第t次迭代的残差rt;当确定第t次迭代的残差rt的二范数与第t-1次迭代的残差rt-1的二范数之差大于所述结束门限值时,设置L=L+1,t=t+1,执行下一次迭代;当确定第t次迭代的残差rt的二范数与第t-1次迭代的残差rt-1的二范数之差小于所述结束门限值时,将求得的稀疏信道的估计值作为所述优化问题模型的最优解;
基于所述优化问题模型的最优解,计算得到毫米波信道估计结果。
2.根据权利要求1所述的基于自适应压缩感知的毫米波信道估计方法,其特征在于,所述毫米波信道估计的优化问题模型为:
h=argmin||y-Ah||2
其中,h表示为估计所得信道信息,y表示压缩感知中的观测向量,A表示压缩感知中的感知矩阵,h=vec(H)表示,H表示基站与用户端之间的信道信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910635118.4A CN110380997B (zh) | 2019-07-15 | 2019-07-15 | 基于自适应压缩感知的毫米波信道估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910635118.4A CN110380997B (zh) | 2019-07-15 | 2019-07-15 | 基于自适应压缩感知的毫米波信道估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110380997A CN110380997A (zh) | 2019-10-25 |
CN110380997B true CN110380997B (zh) | 2022-03-22 |
Family
ID=68253263
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910635118.4A Active CN110380997B (zh) | 2019-07-15 | 2019-07-15 | 基于自适应压缩感知的毫米波信道估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110380997B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112769461B (zh) * | 2020-12-11 | 2022-04-22 | 华南理工大学 | 基于毫米波智能反射面通信的大规模天线信道估计方法 |
CN112929305B (zh) * | 2021-02-08 | 2021-11-16 | 武汉力通通信有限公司 | 基于iq失衡毫米波通信***的多径信道增益估计方法 |
CN113242042B (zh) * | 2021-04-08 | 2022-07-05 | 浙江大学 | 基于分段列相关矩阵的分组并行化的稀疏信道估计方法 |
CN113595944A (zh) * | 2021-07-30 | 2021-11-02 | 电子科技大学 | 一种用于毫米波mimo混合预编码***的信道估计方法 |
CN113630350B (zh) * | 2021-08-13 | 2022-08-26 | 北京邮电大学 | 一种毫米波多天线信道估计方法 |
CN114338301B (zh) * | 2021-12-13 | 2023-10-31 | 上海电机学院 | 一种基于压缩感知的ris辅助毫米波***的信道估计方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9614700B2 (en) * | 2013-01-02 | 2017-04-04 | Qualcomm Incorporated | Techniques for channel estimation in millimeter-wave communication systems |
CN106487725A (zh) * | 2016-10-20 | 2017-03-08 | 电子科技大学 | 一种多用户mimo***毫米波信道估计方法 |
CN107040296B (zh) * | 2017-02-28 | 2020-05-26 | 北京航空航天大学 | 毫米波通信中的信道估计方法 |
CN107018099B (zh) * | 2017-03-07 | 2020-06-09 | 上海交通大学 | 一种针对毫米波多用户mimo***的时变信道估计方法 |
CN107566305B (zh) * | 2017-08-15 | 2020-03-31 | 南京邮电大学 | 一种低复杂度的毫米波***信道估计方法 |
CN107508774B (zh) * | 2017-08-21 | 2020-11-03 | 安徽师范大学 | 联合信道表示和波束设计的毫米波mimo信道估计方法 |
CN108322409B (zh) * | 2018-01-25 | 2020-07-31 | 杭州电子科技大学 | 基于广义正交匹配追踪算法的稀疏ofdm信道估计方法 |
CN108599820B (zh) * | 2018-05-07 | 2021-10-15 | 东北大学 | 基于块结构自适应压缩采样匹配追踪算法的大规模mimo***信道估计方法 |
CN109039960A (zh) * | 2018-06-26 | 2018-12-18 | 东南大学 | 一种水声稀疏信道估计变步长稀疏度自适应匹配追踪方法 |
CN109462559B (zh) * | 2018-11-26 | 2020-12-29 | 东南大学 | 一种存在互耦时的稀疏毫米波信道估计方法 |
CN109617850A (zh) * | 2019-01-07 | 2019-04-12 | 南京邮电大学 | 基于自适应压缩感知的ofdm稀疏信道估计方法 |
-
2019
- 2019-07-15 CN CN201910635118.4A patent/CN110380997B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110380997A (zh) | 2019-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110380997B (zh) | 基于自适应压缩感知的毫米波信道估计方法 | |
EP3406059B1 (en) | Method for decoding symbol and receiver for receiving and decoding symbol | |
JP6685623B2 (ja) | 圧縮センシングに基づく効率的なスパースチャネル推定 | |
EP3348001B1 (en) | Systems and methods for beam selection for hybrid beamforming | |
US8649457B2 (en) | Precoding process for a transmitter of a MU-MIMO communication system | |
CN107566305B (zh) | 一种低复杂度的毫米波***信道估计方法 | |
CN107171985B (zh) | 一种信道估计方法及*** | |
CN110870216B (zh) | 一种波束成形方法及设备 | |
CN104779985A (zh) | 一种基于信道空间稀疏特性的迭代波束成形方法 | |
CN104779988A (zh) | 一种快速迭代波束成形的方法 | |
CN113315560A (zh) | 一种紧凑型平面阵列Massive MIMO***的波束赋形方法 | |
CN114520699B (zh) | 可重构智能表面辅助毫米波通信的低复杂度信道估计方法 | |
CN106130938B (zh) | Tdd大规模mimo***多用户联合信道估计方法 | |
KR101669857B1 (ko) | 대규모 다중-입력 다중-출력 통신 시스템에서 채널 추정 및 피드백 방법 | |
CN114567525A (zh) | 一种信道估计方法及装置 | |
CN114143896A (zh) | 一种大规模mimo跨频协作鲁棒传输方法 | |
CN107181705B (zh) | 一种半盲信道估计方法及*** | |
CN106878225B (zh) | 一种设备指纹与信道分离的方法及装置 | |
CN104994041A (zh) | 三维大规模天线***基于导频的波达波离角联合估计方法 | |
CN114221838B (zh) | 大规模mimo***中利用信道共轭数据的信道估计方法及*** | |
KR101447027B1 (ko) | 다중 안테나 시스템의 빔포밍 장치 및 방법 | |
US11962436B2 (en) | User detection technique, and method and apparatus for channel estimation in wireless communication system supporting massive multiple-input multiple-output | |
Ardah et al. | Low-complexity millimeter wave CSI estimation in MIMO-OFDM hybrid beamforming systems | |
CN114884776B (zh) | 信道估计方法、装置、电子设备及存储介质 | |
CN115278869B (zh) | 幂迭代辅助Root-MUSIC的快速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 |