寧小玲，張林森，梁玥

?

一種改進(jìn)LS信道估計算法在稀疏多徑水聲信道中的應(yīng)用

寧小玲1，張林森2，梁玥3

(1. 海軍工程大學(xué)電子工程學(xué)院，湖北武漢 430033； 2. 海軍工程大學(xué)兵器工程系，湖北武漢 430033；3. 海軍工程大學(xué)教務(wù)處，湖北武漢 430033)

為了消除水聲正交頻分復(fù)用調(diào)制(Orthogonal Frequency Division Multiplexing，OFDM)系統(tǒng)中噪聲對稀疏多徑信道估計的影響，提出了一種改進(jìn)的最小二乘(Least-Square，LS)信道估計算法。該方法在傳統(tǒng)基于離散傅里葉變換(Discrete Fourier Transform，DFT)插值的信道估計結(jié)構(gòu)上進(jìn)行了改進(jìn)，得到了基于閾值探測的DFT插值信道估計方法，該方法將小于閾值的時域信道響應(yīng)置零，探測最有效信道抽頭，有效消除噪聲干擾的影響。仿真結(jié)果驗證了DFT插值在稀疏多徑水聲信道估計中的實用性；得到了當(dāng)循環(huán)前綴長度與信道長度越接近時信道估計性能(Bit Error Rate，BER)越好的結(jié)論；確定了在兩種調(diào)制方式下算法的閾值系數(shù)范圍；且該新算法與已有算法所需的信噪比可低約2 dB，解決了用循環(huán)前綴長度來近似信道真實長度的實際問題。

LS信道估計；正交頻分復(fù)用；稀疏多徑水聲信道；DFT插值；

0 引言

大量研究表明，多徑效應(yīng)是水聲信號在海洋聲信道中傳播時發(fā)生畸變的重要原因，它會導(dǎo)致信號嚴(yán)重衰落和碼間干擾。為了保證水聲通信系統(tǒng)的性能不受信道多徑和衰落的影響，實現(xiàn)高速率的水下數(shù)據(jù)傳輸，就需要對水聲信道進(jìn)行實時估計與均衡，及時調(diào)整接收處理方式，抑制稀疏多徑衰落信道的不利影響。OFDM(正交頻分復(fù)用調(diào)制)具有抗多徑干擾能力強(qiáng)、帶寬利用率高以及有效對抗符號間干擾等特性，已廣泛用于無線局域網(wǎng)、無線廣域網(wǎng)等。近年來OFDM技術(shù)已成為水聲通信領(lǐng)域的研究熱點。其中信道估計是OFDM系統(tǒng)中補(bǔ)償水下聲信道畸變的一項非常重要的技術(shù)[1-3]。

目前，基于導(dǎo)頻的信道估計一般包括三個步驟：(1) 設(shè)計導(dǎo)頻結(jié)構(gòu)；(2) 獲得導(dǎo)頻位置處的信道響應(yīng)；(3) 根據(jù)導(dǎo)頻位置處的信道響應(yīng)恢復(fù)出所有位置的信道響應(yīng)。三個步驟分別涉及到不同的技術(shù)，目前，導(dǎo)頻結(jié)構(gòu)應(yīng)用相對較多的是梳狀導(dǎo)頻結(jié)構(gòu)[4-5]。而在對水聲信道基于導(dǎo)頻的信道估計的研究中，主要有LS估計和線性最小均方誤差(Linear Minimum Mean Square Error，LMMSE)估計兩大類，LS信道估計算法相比LMMSE估計算法簡單易行，不需要知道信道統(tǒng)計特性和系統(tǒng)噪聲統(tǒng)計特性兩個先驗信息，其在工程中應(yīng)用相對較多[6-8]。由導(dǎo)頻位置的信道響應(yīng)恢復(fù)全部信道響應(yīng)，主要依靠內(nèi)插算法?；贒FT的內(nèi)插器最大的優(yōu)點就是對輸入信號采樣率不敏感[9]，只要輸入信號的采樣率滿足Nyquist采樣定理，不論采樣率具體為多少，內(nèi)插器的性能都比較穩(wěn)定[10]，且只要獲得信道的最大時延長度，基于DFT的信道估計能夠得到很好的性能[11]?；贒FT內(nèi)插的估計方法中需要確認(rèn)的是信道長度，文獻(xiàn)[12]把循環(huán)前綴長度L近似為信道長度，這種方法雖然從一定程度上提高了信道的估計精度，但估計性能還有待進(jìn)一步改善。在關(guān)于稀疏多徑信道估計的研究中，文獻(xiàn)[13]引入多徑數(shù)建立由多徑時延、幅度系數(shù)及徑數(shù)組成的水聲信道多徑參數(shù)模型，提出一種針對單頻水聲信號的多徑時延估計算法。文獻(xiàn)[14]利用基于廣義Akaike信息論(Generalized Akaike Information Criterion，GAIC)準(zhǔn)則確定稀疏信道模型的階次及抽頭位置。

