您的当前位置:首页正文

利用频谱特性的速度自适应信道估计方法

2023-12-30 来源:年旅网
维普资讯 http://www.cqvip.com 第23卷第5期 信号处理 Vo1.23. No.5 2007年1O月 SIGNAL PROCESSING 0ct.2oo7 利用频谱特性的速度自适应信道估计方法 何婷玉 邱玲朱近康 (中国科学技术大学电子工程与信息科学系 个人通信与扩频实验室,合肥230027) 摘要:本文利用信道在移动环境下频谱带限和对称的特征,提出一种新的对移动速度自适应的信道估计方法。该方 法通过寻找频谱自相关的峰值来确定频域滤波器的长度。相对传统速度自适应方法,该方法不需要估计多普勒频率或移动 速度,降低了系统实现的复杂度;也不需要信道的长时统计量或信噪比信息,实时性较高。信号在频域的处理必然会有频 谱泄漏问题,传统方法滤除噪声的同时也去掉了部分信号,对此本文同时还提出一种恢复泄漏的频域信号的方法,使估计 性能有所提高。仿真结果显示,从步行到高速移动的环境下,本文提出的信道估计方法的均方误差(MSE)都保持在较低的 水平,具有良好的对速度自适应的估计性能。 关键词:信道估计;多普勒频移;自相关函数;频谱泄漏 Spectrum—Based Velocity-Adaptive Channel Estimation HE Ting・-yu QIU Ling ZHU Jin・-kang (PCN&SS Lab,Dept.of EEIS,University of Science and Technology of China,Hefei,230027) Abstract: A velocity adaptive channel estimation method using a band—limited filter is proposed in this paper.The pass—band of the filter is determined by searching the second highest peak location of the spectrum autocorrelation,which changes with velocity.Corn— pared to the usual velocity adaptive estimator,no velocity or Doppler frequency estimator is needed,which greatly lowers the complexity, and no channel’S higher order statistics or signal—to—noise ratio(SNR)information is needed,which reduces delay.The components which are hi gher than the maximum Doppler frequency ale forced to zeros because they are just noise if there is no spectrum leakage, which is in fact not the real case due to the truncation in the time domain.Therefore,a compensation for the removed leaking signal is made after zero forcing.Computer simulation shows that,in a wide range of mobile speed,the performance of the proposed correlation method in terms of the mean square error(MSE)of estimation,only degrades a little compared to the idela case in which Doppler fre— quency is perfecdy known,and the compensation for the leaking signals improves the accuracy of channel estimation. Key words: channel estimation;Doppler frequency;autocorrelation function;specturm leakage 1 引言 和[2]的自适应思想都是从一组预先确定的滤波器中选出一 个适合于当前移动速度的。文献[3]在频域采用迫零滤波 未来移动通信系统需要支持在各种移动速度下的良好 器,其带宽自适应于速度,但需要假设已知或采用迭代的方 通信性能。但是当移动台作高速移动时,时变衰落信道将使 法逼近多普勒频率,并且迭代的周期较长,次数较多。如果 通信系统中的信道估计技术面临诸多困难。解调、解码,以 只是为了信道估计而引入速度或多普勒频率估计模块会使 及链路自适应处理等模块都需要用到信道估计提供的信息, 系统的开销和复杂度大大增加,因为成熟的速度估计或多普 因而系统的性能受移动性的影响较大。所以速度自适应信 勒频率估计方法通常需要信道的高阶统计量[4],或者长时 道估计方法的研究成为移动通信中非常重要的课题。 统计量如电平交叉率(LCR)[5]一[7],并且需要精确的信噪 目前已有的对速度自适应的信道估计方法均需要在信 比信息。在移动环境下为了使信道估计的更新能够及时跟 道估计中添加速度估计或者多普勒频率估计的模块[1]一 上信道的变化,实用的估计算法的计算不能过于复杂,延时 [3]。文献[1]采用一阶IIR滤波器,其系数自适应于速度。 不能太大。 文献[2]采用FIR滤波器,其阶数自适应于多普勒频率。[1] 本文提出一种在频域估计信道的新方法,通过寻找频谱 收稿日期:2006年2月10日;修回日期:2006年6月22日 ’本项研究得到国家自然科学基金(6O572O66)和863项目(2005AA123920)的支持。 维普资讯 http://www.cqvip.com

