CN110519189B - 高度移动场景下基于压缩感知的毫米波信道估计方法 - Google Patents
高度移动场景下基于压缩感知的毫米波信道估计方法 Download PDFInfo
- Publication number
- CN110519189B CN110519189B CN201910822199.9A CN201910822199A CN110519189B CN 110519189 B CN110519189 B CN 110519189B CN 201910822199 A CN201910822199 A CN 201910822199A CN 110519189 B CN110519189 B CN 110519189B
- Authority
- CN
- China
- Prior art keywords
- matrix
- channel
- millimeter wave
- algorithm
- iteration
- 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
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
-
- 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
- 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
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
-
- 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
-
- 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
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Mathematical Physics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供了一种高速移动场景下基于压缩感知的毫米波信道估计方法。包括:步骤1:分析高速移动场景下的下行链路通信,考虑多普勒频移的影响,基于毫米波通信的低秩特性以及空间相关特性,进行信道建模,写成矩阵的形式;步骤2:对步骤1中所述的信道矩阵进行量化处理后,此时信道矩阵会呈现出稀疏特性,设计感知矩阵以及测量矩阵;步骤3:对步骤2中呈现出稀疏特性的信道矩阵进行向量化处理,基于压缩感知算法以及步骤2中所述的感知矩阵与测量矩阵对向量化处理后的信道模型进行重构,计算归一化均方误差。本发明降低计算成本,提升了性能。
Description
技术领域
本发明属于无线通信技术领域,涉及在高度移动场景下基于压缩感知的毫米波信道估计方法。
背景技术
近年来,高速列车(HST)系统迅猛发展,作为一种快速、便捷、绿色的公共交通系统,正在成为世界铁路运输的未来趋势。如何为乘坐高铁的乘客提供快速无缝的无线通信服务,在超高流量密度、超高连接密度和超高移动性等方面对未来的无线通信网络提出了挑战。在传统的HST场景中,存在一些问题:快速移动、超密集的用户分布以及屏蔽良好的车厢,这些会导致多普勒频移严重、穿透损耗高、频繁且大规模的切换过程、移动中断的功耗增加、数据率降低。为了满足未来高速移动通信对数据传输速率的高要求,一个最直接有效的方法是利用更高的频段——毫米波(mmWave)波段。
毫米波(mmWave)通信是未来室外蜂窝系统一种有发展前途的技术,数据速率可达到千兆每秒,这要归功于毫米波频段中可用的高带宽。为了实现足够的链路余量,毫米波系统将在收发端采用具有大规模天线阵列的定向波束成形。高采样率的混合信号设备成本和功耗较高,毫米波预编码可以划分为模拟域和数字域。随着天线数量的增加以及模拟波束成形技术的出现,人们需要构建毫米波系统专用的信道估计算法以及预编码算法。定向预编码采用大规模天线阵列,具有较高的预编码增益,而可以保证室外链路更长,接收信号的功率足够大。通常基于完善的信道状态信息来设计预编码矩阵,必须开发适合于毫米波蜂窝系统的新信道估计和预编码算法。
为了克服RF硬件的局限性,Optimized Microstrip Antenna Arrays forEmerging Millimeter-Wave Wireless Applications中提供了解决方案:模拟波束。其主要理念是基于一个模拟移相器网络以达到控制传输信号相位的目的。提出了解决方案波束训练算法,用于在发射端未知信道信息时,迭代地设计系统中模拟波束成形系数。旨在达到更高的预编码增益,并实现对多个数据流的预编码,Variable-phase-shift-based RF-baseband codesign for MIMO antenna selection中建议预编码操作在模拟域和数字域之间同时进行。空间分集和复用系统都考虑了联合模拟数字预编码设计问题。第一,最优无约束射频预处理信号变换,接着提出了基带预编码矩阵,然后给出了可变移相器约束射频处理时的闭式次最优解。开发了混合模拟、数字预编码算法,当在只有量化相位的移相器可用时,在存在干扰的情况下最小化接收信号的均方误差。然而,工作并不是专门针对毫米波系统的,也没有考虑到毫米波信道的特性。在Analog Beamforming in MIMOCommunications With Phase Shift Networks and Online Channel Estimation中,利用毫米波信道的稀疏多路径结构和基追踪算法,设计用于低复杂度混合预编码器,假设接收端获取的信道信息是完美的,该预编码器试图接近容量。另外一些混合预编码设计问题被考虑到系统中,其中在发送机处已知部分状态信息。开发的混合预编码算法克服了射频硬件的局限性,可以支持多流传输。然而,在设计预编码矩阵之前,实现这些增益需要知晓发射端信道的部分信息。这推动了多径毫米波信道估计算法的发展,使混合预编码的性能能够与数字预编码算法比较接近。
毫米波(mmWave)通信重新定义了无线通信,开启了无线通信的新时代。信号处理在实现下一代毫米波通信至关重要。由于收发机两端均采用大型天线阵列,结合无线频率以及混合信号功率限制,促使了新的多输入多输出(MIMO)通信信号处理技术的诞生。由于宽带的影响,低复杂度收发器算法变得尤为重要。可以利用压缩感知(CS)技术进行信道估计以及波束成形处理。本小节概述了在毫米波无线系统中信号处理的挑战,着重介绍了在高载波频率下MIMO通信所面临的挑战。
毫米波信道模型是不同的,原因在于传播环境对小波信号具有不同的影响。毫米波信道使用低频段系统的常见特性,具有一些不同的参数(例如少数几个聚集的路径可以使信道呈现出更多的稀疏特性)。此外,同时一些新特性被引入用于解释对阻塞物所表现出的高度敏感以及不同传播条件的强烈差异。我们可以利用信道稀疏性的数学特性来进行信道估计。
发明内容
发明目的:为了克服现有技术的不足,本发明提供一种高度移动场景下基于压缩感知的毫米波信道估计方法,降低计算成本,提升了性能。
技术方案:一种高速移动场景下基于压缩感知的毫米波信道估计方法,包括以下步骤:
步骤1:分析高速移动场景下的下行链路通信,考虑多普勒频移的影响,基于毫米波通信的低秩特性以及空间相关特性,进行信道建模,写成矩阵的形式;
步骤2:对步骤1中所述的信道矩阵进行量化处理后,此时信道矩阵会呈现出稀疏特性,设计感知矩阵以及测量矩阵;
步骤3:对步骤2中呈现出稀疏特性的信道矩阵进行向量化处理,基于压缩感知算法以及步骤2中所述的感知矩阵与测量矩阵对向量化处理后的信道模型进行重构,计算归一化均方误差。
所述的高速移动场景下基于压缩感知的毫米波信道估计方法,步骤1中所述的信道建模具体包括如下步骤:
步骤1.1:建立HST毫米波MIMO信号模型:
发射端配置Nt根天线,接收端的天线数目为Nr;
假设发送端用向量fp进行预编码,接收端用向量wq对训练序列进行合成处理,则接收端的信号可以表示成
其中sp为波束成形向量fp上发送的训练序列,P为每个发送训练符号的平均能量,假设在Mr个连续时隙,接收端用Mr个合成向量wq,q=1,2,…,Mr对经预编码处理后的训练符号进行合成处理,则接收端输出的信号矩阵的表达式为
其中假设发送端在Mt个连续时隙用Mt个预编码向量fp,p=1,2,...,Mt进行预编码处理,接收端用矩阵W合成接收的训练序列,则接收端的输出信号可以表示成Mt个向量yp的组合,p=1,2,...,Mt:
Y=WHHFS+Q (3)
步骤1.2:建立HST毫米波信道模型:
基于毫米波通信的低秩特性以及空间相关特性,Saleh-Vjenzuela信道建模为:
L是对应于有限数量散射体的信道路径,是第l径的增益,不同路径的增益是独立同分布的,均服从复高斯分布 和分别是AoAs、AoDs;和分别代表在特定AoA和AoD的接收和发射天线增益,为了简化且不失一般性,和在AoAs/AoDs的范围内设置成1;和分别是收发两端阵列响应向量,向量的表达式取决于天线阵列结构;
对于U元均匀线性阵列,阵列响应表达式表示为:
其中λ是信号波长,d是相邻天线元之间的距离,通常假设d=λ/2,在(6)中没有写上标{R,S},由于ULA响应向量与仰角无关,也不包含θ;
在高速列车环境下毫米波信道模型表达式为:
其中vl代表l径的多普勒频移。与典型的城市环境不同的是,视距(LoS)路径总是存在于快速时变的高速列车(HST)环境中l=1,其他路径是非视距(NLoS)/反射路径;
应用ULA,则(7)改写成
在公式(8)中AoA和AoD服从拉普拉斯分布,其均值均匀分布在[-π,π],标准差为σAS。
所述的一种高速移动场景下基于压缩感知的毫米波信道估计方法,步骤2中所述信道矩阵进行量化处理后,设计感知矩阵以及测量矩阵,具体是指:对所述的HST毫米波信道模型中的信道矩阵对信道矩阵进行量化处理后,设计感知矩阵以及测量矩阵,包括如下步骤:
步骤2.1:将(8)换算写成如下的矩阵形式
H=ARHaAT (9)
矩阵是NtNr×L的矩阵,每一列都具有的形式,也就是说每一列代表与第l径的AoA、AoD相关的发送端和接收端的阵列响应的克罗内克积;实际上,AoAs和AoDs是连续的,将其量化在Gr和Gt点的均匀网格中,GrGt>>L,G∈{Gr,Gt},忽略网格量化误差,信道矩阵可以表示为
当Nt=Gr时,AR,d为IDFT矩阵;
将(10)中的yv近似成
令M=MrMt,N=GrGt,则za是包含L个非零元素的稀疏向量,稀疏度L<<N,将(14)看成是接收信号的稀疏处理,利用CS理论设计信道估计算法,即M个测量向量来估计za的非零元素(M<<N),定义测量矩阵感知矩阵则(14)可以简写成
yv=Qza+nQ (15)。
所述的一种高速移动场景下基于压缩感知的毫米波信道估计方法,步骤3中所述对呈现出稀疏特性的信道矩阵进行向量化处理,基于压缩感知算法以及步骤2中所述的感知矩阵与测量矩阵对向量化处理后的信道模型进行重构,计算归一化均方误差的具体方法为:
接收端使用yv中M个元素重构za,重构就是对(15)求解za中的N个元素,方程式未知数的个数多于方程组的个数,因此za的解不唯一,依据一些准则可以搜索到最优解,从而重建za,因此,za的重建可以转换成求解最优解的问题:
在衡量HST毫米波信道估计算法的性能有两种度量:一种是归一化均方误差,即NMSE,NMSE的表达式为
其中是对实际HST毫米波信道H的估计,另一种度量是恢复的成功率,在使用贪婪算法对HST毫米波信道估时,成功率可以由成功恢复观测信号的次数除以运行的总次数的次数求得,运行过程中,如果残差小于10-6则认为恢复成功;
在高速移动场景下,建模的信道呈现出稀疏性,因此利用信道低秩的特性,基于OMP算法对信道进行估计,算法总结如下:
B.找出最佳的的AoA/AoD对:找出r和Q的列qj的内积最大值所对应的下标λj;
C.更新AoA/AoD的集合:
(a)更新索引集:
Λt=Λt-1∪{λj} (20)
(b)重建原子集合:
D.求解最小二乘解:
E.更新残差:
F.更新迭代的索引号:t=t+1,判断||rt-1-rt||2>δ,如果成立,返回B进行迭代,否则进入步骤G;
G.重建信道矩阵H。
所述的一种高速移动场景下基于中压缩感知的毫米波信道估计方法,步骤B中所述的找出最佳的的AoA/AoD对,采用内积准则
gt=argmax|<rt-1,ψ>| (24)
(24)对残差值和字典矩阵元素进行相似度匹配计算,g的含义是此次迭代中选择与残差值相似程度最高的元素,ψ代表字集合,绝对值越大意味着残差值越接近于挑选出来的元素,本质上,它是求解冗余字典中的匹配因子与残差值夹角余弦值。当余弦值变大时,向量之间夹角变小,相似程度也会相应地变高,内积准则在度量相似度时,初始信号的部分重要信息会在匹配过程中被丢弃。为了更高的匹配度,改进步骤B,引入广义Dice系数匹配准则:
假设x和y是两个任意向量,x=(x1,x2,…,xn),y=(y1,y2,...,yn),则(25)给出内积法的定义:
在(25)中,sim()代表两个向量的相似程度,由x和y的夹角余弦值求得,OMP算法基于内积法进行迭代计算找出与残差匹配的最优元素,直到字典库匹配计算的完成,但是迭代过程中会造成部分重要信息的缺失,利用算数平均的方法可以有效解决这个问题,从而可以替代内积匹配准则,广义Dice系数可以定义为
广义Dice系数的分母是对向量长度的平方求算数平均:1/2(d2(x)+d2(y)),可以保留重要的信号,有效地区分相似原子,挑选出最佳的匹配原子,提高了信号重构的性能,运用Dice系数因子匹配的OMP算法与OMP算法的改进之处是将匹配准则引入了Dice系数,算法更新如下:
Λt=Λt-1∪{λj} (28)
将该算法记为Dice_OMP算法。
所述的一种高速移动场景下基于中压缩感知的毫米波信道估计方法,步骤B中所述的找出最佳的的AoA/AoD对,采用如下步骤:找出最佳的S个AoA/AoD对:计算u=abs(D(rt-1,Qt)),选择出u中最大的S个值,将这些值对应Q的列序号λj构成列序号结合J0。
所述的一种高速移动场景下基于中压缩感知的毫米波信道估计方法,所述在高速移动场景下,建模的信道呈现出稀疏性,因此利用信道低秩的特性,基于SAMP的算法对HST毫米波信道估计算法归纳如下:
输入观测矩阵Q,观测信号za,步长s;在信道估计的过程中,rt代表第t次迭代的残差,Λt是第t次迭代选出的列索引号的集合(元素个数为LF,其中LF是步长s的整数倍),qj代表Q的第j列,候选集Ck,表示按照Ck选出的Q的列的集合(Lt列),每次迭代计算的za的最小二乘解为Lt×1的列向量,初始化:r0=za,t=1,LF=s;计算u=abs|QTrt-1|,并从u中选择出LF个最大值,将这些值对应Q的列序号j组成列序号的集合Sk;Ck=Λt-1∪Sk,Q={qj},all j∈Ck;求解再从中选出绝对值最大的LF项放在Qt中对应的LF列记为QtL,相应的列索引号记为ΛtL,集合F=ΛtL;更新残差:如果||rnew||2≥||rt-1||2,更新步长,返回继续迭代(LF=LF+s),如果rnew=0,则停止迭代。
所述的一种高速移动场景下基于压缩感知的毫米波信道估计方法,所述基于SAMP的算法对HST毫米波信道估计算法进一步改进如下:
首先广义Dice系数;其次,在预选阶段,F的大小等于步长,固定步长s=1,在HST毫米波信道估计的初始阶段,采用一个步长,当估计的稀疏度逐渐增加,接近于真实的稀疏度,通过逐渐减小步长的方法来改善信道矩阵的恢复精度,考虑指数函数:
g(x)=ax (29)
其中a∈(0,1),是一个固定的常数,当x≥0时,函数是值域为(0,1]的单调递减函数,(29)的导数为g′(x)=axlna,当x增加时,g′(x)会减小,因此斜率单调递减,根据g(x)的变化趋势对步长进行调整,恢复信道向量的初始阶段选取一个较大的s,F的长度LF=LF+s,用来提高信道恢复的效率;第二阶段逐步调整步长,F的长度为其中stage代表迭代阶段,为向上取整,因此估计稀疏度与实际的稀疏度越接近时,F长度LF的改变随着stage的增加减小,直到LF=LF+1,记为Dice_LSAMP算法,具体步骤如下:输入:M×N观测矩阵Q,M维观测信号za,初始步长s;输出:稀疏信号的估计
B.预选出L对最佳的AoA/AoD:计算u=abs(D(rt-1,Qt)),从u中选择出LF个最大值,将这些值对应Q的列序号j组成列序号的集合St;
C.候选AoA/AoD集合:
Ct=Ft-1∪St,Q={qj},all j∈Ct (30)
D.求解最小二乘解:
E.计算残差:
F.如果||r||2≤ε1,转到G,否则转到H;
H.||r||2≥||rt-1||2,更新支撑集长度LF=LF+s和迭代阶段stage=stage+1,转到B;否则更新Ft=F,rt=r,t=t+1,转到B进行下一次迭代;
有益效果:与现有技术相比,本发明的优点在于:改进的压缩感知算法用于毫米波信道估计时可以找出更加匹配的原子,改善重构性能;在此基础上,每次迭代选择多个原子,有效地降低了累积误差的影响;在SAMP的基础上,引入广义Dice系数,改进固定步长,采用变步长的思想,降低计算成本,提升了性能。
附图说明
图1为OMP改进算法的重构性能图,其中图1(a)为采用Dice_OMP算法的重构性能图,
图1(b)为采用Dice_MOMP算法的重构性能图。
图2为Dice_LSAMP与SAMP算法重构性能对比曲线图。
图3为参数a对Dice_LSAMP算法重构成功率影响曲线。
图4为不同重构算法性能对比曲线图,其中图4(a)是不同算法的观测矩阵维度对比图;图4(b)是不同算法的稀疏度对比图。
图5为不同估计算法随信噪比变化的性能曲线图。
图6是本发明的方法流程图。
具体实施方式
一种高速移动场景下基于压缩感知的毫米波信道估计方法,包括以下步骤:
步骤1:分析高速移动场景下的下行链路通信,考虑多普勒频移的影响,基于毫米波通信的低秩特性以及空间相关特性,进行信道建模,写成矩阵的形式;
步骤2:对步骤1中所述的信道矩阵进行量化处理后,此时信道矩阵会呈现出稀疏特性,设计感知矩阵以及测量矩阵;
步骤3:对步骤2中呈现出稀疏特性的信道矩阵进行向量化处理,基于压缩感知算法以及步骤2中所述的感知矩阵与测量矩阵对向量化处理后的信道模型进行重构,计算归一化均方误差。
所述的高速移动场景下基于压缩感知的毫米波信道估计方法,步骤1中所述的信道建模具体包括如下步骤:
步骤1.1:建立HST毫米波MIMO信号模型:
发射端配置Nt根天线,接收端的天线数目为Nr;
假设发送端用向量fp进行预编码,接收端用向量wq对训练序列进行合成处理,则接收端的信号可以表示成
其中sp为波束成形向量fp上发送的训练序列,P为每个发送训练符号的平均能量,假设在Mr个连续时隙,接收端用Mr个合成向量wq,q=1,2,…,Mr对经预编码处理后的训练符号进行合成处理,则接收端输出的信号矩阵的表达式为
Y=WHHFS+Q (3)
步骤1.2:建立HST毫米波信道模型:
基于毫米波通信的低秩特性以及空间相关特性,Saleh-Vjenzuela信道建模为:
L是对应于有限数量散射体的信道路径,是第l径的增益,不同路径的增益是独立同分布的,均服从复高斯分布 和分别是AoAs、AoDs;和分别代表在特定AoA和AoD的接收和发射天线增益,为了简化且不失一般性,和在AoAs/AoDs的范围内设置成1;和分别是收发两端阵列响应向量,向量的表达式取决于天线阵列结构;
对于U元均匀线性阵列,阵列响应表达式表示为:
其中λ是信号波长,d是相邻天线元之间的距离,通常假设d=λ/2,在(6)中没有写上标{R,S},由于ULA响应向量与仰角无关,也不包含θ;
在高速列车环境下毫米波信道模型表达式为:
其中vl代表l径的多普勒频移。与典型的城市环境不同的是,视距(LoS)路径总是存在于快速时变的高速列车(HST)环境中l=1,其他路径是非视距(NLoS)/反射路径;
应用ULA,则(7)改写成
在公式(8)中AoA和AoD服从拉普拉斯分布,其均值均匀分布在[-π,π],标准差为σAS。
所述的一种高速移动场景下基于压缩感知的毫米波信道估计方法,步骤2中所述信道矩阵进行量化处理后,设计感知矩阵以及测量矩阵,具体是指:对所述的HST毫米波信道模型中的信道矩阵对信道矩阵进行量化处理后,设计感知矩阵以及测量矩阵,包括如下步骤:
步骤2.1:将(8)换算写成如下的矩阵形式
H=ARHaAT (9)
矩阵是NtNr×L的矩阵,每一列都具有的形式,也就是说每一列代表与第l径的AoA、AoD相关的发送端和接收端的阵列响应的克罗内克积;实际上,AoAs和AoDs是连续的,将其量化在Gr和Gt点的均匀网格中,GrGt>>L,G∈{Gr,Gt},忽略网格量化误差,信道矩阵可以表示为
当Nt=Gr时,AR,d为IDFT矩阵;
将(10)中的yv近似成
令M=MrMt,N=GrGt,则za是包含L个非零元素的稀疏向量,稀疏度L<<N,将(14)看成是接收信号的稀疏处理,利用CS理论设计信道估计算法,即M个测量向量来估计za的非零元素(M<<N),定义测量矩阵感知矩阵则(14)可以简写成
yv=Qza+nQ (15)。
所述的一种高速移动场景下基于压缩感知的毫米波信道估计方法,步骤3中所述对呈现出稀疏特性的信道矩阵进行向量化处理,基于压缩感知算法以及步骤2中所述的感知矩阵与测量矩阵对向量化处理后的信道模型进行重构,计算归一化均方误差的具体方法为:
接收端使用yv中M个元素重构za,重构就是对(15)求解za中的N个元素,方程式未知数的个数多于方程组的个数,因此za的解不唯一,依据一些准则可以搜索到最优解,从而重建za,因此,za的重建可以转换成求解最优解的问题:
在衡量HST毫米波信道估计算法的性能有两种度量:一种是归一化均方误差,即NMSE,NMSE的表达式为
其中是对实际HST毫米波信道H的估计,另一种度量是恢复的成功率,在使用贪婪算法对HST毫米波信道估时,成功率可以由成功恢复观测信号的次数除以运行的总次数的次数求得,运行过程中,如果残差小于10-6则认为恢复成功;
在高速移动场景下,建模的信道呈现出稀疏性,因此利用信道低秩的特性,基于OMP算法对信道进行估计,算法总结如下:
B.找出最佳的的AoA/AoD对:找出r和Q的列qj的内积最大值所对应的下标λj;
C.更新AoA/AoD的集合:
(a)更新索引集:
Λt=Λt-1∪{λj} (20)
(b)重建原子集合:
D.求解最小二乘解:
E.更新残差:
F.更新迭代的索引号:t=t+1,判断||rt-1-rt||2>δ,如果成立,返回B进行迭代,否则进入步骤G;
G.重建信道矩阵H。
所述的一种高速移动场景下基于中压缩感知的毫米波信道估计方法,步骤B中所述的找出最佳的的AoA/AoD对,采用内积准则
gt=argmax|<rt-1,ψ>| (24)
(24)对残差值和字典矩阵元素进行相似度匹配计算,g的含义是此次迭代中选择与残差值相似程度最高的元素,ψ代表字集合,绝对值越大意味着残差值越接近于挑选出来的元素,本质上,它是求解冗余字典中的匹配因子与残差值夹角余弦值。当余弦值变大时,向量之间夹角变小,相似程度也会相应地变高,内积准则在度量相似度时,初始信号的部分重要信息会在匹配过程中被丢弃。为了更高的匹配度,改进步骤B,引入广义Dice系数匹配准则:
假设x和y是两个任意向量,x=(x1,x2,...,xn),y=(y1,y2,...,yn),则(25)给出内积法的定义:
在(25)中,sim()代表两个向量的相似程度,由x和y的夹角余弦值求得,OMP算法基于内积法进行迭代计算找出与残差匹配的最优元素,直到字典库匹配计算的完成,但是迭代过程中会造成部分重要信息的缺失,利用算数平均的方法可以有效解决这个问题,从而可以替代内积匹配准则,广义Dice系数可以定义为
广义Dice系数的分母是对向量长度的平方求算数平均:1/2(d2(x)+d2(y)),可以保留重要的信号,有效地区分相似原子,挑选出最佳的匹配原子,提高了信号重构的性能,运用Dice系数因子匹配的OMP算法与OMP算法的改进之处是将匹配准则引入了Dice系数,算法更新如下:
Λt=Λt-1∪{λj} (28)
将该算法记为Dice_OMP算法。
所述的一种高速移动场景下基于中压缩感知的毫米波信道估计方法,步骤B中所述的找出最佳的的AoA/AoD对,采用如下步骤:找出最佳的S个AoA/AoD对:计算u=abs(D(rt-1,Qt)),选择出u中最大的S个值,将这些值对应Q的列序号λj构成列序号结合J0。
所述的一种高速移动场景下基于中压缩感知的毫米波信道估计方法,所述在高速移动场景下,建模的信道呈现出稀疏性,因此利用信道低秩的特性,基于SAMP的算法对HST毫米波信道估计算法归纳如下:
输入观测矩阵Q,观测信号za,步长s;在信道估计的过程中,rt代表第t次迭代的残差,Λt是第t次迭代选出的列索引号的集合(元素个数为LF,其中LF是步长s的整数倍),qj代表Q的第j列,候选集Ck,表示按照Ck选出的Q的列的集合(Lt列),每次迭代计算的za的最小二乘解为Lt×1的列向量,初始化:r0=za,t=1,LF=s;计算u=abs|QTrt-1|,并从u中选择出LF个最大值,将这些值对应Q的列序号j组成列序号的集合Sk;Ck=Λt-1∪Sk,Q={qj},all j∈Ck;求解再从中选出绝对值最大的LF项放在Qt中对应的LF列记为QtL,相应的列索引号记为ΛtL,集合F=ΛtL;更新残差:如果||rnew||2≥||rt-1||2,更新步长,返回继续迭代(LF=LF+s),如果rnew=0,则停止迭代。
所述的一种高速移动场景下基于压缩感知的毫米波信道估计方法,所述基于SAMP的算法对HST毫米波信道估计算法进一步改进如下:
首先广义Dice系数;其次,在预选阶段,F的大小等于步长,固定步长s=1,在HST毫米波信道估计的初始阶段,采用一个步长,当估计的稀疏度逐渐增加,接近于真实的稀疏度,通过逐渐减小步长的方法来改善信道矩阵的恢复精度,考虑指数函数:
g(x)=ax (29)
其中a∈(0,1),是一个固定的常数,当x≥0时,函数是值域为(0,1]的单调递减函数,(29)的导数为g′(x)=axlna,当x增加时,g′(x)会减小,因此斜率单调递减,根据g(x)的变化趋势对步长进行调整,恢复信道向量的初始阶段选取一个较大的s,F的长度LF=LF+s,用来提高信道恢复的效率;第二阶段逐步调整步长,F的长度为其中stage代表迭代阶段,为向上取整,因此估计稀疏度与实际的稀疏度越接近时,F长度LF的改变随着stage的增加减小,直到LF=LF+1,记为Dice_LSAMP算法,具体步骤如下:输入:M×N观测矩阵Q,M维观测信号za,初始步长s;输出:稀疏信号的估计
B.预选出L对最佳的AoA/AoD:计算u=abs(D(rt-1,Qt)),从u中选择出LF个最大值,将这些值对应Q的列序号j组成列序号的集合St;
C.候选AoA/AoD集合:
Ct=Ft-1∪St,Q={qj},all j∈Ct (30)
D.求解最小二乘解:
E.计算残差:
F.如果||r||2≤ε1,转到G,否则转到H;
H.||r||2≥||rt-1||2,更新支撑集长度LF=LF+s和迭代阶段stage=stage+1,转到B;否则更新Ft=F,rt=r,t=t+1,转到B进行下一次迭代;
下面结合附图对本发明的技术方案作进一步说明。
为了验证OMP算法的重构恢复性能,随机生成稀疏度为10,长度为256的信号x,采用高斯矩阵作为测量矩阵。观测向量数M≥Ο(Llog(N/L)),在这里取值为64。OMP算法单次重构初始信号,重构时间为0.056089秒,输出残差为6.7532e-15。
验证Dice_OMP以及Dice_MOMP算法的重构效率,在图1(a)中给出了两种算法重构性能的曲线图,选择长度为256的信号,使用Gauss随机稀疏信号矩阵重构信号,在不同维度的测量矩阵中反复迭代1000次,分别得到在稀疏度不同的情况下计算信号的重构成功率。我们发现在观测维度一定时,两种重构算法重构稀疏度较小的信号的成功率较大。稀疏度较大的信号为了获得较大的成功率需要更大的测量维度。此外,Dice_OMP算法在M和L都相同的情况下重构性能要优于OMP算法,因此具有一定的可行性。
图1(b)中验证了Dice_MOMP算法的重构性能,选择不同原子数信号重构的曲线图。Dice_MOMP在Dice_OMP算法的基础上进行了修正,每次迭代匹配搜索筛选出多个原子放入到候选集中。每次迭代筛选出来的多原子集合在一定程度上可以降低OMP算法的误差累积效应的影响。进行验证时,L为10。当S=3时的Dice_MOMP重构曲线图要优于其他算法。当M较小时,OMP算法以及Dice_OMP算法要优于S=6和S=9时的Dice_MOMP算法的重构性能。但是对重构精度要求较高时,所需观测矩阵维度较高,此时Dice_MOMP算法的重构成功率高于OMP算法以及Dice_OMP算法。但是呢,我们也可以发现,当S取值变大时,Dice_MOMP算法的重构成功率反而变小了。分析原因如下:当S=1时,Dice_MOMP算法就变成了Dice_OMP算法,此时的重构性能是低于S取值为3的时候,但是当S取值较大时,每次迭代中挑选的原子数目就会比较多,计算效率降低,此外可能会选出相关性较低的原子,降低了重构效率。
分别比较两种SAMP算法的重构性能,在实验中采用长度N=256,L=20的高斯随机信号,M×N阶高斯随机矩阵作为测量矩阵,固定步长s=5,重复做1000次实验。图2中的横坐标代表观测矩阵维度M,M的取值是30-80,每个间隔为5,纵坐标代表两种算法重构信号的成功率。可以看出,Dice_LSAMP算法的重构性能比SAMP算法要好些。Dice_LSAMP算法在M<40时就可以成功地重构信号了;而只有当M>40时,SAMP算法才能正确重构信号;在重构性能达到最高之前,Dice_LSAMP算法在信号的重建成功率是优于SAMP算法的,这说明了Dice_LSAMP算法在SAMP的基础上又优化了性能。
考虑参数a对Dice_LSAMP算法重构性能的影响。步长s取值为5,a的取值范围在0.1到1之间,间隔为0.1,图3给出了重构性能曲线图,可以看出实验重构性能随着a取值的增加而恶化,当a取值为1时,等价于算法重构过程的采用的步长是固定的,此时重构成功率最低,效果最差。多次迭代证明了Dice_LSAMP算法在a∈[0.3,0.5]时有着较好的重构效果。因此,在后面的仿真中a的取值为0.5。
为了对比不同算法在不同稀疏度或者观测维度下重构信号的成功率,我们在图4中给出了性能曲线图。图4表明:增加测量数可以提升各算法的性能,Dice_LSAMP算法的重构效率是远远超过OMP算法、StOMP算法以及CoSaMP算法。我们可以发现当步长变大时,Dice_LSAMP算法和SAMP算法的重构性能降低。由于Dice_LSAMP算法中采用了指数变步长,其对重构性能影响较小。在步长相同时,Dice_LSAMP算法的重构性能要比SAMP算法要好。图4(b)中比较了算法在不同稀疏度下重构信号的正确率。测量值维度为M=256,稀疏度度的集合L={10,15,20,25,30,35,40},发现各个算法的性能都随着稀疏度的增加而减少。当稀疏度为20时,其他算法的成功率已经降为0,但是Dice_LSAMP和SAMP算法重构率可以仍旧保持在20%。
图5描绘了不同估计算法随信噪比变化的性能曲线图。发送端天线数量Nt=16,接收端天线数量Nr=16,L取值为8,即za中有L个非零元素。Mr和Mt都设置为16,Gr和Gt都设置为64,则M=256,N=4096,满足L<<M<<N。各种信道估计算法的NMSE随着信噪比增加而降低,正确重构的概率越大。估计算法性能从低往高排序:LS,OMP,Dice_OMP,Dice_MOMP,Dice_LSAMP,其中Dice_MOMP算法中S=3,Dice_LSAMP中采用的固定步长为s=5。
应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (8)
1.一种高速移动场景下基于压缩感知的毫米波信道估计方法,其特征在于:包括以下步骤:
步骤3:对步骤2中呈现出稀疏特性的信道矩阵进行向量化处理,基于压缩感知算法以及步骤2中所述的感知矩阵与测量矩阵对向量化处理后的信道模型进行重构,计算归一化均方误差;
步骤1中所述的信道建模具体包括如下步骤:
是对应于有限数量散射体的信道路径,是第径的增益,不同路径的增益是独立同分布的,均服从复高斯分布 和分别是 和分别代表在特定和的接收和发射天线增益,为了简化且不失一般性,和在的范围内设置成1;和分别是收发两端阵列响应向量,向量的表达式取决于天线阵列结构;
在高速列车环境下毫米波信道模型表达式为:
2.根据权利要求1所述的一种高速移动场景下基于压缩感知的毫米波信道估计方法,其特征在于:步骤2中所述信道矩阵进行量化处理后,设计感知矩阵以及测量矩阵,具体是指:对所述的毫米波信道模型中的信道矩阵对信道矩阵进行量化处理后,设计感知矩阵以及测量矩阵,包括如下步骤:
步骤2.1:将(8)换算写成如下的矩阵形式
3.根据权利要求2所述的一种高速移动场景下基于压缩感知的毫米波信道估计方法,其特征在于:步骤3中所述对呈现出稀疏特性的信道矩阵进行向量化处理,基于压缩感知算法以及步骤2中所述的感知矩阵与测量矩阵对向量化处理后的信道模型进行重构,计算归一化均方误差的具体方法为:
(24)对残差值和字典矩阵元素进行相似度匹配计算,的含义是此次迭代中选择与残差值相似程度最高的元素,代表字集合,绝对值越大意味着残差值越接近于挑选出来的元素,本质上,它是求解冗余字典中的匹配因子与残差值夹角余弦值;当余弦值变大时,向量之间夹角变小,相似程度也会相应地变高,内积准则在度量相似度时,初始信号的部分重要信息会在匹配过程中被丢弃;为了更高的匹配度,改进步骤引入广义系数匹配准则:
在(25)中,代表两个向量的相似程度,由和的夹角余弦值求得,算法基于内积法进行迭代计算找出与残差匹配的最优元素,直到字典库匹配计算的完成,但是迭代过程中会造成部分重要信息的缺失,利用算数平均的方法解决这个问题,从而替代内积匹配准则,广义系数定义为
广义系数的分母是对向量长度的平方求算数平均:保留重要的信号,有效地区分相似原子,挑选出最佳的匹配原子,提高了信号重构的性能,运用系数因子匹配的算法与算法的改进之处是将匹配准则引入了系数,算法更新如下:
首先广义系数;其次,在预选阶段,的大小等于步长,固定步长在毫米波信道估计的初始阶段,采用一个步长,当估计的稀疏度逐渐增加,接近于真实的稀疏度,通过逐渐减小步长的方法来改善信道矩阵的恢复精度,考虑指数函数:
其中是一个固定的常数,当时,函数是值域为(0,1]的单调递减函数,(29)的导数为当增加时,会减小,因此斜率单调递减,根据的变化趋势对步长进行调整,恢复信道向量的初始阶段选取一个较大的的长度用来提高信道恢复的效率;第二阶段逐步调整步长,的长度为其中代表迭代阶段,为向上取整,因此估计稀疏度与实际的稀疏度越接近时,长度的改变随着的增加减小,直到记为算法,具体步骤如下:输入:观测矩阵维观测信号初始步长输出:稀疏信号的估计
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910822199.9A CN110519189B (zh) | 2019-08-30 | 2019-08-30 | 高度移动场景下基于压缩感知的毫米波信道估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910822199.9A CN110519189B (zh) | 2019-08-30 | 2019-08-30 | 高度移动场景下基于压缩感知的毫米波信道估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110519189A CN110519189A (zh) | 2019-11-29 |
CN110519189B true CN110519189B (zh) | 2022-12-09 |
Family
ID=68630298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910822199.9A Active CN110519189B (zh) | 2019-08-30 | 2019-08-30 | 高度移动场景下基于压缩感知的毫米波信道估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110519189B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110943768B (zh) * | 2019-12-09 | 2021-09-28 | 中南大学 | 一种毫米波大规模mimo系统的混合预编码码本联合设计方法 |
CN110928328B (zh) * | 2019-12-17 | 2020-08-28 | 中国人民解放军火箭军工程大学 | 一种波束倾斜效应下无人机系统时变信道跟踪方法 |
CN112202481A (zh) * | 2020-06-24 | 2021-01-08 | 北京邮电大学 | 一种基于自适应感知矩阵的压缩感知信道估计算法及实现装置 |
CN114079491A (zh) * | 2020-08-21 | 2022-02-22 | 华为技术有限公司 | 数据传输方法及装置、无线通信系统、存储介质 |
EP3962006A1 (en) * | 2020-09-01 | 2022-03-02 | Vestel Elektronik Sanayi ve Ticaret A.S. | Channel estimation for configurable surfaces |
CN112769461B (zh) * | 2020-12-11 | 2022-04-22 | 华南理工大学 | 基于毫米波智能反射面通信的大规模天线信道估计方法 |
CN112671435B (zh) * | 2020-12-22 | 2022-07-22 | 杭州电子科技大学 | 一种fdd大规模天线中基于字典迭代的信道估计方法 |
CN112737649B (zh) * | 2020-12-25 | 2022-04-01 | 杭州电子科技大学 | 一种基于角度网格优化与范数约束的毫米波信道估计方法 |
CN113630350B (zh) * | 2021-08-13 | 2022-08-26 | 北京邮电大学 | 一种毫米波多天线信道估计方法 |
CN114422302B (zh) * | 2021-11-29 | 2024-03-22 | 新华三技术有限公司成都分公司 | 一种信道估计方法、装置及设备 |
CN115987739A (zh) * | 2023-01-30 | 2023-04-18 | 东南大学 | 一种基于双时间尺度的通信感知一体化信号处理方法 |
CN116360603A (zh) * | 2023-05-29 | 2023-06-30 | 中数元宇数字科技(上海)有限公司 | 基于时序信号匹配的交互方法、设备、介质及程序产品 |
CN117278367B (zh) * | 2023-11-23 | 2024-01-23 | 北京中关村实验室 | 一种分布式压缩感知稀疏时变信道估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106059731A (zh) * | 2016-05-19 | 2016-10-26 | 重庆大学 | 一种适用于快时变稀疏估计的最优导频图样的设计方法 |
CN108322409A (zh) * | 2018-01-25 | 2018-07-24 | 杭州电子科技大学 | 基于广义正交匹配追踪算法的稀疏ofdm信道估计方法 |
EP3406059A1 (en) * | 2017-03-31 | 2018-11-28 | Mitsubishi Electric Corporation | Method for decoding symbol and receiver for receiving and decoding symbol |
CN109560841A (zh) * | 2018-12-13 | 2019-04-02 | 东北大学 | 基于改进的分布式压缩感知算法的大规模mimo系统信道估计方法 |
-
2019
- 2019-08-30 CN CN201910822199.9A patent/CN110519189B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106059731A (zh) * | 2016-05-19 | 2016-10-26 | 重庆大学 | 一种适用于快时变稀疏估计的最优导频图样的设计方法 |
EP3406059A1 (en) * | 2017-03-31 | 2018-11-28 | Mitsubishi Electric Corporation | Method for decoding symbol and receiver for receiving and decoding symbol |
CN108322409A (zh) * | 2018-01-25 | 2018-07-24 | 杭州电子科技大学 | 基于广义正交匹配追踪算法的稀疏ofdm信道估计方法 |
CN109560841A (zh) * | 2018-12-13 | 2019-04-02 | 东北大学 | 基于改进的分布式压缩感知算法的大规模mimo系统信道估计方法 |
Non-Patent Citations (2)
Title |
---|
Compressed Channel Sensing:A New Approach to Estimating Sparse Multipath Channels;Waheed U. Bajwa 等;《Proceedings of the IEEE》;20100426;第98卷(第6期);第1058-1076页 * |
基于自适应压缩感知的OFDM 稀疏信道估计研究;费洪涛 等;《南京邮电大学学报(自然科学版)》;20190430;第39卷(第2期);第49-54页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110519189A (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110519189B (zh) | 高度移动场景下基于压缩感知的毫米波信道估计方法 | |
CN110463147B (zh) | 解码符号的方法以及接收并解码符号的接收器 | |
US10298301B2 (en) | Channel estimation for millimeter-wave communication/data link and the corresponding codebook design | |
CN107135024B (zh) | 一种低复杂度的混合波束赋形迭代设计方法 | |
CN104935366B (zh) | 毫米波通信中的波束搜索方法 | |
CN108650200B (zh) | 高低频混合组网系统的低频辅助信道估计方法 | |
CN107566305B (zh) | 一种低复杂度的毫米波系统信道估计方法 | |
CN109379122B (zh) | 一种毫米波通信多径信道动态波束训练方法 | |
Heng et al. | Learning site-specific probing beams for fast mmWave beam alignment | |
CN106330284B (zh) | 一种低复杂度大规模mimo信道估计方法 | |
KR102228091B1 (ko) | 밀리미터파 기반 대규모 mimo 시스템의 하이브리드 빔포밍 장치 및 방법 | |
CN112468202B (zh) | 低复杂度毫米波大规模mimo混合预编码方法 | |
CN114567525B (zh) | 一种信道估计方法及装置 | |
CN112953865A (zh) | 一种用于大规模多输入多输出系统的信道估计方法 | |
CN110719127B (zh) | 一种具有恒模约束的毫米波mimo系统波束成形方法 | |
CN109861933A (zh) | 一种基于music算法和预编码的毫米波mimo信道估计方法 | |
CN117544205A (zh) | 超大规模mimo系统的近场全维非正交码本设计方法 | |
Li et al. | DeFe: Indoor localization based on channel state information feature using deep learning | |
Lu et al. | An adaptive channel estimation algorithm for millimeter wave cellular systems | |
Shaham et al. | Raf: Robust adaptive multi-feedback channel estimation for millimeter wave mimo systems | |
CN112910524A (zh) | 基于概率模型的毫米波波束跟踪方法 | |
You et al. | Exploiting angular spread in channel estimation of millimeter wave MIMO system | |
CN114884776B (zh) | 信道估计方法、装置、电子设备及存储介质 | |
Kim et al. | Channel estimation for reconfigurable intelligent surface-aided mmWave communications | |
CN105681236B (zh) | 一种球形译码算法的初值选取方法 |
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 |