本文采用基于導(dǎo)頻的LS信道估計，與上面三步驟不同的是，本文在第三個步驟前增加了一個步驟：首先確定時域信道響應(yīng)中有效抽頭數(shù)的“有效閾值”，將小于閾值的信道響應(yīng)置零，然后對新得到的有用信道響應(yīng)進(jìn)行DFT插值，恢復(fù)信道的所有響應(yīng)。這樣可探測最有效信道抽頭，大大減小了噪聲對估計性能的影響，顯著改善LS 算法的信道估計效果。

1 OFDM系統(tǒng)模型和稀疏水聲信道模型

水聲OFDM系統(tǒng)的等效基帶信號模型如圖1所示。

(1)

在高速水下數(shù)據(jù)傳輸中通常面臨的都是稀疏的多徑衰落信道，并且多徑具有可分辨性，如果能夠充分利用這種信道的稀疏性，可以提高信道估計的精度。Stojanovic[15]提出了圖2中的稀疏水聲信道模型。

那么，稀疏信道也可表示為

在頻域內(nèi)的接收信號為

(3)

2 基于LS準(zhǔn)則的DFT插值信道估計算法

2.1 經(jīng)驗的DFT插值LS估計方法

基于DFT內(nèi)插的LS信道估計原理圖如圖3所示[16]。

基于DFT內(nèi)插的信道估計方法基本步驟是：對導(dǎo)頻位置上的信道響應(yīng)作離散傅里葉逆變換(Inverse Discrete Fourier Transform，IDFT)變換，得到信道的時域沖激響應(yīng)。如果導(dǎo)頻的個數(shù)大于信道沖激響應(yīng)的采樣長度，在的高頻部分插零，然后再通過一個點的DFT，就可以內(nèi)插得到全部信道頻域響應(yīng)?；贒FT內(nèi)插的信道估計方法的前提是導(dǎo)頻的結(jié)構(gòu)必須是等間隔的。

信道估計就是利用信號的確知信息來估計出實際信道的徑數(shù)和徑的系數(shù)，就是真實徑或有用徑。由于噪聲的存在和部分導(dǎo)頻對應(yīng)的偽隨機(jī)噪聲碼相位不完全正交就會引起非真實徑，也就是噪聲徑。

圖4 文獻(xiàn)[12]信道估計結(jié)構(gòu)

Fig.4 Channel estimation structure in Literature [12]

最后，可得到基于DFT內(nèi)插的信道估計為

(6)

式(5)、式(6)從一定程度上改善了信道的估計精度，但從式(5)可以看出，離真實的信道長度有一定的差距，這是導(dǎo)致估計性能還需進(jìn)一步改善的原因之一。其次，我們也看出，在進(jìn)行DFT內(nèi)插之前的時域信道響應(yīng)是真實的時域信道響應(yīng)和噪聲的疊加，這也是估計性能需要改善的另一個原因。

2.2 改進(jìn)的LS方法

(1) 根據(jù)LS估計準(zhǔn)則得到導(dǎo)頻位置處的信道估計；

(2) 把LS信道估計出的結(jié)果通過IDFT變換到時域；

