首页 / 专利库 / 信号处理 / 信号 / 导频信号 / 一种半盲信道估计方法及系统

一种半盲信道估计方法及系统

阅读:478发布:2020-05-08

专利汇可以提供一种半盲信道估计方法及系统专利检索,专利查询,专利分析的服务。并且本 发明 所属无线通信系统技术领域,特别涉及一种半盲信道估计方法及系统。所述半盲信道估计方法包括:步骤a:在第一个时间 块 中,采用稀疏 信号 重建 算法 建立基于训练 导频信号 的稀疏信道重构模型,通过所述稀疏信道重构模型估计所述第一个时间块的信道信息;步骤b:在第t个时间块中,基于第t‑1个时间块的信道信息,采用最大似然信号检测器检测所述第t个时间块的传输信号;其中,2≤t≤T,T为所述时间块总数;步骤c:基于所述第t个时间块的传输信号,采用半盲信道估计建立半盲稀疏信道重构模型,通过所述半盲稀疏信道重构模型估计所述第t个时间块的信道信息。本发明能够在连续域 角 度空间域高度紧密逼近信道参数的稀疏性,能够大幅度提高系统的 频谱 效率。,下面是一种半盲信道估计方法及系统专利的具体信息内容。

1.一种半盲信道估计方法,其特征在于,包括以下步骤:
步骤a:在第一个时间中,采用无格点压缩感知技术建立基于训练导频信号的稀疏信道重构模型,通过所述稀疏信道重构模型估计所述第一个时间块的信道信息;
步骤b:在第t个时间块中,基于第t-1个时间块的信道信息,采用最大似然信号检测器检测所述第t个时间块的传输信号;其中,2≤t≤T,T为所述时间块总数;
步骤c:基于所述第t个时间块的传输信号,建立具有解码误差纠正功能的半盲稀疏信道重构模型,通过所述半盲稀疏信道重构模型估计所述第t个时间块的信道信息;
步骤d:迭代执行步骤b和步骤c,直到完成T个时间块的信道信息估计;
其中,所述在第一个时间块中,采用无格点压缩感知技术建立基于训练导频信号的稀疏信道重构模型,具体包括:
假设总共T个时间块,在第一个时间块的K个时隙中,基于发射机发射的训练导频信号S=[s1,...,sK],采用基于无格点压缩感知的稀疏信号重建算法建立连续度空间域的基于训练导频信号的稀疏信道重构模型,通过该稀疏信道重构模型估计第一个时间块的信道信息
其中,连续角度空间域的基于训练导频信号的稀疏信道重构模型由以下过程推导而得:
第一个时间块的第k个时隙,发射端发射训练导频信号sk,则接收端接收到的信号为:
y1,k=W1HH1F1sk+W1Hn1,k
y1,k为第一个时间块的第k个时隙接收的信号,F1和W1分别表示发射端的第一预编码矩阵和接收端的第一测量矩阵,S=[s1,...,sK]表示K个时隙的拼接训练导频信号,n1,k是第一个时间块中第k个时隙接收机的加性白高斯噪声且 其中
表示Nr维的全1列向量;
假设接收机在第一个时间块的第k个时隙中保持F1和W1不变,则拼接的第一接收信号矩阵为:
Y1=W1HH1F1S+Q1
Y1为第一个时间块接收的信号,Q1=[W1Hn1,1,...,W1Hn1,K]是第一个时间块接收端的噪声矩阵;
将Y1按列拉伸为新的列向量形式为:
至此,可得连续角度空间域的基于训练导频信号的稀疏信道重构模型为:
表示估计的信道信息 为权重因子,经过按列拉伸
的逆运算由 得到第一个时间块的信道信息
所述基于所述第t个时间块的传输信号,建立具有解码误差纠正功能的半盲稀疏信道重构模型,具体包括:
建立解码误差补偿机制,在第t个时间块的K个时隙中,基于检测的传输信号 采用基于解码误差稀疏逼近机制的半盲信道估计技术建立连续角度空间域的传输数据辅助的半盲稀疏信道重构模型,并估计当前时间块的信道信息
其中,连续角度空间域的传输数据辅助的半盲稀疏信道重构模型由以下过程推导而得:
第t个时间块中,其中2≤t≤T,将Yt按列拉伸为列向量形式为:
考虑检测的数据信号存在误差,则定义et,k为:
则接收端拼接的第二接收信号矩阵Yt可重新定义为:
是拼接的接收端检测的传输信号矩阵,Et=[et,1,...,et,K]是拼
接的由于解码误差造成的接收信号误差矩阵,得到:
其中,
采用 范数对信号检测误差进行稀疏重建,连续角度空间域的传输数据辅助的半盲稀疏信道重构模型 其中2≤t≤T,为:
上述公式中, 表示估计的 Et=[et,1,...,et,K]是拼接的由于
解码误差造成的接收信号误差矩阵, 为权重因子,||·||1为 范数
操作;经过按列拉伸的逆运算,由 得到所述第t个时间块的信道信息
2.根据权利要求1所述的半盲信道估计方法,其特征在于,所述步骤a前还包括:将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式,具体为:
将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式
毫米波参数信道信息模型为:
其中,Ht为第t个时间块的信道信息,Lt表示第t个时间块的路径数,αt,l是第t个时间块第l条路径的增益, 和 分别表示第t个时间块第l条路径的Nt
个发射天线阵元和Nr个接收天线阵元的均匀线性天线阵列响应,其中,
其中, 和 分别为第t个时间块中第l条路径的离开角和到达角,Nt和Nr分别表示发射机和接收机上的天线阵元数,λ表示毫米波的波长,d是相邻天线阵元间距;
进一步将毫米波参数信道信息模型的公式转化成紧凑的矩阵形式:
其中, 是Lt条路径的增益向量,ht,j表示矩阵Ht的第j列,j=
1,...,Nt;
收发天线阵列响应矩阵定义为:
信道信息的二维原子范式由以下过程推导而得:
将信道信息Ht进行按列拉伸运算,得到信道信息新的列向量表达形式
其中, 为一个二
维复正弦原子,且满足 的限制,e和 分别表示Khatri-Rao积和克罗内克积,vec(·)表示矩阵的按列拉伸成列向量运算;
定义A={c(fl)||fl∈[0,1)×[0,1)}是第l条路径归一化的二维复正弦原子的集合,可得信道信息的二维原子范式表示为:
其中, 为第t个时间块的信道信息的二维原子范式表示,inf为取下确界操作,表示复数域,|·|表示取模操作。
3.一种半盲信道估计系统,其特征在于,包括:
第一信道信息估计模块:用于在第一个时间块中,采用稀疏信号重建算法建立基于训练导频信号的稀疏信道重构模型,通过所述稀疏信道重构模型估计所述第一个时间块的信道信息;
传输信号估计模块:用于在第t个时间块中,基于第t-1个时间块的信道信息,采用最大似然信号检测器检测所述第t个时间块的传输信号;其中,2≤t≤T,T为所述时间块总数;
第二信道信息估计模块:用于基于所述第t个时间块的传输信号,采用半盲信道估计建立半盲稀疏信道重构模型,通过所述半盲稀疏信道重构模型估计所述第t个时间块的信道信息;
迭代判断模块:用于判断是否完成T个时间块的信道信息估计,如果没有完成T时间块的信道信息估计,通过所述传输信号估计模块和第二信道信息估计模块进行迭代,直到完成T个时间块的信道信息估计;
其中,所述在第一个时间块中,采用无格点压缩感知技术建立基于训练导频信号的稀疏信道重构模型,具体包括:
假设总共T个时间块,在第一个时间块的K个时隙中,基于发射机发射的训练导频信号S=[s1,...,sK],采用基于无格点压缩感知的稀疏信号重建算法建立连续角度空间域的基于训练导频信号的稀疏信道重构模型,通过该稀疏信道重构模型估计第一个时间块的信道信息
其中,连续角度空间域的基于训练导频信号的稀疏信道重构模型由以下过程推导而得:
第一个时间块的第k个时隙,发射端发射训练导频信号sk,则接收端接收到的信号为:
y1,k=W1HH1F1sk+W1Hn1,k
y1,k为第一个时间块的第k个时隙接收的信号,F1和W1分别表示发射端的第一预编码矩阵和接收端的第一测量矩阵,S=[s1,...,sK]表示K个时隙的拼接训练导频信号,n1,k是第一个时间块中第k个时隙接收机的加性白高斯噪声且 其中
表示Nr维的全1列向量;
假设接收机在第一个时间块的第k个时隙中保持F1和W1不变,则拼接的第一接收信号矩阵为:
Y1=W1HH1F1S+Q1
Y1为第一个时间块接收的信号,Q1=[W1Hn1,1,...,W1Hn1,K]是第一个时间块接收端的噪声矩阵;
将Y1按列拉伸为新的列向量形式为:
至此,可得连续角度空间域的基于训练导频信号的稀疏信道重构模型为:
表示估计的信道信息 为权重因子,经过按列拉伸
的逆运算由 得到第一个时间块的信道信息
所述基于所述第t个时间块的传输信号,建立具有解码误差纠正功能的半盲稀疏信道重构模型,具体包括:
建立解码误差补偿机制,在第t个时间块的K个时隙中,基于检测的传输信号 采用基于解码误差稀疏逼近机制的半盲信道估计技术建立连续角度空间域的传输数据辅助的半盲稀疏信道重构模型,并估计当前时间块的信道信息
其中,连续角度空间域的传输数据辅助的半盲稀疏信道重构模型由以下过程推导而得:
第t个时间块中,其中2≤t≤T,将Yt按列拉伸为列向量形式为:
考虑检测的数据信号存在误差,则定义et,k为:
则接收端拼接的第二接收信号矩阵Yt可重新定义为:
是拼接的接收端检测的传输信号矩阵,Et=[et,1,...,et,K]是拼
接的由于解码误差造成的接收信号误差矩阵,得到:
其中,
采用 范数对信号检测误差进行稀疏重建,连续角度空间域的传输数据辅助的半盲稀疏信道重构模型 其中2≤t≤T,为:
上述公式中, 表示估计的 Et=[et,1,...,et,K]是拼接的由于
解码误差造成的接收信号误差矩阵, 为权重因子,||·||1为 范数
操作;经过按列拉伸的逆运算,由 得到所述第t个时间块的信道信息
4.根据权利要求3所述的半盲信道估计系统,其特征在于,还包括模型转换模块,所述模型转换模块用于将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式,具体为:
将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式
毫米波参数信道信息模型为:
其中,Ht为第t个时间块的信道信息,Lt表示第t个时间块的路径数,αt,l是第t个时间块第l条路径的增益, 和 分别表示第t个时间块第l条路径的Nt
个发射天线阵元和Nr个接收天线阵元的均匀线性天线阵列响应,其中,
其中, 和 分别为第t个时间块中第l条路径的离开角和到达角,Nt和Nr分别表示发射机和接收机上的天线阵元数,λ表示毫米波的波长,d是相邻天线阵元间距;
进一步将毫米波参数信道信息模型的公式转化成紧凑的矩阵形式:
其中, 是Lt条路径的增益向量,ht,j表示矩阵Ht的第j列,j=
1,...,Nt;
收发天线阵列响应矩阵定义为:
信道信息的二维原子范式由以下过程推导而得:
将信道信息Ht进行按列拉伸运算,得到信道信息新的列向量表达形式
其中, 为一个二
维复正弦原子,且满足 的限制,e和 分别表示Khatri-Rao积和克罗内克积,vec(·)表示矩阵的按列拉伸成列向量运算;
定义A={c(fl)||fl∈[0,1)×[0,1)}是第l条路径归一化的二维复正弦原子的集合,可得信道信息的二维原子范式表示为:
其中, 为第t个时间块的信道信息的二维原子范式表示,inf为取下确界操作,表示复数域,|·|表示取模操作。

