張 曄,吉 磊,熊 剛
(中國(guó)電子科技集團(tuán)公司第三十研究所,四川 成都 610041)
隨著通信技術(shù)的迅猛發(fā)展,電磁空間環(huán)境日益復(fù)雜,頻譜資源更加緊張。正交頻分復(fù)用(Orthogonal Frequency Division Multiplexing,OFDM)信號(hào)由于具有頻譜利用率高、抗多徑能力好、可有效抑制符號(hào)間干擾以及分配信道資源靈活等優(yōu)點(diǎn),已經(jīng)在民用和軍事通信系統(tǒng)被廣泛采用[1-2],如數(shù)字音視頻廣播業(yè)務(wù)系統(tǒng)、無(wú)線局域網(wǎng)系統(tǒng)、4G 移動(dòng)通信系統(tǒng)、無(wú)人機(jī)測(cè)控通信系統(tǒng)和美國(guó)陸軍戰(zhàn)術(shù)作戰(zhàn)人員信息網(wǎng)(Warfighter Information Network-Tactical,WIN-T)等,且被運(yùn)用于新一代無(wú)線網(wǎng)絡(luò)通信標(biāo)準(zhǔn)。此外,信號(hào)處理及分析技術(shù)的研究也受到了眾多學(xué)者和機(jī)構(gòu)的關(guān)注。在非協(xié)作通信情況下,對(duì)參數(shù)的估計(jì)識(shí)別是實(shí)現(xiàn)后續(xù)解析的前提,能更好地提取到有用參數(shù)信息。OFDM 信號(hào)具有多種重要參數(shù),包括循環(huán)前綴(Cyclic Prefix)長(zhǎng)度、載波頻偏、定時(shí)誤差和過(guò)采樣因子等。通過(guò)估計(jì)算法獲得這些OFDM 參數(shù),可為實(shí)現(xiàn)同步、均衡和解調(diào)等提供依據(jù)。
對(duì)OFDM 信號(hào)的參數(shù)估計(jì)識(shí)別方法一般包括有預(yù)處理、特征提取與參數(shù)計(jì)算估計(jì)等步驟。其中,預(yù)處理主要是對(duì)采集的信號(hào)樣本進(jìn)行正交下變頻、濾波和平滑等;特征提取是指在時(shí)域、頻域或其他一些變換后分析OFDM 信號(hào)性質(zhì),可提取出多種有用特征量,從而實(shí)現(xiàn)對(duì)信號(hào)參數(shù)的計(jì)算估計(jì)。一些研究人員和學(xué)者對(duì)OFDM 信號(hào)的參數(shù)估計(jì)識(shí)別開(kāi)展了研究。其中:Ishii 和Wornell 提出了一種參數(shù)估計(jì)算法,主要通過(guò)OFDM 信號(hào)幀中的大量符號(hào)累積實(shí)現(xiàn)計(jì)算,且需要導(dǎo)頻等先驗(yàn)信息[3];文獻(xiàn)[3-4]是在已知采樣頻率的假設(shè)下,采用相關(guān)法來(lái)研究OFDM 系統(tǒng)的系統(tǒng)參數(shù),其中文獻(xiàn)[4]提出的方法在符號(hào)定時(shí)未同步的情況下性能并不理想。
針對(duì)上述問(wèn)題,本文在信號(hào)循環(huán)平穩(wěn)性分析的基礎(chǔ)上,提出了一種基于循環(huán)譜的OFDM 信號(hào)參數(shù)識(shí)別改進(jìn)思路,抗噪性較強(qiáng),利于工程實(shí)現(xiàn)。該方法通過(guò)構(gòu)建OFDM 信號(hào)循環(huán)平穩(wěn)特性分析對(duì)應(yīng)的計(jì)算模型,準(zhǔn)確提取了包括OFDM 循環(huán)前綴長(zhǎng)度、子載波數(shù)目、過(guò)采樣率、載波頻率以及定時(shí)偏差等多個(gè)重要參數(shù),在對(duì)子載波識(shí)別運(yùn)算過(guò)程中利用循環(huán)相關(guān)性實(shí)現(xiàn)了優(yōu)化處理,并進(jìn)行OFDM 信號(hào)的循環(huán)前綴與時(shí)頻域特征融合分析,提高了算法在低信噪比條件下的穩(wěn)健性能。
圖1 表示了OFDM 信號(hào)的基本傳輸原理。將信號(hào)發(fā)送端的串行數(shù)據(jù)流經(jīng)過(guò)串-并轉(zhuǎn)化變?yōu)椴⑿袛?shù)據(jù),然后把對(duì)應(yīng)于各子載波上的數(shù)據(jù)符號(hào)分別進(jìn)行星座映射調(diào)制后再進(jìn)行傅里葉逆變換IFFT,并加上循環(huán)前綴(Cyclic Prefix,CP)形成OFDM 調(diào)制信號(hào)。對(duì)于接收方,把收到的OFDM 信號(hào)進(jìn)行去循環(huán)前綴、同步解調(diào)以及并-串轉(zhuǎn)換等處理,還原為比特流信息數(shù)據(jù)。
圖1 OFDM 信號(hào)傳輸原理
信號(hào)x(t)的二階循環(huán)平穩(wěn)特性可以利用時(shí)變相關(guān)函數(shù)式表示:
式中,R(·)表示相關(guān)函數(shù),E(·)表示數(shù)學(xué)期望,t表示周期性時(shí)間變量,τ表示時(shí)延。
設(shè)M為信號(hào)數(shù)據(jù)長(zhǎng)度,在循環(huán)頻率α位置處的循環(huán)自相關(guān)函數(shù)即為Cxx(α,τ)。它實(shí)質(zhì)上是Rxx(t,τ)關(guān)于時(shí)間t的傅里葉級(jí)數(shù),數(shù)學(xué)表達(dá)式為:
對(duì)于OFDM 信號(hào),設(shè)x(t)表示接收機(jī)收到的信號(hào)。經(jīng)過(guò)采樣,下變頻和濾波處理后可表示為:
式中:s(m)表示循環(huán)前綴為CP 的OFDM 基帶信號(hào);Tb表示信號(hào)在時(shí)域上每個(gè)采樣點(diǎn)的持續(xù)時(shí)間長(zhǎng)度,即采樣時(shí)間間隔;h(t)表示系統(tǒng)傳輸函數(shù),主要由接收濾波器響應(yīng)和信道衰落函數(shù)組成;w(t)表示方差為的加性高斯白噪聲。
為簡(jiǎn)化推導(dǎo)過(guò)程,可假設(shè)基帶傳輸信號(hào)包絡(luò)為矩形脈沖且信道為高斯信道。設(shè)過(guò)采樣因子可表示成q取為足夠大的數(shù)值,Tb和分別為變換前后的采樣間隔。已接收到的采樣信號(hào)序列x(n)是通過(guò)采樣頻率為的過(guò)采樣接收信號(hào)得到的,即:
式中,s(m)表示信號(hào)序列,可歸一化為1,g(l)為矩形脈沖,可表示為:
當(dāng)q為整數(shù)時(shí),在τ=11處的離散時(shí)變相關(guān)函數(shù)Rxx(n,τ)是周期函數(shù),周期為q;當(dāng)q為分?jǐn)?shù)時(shí),Rxx(n,τ)是近似周期函數(shù)。于是,有:
式中,l為大于1 的任意整數(shù)。接收到的信號(hào)經(jīng)過(guò)采樣后轉(zhuǎn)換為離散序列,因此可使用離散數(shù)據(jù)表示時(shí)間延遲。當(dāng)q為整數(shù)時(shí),OFDM 信號(hào)序列與噪聲是相互獨(dú)立的,信號(hào)功率為,Rxx(n,τ)是一個(gè)矩形脈沖序列,可表示為:
因?yàn)槭剑?)是一個(gè)以q為周期的函數(shù),所以當(dāng)延遲q>τ=1 時(shí),離散時(shí)變相關(guān)函數(shù)Rxx(n,τ)即為一個(gè)周期為q的函數(shù)。
實(shí)際中由于引入了循環(huán)前綴CP,OFDM 信號(hào)具有循環(huán)周期性。OFDM 信號(hào)時(shí)變相關(guān)函數(shù)s(m)包絡(luò)為周期性的矩形脈沖,其周期Tall=Ts+Tcp,其中Tcp=GTb代表循環(huán)前綴CP 的持續(xù)時(shí)間,G表示循環(huán)前綴的長(zhǎng)度,且有Ts=NTb表示OFDM 有用符號(hào)的持續(xù)時(shí)間。
若OFDM 信號(hào)以α為循環(huán)頻率,時(shí)延為τ,其對(duì)參數(shù)的估計(jì)識(shí)別數(shù)學(xué)表達(dá)式為:
2.1.1 OFDM 過(guò)采樣因子估計(jì)識(shí)別
OFDM 參數(shù)識(shí)別的第一步是提取出OFDM 過(guò)采樣因子特征。對(duì)于OFDM 信號(hào),一般都采用高采樣頻率生成波形,導(dǎo)致有效符號(hào)長(zhǎng)度中含有的采樣點(diǎn)數(shù)量大于其子載波數(shù)量即OFDM 調(diào)制解調(diào)時(shí)IFFT/FFT 點(diǎn)數(shù),原因是相比帶寬使用了比最小不失真采樣率更大的采樣率,由此引入了過(guò)采樣因子參數(shù)。
過(guò)采樣因子的大小由接收機(jī)和發(fā)射機(jī)使用的采樣頻率之間的比值決定。在OFDM 信號(hào)接收機(jī)端進(jìn)行FFT 處理前,須使得有效符號(hào)長(zhǎng)度內(nèi)采樣點(diǎn)數(shù)目和子載波數(shù)目相等。對(duì)過(guò)采樣因子參數(shù)識(shí)別的作用是為了消除OFDM 解調(diào)過(guò)程中由于過(guò)采樣帶來(lái)的影響。
當(dāng)τ<q時(shí),接收到的過(guò)采樣OFDM 信號(hào)具有循環(huán)平穩(wěn)特性。設(shè)τ=1,時(shí)延τ是固定值,在循環(huán)頻率α處的是一個(gè)復(fù)數(shù)隨機(jī)過(guò)程。
過(guò)采樣因子可表示為:
通過(guò)采用循環(huán)譜變換進(jìn)行分析,式(11)還可表示為:
式中,P表示FFT 長(zhǎng)度,表示循環(huán)譜頻率的峰值之間對(duì)應(yīng)的距離。
2.1.2 OFDM 子載波數(shù)量估計(jì)識(shí)別
OFDM 子載波數(shù)量是OFDM 重要的時(shí)域參數(shù)之一[5],傳統(tǒng)方法是通過(guò)對(duì)OFDM 信號(hào)帶寬Bw和有用符號(hào)持續(xù)時(shí)間分別進(jìn)行估計(jì)計(jì)算[6],再采用式(13)求得結(jié)果[7]。對(duì)OFDM 有用符號(hào)持續(xù)時(shí)間的計(jì)算思路可結(jié)合前述的過(guò)采樣因子,然后得到子載波數(shù)目:
為簡(jiǎn)化運(yùn)算,這里采取一種基于循環(huán)譜的改進(jìn)思路實(shí)現(xiàn)子載波數(shù)目的估計(jì)識(shí)別。
根據(jù)OFDM 信號(hào)的循環(huán)平穩(wěn)特性,若對(duì)信號(hào)的時(shí)變自相關(guān)函數(shù)Rxx(n,τ)作FFT 變換,將出現(xiàn)離散的譜線峰值。如果譜線峰值間距離為D,F(xiàn)FT 點(diǎn)數(shù)是M,則可得周期q=M/D。根據(jù)對(duì)Rxx(n,τ)在FFT變換后峰值譜線間距離的搜索計(jì)算,可求出過(guò)采樣率,然后進(jìn)一步估計(jì)出OFDM 信號(hào)的有效符號(hào)長(zhǎng)度,即:
式中:
L表示偏移長(zhǎng)度量的最大值,Ci(τ)表示相關(guān)函數(shù)值(τ∈[1,L]),Ei(τ)表示信號(hào)能量均值,Mi表示對(duì)參數(shù)估計(jì)的數(shù)目。計(jì)算有效的條件是偏移長(zhǎng)度量需要大于有效符號(hào)長(zhǎng)度。
OFDM 信號(hào)子載波數(shù)目為:
基于循環(huán)譜分析算法可實(shí)現(xiàn)對(duì)OFDM 信號(hào)多參數(shù)的識(shí)別,包括循環(huán)前綴長(zhǎng)度、載波頻偏以及定時(shí)誤差等。
因?yàn)镺FDM 信號(hào)是具有循環(huán)平穩(wěn)性的信號(hào),所以可將τ=Ts代入式(8)來(lái)估計(jì)識(shí)別相應(yīng)的循環(huán)頻率或符號(hào)持續(xù)時(shí)間,并可得出:
傳統(tǒng)的OFDM 參數(shù)估計(jì)器是通過(guò)FFT 變換后的序列系數(shù)計(jì)算來(lái)進(jìn)一步求取結(jié)果[8],且需要大量的采樣樣本點(diǎn)數(shù)以提高精度[9-10]。本文根據(jù)循環(huán)平穩(wěn)統(tǒng)計(jì)量得到了對(duì)頻率偏移和定時(shí)誤差等參數(shù)的聯(lián)合估計(jì)識(shí)別改進(jìn)方法。
首先對(duì)信號(hào)進(jìn)行降采樣,設(shè)信號(hào)的頻率偏移為v,定時(shí)誤差為ε,OFDM 信號(hào)可表示為:
由于OFDM 信號(hào)中CP 具有相關(guān)性,因此對(duì)OFDM 有用符號(hào)持續(xù)時(shí)間Ts的估計(jì)器可以表示為:
其次,將式(19)代入式(20),其中的相關(guān)部分對(duì)應(yīng)為循環(huán)頻率為0 時(shí)的循環(huán)譜值。由于引入了CP,在τ=N時(shí),相關(guān)部分轉(zhuǎn)化為:
式中,I表示CP 中k值的集合,且有:
因?yàn)樵肼暸cOFDM 信號(hào)相互獨(dú)立,且CP 部分與除去CP 的OFDM 信號(hào)部分也相互獨(dú)立,所以根據(jù)中心極限定理,Z(·)可以近似為高斯分布的噪聲。當(dāng)用于估計(jì)的OFDM 符號(hào)數(shù)目較大時(shí),ε可表示為:
頻率偏移的估計(jì)為:
以N=64、G=16、FFT 長(zhǎng)度為1 024 的OFDM信號(hào)為例,設(shè)該接收信號(hào)的有用符號(hào)共有100 個(gè),OFDM 循環(huán)前綴CP 的最大長(zhǎng)度為有用符號(hào)持續(xù)時(shí)間的1/4。圖1 表示在信噪比為5 dB 的情況下對(duì)OFDM 信號(hào)進(jìn)行循環(huán)譜變換(τ=Ts),循環(huán)譜分析中FFT 點(diǎn)數(shù)取為4 096,其中橫軸表示循環(huán)頻率,縱軸表示循環(huán)相關(guān)函數(shù)值。
圖2 改進(jìn)方法與傳統(tǒng)方法的性能比較
圖1 中可看出最左側(cè)位置有3 根較大值的譜線聚集在一起,對(duì)應(yīng)于循環(huán)頻率及其倍數(shù)。在估計(jì)出CP 長(zhǎng)度的參數(shù)信息后,還可進(jìn)一步改善估計(jì)器的性能。根據(jù)前述的式(8)、式(20)、式(24)和式(25),可得=36。
為了驗(yàn)證本文識(shí)別方法的有效性,采用MATLAB 軟件進(jìn)行仿真。仿真參數(shù):OFDM 信號(hào)的采樣率設(shè)為100 MHz,接收信號(hào)共有50 個(gè)有用符號(hào),子載波數(shù)目為N=64,循環(huán)前綴長(zhǎng)度為G=16,過(guò)采樣因子為q=5/2,定時(shí)誤差為ε=0.3,頻率偏移為v=20Tb,Monte-Carlo 仿真次數(shù)為1 000 次,噪聲為加性高斯白噪聲。
采用基于循環(huán)譜分析的方法進(jìn)行OFDM 信號(hào)參數(shù)估計(jì)識(shí)別,得出對(duì)過(guò)采樣因子、子載波數(shù)目、循環(huán)前綴長(zhǎng)度和頻率偏移、定時(shí)誤差參數(shù)估計(jì)的歸一化最小均方誤差(Normalized Mean Square Error,NMSE)。仿真結(jié)果如圖3 所示,其中橫軸表示信噪比(單位dB),縱軸表示估計(jì)誤差。
圖3 本文方法對(duì)OFDM 多參數(shù)的估計(jì)誤差曲線
從圖3 可以看出,本文的改進(jìn)方法是一種有效的OFDM 參數(shù)估計(jì)識(shí)別方法。在信噪比大于5 dB 時(shí),對(duì)OFDM 信號(hào)的過(guò)采樣因子、頻率偏移和定時(shí)誤差等參數(shù)估計(jì)的歸一化最小均方誤差都小于10-2。
分別對(duì)傳統(tǒng)基于最大似然判決的方法和本文改進(jìn)方法進(jìn)行估計(jì)識(shí)別性能的仿真對(duì)比分析,驗(yàn)證對(duì)OFDM 信號(hào)多參數(shù)的平均識(shí)別正確概率結(jié)果。條件與前述的仿真實(shí)驗(yàn)1 相同,仿真結(jié)果如圖4 所示,橫軸表示信噪比,縱軸表示不同估計(jì)識(shí)別方法的平均正確率。
由仿真結(jié)果可以得出,新方法比傳統(tǒng)的估計(jì)識(shí)別方法具有更好的性能,適用于認(rèn)知無(wú)線電系統(tǒng)和其他非協(xié)作通信等應(yīng)用對(duì)OFDM信號(hào)進(jìn)行分析處理。
圖4 改進(jìn)方法與傳統(tǒng)方法的性能比較
隨著通信系統(tǒng)的不斷發(fā)展,OFDM 信號(hào)的應(yīng)用越來(lái)越廣泛,作用越來(lái)越重要,使得所處的無(wú)線環(huán)境日益復(fù)雜,給OFDM 參數(shù)估計(jì)識(shí)別算法的研究提出了更高要求。本文提出了一種基于循環(huán)譜分析的OFDM 參數(shù)估計(jì)識(shí)別改進(jìn)方法。仿真試驗(yàn)表明,該方法比傳統(tǒng)識(shí)別算法性能更優(yōu),穩(wěn)定性更好。今后將進(jìn)一步拓展,為解決OFDM 信號(hào)參數(shù)估計(jì)識(shí)別問(wèn)題提供更有效的手段,為無(wú)線電監(jiān)測(cè)、智能化通信等領(lǐng)域的研究貢獻(xiàn)力量。