第5期 利用频谱特性的速度自适应信道估计方法 797 自相关峰值出现的位置来确定频域去噪滤波器的通带长度, 去掉带外的噪声分量。相对于传统速度自适应方法,该方法 度而性能接近多普勒频率已知的理想情况,并且本方法不需 要长时统计量,实时陛较高。信号在频域的处理必然会有频 低频部分,理想情况下高于最大多普勒频率的成分全部为噪 声,可用低通滤波器滤除。移动速度越快,信号带宽的多普 中间最高峰的位置对应作相关运算的过程中频谱自身完全 首尾对齐的时刻,而左右两个次高峰的位置则对应信号频带 完全重合的时刻,即是所要设计的滤波器通带的边界。移动 速度越快,两个次高峰越向中心主峰靠拢,通过搜索第一个 次高峰出现的位置就可以让滤波器的带宽自适应于速度。 fd/fs=0.01的频谱 fd/fs=0,01的频谱自相关 不需要估计移动速度或多普勒频率,降低了系统实现的复杂 勒扩展越严重。从频谱自相关图上可以看到出现三个峰值, 谱泄漏问题,传统频域估计方法抑制噪声的同时也损失了部 分信号,因此本文同时还提出一种能够恢复泄漏的频域信号 的方案,使估计性能有所提高。 本文余下部分安排如下:第二部分是系统模型和本文提 出的速度自适应的原理。第三部分详细分析了速度自适应 信道估计的方法和解决频谱泄漏的方法。第四部分给出仿 真结果。第五部分是结论。 2系统模型和速度自适应原理 100 200 300 0 200 I ” 4OO 600 本文采用的系统在发送端发送数据的过程中周期性地 fd/fs=O,1的频谱 fd/fs=0.1的频谱自相关 插入导频信号。接收端假设时偏和频偏均已准确校正,则经 过平坦衰落信道后的基带信号采样后表示成: r(i)=h(i) (i)+n(i), (1) 其中 (i)表示发送的第i个符号,n(i)为方差为 :的加性高 斯白噪声,h(i)表示衰落信道的复增益系数[8],可以建模成 k…— , … 0 200 4OO 600 复高斯随机变量,方差为 。假设该移动无线信道的功率谱 符合Jakes模型,那么间隔为 的信道系数的自相关函数表示 为: R (T): (2 T), 图1一次256点FfTr实现的信号频谱及其自相关 3速度自适应信道估计方法 瑚 抛 啪 m∞0 (2) 枷 瑚 瑚 m 0 本文提出的利用频谱特性的信道估计方法实现框图如 3,1利用频谱特性的信道估计方法 其中 (・)是第一类零阶贝塞尔函数 表示最大多普勒频 率,它是移动速度v和载波波长九的函数: = 图2所示,其中虚线表示提出的解决频谱泄漏的方案,将在 (3) 九, ,,一 3.2节中论述。以Ⅳ 点为一次处理周期。首先,我们将序 列r口(n)截取Ⅳ 点: (m),m=0,1…Nr, 一1,通过FFr变 信道的功率谱密度是带限的,在带宽内表达式为: s (f): /《 √1一( ) ),l/l . (4) 换至频域: ( )=FFTt rp(m)}= rp(m)e啊 , =0,l…・ 一1. (6) 在以下的信道估计算法描述中,只是针对导频信号,数 据部分的信道估计可以通过各种内插方式得到。假设导频 的数值为1,插入的频率为 ,在发送序列中每间隔 点插入 一次。因此将接收到的导频序列重新写为: (n)=h(nNs) (nN,)+n(nN,)=h(n )+n(n ) =再将频域信号进行自相关运算,用以下计算方法作为对自相 关的逼近: h (n)+n (n),n=0,1,2.... (5) R (△) 磊 ( ) ( +△一(NFFr一1)), (7) 则 (n)又是信道在导频处的最小均方(Ls)估计[9]。 信道估计对速度自适应,即是让估计器的某个参量随速 度变化而改变,使信道估计的误差在较大移动速度范围内均 保持较低。传统的自适应方案先要估计移动速度然后选择 合适的信道估计滤波器。本文提出的算法基于移动环境下 其中・ 表示取共轭。在[0,Nr, ]范围内搜索l (△)l的第 二极大值点的位置(I.I表示求绝对值),或者等效地,搜索第 一个峰值出现的位置,设为第Ⅳ,点,序号为Ⅳ,一1。于是可 以采用频域迫零的方法[3],将信号带外的频谱迫零从而抑 频谱带限和对称的特征,能够减小信道估计模块的开销和时 延。用归一化多普勒频率( )表征移动的快慢,图1给出 了在10dB的信噪比下两种移动环境( /f,=0.01和0.1)一 制噪声: (k)= //p_zp’ 1] u[N rrr -N ̄, . 一次256点快速傅立叶变换(FFr)实现的信号频谱幅度及其自 相关示意图。从频谱图中可以看到,信号的频谱能量集中在 (8) 最后再用IFFr把 ,(k)变换成时域信道估计值: 维普资讯 http://www.cqvip.com 798 信号处理 第23卷 (m)=IFFTt z,( ) ’FFr… 磊 m=0,1…Ⅳ r一1. ( )e 一, (9) r1一 ) (O) Wt—N +N} (一 r+2N,一1) (0) ( ."一2N,+1) ( ,, 一 ,) 由于对时域信号的截断操作,估计出的h(m)两端误差较大, 我们对h(m)只取中间的 点作为最终的估计值,即输出 ( 一1) [.] 表示转置。用LS准则估计c : C。=W H h (m)=h(m),m∈[Ⅳ r/2一 /2,Ⅳ r/2+ /2]。余下 数据的处理,即从 (Ⅳ /2+Ⅳr/2+1)开始,又取Ⅳ 点进 其中[・] 表示对方阵求逆。得到c。后,代入(12)式就可以 行完全同上所述一个周期的操作,直至导频序列最后一点的 信道估计值输出。 图2速度自适应信道估计方法一次处理周期的实现框图 3.2频谱泄漏问题的解决 由于这种频域估计方法需要在Frr前先对时域信号进 行截断,必然会造成频谱能量的泄漏[10]。实际得到的真实 频域信号不再是带宽受限于最大多普勒频率,而是随旁瓣弥 散在整个频带范围。虽然可以对截取的时域信号加窗函数 (如升余弦窗)来降低旁瓣高度,使能量尽量集中在频谱的主 瓣内,但是泄漏还是不可避免。因此本文提出一种在频域迫 零之后对泄漏的频谱进行补偿的方法。时域信号 (n)截取 后加窗函数 (n)的表达式为: (n)=r(n)・ (n)=[h (n)+n (n)]・ (n).(10) (n)变换到频域,令日 ( ): { (n)}, ( ): { (n)}, ( ):FF { (n)},Cp( )=FF {h (n)}, ( ) = { (n)},o表示卷积运算,则在 e[0,Ⅳ胛一1]上有: 1 ( )=。 ( )o ( ) 、FFT 1 =。 1[ ( )o ( )+ ( )o ( )]. (11) 、F●1 瑞利信号h (n)是带限的,因此c ( )只在区间 e[0,N,一 1]u[ , 一Ⅳ,, 一1]内有非零值。但带限的 ( )卷积 了分布于整个频带的 ( ),使信号能量泄漏于最大多普勒 频率之外,我们的目的就是恢复出在 e[Ⅳ,,Ⅳ 一Ⅳ,一1] 被迫零的信号: 日。 ( )= ( )o ( ). (12) 在 ∈[0,N,一1]u[Ⅳ 一Ⅳ,,Ⅳ 一1]内,记带内噪声 Ⅳ ( )= ( )ow(k),于是: H. =W…C+N (13) 其中: H. =[Hi (0)… (N/一1), (Ⅳ,盯一Ⅳ,)… ( 盯一1)] , C =[ (0)… (N,一1), (Ⅳ r一Ⅳ,)… (N, 一1)] , N =[ (0)… (N/一1), (Ⅳ 一Ⅳ『)... (Ⅳ一一1)] , 计算出 .( ),写成矩阵的形式即: H :w C。, (15) 其中,H一。 .=[屁 ,(Ⅳ,)… (_ ̄即 一N/一1)] , W ==[IL ( N二r ,一; ,一1 ) :‘ :( m_I一)2 ^,) 一 (1一; 1 ) :: ‘一 ’(一 ,) ]Jl ・ 在频域迫零后加上恢复出的 ( )即可估计出更准确的信 号频谱: H ((  ):)=J』l  ‘H (  )、' E ’Ⅳ,一 u Ⅳ =-, E[,v,,,v,,r—N,一1] Ⅳ”Ⅳ”r一 .(.(、 16) 4仿真结果 用计算机进行模特卡罗仿真,参数设置如下:载波频率 2GHz,采用Jakes模型产生瑞利衰落信道,导频信号的传送率 为 =2K//z, 取值范围在[0.003,0.01,0.05,0.1,0.2, 0.3],一次Frr点数为Ⅳ =512,一次估计周期输出Ⅳr=1/2 Ⅳ 点,仿真数据总长度102400点。 仿真实验一 在三种移动环境和三种信噪比(SNlt)条 件下,将本文提出的通过搜索频谱自相关图的次高峰出现位 置来确定频域迫零的边界Ⅳ『与文献[3]中提出的在已知 的理想情况下计算出的Ⅳ,值(N,= Ⅳ + 1,其中过 渡带F :(3 ̄g一0.5)J3 Ⅳ )相比较。整个仿真过程中 找出的位置序号的平均值如表1所示。可见本方法完全不需 要知道速度或多普勒频率就能比较准确地找到Ⅳ,值,尤其在 高信噪比下已经非常接近已知 时的理论计算值,体现出良 好的速度自适应性能。文献[3]采用迭代和长时统计量LCR 来估计. 以计算Ⅳ,,在初始化时要求信道至少在50次处理周 期中不变(一次迭代用10 点统计LCR,迭代5次达到收敛), 收敛后仍需要信道至少在10次处理周期中不变,而本文提出 的方法只要求信道在1次处理周期中不变,实时性大大增强。 {a,j 0.oo1 0.01 0.1 已知 ,计算出的Ⅳ, 2 6 57 SNR=0dB,相关找出的Ⅳ, 3.8 8.5 74.3 SNR:10dB,相关找出的Ⅳ, 2.3 6.6 64.1 SNR=20dB,相关找出的Ⅳ, 2.1 6.2 60.0 表1 512点FFTr搜索次峰的平均位置 仿真实验二在SNR=10dB和20dB下,按照图2中实 线的流程仿真了本文提出的频谱相关法信道估计的MSE性 能,结果如图3所示(图中注为“未知 ”)。为了比较,还做 维普资讯 http://www.cqvip.com 第5期 利用频谱特性的速度自适应信道估计方法 799 出已知, 的理想情况下的结果(图中注为“已知 ”)和不进 行任何处理直接最小均方估计(Ls)的结果(图中注为 “Ls”)。从仿真结果可以看出,在较大范围内改变多普勒频 率值,Ls估计性能主要受噪声影响,因此相同SNR时性能几 乎都处于很差的水平上不变;本文提出的方法估计性均能保 持比较接近已知. 时可以达到的性能。无需专门的多普勒 频率或速度估计模块和长时统计量是本方法的优势,在降低 复杂度的同时MSE性能比较接近多普勒频率精确知道的理 保持在较低的水平。该方法在频域寻找频谱自相关峰值出 现的位置来确定频域滤波器的通带长度,去掉带外的噪声分 量。从仿真结果可以看出,该信道估计方法性能接近已知多 普勒频率的理想情况。相对于传统速度自适应方法,该方法 的特点是不需要预先估计移动速度或多普勒频率,降低了系 统实现的复杂度;也不需要长时统计量,实时性较高。在频 域的信道估计会有频谱泄漏的问题,采用本文提出的一种恢 复泄漏信号的简单方案后,无论是对于多普勒频率先验已知 还是对于本文提出的在未知多普勒频率下的估计,性能都有 想情况。本方法的特点是仅需要作相关运算,许多DSP芯片 都有快速相关运算指令,因此本文方法易于实际实现。 仿真实验三在实验二的基础上加上本文提出的解决 频谱泄漏的方案,在SNR=10dB时的结果如图4所示。在已 知多普勒频率下,加上谱泄漏的补偿(注为“已知 +补偿”) 比不加时性能明显提高,同样效果也体现在不知道多普勒频 率下,采用本文提出的相关法确定频域带宽,加上谱泄漏的 补偿后(注为“未知, +补偿”)性能也比不加时有所提高,并 且几乎和已知, 的理想情况性能相同。可见采用本文提出 的频谱泄漏补偿方法后,由于减小了频域迫零对信号造成的 失真,因此能改善估计的性能。 图3 频谱相关法与理想情况的MSE性能比较 图4补偿频谱泄漏与否的MSE性能比较 5结论 本文提出一种新的对速度自适应、低复杂度的信道估计 方法,从步行到高速移动的速度范围内,估计的均方误差都 所提高。本文提出的信道估方法适用面较广,可以直接应用 于各种窄带系统;对于宽带CDMA系统,可以应用在每个 Rake支路上;对于宽带OFDM系统,则可以应用于各个平坦 的导频子载波上,或者将其IFrr转换到多径域后用于能量最 大的几个径[91上。 参考文献 『11 Hyuk Jun Oh and John M.Ciofif,“An adaptive channel es— timation scheme for DS—CDMA systems,”IEEE VTC’00 [C],vo1.6,PP.2839—2843,Sept.2000. [2] Zhou Chunhui,Meng Lin,Li Gang,Wang Jing,Li Wei— dong.“An order—adaptive filter for channel estimation in ds—cdma mobile systems,”ICCT’00[C],vo1.1,PP.355— 358,Aug.2000. [3]Takki Yu,Gooh ̄run Park,Daesik Hong and Changeon Kang, “A velocity—adaptive channel estimation scheme using the simple zero—forcing technique in the ̄equency domain,” IEEE VTC 2005一Spring[C],vo1.30,PP.557-561,30 May一1 June 2005. [4] K.D Anim—Appiah,“On generlaized covariance—based ve— locity estimation,”IEEE Trans.Veh.Tech.[J],vo1.48, PP.1546—1557,Sept.1999. 5 l A.Sampath and J.M.Holtzman,“Estimation of maximum doppler ̄equency ofr handof decisions,”IEEE VTC[C], PP.859—862,May 1993. [6]Goohyun Park,Daesik Hong,Changeon Kang,“Level cross— ing rate estimation with doppler adaptive noise suppression technique in ̄equency domain,”IEEE VTC’03[C].vo1. 2,PP.i192—1195,Oct.2003. [7]Chengshan Xiao,Karl D.Mann,and Jan C.Oliver,“Mobile speed estimation for TDMA--Based Hierarchical cellular sys-- tems,”IEEE Trans.Veh.Tech.[J],vo1.50,No.4,July 2001. [8] John G.Proakis,Digitla communications[M],Pub}ishing House of Electronics Industry,2001. 9 l J.一J.Van de Beck,O.Edfors,M.Sandell,S.K.Wilson,and P.O.Borjesson,“On channel estimation in OFDM sys— tems.”in Proc.45 IEEE Vehicular Technology corff. [C],Chicago,IL,July 1995,PP.815—819. [1O]Alan V.Oppenheim and Ronald W.Schafer,Discrete—time signal processing[M],Prentice Hall,1989. 作者简介 何婷玉:女,汉族,1982年生,四川成都人,硕士研究生。 研究方向为无线通信中MIMO,OFDM,CDMA系统的信道估 计与均衡技术。 

因篇幅问题不能全部显示,请点此查看更多更全内容