OFDM水聲通信信道估計(jì)技術(shù)研究

　　水聲信道是一個(gè)十分復(fù)雜的時(shí)-空-頻變信道，其主要特征是復(fù)雜性、多變性、強(qiáng)多途和有限帶寬。聲傳播損失和海水吸收損失使得水聲信道帶寬受到極大限制，海洋水聲信道中多徑效應(yīng)的存在造成接收信號(hào)的畸變和嚴(yán)重的碼間干擾，給水聲通信系統(tǒng)的設(shè)計(jì)帶來了巨大的困難，信道中的相位起伏使得載波恢復(fù)和相干解調(diào)變得十分困難。在常用的高速水聲通信技術(shù)中，采用相位相干(PSK/QAM)調(diào)制要面對(duì)信道起伏時(shí)的相干解調(diào)問題，而且要適應(yīng)收發(fā)端相對(duì)運(yùn)動(dòng)所帶來的多普勒頻移。OFDM作為一種可有效對(duì)抗碼間干擾、頻譜利用率高的高速傳輸系統(tǒng)，引起人們廣泛關(guān)注。作為OFDM技術(shù)的關(guān)鍵之一，信道估計(jì)的好壞直接影響整個(gè)系統(tǒng)的性能。目前常用的信道估計(jì)算法主要有導(dǎo)頻信道估計(jì)和盲信道估計(jì)兩種。本文主要研究在具有導(dǎo)頻插入情況下常用的OFDM信道估計(jì)方法，并基于相同的導(dǎo)頻圖案下，對(duì)不同的方法進(jìn)行分析比較。

　　1 OFDM水聲通信系統(tǒng)模型

　　根據(jù)水聲信道的特點(diǎn)，以及所傳數(shù)據(jù)的一些參數(shù)要求，給出了OFDM水聲通信系統(tǒng)的模型，如圖1所示。

　　這里有如下假設(shè)：已使用了循環(huán)前綴;信道沖擊響應(yīng)時(shí)間小于循環(huán)前綴;發(fā)送機(jī)和接收機(jī)完全同步;信道噪聲是復(fù)的加性高斯白噪聲。輸入的比特序列經(jīng)過分組，根據(jù)采用的調(diào)制方式，完成相應(yīng)的調(diào)制映射，在通過串并變化和插入導(dǎo)頻信息后，形成信息序X(k)，對(duì)X(k)進(jìn)行IFFT，計(jì)算出OFDM已調(diào)信號(hào)的時(shí)域抽樣序列，加上循環(huán)前綴CP(循環(huán)前綴可以使OFDM系統(tǒng)消除信號(hào)的多徑時(shí)延造成符號(hào)問干擾(ISI)和載波間干擾(ICI)，再作D/A變換，得到OFDM已調(diào)信號(hào)的時(shí)域波形。接收端先對(duì)接收信號(hào)進(jìn)行A/D變換，去掉循環(huán)前綴CP，得到OFDM已調(diào)信號(hào)的抽樣序列，對(duì)該抽樣序列作DFT得到原調(diào)制信息序列X(k)。

　　從接收方的角度看，當(dāng)循環(huán)前綴的時(shí)間大于信道沖擊響應(yīng)時(shí)間，可以將與信道線性卷積轉(zhuǎn)化為圓周卷積，可以得出OFDM系統(tǒng)的表達(dá)式：

　　其中，Y(k)長度為N的接收序列，X(k)是長度為N的發(fā)送序列，h(n)為信道沖擊響應(yīng)(不足N長的部分添零補(bǔ)足)，ω(n)為信道噪聲，H(k)為信道的頻域響應(yīng)。

　　2 OFDM水聲通信信道估計(jì)

　　OFDM信道估計(jì)方法大致分為盲信道估計(jì)和非盲(基于導(dǎo)頻的)信道估計(jì)2類，盲信道估計(jì)不需要傳送導(dǎo)頻信息，因此可以提高有效數(shù)據(jù)的傳送效率，缺點(diǎn)是需要收集大量的數(shù)據(jù)作為運(yùn)算的依據(jù)，運(yùn)算量大，因此在時(shí)變信道中較難實(shí)現(xiàn)。另一類是基于導(dǎo)頻方式的估計(jì)，此類信道估計(jì)常用的方法有兩種：基于導(dǎo)頻信道的估計(jì)和基于導(dǎo)頻符號(hào)的估計(jì)?；趯?dǎo)頻信道的方法是在系統(tǒng)中設(shè)置專用導(dǎo)頻信道來發(fā)送導(dǎo)頻信號(hào)。由于OFDM系統(tǒng)具有時(shí)頻二維結(jié)構(gòu)，因此采用導(dǎo)頻符號(hào)輔助信道估計(jì)更加靈活。所謂的基于導(dǎo)頻符號(hào)的信道估計(jì)是指在發(fā)送端的信號(hào)中的某些位置插入一些接收端已知的符號(hào)或序列，接收端根據(jù)這些信號(hào)或序列受傳輸衰落影響的程度利用某些算法來估計(jì)信道的衰落性能。這一技術(shù)叫做導(dǎo)頻信號(hào)輔助調(diào)制(PSAM)。