说明书全文

一种半盲信道估计方法及系统

技术领域

[0001] 本发明所属无线通信系统技术领域,特别涉及一种半盲信道估计方法及系统。

背景技术

[0002] 毫米波大规模多天线系统较传统多天线系统有望达到吉比特每秒的数据传输速率。为了获得更高的通信系统吞吐量,收发机获取信道信息是非常关键的一步。实际上,信道部分是有效传送信息必须经过的媒质,准确、实时的信道信息是毫米波通信系统中实现波束成形、多用户调度、速率分配、自适应、信道均衡、信号检测和解码等几乎所有操作的基本依据。然而,在大规模多天线系统中,由于收发机天线数目的增加使得整个信道参数也成指数倍的增加,这样导致为了准确获取信道信息,在信道估计过程中系统的训练需要被分配更多的时频资源。随着毫米波参数化信道模型的提出,使得将信道估计问题转化为稀疏信号重构问题变得可行。
[0003] 对现有技术进行检索发现,A.Alkhateeb等在2014年IEEE Journal of Selected Topics on Signal Processing上发表的 Channel estimation and hybrid precoding for millimeter wave cellular systems[1](毫米波蜂窝系统信道估计和混合预编码) 中基于毫米波信道信息的稀疏表示模型:参数化信道,利用信道的空间稀疏性,提出了一种基于收发机多阶段训练波束交互的闭环方式估计有限路径的离开波束和到达波束对,以及相应的波束方向增益。然而在实际使用中,该方法将严重受限于训练波束图的设计,波束方向估计精度不高。此外,随着收发机数量的增加,训练波束开销成倍增加。而基于压缩感知网格匹配基的开环信道估计方法易造成基误匹配,不符合离开和到达角连续取值的事实,严重降低了信道估计的精度。
[0004] Junho Lee等在2016年IEEE Transactions on Communications 上发表的Channel Estimation via Orthogonal Matching Pursuit for Hybrid MIMO Systems in Millimeter Wave Communications [2](基于正交匹配追踪的毫米波通信混合多入多处系统信道估计) 中,基于收发角度量化的参数信道模型将信道估计问题转化为稀疏信号恢复问题,采用基于[0,π)区间内非均匀角度格分布的正交匹配追踪方法估计信道。然而,正交匹配追踪本质上是一种基于穷尽搜索的贪婪算法,会造成不能容忍的训练时延、能耗和计算复杂度。
[0005] Pejoski Slavche等在2015年IEEE Wireless Communications Letters上发表的Estimation of sparse time dispersive channels in pilot aided OFDM using atomic norm[3](基于原子范数的导频辅助正交频分复用系统的稀疏时间弥散信道估计)和 Zhang Peng等在2015年的IEEE International Conference on Communications上发表的Atomic norm denoising-based channel estimation for massive multiuser MIMO systems[4](基于原子范数解噪的大规模多用户多入多出系统的信道估计)中基于导频序列,采用一维原子范数最小化方法分别在[0,2π)和[-π/2,π/2]区间内逼近信道发射或者接收角度,估计信道信息矩阵。而由于毫米波信道模型中收发天线阵列响应与现有方法中单天线阵列响应矩阵表示形式的差异,现有基于一维原子范数最小化的信道估计方法无法直接扩展到同时具有收发双天线阵列响应的毫米波信道信息估计中。
[0006] 同时,综上所述的信道估计方法都仅限于基于导频的信道估计机制,产生大量的导频训练负载,严重降低了系统的频谱效率。