(3) 根據(jù)時域信道估計值，確定有用徑閾值，再根據(jù)閾值確定有用徑。

(4) 對得到的有用信道估計進(jìn)行DFT插值，恢復(fù)信道的所有響應(yīng)。

本文采用基于導(dǎo)頻的信號結(jié)構(gòu)，方法結(jié)構(gòu)圖如圖5所示。

改進(jìn)的LS方法整個流程如下：

第一步，基于梳妝導(dǎo)頻的LS信道估計表達(dá)式為

(8)

第三步，為了從式(8)中將信道的有用徑區(qū)別開來，有效減少噪聲的影響，本文提出基于閾值的新的LS估計方法：

式中，

(11)

第四步，被估計信道的頻率響應(yīng)為

(14)

3 仿真結(jié)果與分析

稀疏水聲信道模型2為負(fù)聲速梯度水聲信道[18]，其響應(yīng)為：

(16)

為了得到信道估計性能，假設(shè)接收機(jī)實現(xiàn)完全同步。同時，為了避免碼間干擾，選擇循環(huán)前綴大于信道的最大延遲。

3.1 不同插值方法比較

從圖6(a)和圖6(b)可以看出，在兩種信道條件下，四種插值方法的估計效果相差明顯：DFT插值方法在兩種不同信道條件下的均方誤差(Mean Square Error，MSE)性能保持穩(wěn)定，這是因為基于DFT的內(nèi)插器最大的優(yōu)點就是對輸入信號采樣率不敏感，在采樣率滿足Nyquist采樣定理的前提下，不論采樣率具體為多少，內(nèi)插器的性能都比較穩(wěn)定；Cubic插值方法受信道的影響較大，在復(fù)雜信道條件下幾乎不能完成信道估計；Spline插值和Linear插值方法受信道的影響也較大，然而兩種插值方法在不同信道條件下的MSE性能相當(dāng)。所以，仿真驗證了在實際稀疏水聲信道估計應(yīng)用中使用DFT插值方法的實用性和優(yōu)勢。

(a) 信道1, QPSK調(diào)制，L=20

(b) 信道2, QPSK調(diào)制，L=68

圖 6 兩種稀疏多徑水聲信道的LS信道估計的BER性能比較

Fig.6 BER comparison between LS channel estimations for two sparse multi-path underwater acoustic channels

3.2 本文算法和文獻(xiàn)[12]算法比較

由于不能事先知道信道長度的先驗信息，文獻(xiàn)[12]算法把得到的時域沖擊響應(yīng)在CP之外的所有樣值估計點置零，從一定程度上提高了估計精度。圖7通過改變文獻(xiàn)[12]算法中循環(huán)前綴L的值，采用QPSK調(diào)制，分別在信道1(分別取L=60，40，20)和信道2(分別取L=120，100，68)環(huán)境下進(jìn)行仿真，比較其BER性能。

從圖7可以看出，在兩種信道環(huán)境下，隨著循環(huán)前綴L的長度逐漸接近信道的真實長度(信道1：max=20，信道2：max=68)，所得到的BER性能逐漸得到改善，這就指導(dǎo)我們可以找到一種方法，通過這種方法找到信道的實際長度max，并能有效抑制信道的噪聲值，提高信道的估計精度。

圖7(a)和圖7(b)分別比較了在信道1環(huán)境下(取L=max=20)、信道2環(huán)境下(取L=max=68)本文算法和文獻(xiàn)[12]算法兩種方法的BER性能曲線，其中本文算法的閾值系數(shù)在兩種環(huán)境下分別取=4和=3?？梢钥闯?，當(dāng)循環(huán)前綴L的長度和信道真實長度max相等時，文獻(xiàn)[12]算法達(dá)到最優(yōu)估計性能，而且這時，本文算法優(yōu)于文獻(xiàn)[12]算法。

PRRT2基因相關(guān)發(fā)作性疾病的臨床特點與預(yù)后（附1例報告） … …………………… 梁超，王丹，郭虎 469

