一种信道估计方法及系统与流程

本发明涉及无线通信技术领域，特别涉及一种信道估计方法及系统。

背景技术：

毫米波通信自由路损的急剧增加，早期关于毫米波通信的研究大多集中于室内场景。然而，得益于毫米等级的波长，大规模的毫米波通信系统天线阵列可以密集排列在较小体积的印制电路板上，因此基于大规模天线阵列的波束成形技术被广泛应用于实际的毫米波通信系统以对抗剧烈的路径损耗。纽约理工大学根据曼哈顿地区的室外实测数据指出，结合波束成形技术，毫米波通信系统中单个微基站的覆盖范围最高可达200米，与当前微小区在相同功耗条件下的覆盖范围相当。然而，在毫米波多天线系统中，高效的混合模拟数字波束成形设计是以准确的信道状态信息(channelstateinformation，csi)为前提的。即若csi未知，则混合波束成形的设计需要分两个步骤：首先确定模拟域的波束旋转方向，再完成基带部分的设计。相反地，若给定csi，则可通过矩阵分解同时实现模拟和基带部分的波束成形设计，使用较少的频谱资源获得更高的天线阵列增益。因此，要设计新的毫米波无线通信系统，信道估计是首要解决的问题。然而，仅基于导频实现信道估计需要在传输信号中插入大量的训练符号，产生大量的导频训练负载，严重降低了系统的频谱效率。基于传输信号的信道估计能够有效解决这个问题，大幅提高系统平铺效率。而另一方面，由于信号检测性能的好坏将直接影响到信道估计的精度，甚至整个多入多出(multiple-inputmultiple-output，mimo)系统性能的好坏，设计高性能的mimo信号检测算法已成为mimo通信中一项具有重大意义的关键技术。

对现有技术进行检索发现，pejoskislavche等在2015年ieeewirelesscommunicationsletters上发表的estimationofsparsetimedispersivechannelsinpilotaidedofdmusingatomicnorm(基于原子范数的导频辅助正交频分复用系统的稀疏时间弥散信道估计)和zhangpeng等在2015年的ieeeinternationalconferenceoncommunications上发表的atomicnormdenoising-basedchannelestimationformassivemultiusermimosystems(基于原子范数解噪的大规模多用户多入多出系统的信道估计)中基于导频序列，采用一维原子范数最小化方法分别在[0,2π)和[-π/2,π/2]区间内逼近信道发射或者接收角度，估计信道信息矩阵。eunchulyoon在2017年的ieeecommunicationsletters中发表的maximumlikelihooddetectionwithaclosed-formsolutionforthesquareqamconstellation中基于子集分割等手段简化了最大似然检测算法，提出检测信号关于接收信号和估计的信道信息的解析表达式。

上述中，由于毫米波信道模型中收发天线阵列响应与现有方法中单天线阵列响应矩阵表示形式的差异，现有基于一维原子范数最小化的信道估计方法无法直接扩展到同时具有收发双天线阵列响应的毫米波信道信息估计中。而且，上述信道估计方法都仅限于基于导频的信道估计机制，会产生大量的导频训练负载，严重降低了系统的频谱效率。此外，现有的信号检测算法没有考虑利用实际无线通信系统中检测误差较小的特性，导致信号检测误差较高。

技术实现要素：

本发明提供了一种信道估计方法及系统，旨在至少在一定程度上解决现有技术中的上述技术问题之一。

为了解决上述问题，本发明提供了如下技术方案：

一种信道估计方法，包括以下步骤：

步骤a：给定第0个时间块的导频序列；

步骤b：在第t个时间块中，基于t-1个时间块的传输信号，采用稀疏信道估计算法估计出所述第t个时间块的信道信息；其中，1≤t≤t，t为所述时间块总数；且当t＝1时，所述上一个时间块的传输信号为第0个时间块的导频序列；

步骤c：基于所述第t个时间块的信道信息，采用最大似然信号检测器检测出第t个时间块的传输信号；