　　基于導(dǎo)頻OFDM的信道估計(jì)算法的基本過程是：在發(fā)送端適當(dāng)位置插入導(dǎo)頻，接收端利用導(dǎo)頻恢復(fù)出導(dǎo)頻位置的信道信息，然后利用某種處理手段(如內(nèi)插，濾波，變換等)獲得所有時(shí)段的信道信息。這里涉及到3個(gè)主要問題，這也是目前OFDM的非盲估計(jì)算法研究的3個(gè)方向：發(fā)送端導(dǎo)頻的選擇與插入;接收端導(dǎo)頻位置信道信息獲取的方式;通過導(dǎo)頻位置獲取的信道信息如何較好地恢復(fù)出所有時(shí)刻信道的信息。

{{分頁}}

　　由式(1)知：

　　從Y(k)抽取導(dǎo)頻信號(hào)Yp(k)，而發(fā)送的導(dǎo)頻信號(hào)Xp(k)是知的，因此可得到導(dǎo)頻信號(hào)位置的信道響應(yīng)估計(jì)值Hp(k)：

　　在得到導(dǎo)頻信號(hào)位置的信道傳輸函數(shù)的估計(jì)值Hp(k)后，數(shù)據(jù)位置的信道響應(yīng)可通過相鄰導(dǎo)頻信號(hào)信道響應(yīng)內(nèi)插獲得，信道估計(jì)的總體框架如圖2所示。

　　在接收端，輸入倍號(hào)經(jīng)過FFT變換后得到數(shù)據(jù)Y(k)，從Y(k)中抽取導(dǎo)頻位置的數(shù)據(jù)Yp(k)，根據(jù)已知的導(dǎo)頻數(shù)據(jù)Xp(k)，得到導(dǎo)頻位置的信道響應(yīng)估計(jì)值Hp(k)，再通過時(shí)頻方向上的內(nèi)插，得到所有位置的信道響應(yīng)的估計(jì)值H(k)。

　　2.1 導(dǎo)頻插入圖樣的選取

　　基于導(dǎo)頻符號(hào)的信道估計(jì)算法中，導(dǎo)頻插入方案可以分為2類：塊狀導(dǎo)頻分布和梳狀導(dǎo)頻分布。導(dǎo)頻塊狀分布的OFDM系統(tǒng)中，將連續(xù)多個(gè)OFDM符號(hào)分成組，每組中的第一個(gè)OFDM符號(hào)發(fā)送導(dǎo)頻信號(hào)，其余的OFDM符號(hào)傳輸數(shù)據(jù)信息：梳狀導(dǎo)頻的系統(tǒng)中，將N個(gè)子信道均勻地分為M組，在每一組的第一個(gè)子載波中傳輸導(dǎo)頻符號(hào)，其余的M-1個(gè)子載波傳輸數(shù)據(jù)信息。兩種插入方案如圖3和圖4所示。其中，黑色代表導(dǎo)頻符號(hào)，白色代表數(shù)據(jù)符號(hào)。