(a) 信道1, QPSK調(diào)制，=4

(b) 信道2, QPSK 調(diào)制，=3

圖7 在不同循環(huán)前綴長度文獻(xiàn)[12]算法和本文算法的BER性能

Fig.7 BER performance comparison between the proposed algorithm and the algorithm in Literature [12]

以上兩種算法BER性能相差較大的主要原因是因為本文算法有效抑制了噪聲的干擾，比較可以看出，即使文獻(xiàn)[12]采用的循環(huán)前綴長度剛好等于信道的長度L=max，這時

而本文算法由于采用的是基于閾值的信道估計，這時

3.3 參數(shù)的選取

由于帶寬資源受限，在高速水聲通信系統(tǒng)中，大多采用信道利用率較高的非常模信號，下面就兩種調(diào)制方式下(分別為QPSK、64QAM)變化不同閾值得到的BER性能曲線進(jìn)行比較。

圖8(a)是在信道1環(huán)境下采用QPSK調(diào)制信號進(jìn)行的仿真。文獻(xiàn)[12]算法取L=max=20，這時其BER性能最優(yōu)。從圖8(a)看出，當(dāng)閾值系數(shù)在區(qū)間[4,7]時，本文算法性能優(yōu)于文獻(xiàn)[12]算法最優(yōu)估計性能，且取=7，BER=1×10-3時，兩算法相差約2.5 dB。

圖8(b)相比圖8(a)改用64QAM調(diào)制信號。從圖8(b)看出，當(dāng)閾值系數(shù)在區(qū)間[4,7]時，本文算法性能優(yōu)于文獻(xiàn)[12]算法的最優(yōu)估計性能，且取=7， BER=1×10-2時，兩者相差約1.5 dB。

圖8(c)是在信道2環(huán)境下采用QPSK調(diào)制信號進(jìn)行的仿真。文獻(xiàn)[12]算法取循環(huán)前綴長度L=max=68，這時其BER性能最優(yōu)。從圖8(c)看出，當(dāng)閾值系數(shù)在區(qū)間[3,7]時，本文算法性能優(yōu)于文獻(xiàn)[12]算法最優(yōu)估計性能。且當(dāng)本文算法取=7，BER=1×10-3時，兩者相差約2 dB。

圖8(d)相比圖8(c)改用64QAM調(diào)制信號。從圖8(d)看出，當(dāng)閾值系數(shù)在區(qū)間[3,7]時，本文算法性能優(yōu)于文獻(xiàn)[12]算法最優(yōu)估計性能。且當(dāng)本文算法取7，BER=1×10-3時，兩者相差約1.5 dB。

綜合以上可以看出，在稀疏多徑水聲信道估計應(yīng)用中使用DFT插值方法的優(yōu)勢和實用性；從不同調(diào)制方式對算法閾值的影響，得出閾值系數(shù)在區(qū)間[4,7]時，本文算法均具有比文獻(xiàn)[12]更好的性能。由于在實際信道估計應(yīng)用中，我們無法確切知道信道的長度，循環(huán)前綴的長度不會和信道長度恰好相等，所以，使用本文算法比文獻(xiàn)[12]具有約2 dB的優(yōu)勢。

(a) 信道1, QPSK調(diào)制

(b) 信道1, 64QAM調(diào)制

(c) 信道2, QPSK調(diào)制

(d) 信道2, 64QAM調(diào)制

圖8 不同閾值系數(shù)時本文算法和文獻(xiàn)[12]算法的BER性能比較

Fig.8 BER performance comparison between the proposed algorithm and the algorithm in Literature [12]

4 結(jié)論

