CN117014258A - 一种毫米波mimo频率选择性信道估计方法 - Google Patents
一种毫米波mimo频率选择性信道估计方法 Download PDFInfo
- Publication number
- CN117014258A CN117014258A CN202311069563.1A CN202311069563A CN117014258A CN 117014258 A CN117014258 A CN 117014258A CN 202311069563 A CN202311069563 A CN 202311069563A CN 117014258 A CN117014258 A CN 117014258A
- Authority
- CN
- China
- Prior art keywords
- millimeter wave
- channel
- matrix
- frequency selective
- channel estimation
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000011159 matrix material Substances 0.000 claims abstract description 62
- 230000005540 biological transmission Effects 0.000 claims description 17
- 238000005562 fading Methods 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 6
- 238000007493 shaping process Methods 0.000 claims description 5
- 238000013139 quantization Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000003491 array Methods 0.000 claims description 3
- 230000003111 delayed effect Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 6
- 239000000969 carrier Substances 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012804 iterative process Methods 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- 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/0204—Channel estimation of multiple channels
-
- 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/0224—Channel estimation using sounding signals
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Transmission System (AREA)
Abstract
本发明公开了一种毫米波MIMO频率选择性信道估计方法,包括构建毫米波MIMO频率选择性信道模型,表示为紧凑的矩阵形式;基于毫米波MIMO频率选择性信道的稀疏结构,将信道估计问题表征为稀疏重构问题;利用压缩感知的思想求取重构问题的最优解,并计算出信道估计值。本发明有效的提高了稀疏信道索引估计的准确度,从而提升了信道估计的性能,使其性能相比传统的基于压缩感知的毫米波信道估计技术得到大幅提升。
Description
技术领域
本发明涉及无线通信技术领域,尤其涉及一种基于载波信道的联合稀疏结构和统计特性的毫米波MIMO频率选择性信道估计方法。
背景技术
无线通信在高数据速率和大连接方面的需求,推动了毫米波多输入多输出(Multiple Input Multiple Output,MIMO)技术的快速发展,并使其成为未来无线通信系统的关键技术之一。信道状态信息(Channel State Information,CSI)的获取对于毫米波MIMO系统无线传输性能有着重要的影响。基站端(Base Station,BS)需要准确的获取下行信道状态信息,以便于其使用波束赋形、预编码等技术对发射信号进行处理实现高效下行传输。与传统的超高频/微波频带相比,毫米波信道的复杂多变、较高传播路径损耗、多径和大多普勒给信道状态信息(Channel State Information,CSI)获取带来了困难,如何高效地估计毫米波MIMO信道,提升上、下行链路性能,成为系统设计的难点。目前毫米波MIMO信道估计方面的研究主要是利用毫米波信道的结构特征从某个维度对毫米波信道进行变换处理,然后根据变换后的信道设计算法完成信道参数的估计和反馈,然而,在某些场景中毫米波信道具有多个维度,且在不同的维度具有不同的结构特征,如果仅将信道矩阵在单个维度进行降维后再进行变换处理,其他维度的结构信息可能会被忽略或被破坏。
发明内容
发明目的:本发明的目的是提供一种基于载波信道的联合稀疏结构和统计特性的毫米波MIMO频率选择性信道估计方法。
技术方案:本发明包括如下步骤:
(1)构建毫米波MIMO频率选择性信道模型,表示为紧凑的矩阵形式;
(2)基于毫米波MIMO频率选择性信道的稀疏结构,将信道估计问题表征为稀疏重构问题;
(3)利用压缩感知的思想求取重构问题的最优解,并计算出信道估计值。
进一步地,所述步骤(1)包括:
(1.1)构建毫米波MIMO下行传输系统的接收信号模型:
(1.2)构建毫米波MIMO频率选择性信道模型。
进一步地,所述步骤(1.1)假设发射端配置的天线数目为Nt,接收端配置的天线数目为Nr,信道为频率选择性信道,系统使用OFDM技术将传输频带划分为Nc个正交子载波进行数据传输,且每个子载波带宽小于信道的相干带宽,每个子载波的信道为平坦衰落信道,发送端对每个子载波上的传输信号进行波束赋形,第k个子载波上发送的第n个信号xk,n可以表示为:
xk,n=fk,nsk,n
其中,sk,n表示第k个子载波上传输的第n个导频符号,fk,n表示第k个子载波上第n个导频符号使用的波束赋形向量。
进一步地,所述步骤(1.1)假设所有子载波使用相同的波束赋形矩阵和合并矩阵,即Fk=F,Wk=W,k=1,...,Nc,每个子载波使用的相同的具有单位功率的导频符号,第k个子载波上的接收信号可以简化为:
Yk=WHHkF+WHNk
其中,F表示波束赋形矩阵。
进一步地,所述步骤(1.2)根据毫米波信道的有限散射特性,延迟为g个采样周期的信道的时域形式可表示为:
其中,L表示散射体的数目。考虑每个散射体只有单条散射路径的情况,αl表示第l个散射体的散射路径的衰落系数;r(τ)表示通信系统中包含的脉冲成形滤波器、匹配滤波器等在τ时刻对信号带来的综合影响,Ts表示采样周期,τl表示第l条路径的延迟;θl和φl分别表示第l条信号传播路径的出发角和到达角,ar(φl)和at(θl)分别表示第l条路径的传播信号在接收端和发射端的阵列流形矢量,
通过对时域信道H(g)进行Nc点的DFT,第k个子载波上的频域信道矩阵可以表示为更加紧凑的矩阵形式
其中,L表示路径数。
进一步地,所述发射端和接收端采用均匀线阵,ar(φl)和at(θl)形式可表示为:
其中,λ为信号的波长,d表示相邻天线的间距。
进一步地,所述步骤(2)假设毫米波信道中传播路径出发角和到达角位于和/>构成的量化角度栅格内,此时,Hk可以近似表示为下列形式:
其中,是一个M×N维的稀疏矩阵,该矩阵具有L个非零元素值,其非零值为第k个子载波在第l条路径的频域衰落系数βl,k,每个非零值的索引为该路径的到达角和方向角在其量化角度区间的位置,/>通常被称为Hk的虚拟信道,信道估计问题可以转化为稀疏矩阵的估计问题,Nc个子载波的信道的联合估计可以表示为如下形式:
supp(Hv(:,1))=,...,=supp(Hv(:,Nc))
其中,Hv(:,k)表示Hv的第k列。
进一步地,所述步骤(3)包括令测量矩阵构建索引集合/>和/>用于存储算法迭代过程中得到的非零索引的估计值,在第i次迭代的过程中,算法首先将前一次迭代中得到的所有子载波残差信号与测量矩阵进行相关操作,得到相关矩阵作为稀疏信号的参考,提取相关矩阵中每个子载波对应的相关向量中绝对值最大的元素的索引并将其存入集合/>中。
进一步地,所述步骤(3)计算集合中不同元素出现的频率,记为/>令
将中出现频率最大元素amax存入非零元素索引集合/>假设存在N个索引元素满足下列不等式
分别计算残差信道相关矩阵中第行的元素绝对值的累加和,记为/>
令
将amax存入非零元素集合假设本次迭代索引集合/>中的元素为/>根据/>中的元素和测量矩阵Φ构建新的测量矩阵/>
基于当前索引集合的信道估计值为/>
进一步地,所述步骤(3)使用和/>计算本次迭代的信号残差R供下次迭代使用
当迭代终止条件满足时,算法结束,得到和集合/>否则,重复上述过程直至迭代收敛条件满足为止,
算法结束后,矩阵包含所有子载波上的路径衰落系数的估计值,其中,/>为第k个子载波的虚拟信道/>的向量化形式/>中非零元素的估计值,/>中的索引为/>中的元素在稀疏向量/>中对应的位置索引估计值,根据/>和/>重构/>的估计值对/>中的元素进行矩阵重排得到/>的估计值,即
根据第k个子载波的信道估计值为/>
有益效果:本发明与现有技术相比,具有如下显著优点:本发明提出的算法在迭代过程中利用了子载波信道之间的联合稀疏结构,使用所有子载波的样本数据联合估计稀疏信道的非零元素的索引,有效的提高了稀疏信道索引估计的准确度,从而提升了信道估计的性能,使其性能相比传统的基于压缩感知的毫米波信道估计技术得到大幅提升。
附图说明
图1为本发明的流程图;
图2为基于稀疏结构求取重构问题最优解的流程图;
图3为本发明与传统算法NMSE性能比较图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1所示,本发明包括如下步骤:
(1)构建毫米波MIMO频率选择性信道模型,表示为紧凑的矩阵形式。
(1.1)构建毫米波MIMO下行传输系统的接收信号模型;
发射端配置的天线数目为Nt,接收端配置的天线数目为Nr,信道为频率选择性信道,为了降低信道的频率选择性对信号传输的影响,系统使用OFDM技术将传输频带划分为Nc个正交子载波进行数据传输,且每个子载波带宽小于信道的相干带宽,此时每个子载波的信道为平坦衰落信道。为了提升下行信号传输的性能,发送端对每个子载波上的传输信号进行波束赋形,第k个子载波上发送的第n个信号xk,n可以表示为:
xk,n=fk,nsk,n (1)
其中,sk,n表示第k个子载波上传输的第n个导频符号,表示第k个子载波上第n个导频符号使用的波束赋形向量。假设系统使用Np个导频符号进行下行信道估计,将所有导频符号联合进行考虑,令
则第k个子载波上发送的信号可以统一表示为:
Xk=FkSk (3)
其中,是第k个子载波上使用的导频符号矩阵,其对角元素为/>根据OFDM的处理流程,系统使用IDFT将每个子载波上的传输信号从频域信号变换为时域信号,为了降低符号间干扰,发射端为时域信号添加循环前缀(Cyclic Prefix,CP),然后将其发送给接收端,接收端通过去CP和DFT处理后将时域信号变换为频域信号。此时,接收端在第k个子载波上的接收信号可以表示为:
Rk=HkFkSk+Nk (4)
其中,为第k个子载波上发送端和接收端的频域信道矩阵,/>是传输过程中的加性噪声,这里,假设Nk中的元素服从均值为0、方差为1的复高斯分布。接收端分别对每个子载波的接收信号进行合并处理,则合并后的接收信号可以表示为:
其中,为接收端对第k个子载波的接收信号进行合并时使用的合并矩阵。由于本章不考虑波束赋形矩阵和合并矩阵的设计问题,因此,为了不失一般性,假设所有子载波使用相同的波束赋形矩阵和合并矩阵,即Fk=F,Wk=W,k=1,...,Nc,每个子载波使用的相同的具有单位功率的导频符号,则第k个子载波上的接收信号可以简化为
Yk=WHHkF+WHNk (6)
(1.2)构建毫米波MIMO频率选择性信道模型:
根据毫米波信道的有限散射特性,延迟为g个采样周期的信道的时域形式可表示为
其中,L表示散射体的数目。考虑每个散射体只有单条散射路径的情况,αl表示第l个散射体的散射路径的衰落系数,假设服从均值为0,方差为1的高斯分布。r(τ)表示通信系统中包含的脉冲成形滤波器、匹配滤波器等在τ时刻对信号带来的综合影响,Ts表示采样周期,τl表示第l条路径的延迟。θl和φl分别表示第l条信号传播路径的出发角和到达角,ar(φl)和at(θl)分别表示第l条路径的传播信号在接收端和发射端的阵列流形矢量。假设发射端和接收端均采用均匀线阵,则ar(φl)和at(θl)具有以下形式
其中,λ为信号的波长,d表示相邻天线的间距。通过对时域信道H(g)进行Nc点的DFT,第k个子载波上的频域信道矩阵可以表示为如下形式
其中,
将L条路径的信道联合考虑,令
根据(9)式,第k个子载波的频域信道矩阵可以表示为更加紧凑的矩阵形式
通过(11)和(12)可以看出,毫米波频率选择性信道具有一定的结构特征,信道估计问题可以表征为信道参数矩阵Ar,和At的估计问题。下一节将基于上述信号和信道模型提出新的信道估计算法。
(2)基于毫米波MIMO频率选择性信道的稀疏结构,将信道估计问题表征为稀疏重构问题。
毫米波信号传播过程中通常具有有限散射特性,导致毫米波信道在角度域具有一定的稀疏性。定义毫米波信号传播路径的AoA和AoD的角度字典分别为:
其中,为预定义的出发角和到达角的量化角度。不失一般性,假设/>和/>在[0,π]之间均匀取值,相邻角度之间的间隔分别为/>和/>此时,量化角度取值分别为
假设毫米波信道中传播路径出发角和到达角位于和/>构成的量化角度栅格内,此时,Hk可以近似表示为下列形式
其中,是一个M×N维的稀疏矩阵,该矩阵具有L个非零元素值,其非零值为第k个子载波在第l条路径的频域衰落系数βl,k,每个非零值的索引为该路径的到达角和方向角在其量化角度区间的位置。/>通常被称为Hk的虚拟信道。根据/>的结构可知,如果/>可以确定,则信道矩阵Hk可以确定,因此,信道估计问题可以转化为稀疏矩阵/>的估计问题。如图2所示,为了方便分析,将(15)带入(6),接收端第k子载波的接收信号可以得到如下的表达式
其中,此时,毫米波MIMO系统的频率选择性信道的估计问题可以表征为Nc个子载波虚拟信道的稀疏重构问题,即
其中,为容忍参数,其值由每个子载波的噪声功率确定。根据压缩感知理论可知,子载波的虚拟信道/>可以通过使用压缩感知中的稀疏恢复算法分别进行估计,然后根据(15)得到/>的估计值。但是,当子载波数目较多时,上述信道估计过程的复杂度是比较高的。
根据信道模型可知,不同的子载波上的信号具有相同的散射路径,因此,其对应传播路径的出发角和到达角相同,但是,由于信道具有频率选择性,不同子载波上经历的信道衰落不同。因此,对比(12)和(15)可以发现,虽然不同子载波的虚拟信道中的非零元素的值不同,但是非零元素的位置相同。为了更好的表征这种结构特征,根据(16)式给出所有子载波接收信号的向量形式:
其中,yk=vec(Yk),(4.17)式中,由于虚拟信道/>仅包含L个非零元素,因此/>为MN×1的稀疏向量,其稀疏度为L。同时,由于/>具有相同的稀疏结构,使得/>中非零元素的位置相同。因此,可以将Nc个子载波上的接收信号联合进行考虑,令/>得到如下形式
其中,
由于为稀疏度为L的稀疏向量,因此,矩阵Hv是一个稀疏度为NcL的稀疏矩阵,而且该矩阵的每一列有且仅有L个非零元素。但是,与传统的稀疏矩阵不同,由于中非零元素的位置相同,所以Hv非零元素集中在对应的L行中。此时,Nc个子载波的信道的联合估计问题可以表示为如下形式:
其中,Hv(:,k)表示Hv的第k列,从(21)式可以看出,毫米波系统中不同子载波信道间具有特殊的稀疏结构,充分利用该结构可以提升信道估计的性能。
(3)利用压缩感知的思想求取重构问题的最优解,并计算出信道估计值。
为了方便表述,令测量矩阵构建索引集合/>和/>用于存储算法迭代过程中得到的非零索引的估计值,在初始化阶段,/>和/>均设置为空集。在第i次迭代的过程中,算法首先将前一次迭代中得到的所有子载波残差信号与测量矩阵进行相关操作,得到相关矩阵作为稀疏信号的参考。提取相关矩阵中每个子载波对应的相关向量中绝对值最大的元素的索引并将其存入集合/>中。注意到,虽然集合/>中包含从Nc个子载波得到的索引值,但是由于部分子载波得到的索引值可能相同,因此/>中不同元素的个数小于或等于子载波的数目,记/>中互异元素的个数为Na,其值分别为/>则Na≤Nc。
其次,计算集合中不同元素出现的频率,记为/>令
将中出现频率最大元素amax存入非零元素索引集合/>需要指出的是,可能会出现N个元素在/>中的频率相等且均大于其他元素出现频率的情况,此时可以对这些索引进行二次选择。假设存在N个索引元素/>满足下列不等式
此时,分别计算残差信道相关矩阵中第行的元素绝对值的累加和,记为
令
将amax存入非零元素集合假设本次迭代索引集合/>中的元素为/>根据/>中的元素和测量矩阵Φ构建新的测量矩阵/>
基于当前索引集合的信道估计值可以通过求解下式得到:
最后,使用和/>计算本次迭代的信号残差R供下次迭代使用
当迭代终止条件满足时,算法结束,得到和集合/>否则,重复上述过程直至迭代收敛条件满足为止。JSD-OMP算法结束后,矩阵/>包含所有子载波上的路径衰落系数的估计值,其中,/>为第k个子载波的虚拟信道/>的向量化形式/>中非零元素的估计值,/>中的索引为/>中的元素在稀疏向量/>和对应的位置索引估计值。因此,根据/>和/>可以重构/>的估计值/>对/>中的元素进行矩阵重排可以得到/>的估计值,即
根据(4.15),第k个子载波的信道估计值可以表示为
本发明采用信道估计的NMSE作为性能指标对算法的信道估计性能进行评估,并结合图3的仿真实验进行分析。NMSE的表达式由下式给出
本发明的仿真实验中,毫米波MIMO系统的收发两端均采用均匀阵列,每条散射路径的出发角和到达角在[0,π]间随机取值,路径衰落系数、预编码矩阵和合并矩阵均随机取值,且服从均值为0,方差为1的高斯分布。
由图3所示,本发明提出的算法信道估计性能均优于传统的正交匹配追踪算法,这是因为本发明提出的算法在迭代过程中利用了子载波信道之间的联合稀疏结构,使用所有子载波的样本数据联合估计稀疏信道的非零元素的索引,有效的提高了稀疏信道索引估计的准确度,从而提升了信道估计的性能。
Claims (10)
1.一种毫米波MIMO频率选择性信道估计方法,其特征在于:包括如下步骤:
(1)构建毫米波MIMO频率选择性信道模型,表示为紧凑的矩阵形式;
(2)基于毫米波MIMO频率选择性信道的稀疏结构,将信道估计问题表征为稀疏重构问题;
(3)利用压缩感知的思想求取重构问题的最优解,并计算出信道估计值。
2.根据权利要求1所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述步骤(1)包括:
(1.1)构建毫米波MIMO下行传输系统的接收信号模型;
(1.2)构建毫米波MIMO频率选择性信道模型。
3.根据权利要求2所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述步骤(1.1)假设发射端配置的天线数目为Nt,接收端配置的天线数目为Nr,信道为频率选择性信道,系统使用OFDM技术将传输频带划分为Nc个正交子载波进行数据传输,且每个子载波带宽小于信道的相干带宽,每个子载波的信道为平坦衰落信道,发送端对每个子载波上的传输信号进行波束赋形,第k个子载波上发送的第n个信号xk,n可以表示为:
xk,n=fk,nsk,n
其中,sk,n表示第k个子载波上传输的第n个导频符号,fk,n表示第k个子载波上第n个导频符号使用的波束赋形向量。
4.根据权利要求2所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述步骤(1.1)假设所有子载波使用相同的波束赋形矩阵和合并矩阵,即Fk=F,Wk=W,k=1,...,Nc,每个子载波使用的相同的具有单位功率的导频符号,第k个子载波上的接收信号可以简化为:
Yk=WHHkF+WHNk
其中,F表示波束赋形矩阵。
5.根据权利要求2所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述步骤(1.2)根据毫米波信道的有限散射特性,延迟为g个采样周期的信道的时域形式可表示为:
其中,L表示散射体的数目。考虑每个散射体只有单条散射路径的情况,αl表示第l个散射体的散射路径的衰落系数;r(τ)表示通信系统中包含的脉冲成形滤波器、匹配滤波器等在τ时刻对信号带来的综合影响,Ts表示采样周期,τl表示第l条路径的延迟;θl和φl分别表示第l条信号传播路径的出发角和到达角,ar(φl)和at(θl)分别表示第l条路径的传播信号在接收端和发射端的阵列流形矢量,
通过对时域信道H(g)进行Nc点的DFT,第k个子载波上的频域信道矩阵可以表示为更加紧凑的矩阵形式
其中,L表示路径数。
6.根据权利要求5所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述发射端和接收端采用均匀线阵,ar(φl)和at(θl)形式可表示为:
其中,λ为信号的波长,d表示相邻天线的间距。
7.根据权利要求1所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述步骤(2)假设毫米波信道中传播路径出发角和到达角位于和/>构成的量化角度栅格内,此时,Hk可以近似表示为下列形式:
其中,是一个M×N维的稀疏矩阵,该矩阵具有L个非零元素值,其非零值为第k个子载波在第l条路径的频域衰落系数βl,k,每个非零值的索引为该路径的到达角和方向角在其量化角度区间的位置,/>通常被称为Hk的虚拟信道,信道估计问题可以转化为稀疏矩阵/>的估计问题,Nc个子载波的信道的联合估计可以表示为如下形式:
supp(Hv(:,1))=,...,=supp(Hv(:,Nc))
其中,Hv(:,k)表示Hv的第k列。
8.根据权利要求1所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述步骤(3)包括令测量矩阵构建索引集台/>和/>用于存储算法迭代过程中得到的非零索引的估计值,在第i次迭代的过程中,算法首先将前一次迭代中得到的所有子载波残差信号与测量矩阵进行相关操作,得到相关矩阵作为稀疏信号的参考,提取相关矩阵中每个子载波对应的相关向量中绝对值最大的元素的索引并将其存入集合/>中。
9.根据权利要求8所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述步骤(3)计算集合中不同元素出现的频率,记为/>令
将中出现频率最大元素amax存入非零元素索引集合/>假设存在N个索引元素/>满足下列不等式
分别计算残差信道相关矩阵中第行的元素绝对值的累加和,记为/>令
将amax存入非零元素集合假设本次迭代索引集合/>中的元素为/>根据/>中的元素和测量矩阵Φ构建新的测量矩阵/>
基于当前索引集合的信道估计值为/>
10.根据权利要求9所述的毫米波MIMO频率选择性信道估计方法,其特征在于:所述步骤(3)使用和/>计算本次迭代的信号残差R供下次迭代使用
当迭代终止条件满足时,算法结束,得到和集合/>否则,重复上述过程直至迭代收敛条件满足为止,
算法结束后,矩阵包含所有子载波上的路径衰落系数的估计值,其中,为第k个子载波的虚拟信道/>的向量化形式/>中非零元素的估计值,/>中的索引为/>和的元素在稀疏向量/>中对应的位置索引估计值,根据/>和/>重构/>的估计值对/>中的元素进行矩阵重排得到/>的估计值,即
根据第k个子载波的信道估计值为/>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311069563.1A CN117014258A (zh) | 2023-08-23 | 2023-08-23 | 一种毫米波mimo频率选择性信道估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311069563.1A CN117014258A (zh) | 2023-08-23 | 2023-08-23 | 一种毫米波mimo频率选择性信道估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117014258A true CN117014258A (zh) | 2023-11-07 |
Family
ID=88565484
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311069563.1A Pending CN117014258A (zh) | 2023-08-23 | 2023-08-23 | 一种毫米波mimo频率选择性信道估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117014258A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117254994A (zh) * | 2023-11-20 | 2023-12-19 | 南京邮电大学 | 基于定秩矩阵流形上近端梯度算法的稀疏信道估计方法 |
-
2023
- 2023-08-23 CN CN202311069563.1A patent/CN117014258A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117254994A (zh) * | 2023-11-20 | 2023-12-19 | 南京邮电大学 | 基于定秩矩阵流形上近端梯度算法的稀疏信道估计方法 |
CN117254994B (zh) * | 2023-11-20 | 2024-03-15 | 南京邮电大学 | 基于定秩矩阵流形上近端梯度算法的稀疏信道估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106559367B (zh) | 基于低秩张量分解的mimo-ofdm系统毫米波信道估计方法 | |
CN108933745B (zh) | 一种基于超分辨率角度和时延估计的宽带信道估计方法 | |
CN109743268B (zh) | 基于深度神经网络的毫米波信道估计和压缩方法 | |
CN104052691B (zh) | 基于压缩感知的mimo‑ofdm系统信道估计方法 | |
US8594221B2 (en) | Model-based channel estimator for correlated fading channels and channel estimation method thereof | |
CN110099016B (zh) | 一种基于深度学习网络的毫米波稀疏阵面信道估计方法 | |
CN105891771B (zh) | 一种提高估计精度的基于连续分布的角度估计方法与设备 | |
CN109005133B (zh) | 双稀疏多径信道模型及基于此模型的信道估计方法 | |
Zhang et al. | Channel estimation and training design for hybrid multi-carrier mmwave massive MIMO systems: The beamspace ESPRIT approach | |
CN101494627A (zh) | 宽带移动通信中利用压缩感知减少导频数的信道估计方法 | |
CN108881074B (zh) | 一种低精度混合架构下宽带毫米波信道估计方法 | |
CN108471325B (zh) | 一种稀疏射频/基带混合预编码方法 | |
CN103634038A (zh) | 基于多天线的联合doa估计和波束形成的多径信号接收方法 | |
CN104698430A (zh) | 一种用于提高精度的基于虚拟天线阵列的角度估计方法 | |
CN117014258A (zh) | 一种毫米波mimo频率选择性信道估计方法 | |
CN112769462B (zh) | 一种基于联合参数学习的毫米波mimo宽带信道估计方法 | |
CN108259397B (zh) | 基于自适应正则化子空间追踪压缩感知算法的大规模mimo系统信道估计方法 | |
CN110034916B (zh) | 基于终端反馈的天线相位同步与信道互易性校准的方法 | |
CN113193895B (zh) | 大规模mimo信道状态信息获取方法、系统及计算机存储介质 | |
CN106972875B (zh) | 一种mimo系统下多维联合估计动态稀疏信道的方法 | |
CN114879173A (zh) | 基于mimo-ofdm信号的感知性能优化方法 | |
Shafin et al. | DoA estimation and RMSE characterization for 3D massive-MIMO/FD-MIMO OFDM system | |
CN105812032B (zh) | 基于波束块结构压缩感知的信道估计方法 | |
CN112803977B (zh) | 波束偏移效应下毫米波通信系统的混合预编码方法 | |
CN114629533A (zh) | 大规模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 |