专利名称：高速单载波频域均衡超宽带系统的信道估计方法
技术领域：
本发明涉及一种高速单载波频域均衡超宽带(SC-FDEUWB)系统的信道估 计方法，具体涉及一种应用于频域均衡器的低复杂度信道估计方法。属于无线通 信系统的高速超宽带技术领域。
技术背景目前，超宽带(UWB)通信系统的物理层传输方案主要分为三类基于脉冲的UWB、多带正交频分复用(OFDM)的UWB以及单载波频域均衡(SC-FDE) 的UWB。其中SC-FDE UWB的传输方案比其它两种方案在整体性能和执行问 题上都具有优势。然而，类似于传统的窄带或宽带通信系统，准确的信道估计 信息对于确保UWB通信环境中可靠的数据传输起着至关重要的作用。一般而言，信道估计方法通常基于最小均方误差(MMSE)或者最小二乘 (LS)准则进行设计。OFDM系统中已有的信道估计算法可以通过简单的修改 推广到SC-FDE UWB系统中。其中的MMSE信道估计方法由于利用了信道的频 域相关特性，所以获得了重要的性能增益，但复杂度较高。尽管LS信道估计方 法执行简单，但信道估计误差较大，具有较高的均方误差(MSE)值。基于时 域最大(ML)似然标准的估计方法，可在一定程度上减少信道估计的MSE值， 但缺点是信道长度(或信道的有限延时扩展)信息需要在ML估计前被准确获 得。再者，OFDM系统中现有的基于时域低通滤波的信道估计方案也需事先获 得信道长度信息。这类信道估计算法涉及了相关信道长度的估计算法，因此， 增加了系统信道估计过程的持续时间。在SC-FDE UWB通信系统中，性能较优的信道估计方法是基于最小均方误 差准卯J的(Y. Wang and X. D. Dong, "Frequency-domain channel estimation for SC-FDE in UWB communications, ，， TSSfi" 7h3"5\3c^z'o/7s朋c。腳〃/w'cs^/ims1, vol. 54， No. 12， pp. 2155-2163， Dec. 2006.)。由于估计过程中利用了信道的统计特性，所以可以获得较优的信道估计性能。然而，该方法未涉及UWB系 统估计过程中所必须的背景噪声功率估计问题，同时该方法不仅需要对信道的 统计特性信息进行估计，而且还需要进行矩阵求逆运算，故复杂度很高。 发明内容本发明的目的在于针对现有技术的不足，提供一种高速单载波频域均衡超 宽带系统的信道估计方法，非但不需要信道统计特性的任何信息，而且还可获 取噪声功率值，运算过程中又能避免矩阵逆的计算，且可以实现较优的估计性能。为实现这一目的，本发明提供的方法中无需任何信道统计特性信息和复杂 矩阵运算，仅涉及了傅立叶变换和简单的乘法运算。接收端首先使用复数值的 频域恒模导频信号进行频域最小二乘信道估计；然后，基于得到的最小二乘信 道估计矢量构造新的数据矢量进行背景噪声功率估计；再将最小二乘信道估计 矢量进行逆傅立叶变换，由频域变换到时域，再送入时域滤波器进行去噪，采 用的判决阈值仅与背景噪声的功率有关；最后将去噪后的信道冲激响应估计通 过傅立叶变换得到最终的信道频响估计矢量。本发明的方法包括如下具体步骤1、 发送端采用单载波块传输方式，在各个传输块的前端插入循环前缀，多个插 入循环前缀后的传输块形成一个数据帧，数据帧内首个传输块用于传输复数 值的频域恒模导频信号，剩余的传输块用于传输数据，导频块及数据块以时 分复用方式传送；然后各传输块经脉冲成型后从天线发射出去，经历超宽带 信道并叠加高斯白噪声后到达接收端。2、 接收端将经过信道衰落和高斯白噪声污染的发射信号匹配滤波，然后对数据 帧内每个传输块去除循环前缀，再对首个去除循环前缀后的传输块进行傅立 叶变换，得到频域数据矢量。3、 将接收端已知的导频块进行傅立叶变换得到导频块的频域数据，由导频块的 频域数据和步骤2得到的频域数据矢量，计算得到最小二乘信道频响估计矢4、利用最小二乘信道频响估计矢量构造一个新数据矢量Afi， Afi中第A维数据的构造方法为A^t-^t-#;_t，A: = l,...,A72-l，其中^t为最小二乘信 道频响估计矢量中第A:个子载波处的频响值，7V为传输块大小；然后根据克 拉美罗界定理，求出Afi矢量的2范数值并进行矢量长度归一化，再乘以0.5得到背景噪声功率的估计值^ 。5、利用式K《)、/2^Mn,1 、计算时域滤波器的判决阈值，其中，参数p取1-Pj介于0和1且接近于1的任意值。6、 对最小二乘信道频响估计矢量作逆傅立叶变换，由频域转换到时域，得到信 道冲激响应估计矢量，然后将信道冲激响应估计矢量送入时域滤波器进行滤波，滤波准则为将信道冲激响应估计矢量中每一维系数的幅度与步骤5所计算的判决阈值K^J比较，保留不小于该阈值的信道冲激响应系数，并使其余元素为零，得到去噪后的信道冲激响应估计矢量。7、 对去噪后的信道冲激响应估计矢量作傅立叶变换，由时域转到频域，得到最 终的信道频响估计矢量，此时信道频响估计矢量的每一维具有更低的噪声水 平。本发明方法使用复数值的频域恒模导频信号进行频域最小二乘信道估计， 得到信道频响估计矢量，通过逆傅立叶变换得到信道冲激响应估计矢量，再使 用阈值滤波对信道冲激响应估计矢量进行去噪，滤波器的阈值计算仅需背景噪 声的功率值，而背景噪声的功率值可直接由最小二乘信道频响估计矢量构造的 新数据矢量计算得到。整个过程无需任何的信道统计特性信息，因此它的复杂 度远小于需要矩阵求逆的MMSE估计方法。在高速UWB系统环境下，本发明 信道估计方法的性能非常接近于最小均方误差估计方法，而且发明中同时解决 了背景噪声功率估计的问题，使其更加适用于实际的应用场景。本发明方法的计算量主要源于傅立叶变换，为了分析简单起见，本发明仅比较各算法所需的乘法次数的最高阶。传统的MMSE估计方法由于需要矩阵求 逆，所以复杂度高达O(W)阶。而由于本发明方法主要涉及傅立叶变换的操作， 因此复杂度仅为0(Mog27V)，远远低于0(W3)，所以易于实用化。
图l为本发明信道估计方法的流程框图。图2为CM3信道下传统信道估计算法和本发明的信道估计方法均方误差 性能的比较。
具体实施方式
以下结合附图和实施例对本发明的技术方案作进一步描述。 考虑单天线、单用户的SC-FDEUWB块传输系统，其通信环境为频选较严 重的CM3信道，背景噪声服从零均值的复加性高斯分布。每个长度为256的传 输块的块首插入了长度为64的循环前缀(CP)，此时系统无块间干扰产生。采 用本发明方法进行信道估计的流程框图如图l所示，具体的实施步骤如下1) 发送端采用单载波块传输方式，在各个传输块的前端插入循环前缀CP，多个 插入CP后的传输块形成一个数据帧，每个数据帧内首个传输块用于传输 Chu，s导频信号，以用于信道估计，导频块和剩余的100个数据块以时分复 用方式传送；然后各传输块经脉冲成型后从天线发射出去，经历超宽带信道 并叠加高斯白噪声后到达接收端；2) 接收端将经过信道衰落和高斯白噪声污染的发射信号匹配滤波，然后对数据 帧内每个传输块去除循环前缀，再对首个去除循环前缀后的传输块进行傅立 叶变换，得到频域数据矢量Y-XH + W。其中，矢量H二[i/。,…,/^，…，7/^]T表示阶数为Z的时域信道冲激响应矢量h = [&, ,^一JT的频响，X是所发射的收端已知的Chu's导频信号频域矢量[Io，…，Xt，…,X^]T所形成的对角矩阵，噪声的频域矢量w^『。,…,^,…,『^;r的各维是零均值、方差为^7备=7\^^的复高斯噪声，其中W为时域噪声的功率大小；3) 将接收端已知的导频块进行傅立叶变换得到导频块的频域数据 [A，…，X4，…,X^]T以生成相应的对角矩阵X，由步骤2)得到的频域数据 矢量Y和导频块的频域数据形成的对角矩阵X ，计算最小二乘信道频响估计 HiS=X—、，该值将用于接下来噪声功率的获取；4) 利用fiu构造一个新数据矢量Afi 。首先，以标量的形式表示fi^为 A=i^ = //A+_^,yt = 0，...,iV-l，其中仏为flu中第A个子载波处的频响义yt ^Ct值，iV为传输块大小。然后，根据高速UWB信道的实值性以及傅立叶变换 的共轭对称特性，可以构造Afi中第A维数据为 ^^^,=#4-#;—t，A: = l，...,iV/2-1。那么，由克拉美罗界定理可以得到噪声方差估计值为W =工￡ ;iV-2台15) 利用步骤3)和步骤4)计算时域滤波器的判决阈值。将最小二乘信道估计 矢量表示成fiiS =H + X—。相应的时域冲激相应估计矢量为fi￡S=h + FH(X—'W)二h + b，其中，矢量b二FH(X-^W)代表信道的估计误差或者估计噪声，而且傅立叶变换矩阵满足FHF:iV^的关系。矢量b的协方差矩阵为^bbHH《FH(XHX)—卞。由于系统中使用了复值导频序列，因此矢量b中各维幅度IW将服从瑞利分布，其分布函数为"(|6 |) = 1-6 2《，"=0,".，iV-l。其中，瑞利分布的方差C^ :+^五(blO卜《Zp^，"=0,.--,iV-1。那么，|6 |会按概率; 逼近某一值；>0，艮"-62《=; 。则有，r =、2《ln~!—, " = 0,...,iV-l，则参数p为介于0和1且接近于1 的任意值。利用导频的频域恒模特性并将方差《带入^;中，阈值可简化为r(o^) = r"=」2^1n~!~ ， " = 0，".,W-1。如果取概率/7 = 0.95，则上述阈值的含义为各维噪声元的幅度会以95%的概率逼近值〃(<^)。利用上述特 性，LS信道时域估计矢量fiu中的信道子空间和噪声子空间可以得到有效的分离，而且无需使用任何信道的统计特性； 6)对Hu作逆傅立叶变换，由频域转换到时域，得到信道冲激响应估计矢量fiu，然后将fiu送入时域滤波器进行滤波，滤波准则为将信道冲激响应估计矢量中每一维系数的幅度与所计算的判决阈值r(6。比较，保留不小于该阈值的信道冲激响应系数，并使其余元素为零，得到去噪后的信道冲激响应估计，波后的信道冲激响应估计矢量。此时，估计矢量中所含噪声元减少到原来的 Z /iV，从而实现了原信道冲激响应估计矢量fiu中噪声的滤除。注意到参数丄 是经过滤波的信道径数，它小于信道长度的真实值丄；7)对去噪后的信道冲激响应估计矢量fi-[&，…,^^]T作傅立叶变换，由时域转到频域，得到最终的信道频响估计矢量fi-[》。，…,^vjT，此时所有子载波 的噪声水平更低。上述的信道估计过程，通过时域滤波去噪大大改善了信道的估计性能，同 时，通过使用最小二乘信道估计矢量所构造的新数据矢量解决了背景噪声功率 估计的问题，整个过程未使用任何信道统计特性知识和复杂的矩阵逆操作。因 此，使用本发明的信道估计方法计算频域均衡器系数，不仅复杂度低，而且不A: = 0,...，iV-l。其中，￡i =[《。，...，^—J代表经过滤需要任何信道的统计特性信息，简单实用。图2为本发明方法采用的低复杂度信道估计方法、最小均方误差估计方法 以及最小二乘估计方法均方误差估计性能的比较结果。100次UWB信道是通过 将滚降因子为0.5奈奎斯特脉冲成型滤波器、接收机匹配滤波和正EE 802.15.3a CM3信道模型的任意一次实现进行巻积而形成。其它信道模型下，可得到相同 的算法相对性能比较结果。从图上可以看出，本发明的低复杂度信道估计算法 的均方误差性能明显优于传统的最小二乘信道估计方法，并且十分接近于传统 的MMSE信道估计方法。本发明信道估计方法不需要进行矩阵求逆运算，也无 需事先获取信道的统计特性信息，计算复杂度远小于最小均方误差估计方法， 而它的估计性能十分接近于最小均方误差估计方法，易于实用化。
权利要求
1、一种高速单载波频域均衡超宽带系统的信道估计方法，其特征在于包括如下具体步骤1)发送端采用单载波块传输方式，在各个传输块的前端插入循环前缀，多个插入循环前缀后的传输块形成一个数据帧，数据帧内首个传输块用于传输复数值的频域恒模导频信号，剩余的传输块用于传输数据，导频块及数据块以时分复用方式传送；然后各传输块经脉冲成型后从天线发射出去，经历超宽带信道并叠加高斯白噪声后到达接收端；2)接收端将经过信道衰落和高斯白噪声污染的发射信号匹配滤波，然后对数据帧内每个传输块去除循环前缀，再对首个去除循环前缀后的传输块进行傅立叶变换，得到频域数据矢量；3)将接收端已知的导频块进行傅立叶变换得到导频块的频域数据，由导频块的频域数据和步骤2)得到的频域数据矢量，计算得到最小二乘信道频响估计矢量；4)利用最小二乘信道频响估计矢量构造一个新数据矢量id="icf0001" file="S200810034568XC00011.gif" wi="19" he="4" top= "140" left = "133" img-content="drawing" img-format="tif" orientation="portrait" inline="no"/>中第k维数据的构造方法为id="icf0002" file="S200810034568XC00012.gif" wi="30" he="4" top= "151" left = "60" img-content="drawing" img-format="tif" orientation="portrait" inline="no"/>k＝1，...，N/2-1，其中id="icf0003" file="S200810034568XC00013.gif" wi="4" he="4" top= "150" left = "133" img-content="drawing" img-format="tif" orientation="portrait" inline="no"/>为最小二乘信道频响估计矢量中第k个子载波处的频响值，N为传输块大小；然后求出id="icf0004" file="S200810034568XC00014.gif" wi="5" he="3" top= "162" left = "154" img-content="drawing" img-format="tif" orientation="portrait" inline="no"/>矢量的2范数值并进行矢量长度归一化，再乘以0.5得到背景噪声功率的估计值id="icf0005" file="S200810034568XC00015.gif" wi="6" he="3" top= "173" left = "163" img-content="drawing" img-format="tif" orientation="portrait" inline="no"/>5)利用式id="icf0006" file="S200810034568XC00016.gif" wi="43" he="12" top= "183" left = "48" img-content="drawing" img-format="tif" orientation="portrait" inline="no"/>计算时域滤波器的判决阈值，其中，参数p取介于0和1且接近于1的任意值；6)对最小二乘信道频响估计矢量作逆傅立叶变换，由频域转换到时域，得到信道冲激响应估计矢量，然后将信道冲激响应估计矢量送入时域滤波器进行滤波，滤波准则为将信道冲激响应估计矢量中每一维系数的幅度与所计算的判决阈值id="icf0007" file="S200810034568XC00017.gif" wi="9" he="4" top= "233" left = "35" img-content="drawing" img-format="tif" orientation="portrait" inline="no"/>比较，保留不小于该阈值的信道冲激响应系数，并使其余元素为零，得到去噪后的信道冲激响应估计矢量；7)对去噪后的信道冲激响应估计矢量作傅立叶变换，由时域转到频域，得到最终的信道频响估计矢量。
全文摘要
本发明涉及一种适用于高速单载波频域均衡超宽带系统的信道估计方法，发明中无需任何信道统计特性信息和复杂矩阵运算，仅涉及了傅立叶变换和简单的乘法运算。接收端首先使用复数值的频域恒模导频信号进行频域最小二乘信道估计；然后，基于得到的最小二乘信道估计矢量构造新的数据矢量进行背景噪声功率估计；再将最小二乘信道估计矢量进行逆傅立叶变换，由频域变换到时域，再送入时域滤波器进行去噪，采用的判决阈值仅与背景噪声的功率有关；最后将去噪后的信道冲激响应估计矢量通过傅立叶变换得到最终的信道频响估计矢量。仿真结果表明本发明在保持估计方法低复杂度特点的同时，大大提高了信道估计的精度，甚至可以逼近最小均方误差信道估计方法。
文档编号H04B1/69GK101242388SQ200810034568
公开日2008年8月13日 申请日期2008年3月13日 优先权日2008年3月13日
发明者晨 何, 丹 王, 蒋铃鸽 申请人:上海交通大学