本文在文獻(xiàn)[12]基于DFT插值的信道估計結(jié)構(gòu)上進(jìn)行了改進(jìn)，提出了一種改進(jìn)的基于閾值探測的LS信道估計方法，該方法的基本思路：首先根據(jù)LS估計準(zhǔn)則得到導(dǎo)頻位置處的信道估計；然后把LS信道估計出的結(jié)果通過IDFT變換到時域；再根據(jù)時域信道估計值，確定有用徑閾值；最后根據(jù)閾值確定有用徑。在仿真分析中，首先比較了四種插值方法在稀疏多徑水聲信道中的BER性能，證明了DFT插值的優(yōu)勢和實用性。其次分析了循環(huán)前綴長度的取值對BER的影響，當(dāng)循環(huán)前綴長度與信道長度越接近時，BER性能越好。并在信道1和信道2兩種環(huán)境下采用兩種調(diào)制方式模擬了本文算法對閾值的影響，得出閾值系數(shù)在區(qū)間[4,7]時，本文算法均具有比文獻(xiàn)[12]更好的性能，即具有SNR可以低約2 dB的優(yōu)勢。

[1] Stojanovic M. Low complexity OFDM detector for underwater acoustic Channels[C] //IEEE Oceans Conf, 2006: 1-6.

[2] Berger C R, Zhou S L, Preisig J C, et al．Sparse channel estimation for multicarrier underwater acoustic communication: From subspace methods to compressed sensing[J]. IEEE Transactions on Signal Processing, 2010, 58(3): 1708-1721.

[3] 寧小玲, 劉忠, 張林森. 多載波水聲信道估計技術(shù)及其研究進(jìn)展[J]. 聲學(xué)技術(shù), 2012, 31(1): 77-81． NíNG Xiaoling, LIU Zhong, ZHANG Linsen. Advances in channel estimation for multicarrier underwater acoustic communication[J]. Technical Acoustics, 2012, 31(1): 77-81．

[4] Kyungchul Kwak, Sungeun Lee, Jihyung Kim, et al. A new dft-based channel estimation approach for OFDM with virtual subcarriers by leakage estimation[J]. IEEE Transactions on wireless communications, 2008, 7(6): 2004-2008．

[5] 徐小卡, 桑恩方, 李娜, 等．基于DFT的OFDM水聲通信信道估計算法[J]. 數(shù)據(jù)采集與處理, 2009, 24(1): 44-48． XU Xiaoka, SANG Enfang, LI Na, et al. Channel estimation algorithm in underwater acoustic communication based on DFT and OFDM[J]. Journal of Data Acquisition & Processing, 2009, 24(1): 44-48.

[6] Li B S, Zhou S L, Stojanovic M, et al. Pilot-tone based ZP-OFDM demodulation for an underwater acoustic channel[C] //Proceedings of the MTS/IEEE Oceans Conference, Boston, 2006: 1-5.

[7] Li B S，Zhou S L，Stojanovic M，et al．Multicarrier communication over underwater acoustic channels with nonuniform Doppler shifts [J]. IEEE Journal of Oceanic Engineering, 2008, 33(2): 198-209．

[8] 張永懷, 孫海信, 許芳, 等. 應(yīng)用于水聲OFDM通信系統(tǒng)的基于維納濾波的信道估計[J]. 廈門大學(xué)學(xué)報: 自然科學(xué)版, 2009, 48(1): 51-54. ZHANG Yonghuai, SUN Haixin, XU Fang, et al. Channel estimation based on wiener filter for underwater acoustic OFDM communication systems [J]. Journal of Xiamen University: Natural Science, 2009, 48(1): 51-54.

[9] Fraser D. Interpolation by the FFT revised-an experimental investigation[J]. IEEE Transactions on Acoustics, Speech, and Signal Processing, 1989, 37(5): 1007-1015.

[10] ZHAO Yuping, HUANG Aiping. A novel channel estimation method for OFDM mobile communication systems based on pilot signals and transform-domain processing[C]//IEEE Vehicular Technology Conf.,1998, 46: 931-939.

[11] Dowler A, Doufexi A, Nix A. Performance evaluation of channel estimation techniques for a mobile fourth generation wide area OFDM system[C]//IEEE Vehicular Technology Conf., 2002, 4: 2036-2040.

[12] Qu F Z, Yang L Q. On the estimation of doubly-selective fading channels[J]. IEEE Transactions on Wireless Communications，2008, 9(4): 1261-1265．

