第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系统的信道估 计与均衡技术。
因篇幅问题不能全部显示,请点此查看更多更全内容