　　塊狀導(dǎo)頻插入方式往往只在開始發(fā)送一些訓(xùn)練符號(hào)，估計(jì)出導(dǎo)頻符號(hào)處的信道信息將作為以后所有時(shí)刻的信道信息，直到下一個(gè)含有導(dǎo)頻信息的符號(hào)到來，這就要求信道在相當(dāng)長的時(shí)間內(nèi)變化較小，甚至不變，即所謂準(zhǔn)靜止信道、慢衰落信道(不考慮或者只考慮較小的多普勒頻移)。由此可見，塊狀導(dǎo)頻方式較適合于恒參信道、WLAN信道等。而梳狀導(dǎo)頻分布在不同的OFDM符號(hào)中，因此能夠較好地跟蹤不同符號(hào)下信道狀態(tài)的變化，特別是在信道快變化的條件下這種優(yōu)勢(shì)更加明顯。如果信道中有較大的多普勒頻移，信道變化較快，則不能選用塊狀導(dǎo)頻插入方式，而應(yīng)在整個(gè)信號(hào)的時(shí)頻空間內(nèi)插入梳狀導(dǎo)頻信號(hào)。

　　由于水聲通信信道具有頻率選擇性和時(shí)變特性，也就是說實(shí)際OFDM系統(tǒng)中信道的傳輸特性在時(shí)域和頻域內(nèi)都是時(shí)變的，即使在傳輸一幀OFMD數(shù)據(jù)時(shí)間內(nèi)，信道特性也會(huì)發(fā)生明顯的變化，因此應(yīng)該采用基于間隔式導(dǎo)頻插入，即梳狀導(dǎo)頻的信道估計(jì)方法，跟蹤不同符號(hào)下信道狀態(tài)的變化，對(duì)信道做動(dòng)態(tài)估計(jì)。

{{分頁}}

　　2.2 基于梳狀導(dǎo)頻的信道估計(jì)算法

　　OFDM系統(tǒng)可以被視為多個(gè)高斯信道，如圖5所示。

　　把式(2)寫成矩陣形式：

　　如果信道時(shí)域矢量h是高斯的并且與信道噪聲W不相關(guān)，則對(duì)h的LS估計(jì)可表示為：

　　MMSE信道估計(jì)器雖然在實(shí)際應(yīng)用當(dāng)中的性能較好，但是也存在一定的弱點(diǎn)，即他計(jì)算相當(dāng)復(fù)雜，并依賴于信道的統(tǒng)計(jì)特性，這樣就阻礙了他的應(yīng)用。文獻(xiàn)[2]給出了MMSE估計(jì)器和LS估計(jì)器之間的關(guān)系，即LMMSE估計(jì)器：

　　LMMSE算法可以降低計(jì)算復(fù)雜度，但仍需對(duì)階數(shù)等于子載波個(gè)數(shù)的復(fù)數(shù)矩陣求逆，隨著子載波數(shù)的增加，其計(jì)算量十分可觀，因此在實(shí)際應(yīng)用中，LMMSE算法往往受到限制。采用 SVD(奇異值分解)方法可以對(duì)LMMSE算法進(jìn)一步簡(jiǎn)化，即所謂的LRLMMSE。具體方法如下：

　　對(duì)頻域信道向量自協(xié)方差矩陣RHH進(jìn)行奇異值分解，可以得到：

　　式中U是由歸一化正交向量組成的矩陣，A是一個(gè)對(duì)角矩陣，他的對(duì)角線包含了RHH的奇異值λ0≥λ1≥…≥λN-1≥0，將式(9)代入式(8)，可以得到基于奇異值分解的信道估計(jì)器為：

　　式中，△J為△的左上角J