发明内容

[0007] 本发明提供了一种半盲信道估计方法及系统,旨在至少在一定程度上解决现有技术中的上述技术问题之一。
[0008] 为了解决上述问题,本发明提供了如下技术方案:
[0009] 一种半盲信道估计方法,包括以下步骤:
[0010] 步骤a:在第一个时间中,采用无格点压缩感知技术建立基于训练导频信号的稀疏信道重构模型,通过所述稀疏信道重构模型估计所述第一个时间块的信道信息;
[0011] 步骤b:在第t个时间块中,基于第t-1个时间块的信道信息,采用最大似然信号检测器检测所述第t个时间块的传输信号;其中, 2≤t≤T,T为所述时间块总数;
[0012] 步骤c:基于所述第t个时间块的传输信号,建立具有解码误差纠正功能的半盲稀疏信道重构模型,通过所述半盲稀疏信道重构模型估计所述第t个时间块的信道信息;
[0013] 步骤d:迭代执行步骤b和步骤c,直到完成T个时间块的信道信息估计。
[0014] 本发明实施例采取的技术方案还包括:所述步骤a前还包括:将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式。
[0015] 本发明实施例采取的技术方案还包括:所述二维原子范式表示为:
[0016]
[0017] 上述公式中, 为第t个时间块的信道信息的二维原子范式表示,inf为取下确界操作, 表示复数域,|·|表示取模操作。
[0018] 本发明实施例采取的技术方案还包括:在所述步骤a中,所述通过稀疏信道重构模型估计所述第一个时间块的信道信息具体为:
[0019] 所述稀疏信道重构模型为:
[0020]
[0021] 上述公式中, 表示估计的信道信息 为权重因子;经过按列拉伸的逆运算,由 得到所述第一个时间块的信道信息
[0022] 本发明实施例采取的技术方案还包括:在所述步骤c中,通过所述半盲稀疏信道重构模型估计当前时间块的信道信息具体为:
[0023] 所述半盲稀疏信道重构模型 为:
[0024]
[0025] 上述公式中, 表示估计的  Et=[et,1,...,et,K]是拼接的由于解码误差造成的接收信号误差矩阵, 为权重因子,||·||1为
范数操作;经过按列拉伸的逆运算,由 得到所述第t个时间块的信道信息[0026] 本发明实施例采取的另一技术方案为:一种半盲信道估计系统,包括:
[0027] 第一信道信息估计模块:用于在第一个时间块中,采用稀疏信号重建算法建立基于训练导频信号的稀疏信道重构模型,通过所述稀疏信道重构模型估计所述第一个时间块的信道信息;
[0028] 传输信号估计模块:用于在第t个时间块中,基于第t-1个时间块的信道信息,采用最大似然信号检测器检测所述第t个时间块的传输信号;其中,2≤t≤T,T为所述时间块总数;
[0029] 第二信道信息估计模块:用于基于所述第t个时间块的传输信号,采用半盲信道估计建立半盲稀疏信道重构模型,通过所述半盲稀疏信道重构模型估计所述第t个时间块的信道信息;
[0030] 迭代判断模块:用于判断是否完成T个时间块的信道信息估计,如果没有完成T时间块的信道信息估计,通过所述传输信号估计模块和第二信道信息估计模块进行迭代,直到完成T个时间块的信道信息估计。
[0031] 本发明实施例采取的技术方案还包括模型转换模块,所述模型转换模块用于将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式。
[0032] 本发明实施例采取的技术方案还包括:所述二维原子范式表示为:
[0033]
[0034] 上述公式中, 为第t个时间块的信道信息的二维原子范式表示,inf为取下确界操作, 表示复数域,|·|表示取模操作。
[0035] 本发明实施例采取的技术方案还包括:所述第一信道信息估计模块通过稀疏信道重构模型估计所述第一个时间块的信道信息具体为:
[0036] 所述稀疏信道重构模型为:
[0037]
[0038] 上述公式中, 表示估计的信道信息 为权重因子;经过按列拉伸的逆运算,由 得到所述第一个时间块的信道信息
[0039] 本发明实施例采取的技术方案还包括:所述第二信道信息估计模块通过半盲稀疏信道重构模型估计第t个时间块的信道信息具体为:
[0040] 所述半盲稀疏信道重构模型 为:
[0041]
[0042] 上述公式中, 表示估计的  Et=[et,1,...,et,K]是拼接的由于解码误差造成的接收信号误差矩阵, 为权重因子,||·||1为
范数操作;经过按列拉伸的逆运算,由 得到所述第t个时间块的信道信息[0043] 相对于现有技术,本发明实施例产生的有益效果在于:本发明实施例的半盲信道估计方法及系统采用基于传输信号的半盲信道估计技术,能够在连续域角度空间域高度紧密逼近信道参数的稀疏性,避免了基误匹配问题,相比基于导频的估计技术,降低系统训练负载,能够大幅度提高系统的频谱效率。同时为保证信道估计的准确性,同时引入解码误差补偿机制,提高系统频谱效率的同时改善系统信道解码误差,大幅提高信道估计算法的精确度。
附图说明
[0044] 图1是本发明实施例的半盲信道估计方法的流程图
[0045] 图2是本发明实施例的半盲信道估计系统的结构示意图;
[0046] 图3(a)至图3(c)是本发明实施例的仿真结果对比示意图。

