本发明涉及移动通信中信道估计领域,尤其涉及毫米波大规模mimo系统在混合模-数预编码架构下基于超分辨率角度和时延估计的宽带信道估计。
背景技术
随着移动流量业务呈现出爆炸式增长的趋势,可预见当前基于lte/lte-a的4g移动蜂窝网络将愈加难以满足这种几何式增长的容量需求。毫米波大规模mimo(多输入多输出,multiple-inputmultiple-output)技术被广泛认为是实现通信网络千倍容量提升这一宏伟愿景的物理层关键技术之一。首先,毫米波频段(30-300ghz)具有丰富的频谱资源,能有效地缓解当前频谱资源短缺的问题。其次,毫米波频段波长较短,更有利于大规模天线阵列的紧凑封装,同时,大规模天线阵列又可为毫米波通信提供足够大的阵列增益。最后,毫米波信号的高路损、易遮挡的传输特性能有效避免小区间的干扰,以便形成超密集组网。为了充分挖掘和利用毫米波大规模mimo技术以上诸多优势,可靠的信道估计是一个先决条件。然而,在毫米波大规模mimo系统中,为了降低系统所需的硬件成本和功耗,收发机通常采用混合模-数预编码架构,这将会使得该系统中的信道估计问题研究更具有挑战性。
目前,很多关于毫米波混合大规模mimo系统中信道估计研究的重点都放在了窄带信道中。但是,在实际的情况中,毫米波大规模mimo的系统带宽往往能高达数ghz,这时,系统在进行信道估计时必将面临宽带与频率选择性的问题。同时,由于毫米波信道中信号在非直射径情形下的路径损耗非常严重,使得接收机只能接收到对应于有限几个显著散射体的多径信号,此时的无线信道呈现出稀疏性。如何同时准确地估计出宽带信道中稀疏多径成分所对应的角度以及多径时延是一个亟待解决的难题。
在毫米波混合大规模mimo系统中,目前已有多种信道估计方法被提出,用以解决宽带频率选择性衰落信道的估计问题。为了对抗宽带多径效应的影响,系统采用ofdm(正交频分复用,orthogonalfrequencydivisionmultiplexing)传输技术。具体地,venugopal等人提出了一种频域稀疏方法来估计频率选择性衰落信道,利用了omp(正交匹配追踪,orthogonalmatchingpursuit)算法来找到包含了信道增益的稀疏向量中的最优稀疏近似值,但是,该方法在计算时要进行与子载波数目相同次数的omp算法循环,这将导致相当高的计算复杂度。fernández等人根据系统带宽内毫米波信道的角度域稀疏性是不变的性质,通过利用不同子载波信道之间的空间共同稀疏性,提出了一种sw(同时加权,simultaneousweighted)-omp算法来降低计算复杂度。这些宽带信道估计方法基本上是利用毫米波mimo信道在角度域上的稀疏性,通过各种压缩感知理论来解决近优的稀疏重建问题,即估计稀疏向量的支撑集以及支撑集上元素所对应的值,以降低系统在信道估计时所需的导频开销。
然而,对于以上传统基于压缩感知理论的宽带信道估计方法,一方面,在信道估计问题表达成稀疏重建问题时,其中假定的稀疏向量并非是真正意义上完全稀疏的。另一方面,基于压缩感知理论的宽带信道估计方法均会对连续的角度量化为离散化的网格点,这样有限的量化角度网格分辨率以及码本大小,就不可避免地引入了量化误差,在高信噪比下,这种量化误差会愈发明显。
技术实现要素:
有鉴于此,本发明提供了一种基于超分辨率角度和时延估计的宽带信道估计方法,用以解决传统基于压缩感知的宽带信道估计方法中由有限的量化角度网格分辨率所造成的量化误差问题。
为了解决上述技术问题,本发明是这样实现的:
一种基于超分辨率角度和时延估计的宽带信道估计方法,包括:
步骤1、针对毫米波大规模mimo-ofdm系统在混合模-数预编码架构下的信道估计过程,设计了收发端的训练信号,从而得到保留有阵列响应移不变性的低维等效信道;
所述训练信号包括发射端的组合预编码器
其中,
第j个混合预编码器fj由模拟预编码器frf,j和数字预编码器fbb,j构成,即fj=frf,jfbb,j,第i个混合合并器wi由模拟合并器wrf,i和数字合并器wbb,i构成,即wi=wrf,iwbb,i;
将酉矩阵
发射端的第j个模拟预编码器
接收端的第i个模拟合并器
以上设计好的
步骤2、利用步骤1设计好的训练信号,在接收端获得第k个子载波对应的低维等效信道矩阵
步骤3、利用步骤2估计到的角度和时延信息,结合等效信道向量
优选地,所述步骤1中,模拟预编码器
其中,nq为设定的角度量化比特值;当nq=1时,所述酉矩阵
优选地,所述步骤2中,低维等效信道矩阵到等效信道向量
定义第k个子载波对应的等效信道向量为:
其中,
同时考虑所有k个子载波的等效信道向量,并将其组合成一个矩阵
其中aτ=[τ[0],τ[1],…,τ[k-1]]t,
那么,再对矩阵
其中,
优选地,所述步骤3为:
步骤3.1、利用估计出的所述到达角、离开角以及多径时延的超分辨率估计值,重建出发射端和接收端对应的导向矢量矩阵
步骤3.2、结合设计好的训练信号以及等效信道向量
步骤3.3、利用所述到达角、离开角、多径时延以及路径复增益
有益效果:
(1)本发明通过在收发端设计出合适的训练信号,能大幅度地降低信道估计时所需的导频开销,得到保留有阵列响应移不变性的低维等效信道。
(2)本发明将系统的相移网络中移相器的相位值量化为有限的比特数,即利用哈达玛矩阵或者dft矩阵作为训练信号的取值集合,来设计混合预编码器和混合合并器,以降低实际系统的实现难度。
(3)本发明联合处理所有子载波对应的低维等效信道矩阵,变换为最终的等效信道向量,并对该向量利用经典空间谱估计理论中的三维酉esprit算法,可同时估计出宽带毫米波信道中稀疏多径成分对应的到达角、离开角以及多径时延的超分辨率估计值。
(4)本发明利用估计到的到达角、离开角和多径时延,同时结合设计好的训练信号以及最终的等效信道向量,可计算信道中路径复增益的最小二乘解,并重建出原始频域信道矩阵。
附图说明
图1为典型混合模-数预编码架构下毫米波大规模mimo-ofdm系统框图。
图2为本发明基于超分辨率角度和时延估计的宽带信道估计方法的流程图。
图3为不同信道估计方法在相同导频开销tpilot=75时随信噪比变化的nmse(归一化均方误差,normalizedmeansquareerror)性能评估对比图。
图4为不同信道估计方法在导频开销tpilot={48,108}时随信噪比变化的nmse误差性能评估对比图。
图5为不同信道估计方法在信噪比为{0,10}db下随路径数l变化的nmse性能评估对比图。
图6为不同信道估计方法随信噪比变化的平均频谱效率性能评估对比图。
具体实施方式
本发明的基本思想是利用毫米波信道的稀疏性,通过设计收发端的训练信号,能大幅度地降低信道估计时所需的导频开销,得到保留有阵列响应移不变性的低维等效信道,之后通过联合利用ofdm系统中多个子载波的信息来对接收信号进行处理,便可以通过经典空间谱估计中的三维酉esprit(借助旋转不变技术估计信号参数,estimatingsignalparametersviarotationalinvariancetechniques)算法来高精确地估计出宽带毫米波mimo信道中稀疏多径成分的到达角,离开角以及相应时延的超分辨率估计值,并以估计到的角度和时延信息重建出原始频域信道,从而显著提高信道估计的准确性。此外,本发明还将系统的相移网络中移相器的相位值量化为有限的比特数,以便降低实际系统的实现难度。
下面结合附图并举实施例,对本发明进行详细描述。
本发明考虑典型混合模-数预编码架构下毫米波混合大规模mimo-ofdm系统,如图1所示。在该系统中,发射端和接收端分别装备有nt和nr根天线,以及
这里(·)h为共轭转置符号,接收端的混合合并器w[k]=wrfwbb[k]是由模拟合并器wrf和基带数字合并器wbb[k]级联而成,h[k]是第k个子载波对应的毫米波频率选择性衰落信道,发射端的混合合并器f[k]=frffbb[k]则可看做是模拟预编码器frf和基带数字预编码器fbb[k]的级联。s[k]和n[k]分别是发送信号向量和复高斯白噪声。注意,因为这里的模拟预编码器frf和模拟合并器wrf均是由模拟相移网络实现的,故它们的每一项都需要满足一个恒定模值的约束条件,即对于frf和wrf的第(m,n)项,有
由于毫米波信号具有高路损、易遮挡的传输特性,使得毫米波系统中的通信链路往往是由有限多个显著的散射体产生的。因此,本发明考虑接收端和发射端之间只有l个不同的散射体对应l个主要稀疏多径成分的几何毫米波宽带频率选择性衰落信道模型,其中,时域的频率选择性衰落信道矩阵可表示为
这里δ(·)表示狄拉克δ函数。对于第l条路径而言,αl是服从复高斯分布
分别表示对应于接收端和发射端的空间频率,其中λ和d分别为波长和相邻天线的间隔,θr,l和
那么,变换到频域上,第k个子载波对应的频域信道矩阵则为
这里fs表示系统带宽,也是系统的采样速率。于是,公式(5)中的频域信道矩阵h[k]可进一步写成如下更紧凑的形式
其中ar=[a(μr,1),…,a(μr,l)]和at=[a(μt,1),…,a(μt,l)]分别为接收端和发射端对应的导向矢量矩阵。d[k]=diag(d[k])是一个对角矩阵,其对角元素组成的向量d[k]可分解为d[k]=diag(α)τ[k],且复增益向量为
为了能有效地利用三维酉esprit算法来估计宽带毫米波mimo信道,本发明考虑在信道估计阶段用时块来传输信号,且每个时块包含了ns个ofdm符号。同时,为了有助于后续对信道估计问题进行数学建模,这里不同子载波均使用相同的数字预编码器/合并器和导频信号,也就是说,对于
这里
其中对于
分别是需要设计的组合预编码器和组合合并器。块对角矩阵
这里
基于上述分析,下面结合图2分步骤对本发明的信道估计过程进行详细描述。
步骤1、设计训练信号并量化处理
对于毫米波混合大规模mimo系统,每个基带观测值是不同天线通过模拟相移网络混合后形成的,这就意味着相移网络破坏了阵列响应的移不变性,因此不能直接使用传统esprit之类的算法来估计毫米波信道中的到达角/离开角。为了解决该难题,本发明进行了以下步骤。
步骤1.1、设计收发端的训练信号
本发明需要在系统的收发端设计出合适的训练信号,也就是组合预编码器
这里αf和αw是为了保证
其中hm,n[k]代表信道
现在需要设计收发端的训练信号,也就是组合预编码器
分别由(j-1)ns和nt-jns个相同的
类似地,对于
步骤1.2、对相移网络中移相器的相位值量化处理
为了便于实际系统的实现,本发明将系统的相移网络中移相器的相位值量化为有限的比特数,即对上述设计好的模拟预编码器
其中,nq为设定的角度量化的比特值。为了防止相位量化过程影响到酉矩阵
其他的形式以此类推。而当
其他的形式也以此类推。公式(15)中的j表示数学中的虚数单位,即
步骤2、联合处理与三维酉esprit算法
接下来,本发明需要联合处理所有子载波对应的低维等效信道矩阵,并将其变换为最终的等效信道向量。之后,对该向量利用经典空间谱估计理论中的三维酉esprit算法,可同时估计出宽带毫米波信道中稀疏多径成分对应的到达角、离开角以及多径时延的超分辨率估计值。
步骤2.1、联合处理所有子载波的信号
由公式(10)以及设计好的组合预编码器
这里定义了
其中
这里定义
步骤2.2、三维酉esprit算法
为了简化后续的表达式,这里重新定义与接收端和发射端相关联的等价导向矢量矩阵
1、3d空间平滑处理
首先定义三个空间平滑参量mr,mt以及mτ,且它们分别满足
于是,通过这三个1d选择矩阵可以定义m=mrmtmτ个3d空间平滑选择矩阵,其中对于第(mr,mt,mτ)个3d空间平滑选择矩阵可表示为
那么,进行3d空间平滑后的信道矩阵为
2、实值处理
为了能在降低算法中计算复杂度的同时充分地利用平滑后的数据,可以通过前后向平均将3d空间平滑后的信道矩阵
这里矩阵πn表示维度为n×n的交换矩阵,且左实转换矩阵qn的定义如下
这里的j与公式(15)中的一样,表示数学中的虚数单位。
3、信号子空间近似
通过对实值矩阵
4、求解移不变等式
对于某一个非奇异矩阵t,可以得到以下三个实值的移不变等式
这里,θ,φ以及ω是三个对角矩阵,可分别表示为
而
且其中的三个3d选择矩阵jr,jt以及jτ定义为
那么,通过利用最小二乘估计或者是总体最小二乘估计来求解上面的三个实值的移不变等式为
5、利用ssd算法联合对角化
本发明考虑采取一种改进的实值schur分解算法,ssd(同时schur分解,simultaneousschurdecomposition)算法,来对tθt-1,tφt-1以及tωt-1这三个实值矩阵进行联合对角化(对于ssd算法,具体参见文献“译名:利用多个非对称矩阵的同时schur分解来实现多维谐波恢复的自动配对”,其作者、英文名称及出处为“haardtm,nossekja.simultaneousschurdecompositionofseveralnonsymmetricmatricestoachieveautomaticpairinginmultidimensionalharmonicretrievalproblems[j].ieeetransactionsonsignalprocessing,1998,46(1):161-169.”)。具体地,定义一个代价函数ψ(ξ)为
这里
步骤3、重建原始频域信道
本发明利用以上等效信道向量
步骤3.1、重建导向矢量矩阵
由公式(19)看出到达角,离开角以及多径时延与相应的空间频率呈一一对应的关系。那么,通过计算,将相应的空间频率的估计
利用公式(20)中得到的导向矢量,可以分别重建出对应于接收端的导向矢量矩阵
步骤3.2、计算复增益
已知设计好的训练信号
之后,再利用ls(最小二乘,leastsquare)估计器就可求得如下路径复增益α的ls解
步骤3.3、重建频域信道矩阵
根据以上获得的配对好的估计值,即
以上便是本发明公开的基于超分辨率角度和时延估计的宽带信道估计方法。
为了说明本发明与传统基于压缩感知理论的宽带信道估计方法(这里考虑基于omp算法和sw-omp算法的宽带信道估计方法作为对比方法)在提高信道估计性能上的优势,这里用图3~图6来说明本发明的效果。具体地,
图3对比了基于超分辨率角度和时延估计的宽带信道估计方法与基于omp算法和sw-omp算法的宽带信道估计方法在相同导频开销tpilot=75时随信噪比变化的nmse(归一化均方误差,normalizedmeansquareerror)性能。从图3可以看出,本发明中信道估计方法的nmse性能要显著地优于其余两种对比方法。这是因为本发明中的信道估计方法能高精确地获得到达角,离开角以及多径时延的超分辨率估计值,而相比之下,基于omp算法和sw-omp算法的这些压缩感知类信道估计方法会因为有限的码本大小和量化角度网格分辨率而在高信噪比时信道估计的性能趋于平缓。
图4比较了图3中所述三种信道估计方法在不同导频开销tpilot={48,108}下随信噪比变化的nmse性能。从图4中看出,当导频开销从tpilot=48增大到tpilot=108时,本发明中信道估计方法的nmse性能会有较为显著的提升,对应的nmse性能曲线呈两条近乎平行的直线,相差大约8db左右。这是因为导频开销越多,接收端能获得的等效信道矩阵的维度也会越大,从而使得信道估计的nmse性能越好。相比之下,两种对比方法在导频开销增大时提升的nmse性能十分有限,大约在2~3db左右,且它们在低导频开销下会出现明显的nmse性能平台。
图5对比了所述三种信道估计方法在不同信噪比{0,10}db下随路径数变化的nmse性能。从图5可以看出,随着路径数变多,所述三种信道估计方法的nmse性能均会有一定程度的下降。然而,相比于两种对比方法,本发明中信道估计方法的nmse性能曲线更陡峭一些,也就是说,基于超分辨率角度和时延估计的宽带信道估计方法在路径数越少的情况下,到达角/离开角的分辨能力就会越强,其nmse性能也就越好。此外,从图5中不同信噪比0db和10db的nmse性能曲线可以看出,本发明中信道估计方法的nmse性能的间距更大,也就说明了本发明在高信噪比下能获得更好的信道估计性能。
图6对比了所述三种信道估计方法随信噪比变化的平均频谱效率性能,这里考虑以收发端均已知完美信道状态信息情况下的最优性能作为性能比较的上界。从图6可以看出,本发明中信道估计方法的平均频谱效率性能要优于其他两种对比方法,且当信噪比大于-10db时,其性能已经趋近于最优的性能上界。而两种对比方法与最优的性能上界之间还存在着一定的性能差距。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。