步骤d：迭代执行步骤b和步骤c，直到完成t个时间块的信道信息估计和传输信号检测。

本发明实施例采取的技术方案还包括：所述步骤a前还包括：将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式；所述二维原子范式表示为：

上述公式中，为第t个时间块的信道信息的二维原子范式表示，inf为取下确界操作，表示复数域，|·|表示取模操作。

本发明实施例采取的技术方案还包括：在所述步骤b中，所述采用稀疏信道估计算法估计出第t个时间块的信道信息具体为：

基于原子范数最小化的稀疏信道重构模型为：

上述公式中，为权重因子；将上述公式转化为半正定规划的标准表达形式：

上述公式中，tr(·)表示求矩阵的迹操作，vt,g＝[vt,g(-nt+1),vt,g(-nt+2),...,vt,g(nt-1)]^t，是一个块托普利兹矩阵，将上述公式改写成适用于交替方向乘子法的表达形式：

将上述公式改写成增广拉格朗日的形式：

上述公式中，是对偶变量，ρ＞0是罚参数，η∞(·)是指示函数，当括号中的条件为真时，取值为0；否则为无穷；通过上述公式对未知参数求导。

本发明实施例采取的技术方案还包括：所述步骤d后还包括：采用归一化均方误差和误比特率标准分别对所述信道信息估计和传输信号检测结果进行性能评价。

本发明实施例采取的技术方案还包括：所述归一化均方误差标准定义为：

上述公式中，||·||f为范数操作；

所述误比特率标准定义为：

上述公式中，为取均值操作，||·||0为范数操作。

本发明实施例采取的另一技术方案为：一种信道估计系统，包括：

导频序列设定模块：用于给定第0个时间块的导频序列；

信道信息估计模块：用于在第t个时间块中，基于t-1个时间块的传输信号，采用稀疏信道估计算法估计出所述第t个时间块的信道信息；其中，1≤t≤t，t为所述时间块总数；且当t＝1时，所述上一个时间块的传输信号为第0个时间块的导频序列；

传输信号检测模块：用于基于所述第t个时间块的信道信息，采用最大似然信号检测器检测出第t个时间块的传输信号；

迭代判断模块：用于判断是否完成t个时间块的信道信息估计和传输信号检测，如果没有完成t个时间块的信道信息估计和传输信号检测，通过所述信道信息估计模块和传输信号检测模块进行迭代，直到完成t个时间块的信道信息估计和传输信号检测。

本发明实施例采取的技术方案还包括模型转换模块，所述模型转换模块用于将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式；所述二维原子范式表示为：

上述公式中，为第t个时间块的信道信息的二维原子范式表示，inf为取下确界操作，表示复数域，|·|表示取模操作。

本发明实施例采取的技术方案还包括：所述信道信息估计模块采用稀疏信道估计算法估计出第t个时间块的信道信息具体为：

基于原子范数最小化的稀疏信道重构模型为：

上述公式中，为权重因子；将上述公式转化为半正定规划的标准表达形式：

上述公式中，tr(·)表示求矩阵的迹操作，vt,g＝[vt,g(-nt+1),vt,g(-nt+2),...,vt,g(nt-1)]^t，是一个块托普利兹矩阵，将上述公式改写成适用于交替方向乘子法的表达形式：

将上述公式改写成增广拉格朗日的形式：

上述公式中，是对偶变量，ρ＞0是罚参数，η∞(·)是指示函数，当括号中的条件为真时，取值为0；否则为无穷；通过上述公式对未知参数求导。

本发明实施例采取的技术方案还包括性能评价模块，所述性能评价模块用于采用归一化均方误差和误比特率标准分别对所述信道信息估计和传输信号检测结果进行性能评价。

本发明实施例采取的技术方案还包括：所述归一化均方误差标准定义为：

上述公式中，||·||f为范数操作；

所述误比特率标准定义为：

上述公式中，为取均值操作，||·||0为范数操作。