具体实施方式

[0047] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
[0048] 本发明实施例的半盲信道估计方法及系统采用基于无格点压缩感知的毫米波大规模多输入多输出系统传输信号辅助的半盲信道估计方式。首先,在第一个时间块内,接收机根据已知的训练导频信号,采用基于无格点压缩感知的稀疏信号重建算法估计出当前信道信息;然后,在后续时间块,依次执行基于传输信号的半盲信道估计技术:依据上一个时间块中估计的信道信息,采用最大似然法检测出当前时间块内的传输信号;建立解码误差补偿机制,基于检测的数据信号,采用基于解码误差稀疏逼近机制的半盲信道估计技术实现快变信道估计;传输信号的检测和半盲信道估计交替执行,直到完成所有时间块的信道信息估计。
[0049] 具体地,请参阅图1,是本发明实施例的半盲信道估计方法的流程图。本发明实施例的半盲信道估计方法包括以下步骤:
[0050] 步骤100:将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式[0051] 在步骤100中,毫米波参数信道信息模型为:
[0052]
[0053] 公式(1)中,Ht为第t个时间块的信道信息,Lt表示第t个时间块的路径数,αt,l是第t个时间块第l条路径的增益, 和 分别表示第t个时间块第l条路径的Nt个发射天线阵元和Nr个接收天线阵元的均匀线性天线阵列响应,其中,[0054]
[0055]
[0056] 公式(2)和(3)中, 和 分别为第t个时间块中第l条路径的离开角和到达角,Nt和Nr分别表示发射机和接收机上的天线阵元数,λ表示毫米波的波长,d是相邻天线阵元间距。
[0057] 进一步将公式(1)转化成紧凑的矩阵形式:
[0058]
[0059] 公式(4)中, 是Lt条路径的增益向量,ht,j表示矩阵Ht的第j列,j=1,...,Nt。
[0060] 收发天线阵列响应矩阵定义为:
[0061]
[0062]
[0063] 任何求解模型与本发明方法中的式(4)形式等同,例如: yt,k=WtHHtxt,k+WtHnt,k或yt,k=Htxt,k+nt,k,且求解稀疏信号的参数在连续域取值,则视为同本发明方法为同一类问题。
[0064] 上述步骤100中,信道信息的二维原子范式由以下过程推导而得:
[0065] 将信道信息Ht进行按列拉伸运算,得到信道信息新的列向量表达形式[0066]
[0067] 公式(7)中,为一个二维复正弦原子,且满足 的限
制,e和 分别表示Khatri-Rao积和克罗内克积,vec(·)表示矩阵的按列拉伸成列向量运算。
[0068] 定义A={c(fl)‖fl∈[0,1)×[0,1)}是第l条路径归一化的二维复正弦原子的集合。至此,可得信道信息的二维原子范式表示为:
[0069]
[0070] 公式(8)中, 为第t个时间块的信道信息的二维原子范式表示,inf为取下确界操作, 表示复数域,|·|表示取模操作。
[0071] 步骤200:假设总共T个时间块,在第一个时间块的K个时隙中,基于发射机发射的训练导频信号S=[s1,...,sK],采用基于无格点压缩感知的稀疏信号重建算法建立连续角度空间域的基于训练导频信号的稀疏信道重构模型,通过该稀疏信道重构模型估计第一个时间块的信道信息
[0072] 在步骤200中,连续角度空间域的基于训练导频信号的稀疏信道重构模型由以下过程推导而得:
[0073] 第一个时间块的第k个时隙,发射端发射训练导频信号sk,则接收端接收到的信号为:
[0074] y1,k=W1HH1F1sk+W1Hn1,k  (9)
[0075] 公式(9)中,y1,k为第一个时间块的第k个时隙接收的信号, F1和W1分别表示发射端的第一预编码矩阵和接收端的第一测量矩阵,S=[s1,...,sK]表示K个时隙的拼接训练导频信号,n1,k是第一个时间块中第k个时隙接收机的加性白高斯噪声且其中 表示Nr维的全1列向量。
[0076] 假设接收机在第一个时间块的第k个时隙中保持F1和W1不变,则拼接的第一接收信号矩阵为:
[0077] Y1=W1HH1F1S+Q1  (10)
[0078] 公式(10)中,Y1为第一个时间块接收的信号, Q1=[W1Hn1,1,...,W1Hn1,K]是第一个时间块接收端的噪声矩阵。
[0079] 将Y1按列拉伸为新的列向量形式为:
[0080]
[0081] 至此,可得连续角度空间域的基于训练导频信号的稀疏信道重构模型为:
[0082]
[0083] 公式(12)中, 表示估计的信道信息 为权重因子。经过按列拉伸的逆运算即可由 得到第一个时间块的信道信息
[0084] 步骤300:在第t(2≤t≤T)个时间块的K个时隙中,基于上一个时间块中估计的信道信息 采用最大似然信号检测器,检测出当前时间块的传输信号
[0085] 在步骤300中,第t个时间块中(2≤t≤T)的第k个时隙内,发射端发射传输信号xt,k,则接收端接收到的信号为:
[0086] yt,k=WtHHtFtxt,k+WtHnt,k  (13)
[0087] 公式(13)中,yt,k为第t个时间块的第k个时隙接收的传输信号,Ft和Wt分别表示发射端的第二预编码矩阵和接收端的第二测量矩阵,nt,k是第t个时间块中第k个时隙接收端的加性白高斯噪声且
[0088] 假设接收机在第t个时间块的第k个时隙中保持Ft和Wt不变,则拼接的第二接收信号矩阵Yt为:
[0089] Yt=WtHHtFtXt+Qt(14)
[0090] 公式(14)中,Xt=[xt,1,...,xt,K]表示第t个时间块的第k个时隙拼接的发射信号,Qt=[WtHnt,1,WtHnt,2,...,WtHnt,K]是第t个时间块接收端的噪声矩阵。
[0091] 定义最大似然信号检测器为:
[0092]
[0093] 公式(15)中, 表示估计的当前时间块的传输信号xt,k,D 为xt,k可能的码本集合, 为第t-1个时间块估计的信道信息。
[0094] 已知 和Yt,则利用公式(15)即可得到当前时间块的传输信号
[0095] 步骤400:建立解码误差补偿机制,在第t个时间块的K个时隙中,基于检测的传输信号 采用基于解码误差稀疏逼近机制的半盲信道估计技术建立连续角度空间域的传输数据辅助的半盲稀疏信道重构模型,并估计当前时间块的信道信息
[0096] 在步骤400中,连续角度空间域的传输数据辅助的半盲稀疏信道重构模型由以下过程推导而得:
[0097] 第t个时间块中(2≤t≤T),将Yt按列拉伸为列向量形式为:
[0098]
[0099] 考虑检测的数据信号存在误差,则定义et,k为:
[0100]
[0101] 则接收端拼接的第二接收信号矩阵Yt可重新定义为:
[0102]
[0103] 公式(18)中, 是拼接的接收端检测的传输信号矩阵,Et=[et,1,...,et,K]是拼接的由于解码误差造成的接收信号误差矩阵。至此,得到:
[0104]
[0105] 公式(19)中,
[0106] 又因为实用通信系统需达到较高的信号检测精度,所以可以认为信号检测误差是稀疏的,本发明实施例中采用 范数对信号检测误差进行稀疏重建。因此,连续角度空间域的传输数据辅助的半盲稀疏信道重构模型 (2≤t≤T)为:
[0107]
[0108] 公式(20)中, 表示估计的 为权重
[0109] 因子,||·||1为 范数操作;经过按列拉伸的逆运算即可由 得到当前时间块的信道信息
[0110] 步骤500:判断是否完成所有时间块的信道信息估计,如果没有完成所有时间块的信道信息估计,交替执行步骤300和步骤400,直到完成所有时间块的信道信息;如果完成所有时间块的信道信息估计,执行步骤600;
[0111] 在步骤500中,步骤300和步骤400交替执行,直到估计完T个时间块的信道信息,并记作:
[0112] 步骤600:采用归一化均方误差(Normalized Mean Square Error,NMSE)对信道信息估计结果进行性能评价;
[0113] 在步骤600中,采用的NMSE评价标准定义为:
[0114]
[0115] 公式(21)中, 为取均值操作,||·||F为 范数操作。
[0116] 请参阅图2,是本发明实施例的半盲信道估计系统的结构示意图。本发明实施例的半盲信道估计系统包括模型转换模块、第一信道信息估计模块、传输信号估计模块、第二信道信息估计模块、迭代判断模块和性能评价模块。
[0117] 模型转换模块:用于将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式
[0118] 上述中,毫米波参数信道信息模型为:
[0119]
[0120] 公式(1)中,Ht为第t个时间块的信道信息,Lt表示第t个时间块的路径数,αt,l是第t个时间块第l条路径的增益, 和 分别表示第t个时间块第l条路径的Nt个发射天线阵元和Nr个接收天线阵元的均匀线性天线阵列响应,其中,[0121]
[0122]
[0123] 公式(2)和(3)中, 和 分别为第t个时间块中第l条路径的离开角和到达角,Nt和Nr分别表示发射机和接收机上的天线阵元数,λ表示毫米波的波长,d是相邻天线阵元间距。
[0124] 进一步将公式(1)转化成紧凑的矩阵形式:
[0125]
[0126] 公式(4)中, 是Lt条路径的增益向量,ht,j表示矩阵Ht的第j列,j=1,...,Nt。
[0127] 收发天线阵列响应矩阵定义为:
[0128]
[0129]
[0130] 上述中,信道信息的二维原子范式由以下过程推导而得:
[0131] 将信道信息Ht进行按列拉伸运算,得到信道信息新的列向量表达形式[0132]
[0133] 公式(7)中,为一个二维复正弦原子,且满足 的限
制,e和 分别表示Khatri-Rao积和克罗内克积,vec(·)表示矩阵的按列拉伸成列向量运算。
[0134] 定义A={c(fl)‖fl∈[0,1)×[0,1)}是第l条路径归一化的二维复正弦原子的集合。至此,可得信道信息的二维原子范式表示为:
[0135]
[0136] 公式(8)中, 为第t个时间块的信道信息的二维原子范式表示,inf为取下确界操作, 表示复数域,|·|表示取模操作。
[0137] 第一信道信息估计模块:假设总共T个时间块,用于在第一个时间块的K个时隙中,基于发射机发射的训练导频信号 S=[s1,...,sK],采用基于无格点压缩感知的稀疏信号重建算法建立连续角度空间域的基于训练导频信号的稀疏信道重构模型,通过该稀疏信道重构模型估计第一个时间块的信道信息
[0138] 上述中,连续角度空间域的基于训练导频信号的稀疏信道重构模型由以下过程推导而得:
[0139] 第一个时间块的第k个时隙,发射端发射训练导频信号sk,则接收端接收到的信号为:
[0140] y1,k=W1HH1F1sk+W1Hn1,k (9)
[0141] 公式(9)中,y1,k为第一个时间块的第k个时隙接收的信号, F1和W1分别表示发射端的第一预编码矩阵和接收端的第一测量矩阵,S=[s1,...,sK]表示K个时隙的拼接训练导频信号,n1 ,k是 第一个时间块中第k个时隙接收机的加性白 高斯噪 声且其中 表示Nr维的全1列向量。
[0142] 假设接收机在第一个时间块的第k个时隙中保持F1和W1不变,则拼接的第一接收信号矩阵为:
[0143] Y1=W1HH1F1S+Q1  (10)
[0144] 公式(10)中,Y1为第一个时间块接收的信号, Q1=[W1Hn1,1,...,W1Hn1,K]是第一个时间块接收端的噪声矩阵。
[0145] 将Y1按列拉伸为新的列向量形式为:
[0146]
[0147] 至此,可得连续角度空间域的基于训练导频信号的稀疏信道重构模型为:
[0148]
[0149] 公式(12)中, 表示估计的信道信息 为权重因子。经过按列拉伸的逆运算即可由 得到第一个时间块的信道信息
[0150] 传输信号估计模块:用于在第t(2≤t≤T)个时间块的K个时隙中,基于上一个时间块中估计的信道信息 采用最大似然信号检测器,检测出当前时间块的传输信号[0151] 上述中,第t个时间块中(2≤t≤T)的第k个时隙内,发射端发射传输信号xt,k,则接收端接收到的信号为:
[0152] yt,k=WtHHtFtxt,k+WtHnt,k  (13)
[0153] 公式(13)中,yt,k为第t个时间块的第k个时隙接收的传输信号,Ft和Wt分别表示发射端的第二预编码矩阵和接收端的第二测量矩阵,nt,k是第t个时间块中第k个时隙接收端的加性白高斯噪声且
[0154] 假设接收机在第t个时间块的第k个时隙中保持Ft和Wt不变,则拼接的第二接收信号矩阵Yt为:
[0155] Yt=WtHHtFtXt+Qt  (14)
[0156] 公式(14)中,Xt=[xt,1,...,xt,K]表示第t个时间块的第k个时隙拼接的发射信号,Qt=[WtHnt,1,WtHnt,2,...,WtHnt,K]是第t个时间块接收端的噪声矩阵。
[0157] 定义最大似然信号检测器为:
[0158]
[0159] 公式(15)中, 表示估计的当前时间块的传输信号xt,k,D 为xt,k可能的码本集合, 为第t-1个时间块估计的信道信息。
[0160] 已知 和Yt,则利用公式(15)即可得到当前时间块的传输信号
[0161] 第二信道信息估计模块:用于建立解码误差补偿机制,在第t个时间块的K个时隙中,基于检测的传输信号 采用基于解码误差稀疏逼近机制的半盲信道估计技术建立连续角度空间域的传输数据辅助的半盲稀疏信道重构模型,并估计当前时间块的信道信息[0162] 上述中,连续角度空间域的传输数据辅助的半盲稀疏信道重构模型由以下过程推导而得:
[0163] 第t个时间块中(2≤t≤T),将Yt按列拉伸为列向量形式为:
[0164]
[0165] 考虑检测的数据信号存在误差,则定义et,k为:
[0166]
[0167] 则接收端拼接的第二接收信号矩阵Yt可重新定义为:
[0168]
[0169] 公式(18)中, 是拼接的接收端检测的传输信号矩阵,Et=[et,1,...,et,K]是拼接的由于解码误差造成的接收信号误差矩阵。至此,得到:
[0170]
[0171] 公式(19)中,
[0172] 又因为实用通信系统需达到较高的信号检测精度,所以可以认为信号检测误差是稀疏的,本发明实施例中采用 范数对信号检测误差进行稀疏重建。因此,连续角度空间域的传输数据辅助的半盲稀疏信道重构模型 (2≤t≤T)为:
[0173]
[0174] 公式(20)中, 表示估计的 为权重因子,||·||1为范数操作。经过按列拉伸的逆运算即可由 得到当前时间块的信道信息
[0175] 迭代判断模块:用于判断是否估计完所有时间块的信道信息,如果没有估计完所有时间块的信道信息,通过传输信号估计模块和第二信道信息估计模块交替估计下一个时间块的传输信号和信道信息,直到估计完T个时间块的信道信息,并记作:如果完成所有时间块的信道信息估计,通过性能评价模块对信道信息估计结果进行性能评价。
[0176] 性能评价模块:用于在完成所有时间块的信道信息估计后,采用归一化均方误差标准对信道信息估计结果进行性能评价;采用的 NMSE评价标准定义为:
[0177]
[0178] 公式(21)中, 为取均值操作,||·||F为 范数操作。
[0179] 请参阅图3(a)至图3(c),为本发明实施例的仿真结果对比示意图。经在MATLAB平台进行仿真验证,在仿真中,假设T个时间块的路径条数相同,即Lt=L, 以得到NMSE与路径条数之间的关系。从仿真结果可以得出结论,本发明的NMSE评价标准表现随着信噪比的增加而减小(如图3(a)所示),随着路径数增加而增长(如图3(b)所示),随着观察次数(即每个时间块中用于信道估计的时隙数)的增加而减小(如图3(c)所示)。相较于现有主要相关算法,本发明表现出NMSE性能的优越性。
[0180] 本发明实施例的半盲信道估计方法及系统采用基于传输信号的半盲信道估计技术,能够在连续域角度空间域高度紧密逼近信道参数的稀疏性,避免了基误匹配问题,相比基于导频的估计技术,降低系统训练负载,能够大幅度提高系统的频谱效率。同时为保证信道估计的准确性,同时引入解码误差补偿机制,提高系统频谱效率的同时改善系统信道解码误差,大幅提高信道估计算法的精确度。
[0181] 对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
高效检索全球专利

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

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

申请试用

分析报告

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

申请试用

QQ群二维码
意见反馈