[13] 童峰, 許肖梅, 方世良. 一種單頻水聲信號多徑時延估計算法[J]. 聲學(xué)學(xué)報, 2008, 33(1): 62-69. TONG Feng, XU Xiaomei, FANG Shiling. Multipath time-delay estimation of underwater acoustic sinusoidal signals[J]. Acta Acustica, 2008, 9(4): 1261-1265．

[14] 何劍, 王永生, 姚如貴, 等. 基于廣義A IC 和RPC 的O FDM 航空信道估計算法[J]. 西北工業(yè)大學(xué)學(xué)報, 27(2): 250-254. HE Jian, WANG Yongsheng, YAO Rugui, et al. A better algorithm for estimating aeronautical channel OFDM systembased on GAIC and RPC [J]. Journal of Northwestern Polytechnical University, 27(2): 250-254.

[15] Stojanovic M. Retrofocusing techniques for high rate acoustic communications[J]. The Journal of the Acoustical Society of America, 2005, 117(3): 1173-1185．

[16] 宋伯煒. OFDM無線寬帶移動通信系統(tǒng)中信道估計與均衡技術(shù)研究[D]. 上海: 上海交通大學(xué), 2005． SONG Bowei. Study on Channel estimation and equalization technique for OFDM based wireless broadband mobile commnication systems[D]. Shanghai: Shanghai Jiaotong University, 2005.

[17] 王峰. 基于高階統(tǒng)計量的水聲信道盲均衡理論與算法[D]. 西安: 西北工業(yè)大學(xué), 2003. WANG Feng. Blind equalization algorithms using higher-order statistics for underwater acoustic channel[D]. Xi’an: Northwestern Polytechnical University, 2003.

[18] 郭業(yè)才. 自適應(yīng)盲均衡技術(shù)[M]. 合肥: 合肥工業(yè)出版社, 2007. GUO Yecai. Adaptive blind equalization technique[M]. Hefei: Hefei Industry Press, 2007.

Application of an improved LS channel estimation algorithm to sparse multipath underwater acoustic channel

NING Xiao-ling1, ZHANG Lin-sen2, LIANG YUE3

(1. Electronics Engineering College, Naval University of Engineering, Wuhan 430033,Hubei,China;2.Department of Weaponary Engineering, Naval University of Engineering, Wuhan 430033, Hubei, China;3. Dean’s Office of Naval University of Engineering, Wuhan 430033,Hubei, China)

In order to eliminate the effect brought from the noise in underwater orthogonal frequency division multiplexing (OFDM) system, an improved least-square (LS) channel estimation method is proposed for sparse multi-path channel estimation. The conventional discrete Fourier transform (DFT)-based interpolation channel estimation method is improved by using threshold detecting channel estimation, which can decide significant channel taps and erase the effect of noise by comparing the time domain channel estimation value with the predetermined threshold. The proposed method increases estimation precision effectively. Simulation results demonstrate that the DFT-based interpolation is practical in sparse multipath underwater acoustic channel, and the bit error rate (BER) becomes better when the circular prefix approximate actual channel length gradually. Moreover, the threshold coefficient regions of this method for two modulation modes are ascertained. The results show that BER of the proposed LS channel estimation is better than that of the conventional DFT-based channel estimations by about 2 dB, therefore the problem of using circular prefix to substitute actual channel length, which exists in the conventional DFT-based channel estimation method, is solved.

least square channel estimation; orthogonal frequency division multiplexing; sparse multipath underwateracoustic channel; discrete Fourier transform based interpolation

TN911.5

A

1000-3630(2016)-04-0378-07

10.3969/j.issn1000-3630.2016.04.018

2015-10-24;

2016-01-25

國家自然科學(xué)基金(F010202)、海軍工程大學(xué)自然科學(xué)基金(HGDQNJJ3024) 資助項目。

寧小玲(1982－), 女, 湖南邵陽人, 博士, 研究方向為水下高速率數(shù)據(jù)傳輸、信道估計與均衡等。

寧小玲, E-mail: zhang_ning1982b@sina.com.