首页 / 专利库 / 电脑编程 / 算法 / 期望最大化算法 / 具有透镜天线阵的MIMO系统的稀疏信道估计方法

具有透镜天线阵的MIMO系统的稀疏信道估计方法

阅读:260发布:2020-05-18

专利汇可以提供具有透镜天线阵的MIMO系统的稀疏信道估计方法专利检索,专利查询,专利分析的服务。并且本 发明 属于无线通信技术领域,涉及一种具有透镜天线阵的MIMO系统的稀疏信道估计方法。本发明提出了一种具有透镜天线阵列的毫米波大规模MIMO系统进行信道估计方法。基于波束空间中的信道稀疏性,首先将信道估计放到稀疏 信号 恢复的 框架 中,然后利用期望最大化 算法 来学习该框架中的相关参数,利用估计的参数,可以通过简单的替换获得信道估计。相较基于支持检测、SCAMPI的方案,本发明所提出的方案可以充分有效地利用信道响应的特性,因此可以实现更好的性能。,下面是具有透镜天线阵的MIMO系统的稀疏信道估计方法专利的具体信息内容。

1.具有透镜天线阵的MIMO系统的稀疏信道估计方法,所述MIMO系统的基站配备有三维透镜天线阵列,透镜长度和高度分别为Dy和Dz,并将M×N天线阵列置于透镜焦平面上,MN根天线通过Q×MN的选择网络W与Q个射频链相连,即基站利用网络 检测接收信号y;
其特征在于,所述信道估计方法包括以下步骤:
S1、将信道估计问题放在稀疏信号恢复的框架下;令入射平面波方位和入射平面波仰角 在[-1,1]均匀采样P个网格点S={-1+1/P,-1+3/P,...,(P-1)/P},信道响应矢量记为:
其中,s1、s2表示S中离 最近的格子点,
δy、δz表示 的离散化误差;d(s1,s2)=v(s1,s2),dy(s1,s2)和dz(s1,s2)分别是 对求导后在s1、s2处的值;构建3个MN×P2的稀疏使能字典D,Dy,Dz且D(:,p)=d(sp1,sp2),Dy(:,p)=dy(sp1,sp2),Dz(:,p)=dz(sp1,sp2),sn表示S的第n个元素,p=p1+(p2-1)P且1≤p1,p2≤P;获得信道表达式:
其中, 的离散化误差δ1和δ2分别为 且δy=[δy,1,δy,2,...,
δy,P]T, 且δz=[δz,1,δz,2,...,δz,P]T,IP表示P维单位列向量,L为多径数,diag(δ)代表对角元上为矢量δ的对角矩阵;
令 是离 最近的格子点, 表示S的第
(l')个元素, 表示δy(δz)的第l个分支,且x的第p个分支xp=α(l), 由
于信道只有L条多径,δy、δz、x只有L个非零元,若(sp1,sp2)不是离 最近的,则δy,p1=
0、δz,p2=0、xp=0;此时记接收信号为:
测量矩阵 n为复高斯噪声矢量;
S2、通过以下的步骤来实现期望最大化算法迭代
S21、计算信号x的后验分布均值和方差
由于L<<P2,信号x具有稀疏性,令信号x中每个元素服从均值为0,方差为 的相同复高斯分布,且α=[α1,α2,…,αP2]T;噪声n中每个元素服从均值为0,方差为β-1的相同复高斯分布;信号x的先验概率和条件概率p(r|x,β)分别为:
其中|·|和||·||表示分别表示元素的模和向量的二范数;基于期望最大化算法,信号x的后验概率为:
H
其中,<·>q(n)表示相对于高斯噪声后验分布q(n)的期望,(·) 表示共轭转置,const为常数项,因此后验均值μ=βτΦHr,协方差矩阵τ=(βΦHΦ+D)-1,D=diag(α)表示对角线上为矢量α的对角阵;
S22、更新参数α、β、δy、δz
通过最大化完全似然函数F(α,β,δy,δz)=q(x)来实现参数的更新:
对完全似然函数求导并令其一阶导数为0,得到参数αl、β的更新公式:
其中,τ(l,l)表示协方差矩阵τ的第(l,l)个分支,tr(·)表示矩阵的迹,相对参数δy、δz最大化完全似然函数也就是最小化函数F'(δy,δz)=<||r-Φx||2>,经过一系列操作,只保留与δy、δz相关的项,最终得到:
其中,
* T
且 (·) 、(·) 分别代表共轭和转置, 表
示对变量取实部,⊙表示Hadamard积,上述E、F、G、U、V、A、B、C只是作为中间变量方便公式的表示,没有物理含义,将 带入公式(15)并分别对δy、δz求导,得到参
数δy、δz的更新公式:
其中,矩阵E1中的元素满足E1(i,j)=∑E(i:P:P2,j:P:P2),矩阵F1中的元素满足F1(i,j)=∑F((i-1)*P+1:i*P,(j-1)*P+1:j*P),矩阵G1中的元素满足G1(i,j)=∑G(i:P:P2,(j-
1)*P+1:j*P),矩阵G2满足 列矢量U1中的元素满足U1(i)=∑U(i:P:P2),列矢量V1中的元素满足V1(i)=∑V((i-1)*P+1:i*P);用M(i,j)来表示矩阵M的第(i,j)个分支,且(i:
k:j)表示以间隔k从第i个到第j个点的取值;(i:j)表示以单位间隔从第i个到第j个点的取值;由于δy∈[-1/P,1/P]、δz∈[-1/P,1/P],因此对于δy、δz中的任意一个元素δi有以下限制:
当-1/P≤δi≤1/P时,δi=δi;当δi<-1/P时,δi=-1/P;当δi>1/P时,δi=1/P;基于式(16)、(17),最后再次利用 得到δ1、δ2的估计值;
S23、迭代步骤S21、S22直到循环结束条件:迭代次数大于等于N或
其中N表示最大迭代次数,μn代表第n次迭代的后验均值;
S3、信道估计值的获取
根据参数δ1、δ2的估计值和信号x的后验均值,基于公式(7),最终可以得到信道h的估计值hest。

说明书全文

具有透镜天线阵的MIMO系统的稀疏信道估计方法

技术领域

[0001] 本发明属于无线通信技术领域,涉及一种具有透镜天线阵的MIMO系统的稀疏信道估计方法。

背景技术

[0002] 由于巨大的频谱带宽,毫米波通信可以提供Gbps的高速数据速率,因此被认为是未来蜂窝系统的关键技术之一。但是,毫米波上的穿透损耗大于低频波。因此需要制作一个具有高定向发射/接收功能的大型天线阵以补偿高穿透损耗,毫米波信号的短波长使得封装大量具有小形状因子的天线变得可行,从而可以应用大规模MIMO技术来补偿毫米波信道中严重的传播损耗。然而,在硬件成本和功耗方面实际上难以实现为每个天线配备一个专用射频链。为了减少射频链的数量,已经研究了模拟波束形成、混合模拟和数字波束成形以及透镜天线阵列等几种方法。
[0003] 透镜天线阵列相较普通天线阵列有许多优点,特别是透镜可以在非常短的脉冲长度下工作,并可以扫描比任何先前已知设备更宽的波束宽度。此外,透镜能够形成低旁瓣光束。鉴于这些优点,透镜天线阵列作为发射/接收前端引入毫米波通信系统。透镜天线阵列包括两个部分:电磁透镜和位于透镜焦平面的天线阵列。电磁透镜能够让入射路径相移,使得路径能量聚焦在少量天线元件上。由于透镜的度相关能量聚焦能,可以将具有独立角度的不同多径的能量聚焦到不同的天线元件子集。因此,毫米波信道可以由波束空间很好地表征,其中一个波束对应于一个可能的多径角。由于多径的数量是有限的,因此在波束空间中毫米波信道是稀疏的,这也使得几种信道估计方案相继提出,包括基于支持检测、SCAMPI等信道估计算法

发明内容

[0004] 为了进一步提高信道估计性能,本发明提出了一种具有透镜天线阵列的毫米波大规模MIMO系统进行信道估计方法。基于波束空间中的信道稀疏性,首先将信道估计放到稀疏信号恢复的框架中,然后利用期望最大化算法来学习该框架中的相关参数,利用估计的参数,可以通过简单的替换获得信道估计。相较基于支持检测、SCAMPI的方案,本发明所提出的方案可以充分有效地利用信道响应的特性,因此可以实现更好的性能
[0005] 本发明采用了期望最大化算法,期望最大化算法是一种求解未知随机变量的后验分布的算法,通过不断地迭代,得到样本已知的条件下的隐变量的均值。
[0006] 为了便于本领域内技术人员对本发明技术方案的理解,首先对本发明采用的系统模型进行说明。
[0007] 考虑单天线用户向配备有三维(3D)透镜天线阵列的基站(BS)发送信号,基站处的EM透镜长度和高度分别为Dy和Dz,并将M×N天线阵列置于透镜焦平面上,MN根天线通过Q×MN的选择网络W与Q个射频链相连。记天线的索引为(m,n),于是有第(m,n)根天线处的阵列响应为:
[0008]
[0009] 其中, λ、φy、φz分别代表入射平面波波长、入射平面波方位角、入射平面波仰角,且天线数与透镜尺寸满足 且 且
[0010] 记 为第l条多径的虚拟角度,相应整个天线阵列的信道响应矩阵可记为假设共有L条多径,3D透镜天线阵列的信道矩阵为:
[0011]
[0012] 其中,α(l)是第l条路径的复增益,且服从零均值单位方差的复高斯分布。通过矢量化H和A,式(2)另记为:
[0013]
[0014] 表示 的矢量化结果。
[0015] 在上行链路信道估计中,用户向基站发送的训练符号s,此时基站的接收信号为:
[0016]
[0017] 其中, 为服从 分布的复高斯噪声矢量。信噪比定义为不失一般性,假设s=1,因此可以省略符号s以简化表示。
[0018] 由于射频链的数量有限,BS不能直接观察信号 假设有K个射频链,基站利用选择网络 来检测 因此得到 的每一个分支由1位移相器实现。对于q个训练符号,总的接收信号为:
[0019]
[0020] 其中, 且Q=Kq, 且Iq为长度为q的单位矢量,复高斯噪声矢量表示Kronecker积。
[0021] 本发明通过如下步骤实现:
[0022] S1、首先将信道估计问题放在稀疏信号恢复的框架下。令到达角 在[-1,1]均匀采样P个网格点S={-1+1/P,-1+3/P,...,(P-1)/P},信道响应矢量可记为:
[0023]
[0024] 其中,s1、s2表示S中离 最近的格子点,δy、δz表示 的离散化误差。d(s1,s2)=
v(s1,s2),dy(s1,s2)和dz(s1,s2)分别是 对 求导后在s1、s2处的值。基于此可以构建3个MN×P2的稀疏使能字典D,Dy,Dz且D(:,p)=d(sp1,sp2),Dy(:,p)=dy(sp1,sp2),Dz(:,p)=dz(sp1,sp2),sn表示S的第n个元素,p=p1+(p2-1)P且1≤p1,p2≤P。最终(3)中的信道可另记为:
[0025]
[0026] 其中, 的离散化误差δ1和δ2分别为 且δy=[δy,1,δy,2,...,δy,P]T, 且δz=[δz,1,δz,2,...,δz,P]T,IP表示P维单位列向量,L为多径数,diag(δ)代表对角元上为矢量δ的对角矩阵。
[0027] 考虑L条多径,假设 是离 最近的格子点,表示S的第 个元素, 表示δy(δz)的第l个
分支,且x的第p个分支xp=α(l), 由于信道只有L条多径,δy、δz、x只有L个非零元,若(sp1,sp2)不是离 最近的,则δy,p1=0、δz,p2=0、xp=0。此时重记(5)中的接收信号为:
[0028]
[0029] 其中,测量矩阵
[0030] S2、通过以下的步骤来实现期望最大化算法的迭代:
[0031] S21、计算信号x的后验分布均值和方差
[0032] 由于L<<P2,信号x具有稀疏性,假设信号x中每个元素服从均值为0,方差为的相同复高斯分布且下标l表示x的第l个元素,令x中所有元素方差集合为噪声n中每个元素服从均值为0,方差为β-1的相同复高斯分布。因此,信
号x的先验概率和条件概率p(r|x,β)分别为:
[0033]
[0034]
[0035] 其中,|·|和||·||表示分别表示元素的模和向量的二范数,基于期望最大化算法,信号x的后验概率为:
[0036]
[0037] 其中<·>q(n)表示相对于高斯噪声后验分布q(n)的期望,(·)H表示共轭转置,因此H H -1后验均值μ=βτΦr,协方差矩阵τ=(βΦ Φ+D) ,D=diag(α)表示对角线上为矢量α的对角阵。
[0038] S22、更新参数α、β、δy、δz
[0039] 通过最大化完全似然函数F(α,β,δy,δz)=q(x)来实现参数的更新,<·>q(x)表示相对信号x后验分布q(x)的期望:
[0040]
[0041] 对完全似然函数求导并令其一阶导数为0,得到参数αl、β的更新公式:
[0042]
[0043]
[0044] 其中,τ(l,l)表示协方差矩阵τ的第(l,l)个分支,tr(·)表示矩阵的迹,相对参数δy、δz最大化完全似然函数也就是最小化函数F'(δy,δz)=<||r-Φx||2>q(x),经过一系列操作,只保留与δy、δz相关的项,最终得到:
[0045]
[0046] 其中,且 (·)*、(·)T分别代表共轭和转置, 表
示对变量取实部,⊙表示Hadamard积,上述E、F、G、U、V、A、B、C只是作为中间变量方便公式的表示,没有特殊含义。将 带入公式(15)并分别对δy、δz求导,得到参数
δy、δz的更新公式:
[0047]
[0048]
[0049] 其中,矩阵E1中的元素满足E1(i,j)=∑E(i:P:P2,j:P:P2),矩阵F1中的元素满足F1(i,j)=∑F((i-1)*P+1:i*P,(j-1)*P+1:j*P),矩阵G1中的元素满足G1(i,j)=∑G(i:P:P2,(j-1)*P+1:j*P),矩阵G2满足 列矢量U1中的元素满足U1(i)=∑U(i:P:P2),列矢量V1中的元素满足V1(i)=∑V((i-1)*P+1:i*P)。这里用M(i,j)来表示矩阵M的第(i,j)个分支,且(i:k:j)表示以间隔k从第i个到第j个点的取值;(i:j)表示以单位间隔从第i个到第j个点的取值,∑(·)表示求和运算。由于δy∈[-1/P,1/P]、δz∈[-1/P,1/P],因此对于δy、δz中的任意一个元素δi有以下限制:当-1/P≤δi≤1/P时,δi=δi;当δi<-1/P时,δi=-1/P;当δi>1/P时,δi=1/P。基于式(16)、(17),最后再次利用 得到δ1、δ2的估计值。
[0050] S23、迭代步骤S21、S22直到循环结束条件:迭代次数大于等于N或其中N表示最大迭代次数,μn代表第n次迭代的后验均值。
[0051] S3、信道估计值的获取
[0052] 根据参数δ1、δ2的估计值和信号x的后验均值,基于公式(7),最终可以得到信道h的估计值hest。
[0053] 本发明的有益效果为能够对具有透镜天线阵列的毫米波大规模MIMO系统进行有效的信道估计且估计性能相较之前提出的方法有较大提升。附图说明
[0054] 图1是本发明的算法与之前相关文献所提算法的NMSE性能曲线比较;
[0055] 图2是在移相器的不同断开比率下本发明的算法的NMSE性能曲线图。

具体实施方式

[0056] 下面结合附图和仿真示例说明本发明的有效性:
[0057] 仿真基于一个基站配有41×41透镜天线阵列的毫米波系统进行,相应的透镜尺寸为 其中波长λ=1,多径数L为3,同时为了保证字典D的完备性,将P取作41。以下使用归一化均方误差即 来评估每个算法的性能,其中hest表示估
计信道,h表示实际生成的信道,这里考虑100个信道。
[0058] 图1表示在相同的仿真条件下,SD算法、服从均匀分布的SCAMPI算法、本发明所提出的算法在-20至10dB信噪比时的NMSE曲线。通过比较各算法的仿真性能可以发现,本发明所提出算法具有最好的信道估计性能,且随着信噪比增大它的性能不断提升。
[0059] 图2表示当W中的一些移相器断开以节省功率时所提出的方案的性能。令p为总移相器中断开的移相器比率,从图2中可以看出,本发明所提出算法在总的移相器减少50%时的性能接近SCAMPI且10dB处的NMSE增加仅约4×10-3,即该算法对移相器的断开比率具有鲁棒性,相较于网络功耗的降低,性能的下降可忽略。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