OFDM移動通信系統中的最大多普勒估計

孟建剛，陳建春

(西安電子科技大學電子工程學院，陜西西安 710071)

多普勒頻移是影響OFDM系統性能的重要參數之一，OFDM系統的多個環節如編碼、調制、信道估計等都需要多普勒頻移信息，因此多普勒頻移估計將對整個OFDM系統的性能產生較大影響。目前大量的多普勒頻移估計算法被提出［1－3］。文獻［1］提出了一種基于循環前綴的自相關算法，該算法結構簡單，但是在低速和低信噪比情況下，估計精度差。文獻［2］提出了一種基于信道估計的自相關算法，在低速情況下有較好的估計精度，但是其估計精度主要取決于信道估計的精度。文獻［3］中提出了一種基于最大似然函數的估計方法，其估計的精度很高，但實現的復雜度也大，難以應用于實際的系統當中。本文結合文獻［1］與文獻［2］中算法的特點提出了一種改良算法。通過在OFDM符號中插入已知的OVSF序列，利用OVSF序列的特點簡便地求得信道的沖激響應，然后通過信道沖激響應的自相關，獲得多普勒頻移信息。仿真結果表明，本文提出的多普勒頻移估計算法可以有效地估計出多普勒頻移。

1 算法原理描述

Clarke［4］信道模型是用于描述小尺度衰落的一種平坦衰落信道模型，即瑞利衰落信道。其移動臺接收信號強度的統計特性是基于散射的，這正好與市區環境中無直視通路的特點相吻合，因而廣泛應用于市區環境的仿真中。因此，信道具有Jakes［5］功率譜模型

式中，hi是第i路的信道沖激響應;是信號的平均功率;fd是多普勒頻移，它與速度v的公式為fd=fcv/c，這里fc為載波頻率;c是光速，則信道沖激響應的歸一化自相關函數由式(1)可得

這里J0是第一類零階Bessel函數，則

由此可知，通過信道沖激響應的自相關函數，運用式(3)就可得到最大多普勒頻移。

2 系統模型

2.1 多普勒頻移估計算法

假定接收信號r(t)為復基帶信號，經過瑞利衰落信道，在加性高斯白噪聲干擾下，其離散形式可以表示為

式中，a(k)為發送信號;x(k)和y(k)分別為a(k)經過信道以后期望信號的實部和虛部，且E［x2(k)］=E［y2(k)］=η2/2。n(k)=ni(k)+jnq(k)為復高斯加性白噪聲且與信號不相關，其實部和虛部的功率為

接收信號的自相關函數為

J0為第一類零階Bessel函數。

如圖1所示，在OFDM符號前面插入循環前綴，L為循環前綴長度，K為一個OFDM符號的子載波個數，T為采樣時間，則一個包含循環前綴的OFDM符號的持續時間為Ts=(K+L)T。

圖1 OFDM符號結構

圖2是該算法的結構圖，ri(k)和rq(k)分別為接收信號的實部和虛部，ρi(KT，fd)和ρq(KT，fd)分別為實部和虛部所對應的自相關函數。

圖2 算法結構圖

本算法中不需要整個循環前綴GI都參加自相關運算，M的取值根據信道環境的不同而不同，在文獻［1］中有關M取值的討論，在此不作詳細說明。利用式(5)可以求出最大多普勒頻移fd。

2.2 改良的多普勒頻移估計算法

如圖1所示，在原有的循環前綴位置插入OVSF序列，OVSF序列的持續時間為Tg=LT，有效數據持續時間為Td=KT，T為采樣時間，K為一個OFDM符號的子載波數目，L為插入的OVSF序列的長度，完整OFDM符號的持續時間為Ts=(K+L)T。

假設信號經過瑞利衰落信道，在理想的定時同步情況下，第i路徑的接收信號可以表示為［6］

式中，ci(k)是 OVSF序列的元素;k∈［1，L］，hi(k)是信道的沖激響應;w(k)是0均值，方差為σ2的復高斯白噪聲，σ2取值在不同的路徑下可能不同。