相对于现有技术，本发明实施例产生的有益效果在于：本发明实施例的信道估计方法及系统通过在第0个时间块给定接收机已知的导频序列,第一个时间块中，基于给定的导频序列，采用稀疏信道估计算法估计出第一个时间块的信道信息，并依据第一个时间块的信道信息，采用具有解码误差补偿功能的最大似然信号检测器检测出第一个时间块的传输信号，在后续所有时间块内，接收机基于上一个时间块的传输信号，依次交替执行联合信道信息估计和传输信号检测策略，直到完成所有时间块的信道信息估计和传输信号检测，实现低复杂度的在线迭代快速算法，提高了算法的实时性，增强了算法工程实现的可能性。同时，为保证信号检测的准确性，本发明引入解码误差补偿机制，改善系统信道解码误差性能表现，继而大幅提高了信道估计算法的精确度。

附图说明

图1是本发明实施例的信道估计方法的流程图；

图2是本发明实施例的信道估计系统的结构示意图；

图3(a)至图3(d)是本发明实施例的仿真结果对比示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

本发明实施例的信道估计方法及系统提供了一种高精度和低误码率的联合信道估计和信号检测的快速实现方法。在第0个时间块给定接收机已知的导频序列。在后续所有时间块内，接收机依次交替执行联合信道估计和信号检测策略：首先，基于已知的导频序列，采用稀疏信道估计算法估计出第一个时间块的信道信息；然后，依据第一个时间块的信道信息，采用具有解码误差补偿功能的最大似然信号检测器检测出第一个时间块的传输信号；在后续所有时间块内，接收机基于上一个时间块的传输信号，依次交替执行联合信道信息估计和传输信号检测策略，直到完成所有时间块的信道信息估计和传输信号检测。

具体地，请参阅图1，是本发明实施例的信道估计方法的流程图。本发明实施例的信道估计方法包括以下步骤：

步骤100：将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式

在步骤100中，毫米波多输入多输出系统的信号传输过程如下：

第t个时间块中的第k个时隙，发射端发射传输信号xt,k，则接收端接收到的信号为：

yt,k＝wt^hhtftxt,k+wt^hnt,k(1)

公式(1)中，yt,k为第t个时间块中的第k个时隙接收的传输信号，ht为第t个时间块的信道信息，ft和wt分别表示发射端的预编码矩阵和接收端的测量矩阵，nt,k是第t个时间块中第k个时隙接收端的加性白高斯噪声且其中，表示nr维的列向量。

假设接收机在第t个时间块中的第k个时隙保持ft和wt不变，则拼接的第一传输信号矩阵yt为：

yt＝wt^hhtftxt+qt(2)

公式(2)中，xt＝[xt,1,...,xt,k]表示第t个时间块中的第k个时隙拼接的发射信号，qt＝[wt^hnt,1,wt^hnt,2,...,wt^hnt,k]是第t个时间块中接收端的噪声矩阵。

然后，将yt按列拉伸为列向量形式为：

公式(3)中，vec(·)表示矩阵的按列拉伸成列向量运算，和分别表示信道信息矩阵ht和噪声矩阵qt按列拉伸后形成的列向量，表示克罗内克积。

毫米波参数信道信息模型为:

公式(4)中，lt表示第t个时间块中的路径数，αt,l是第t个时间块中第l条路径的增益，和分别表示第t个时间块中第l条路径的nt个发射天线阵元和nr个接收天线阵元的均匀线性天线阵列响应，其中，

公式(5)和(6)中，和分别为第t个时间块中第l条路径的离开角和到达角，nt和nr分别表示发射机和接收机上的天线阵元数，λ表示毫米波的波长，d表示相邻天线阵元间距。

进一步将公式(4)转化成紧凑的矩阵形式：

公式(7)中，是lt条路径的增益向量，ht,j表示矩阵ht的第j列，j＝1,...,nt。收发天线阵列响应矩阵定义为：