由于信號經過瑞利衰落信道多徑傳輸，各個路徑的能量可能不同，那么也就不同。為方便起見，這里假定hi就是由信號能量最強路徑得到的信道沖激響應，根據文獻［2］，由式(6)可以得到hi為

在文獻［1］的算法原理基礎上，進行信道沖激響應的自相關，該算法結構如圖3所示。

圖3 多普勒頻移估計算法結構圖

hi，I，hi，Q是 hi的實部和虛部，N 是參加多普勒頻移估計的符號數量，mTs是相關時間間隔，m是一個常數，它的取值將在下面討論，則歸一化的自相關函數可以表示為

由此，通過式(3)和式(8)就可以求出最大多普勒頻移fd。

3 仿真結果分析

在瑞利衰落信道下進行Matlab仿真，選取6路多徑衰落，延遲的系數分別為 －0.2 μs，0 μs，0.6 μs，1.2 μs，和 5.5 μs，相 對 應 的 衰 落 系 數 分 別 為－3 dB，0 dB，－2 dB，－5 dB，－8 dB和－10 dB，選取數據長度K=3780，OVSF序列長度L=420，進行仿真的OFDM符號數量N=2000，常數P=6.5。符號率為 7.56 Mbit·s－1，那么 T=1/7.56 μs。調制方式為16QAM，載波頻率fc=600 MHz。

算法中用到的第一類零階Bessel函數是一個非線性函數，J0(x)的一個值可能對應著 x的多個值，因此必須把限定在第一個單調區間內，如圖4所示，則

假定高速移動的火車速度達到350 km·h－1，那么多普勒頻移fd=fcv/c將接近200 Hz，Ts=(L+K)T=555.6 μs，因此根據式(9)，m取4最合適。J0(x)的取值可以由查表法求得。

圖4 第一類零階Bessel函數

圖5是不同信噪比下本算法的仿真結果，選取的多普勒頻移范圍在0～200 Hz之間，仿真結果表明，本算法在低速情況下有較高的精度，對信噪比的變化不是特別敏感，特別在40 Hz附近的精確度高于其他部分，在高速情況下也有良好的估計精度。圖6是頻率范圍在0～200 Hz內，不同信噪比情況下多普勒頻移的估計值與實際值平均誤差對比，由圖6可以看出，在不同信噪比下，本算法估計性能比文獻［1］和文獻［2］中算法的估計性能都好，特別在低信噪比情況下，與文獻［1］算法相比，最高時平均誤差將減少約20 Hz，與文獻［2］相比，也有較大優勢。

4 結束語

本文運用文獻［1］中算法的原理，結合文獻［2］中信道估計方法，提出了一種新的改良算法，本算法繼承了兩種方法的優點，發揮了OVSF序列的優勢。仿真結果表明，在瑞利衰落信道下，本算法有著較高的估計精度，能夠適應于較大范圍內多普勒頻移的估計，較文獻［1］與文獻［2］中的算法有著較大的優勢，而且本算法復雜度較低，可以方便地運用到實際的OFDM系統中。

［1］ CAI J P，SONG W T，LI Z.Doppler spread estimation for mobile OFDM systems in rayleigh fading channels［J］.IEEE Trans on Consum Electron，2003，49(4):973 －977.

［2］ LU Qiaoli，CHEN Wei，XIE Tao，et al.A doppler spread estimator design for mobile OFDM syestems［C］.Singapore:2008 11th IEEE Singapore International Conference on Communication Systems，2008.

［3］ CHOI Y S，OZDURAL O C，LIU H P，et al.A maximum likelihood doppler frequency estimator for OFDM systems［J］.IEEE Trans on Consum Electron，2006(10):4572－4576.

［4］ CLARKE R H.Principle of mobile communications［M］.Boston:Luwer Academic Publishers，1996.

［5］ JAKES W C.Microwave mobile communications［M］.1sted.Piscataway:IEEE Press，1993.

［6］ 韋崗，季飛，傅娟.通信系統建模與仿真［M］.北京:電子工業出版社，2007.