任何求解模型与本发明实施例中的公式(4)形式等同，例如：yt,k＝wt^hhtxt,k+wt^hnt,k或yt,k＝htxt,k+nt,k，且求解稀疏信号的参数在连续域取值，则视为同本发明方法为同一类问题。

步骤100中，信道信息的二维原子范式由以下过程推导而得：

将信道信息矩阵ht进行按列拉伸运算，得到信道信息的列向量表达形式

公式(10)中，为一个二维复正弦原子，且满足的限制，⊙表示khatri-rao积。

定义是第l条路径归一化的二维复正弦原子的集合。至此，可得信道信息的二维原子范式表示为：

公式(11)中，为第t个时间块的信道信息的二维原子范式表示，inf为取下确界操作，表示复数域，|·|表示取模操作。

步骤200：假设总共有t个时间块，给定第0个时间块中k个时隙的导频序列

在步骤200中，第0个时间块中k个时隙的导频序列可根据不同的导频设计准则进行设计，在本发明实施例中，将导频序列设计为每列仅一个元素为1，其余元素皆为0。

步骤300：在第t(1≤t≤t)个时间块的k个时隙中，基于上一个时间块的传输信号采用稀疏信道估计算法估计出当前时间块的信道信息

在步骤300中，当t＝1时，上一个时间块的传输信号为步骤200中给定的导频序列当第2≤t≤t时，上一个时间块的传输信号为步骤400中检测出的当前时间块的传输信号

采用稀疏信道估计算法估计出当前时间块的信道信息具体包括：

基于原子范数最小化的稀疏信道重构模型为：

公式(12)中，为权重因子。为便于计算，本发明实施例将公式(12)转化为半正定规划的标准表达形式：

公式(13)中，tr(·)表示求矩阵的迹操作，vt,g＝[vt,g(-nt+1),vt,g(-nt+2),...,vt,g(nt-1)]^t，是一个块托普利兹矩阵，进一步定义为：

公式(14)中，toep(·)表示托普利兹矩阵，且第一列为输入向量，具体定义为：

针对大规模天线阵列配置场景，现有求解半正定规划的工具sedumi和sdpt3的运算速度过慢，不能满足实际无线通信应用的需求。因此，本发明实施例给出一种基于交替方向乘子法的快速稀疏信道估计策略：

首先，将公式(13)改写成以下表达形式以适用于交替方向乘子法：

将公式(16)改写成增广拉格朗日的形式：

公式(17)中，是对偶变量，ρ＞0是罚参数，η∞(·)是指示函数，当括号中的条件为真时，取值为0；否则为无穷。然后，交替方向乘子法包含以下的更新步骤：

为进一步解释公式(18)和(19)，本发明实施例引入如下的分割：

公式(21)和(22)中，和是nrnt×nrnt矩阵，和是nrnt维的列向量，和是常量。令以及正交补矩阵满足令公式(17)分别对各个未知参数求导，可得各未知参数的更新方法：

公式(26)中，i1＝[1,0,0,...,0]^t，表示的伴随矩阵。令其中以及cg＝[cg(-nt+1),...,cg(nt-1)]^t。具体地，

公式(27)中，输出矩阵的第(g,j)个nt×nt的子矩阵a(g,j)，g＝-nr+1,...,nr-1，j＝-nt+1,...,nt-1。

而ψt^l的更新应依据以下半正定核的投影：

可以通过将矩阵进行特征值分解并将负的特征值置为0实现矩阵向半正定核的投影。

经过按列拉伸的逆运算即可由得到当前时间块的信道信息估计矩阵

步骤400：在第t(1≤t≤t)个时间块中的k个时隙中，基于当前时间块的信道信息采用具有解码误差补偿功能的最大似然信号检测器，检测出当前时间块的传输信号

在步骤400中，第t(1≤t≤t)个时间块中，定义的最大似然信号检测器为：

公式(29)中，表示检测的第t个时间块的传输信号xt,k，d为xt,k可能的码本集合，为第t个时间块估计的信道信息。

第t个时间块中(1≤t≤t)，具有解码误差补偿功能的最大似然信号检测器由以下过程推导而得：

考虑采用公式(29)所述的最大似然信号检测器检测的信号存在误差，定义et,k为：

则接收端拼接的第二传输信号矩阵yt可重新定义为：

公式(31)中，表示通过拼接k个时隙的传输信号得到的第t个时间块的传输信号，et＝[et,1,...,et,k]表示拼接的第t个时间块的解码误差矩阵。

因为实用通信系统需达到较高的信号检测精度，所以可以认为信号检测误差是稀疏的，本发明实施例中采用l1范数对信号检测误差进行稀疏重建。至此，得到具有解码误差补偿功能的最大似然信号检测器：

公式(32)中，表示更新的传输信号表示估计的数据信号检测误差et,k，为权重因子，||·||1为范数操作。拼接k个时隙获得的传输信号和估计的数据信号检测误差即可得到第t个时间块拼接的传输信号矩阵和拼接的估计解码误差矩阵

基于上述，任何现有的最大似然信号检测的替代算法与本发明实施例中的检测误差补偿机制的直接加性组合皆视为本发明实施例的替代方案。

步骤500：判断是否完成所有时间块的信道信息估计和传输信号检测，如果没有完成所有时间块的信道信息估计和传输信号检测，交替执行步骤300和步骤400，直到完成所有时间块的信道信息估计和传输信号检测；如果完成所有时间块的信道信息估计和传输信号检测，执行步骤600；

在步骤500中，步骤300和步骤400交替执行，直到估计完t个时间块的信道信息估计和传输信号检测，并分别记作和

步骤600：采用归一化的均方误差(normalizedmeansquareerror，nmse)和误比特率(biterrorratio，ber)标准分别对信道信息估计和传输信号检测结果进行性能评价。

在步骤600中，采用的nmse评价标准定义为：

公式(33)中，||·||f为范数操作。

采用的ber标准定义为：

公式(34)中，为取均值操作，||·||0为范数操作。

请参阅图2，是本发明实施例的信道估计系统的结构示意图。本发明实施例的信道估计系统包括模型转换模块、导频序列设定模块、信道信息估计模块、传输信号检测模块、迭代判断模块和性能评价模块。具体地：

模型转换模块：用于将毫米波参数信道信息模型转换为符合构造二维原子范式表达的形式

其中，毫米波多输入多输出系统的信号传输过程如下：

第t个时间块中的第k个时隙，发射端发射传输信号xt,k，则接收端接收到的信号为：

yt,k＝wt^hhtftxt,k+wt^hnt,k(1)

公式(1)中，yt,k为第t个时间块中的第k个时隙接收的传输信号，ht为第t个时间块的信道信息，ft和wt分别表示发射端的预编码矩阵和接收端的测量矩阵，nt,k是第t个时间块中第k个时隙接收端的加性白高斯噪声且其中，表示nr维的列向量。

假设接收机在第t个时间块中的第k个时隙保持ft和wt不变，则拼接的第一传输信号矩阵yt为：

yt＝wt^hhtftxt+qt(2)

公式(2)中，xt＝[xt,1,...,xt,k]表示第t个时间块中的第k个时隙拼接的发射信号，qt＝[wt^hnt,1,wt^hnt,2,...,wt^hnt,k]是第t个时间块中接收端的噪声矩阵。

然后，将yt按列拉伸为列向量形式为：

公式(3)中，vec(·)表示矩阵的按列拉伸成列向量运算，和分别表示信道信息矩阵ht和噪声矩阵qt按列拉伸后形成的列向量，表示克罗内克积。

毫米波参数信道信息模型为:

公式(4)中，lt表示第t个时间块中的路径数，αt,l是第t个时间块中第l条路径的增益，和分别表示第t个时间块中第l条路径的nt个发射天线阵元和nr个接收天线阵元的均匀线性天线阵列响应，其中，

公式(5)和(6)中，和分别为第t个时间块中第l条路径的离开角和到达角，nt和nr分别表示发射机和接收机上的天线阵元数，λ表示毫米波的波长，d表示相邻天线阵元间距。

进一步将公式(4)转化成紧凑的矩阵形式：

公式(7)中，是lt条路径的增益向量，ht,j表示矩阵ht的第j列，j＝1,...,nt。收发天线阵列响应矩阵定义为：

任何求解模型与本发明实施例中的公式(4)形式等同，例如：yt,k＝wt^hhtxt,k+wt^hnt,k或yt,k＝htxt,k+nt,k，且求解稀疏信号的参数在连续域取值，则视为同本发明方法为同一类问题。

信道信息的二维原子范式由以下过程推导而得：

将信道信息矩阵ht进行按列拉伸运算，得到信道信息的列向量表达形式

公式(10)中，为一个二维复正弦原子，且满足的限制，⊙表示khatri-rao积。

定义是第l条路径归一化的二维复正弦原子的集合。至此，可得信道信息的二维原子范式表示为：

公式(11)中，为第t个时间块的信道信息的二维原子范式表示，inf为取下确界操作，表示复数域，|·|表示取模操作。

导频序列设定模块：假设总共t个时间块，用于给定第0个时间块中k个时隙的导频序列其中，第0个时间块中k个时隙的导频序列可根据不同的导频设计准则进行设计，在本发明实施例中，将导频序列设计为每列仅一个元素为1，其余元素皆为0。

信道信息估计模块：用于在第t(1≤t≤t)个时间块的k个时隙中，基于上一个时间块的传输信号采用稀疏信道估计算法估计出当前时间块的信道信息其中，当t＝1时，上一个时间块的传输信号为给定的第0个时间块中k个时隙的导频序列当第2≤t≤t时，上一个时间块的传输信号为传输信号检测模块检测出的当前时间块的传输信号

基于原子范数最小化的稀疏信道重构模型为：

公式(12)中，为权重因子。为便于计算，本发明实施例将公式(12)转化为半正定规划的标准表达形式：

公式(13)中，tr(·)表示求矩阵的迹操作，vt,g＝[vt,g(-nt+1),vt,g(-nt+2),...,vt,g(nt-1)]^t，是一个块托普利兹矩阵，进一步定义为：

公式(14)中，toep(·)表示托普利兹矩阵，且第一列为输入向量，具体定义为：

针对大规模天线阵列配置场景，现有求解半正定规划的工具sedumi和sdpt3的运算速度过慢，不能满足实际无线通信应用的需求。因此，本发明实施例给出一种基于交替方向乘子法的快速稀疏信道估计策略：

首先，将公式(13)改写成以下表达形式以适用于交替方向乘子法：

将公式(16)改写成增广拉格朗日的形式：

公式(17)中，是对偶变量，ρ＞0是罚参数，η∞(·)是指示函数，当括号中的条件为真时，取值为0；否则为无穷。然后，交替方向乘子法包含以下的更新步骤：

为进一步解释公式(18)和(19)，本发明实施例引入如下的分割：

公式(21)和(22)中，和是nrnt×nrnt矩阵，和是nrnt维的列向量，和是常量。令以及正交补矩阵满足令公式(17)分别对各个未知参数求导，可得各未知参数的更新方法：

公式(26)中，i1＝[1,0,0,...,0]^t，表示的伴随矩阵。令其中以及cg＝[cg(-nt+1),...,cg(nt-1)]^t。具体地，

公式(27)中，输出矩阵的第(g,j)个nt×nt的子矩阵a(g,j)，g＝-nr+1,...,nr-1，j＝-nt+1,...,nt-1。

而ψt^l的更新应依据以下半正定核的投影：

可以通过将矩阵进行特征值分解并将负的特征值置为0实现矩阵向半正定核的投影。

经过按列拉伸的逆运算即可由得到当前时间块的信道信息估计矩阵

传输信号检测模块：用于在第t(1≤t≤t)个时间块中的k个时隙，基于当前时间块中的信道信息采用具有解码误差补偿功能的最大似然信号检测器，检测出当前时间块的传输信号其中，第t个时间块中(1≤t≤t)，定义的最大似然信号检测器为：

公式(29)中，表示检测的第t个时间块的传输信号xt,k，d为xt,k可能的码本集合，为第t个时间块估计的信道信息。

第t个时间块中(1≤t≤t)，具有解码误差补偿功能的最大似然信号检测器由以下过程推导而得：

考虑采用公式(29)所述的最大似然信号检测器检测的信号存在误差，定义et,k为：

则接收端拼接的第二传输信号矩阵yt可重新定义为：

公式(31)中，表示通过拼接k个时隙的传输信号得到的第t个时间块的传输信号，et＝[et,1,...,et,k]表示拼接的第t个时间块的解码误差矩阵。

因为实用通信系统需达到较高的信号检测精度，所以可以认为信号检测误差是稀疏的，本发明实施例中采用范数对信号检测误差进行稀疏重建。至此，得到具有解码误差补偿功能的最大似然信号检测器：

公式(32)中，表示更新的传输信号表示估计的数据信号检测误差et,k，为权重因子，·||1为范数操作。拼接k个时隙获得的传输信号和估计误差即可得第t个时间块拼接的传输信号矩阵和拼接的估计解码误差矩阵

基于上述，任何现有的最大似然信号检测的替代算法与本发明实施例中的检测误差补偿机制的直接加性组合皆视为本发明实施例的替代方案。

迭代判断模块：用于判断是否完成所有时间块的信道信息估计和传输信号检测，如果没有完成所有时间块的信道信息估计和传输信号检测，通过信道信息估计模块和传输信号检测模块交替估计下一个时间块的信道信息和传输信号，直到完成所有时间块的信道信息估计和传输信号检测，并分别记作：和如果完成所有时间块的信道信息估计和传输信号检测，通过性能评价模块对信道信息估计和传输信号检测结果进行性能评价。

性能评价模块：用于采用归一化的均方误差(normalizedmeansquareerror，nmse)和误比特率(biterrorratio，ber)标准分别对信道信息估计和传输信号检测结果进行性能评价。其中，采用的nmse评价标准定义为：

公式(33)中，||·||f为范数操作。

采用的ber标准定义为：

公式(34)中，为取均值操作，||·||0为范数操作。

请参阅图3(a)至图3(d)，为本发明实施例的仿真结果对比示意图。经在matlab平台进行仿真验证，在仿真中，假设t个时间块的路径条数相同，即lt＝l，从仿真结果可以得出结论，一方面，本发明实施例中的解码误差补偿功能能够有效地降低ber，随着snr(信噪比)的增长，ber减小的越明显。另一方面，本发明实施例中的nmse随着信噪比的增加而减小，随着路径数增加而增长，随着观察次数(即每个时间块中用于信道估计的时隙数)的增加而减小。相较于现有相关算法，本发明实施例表现出nmse性能的优越性。

本发明实施例的信道估计方法及系统通过在第0个时间块给定接收机已知的导频序列,第一个时间块中，基于给定的导频序列，采用稀疏信道估计算法估计出第一个时间块的信道信息，并依据第一个时间块的信道信息，采用具有解码误差补偿功能的最大似然信号检测器检测出第一个时间块的传输信号，在后续所有时间块内，接收机基于上一个时间块的传输信号，依次交替执行联合信道信息估计和传输信号检测策略，直到完成所有时间块的信道信息估计和传输信号检测，实现低复杂度的在线迭代快速算法，提高了算法的实时性，增强了算法工程实现的可能性。同时，为保证信号检测的准确性，本发明引入解码误差补偿机制，改善系统信道解码误差性能表现，继而大幅提高了信道估计算法的精确度。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。