孫向前，李晴，范展

?

全相位頻譜校正技術(shù)在水聲通信中的應(yīng)用研究

孫向前1，李晴2,3，范展2,3

(1. 廣東湛江91388部隊(duì)，廣東湛江524022；2. 哈爾濱工程大學(xué)水聲技術(shù)重點(diǎn)實(shí)驗(yàn)室，黑龍江哈爾濱 150001； 3. 哈爾濱工程大學(xué)水聲工程學(xué)院，黑龍江哈爾濱 150001)

多普勒效應(yīng)是影響水下移動(dòng)通信性能的主要因素，準(zhǔn)確估計(jì)多普勒頻移對(duì)提高通信系統(tǒng)可靠性具有重要意義。在進(jìn)行離散頻譜分析時(shí)，時(shí)域非整周期截?cái)鄷?huì)造成頻域能量的泄漏，導(dǎo)致頻譜估計(jì)精度降低。全相位頻譜分析具有良好的抑制頻譜泄漏特性及相位不變性。仿真驗(yàn)證了全相位頻譜校正技術(shù)相對(duì)于傳統(tǒng)頻譜校正技術(shù)在估計(jì)性能上的優(yōu)勢(shì)，并在此基礎(chǔ)上探討了全相位頻譜校正技術(shù)在水聲通信中的應(yīng)用。采用全相位頻譜校正技術(shù)進(jìn)行多普勒頻移估計(jì)，進(jìn)而進(jìn)行多普勒補(bǔ)償以降低通信系統(tǒng)誤碼率。仿真結(jié)果表明，全相位頻譜校正技術(shù)能夠?qū)崿F(xiàn)高精度多普勒頻移估計(jì)，從而提高水下移動(dòng)通信系統(tǒng)的可靠性。

全相位頻譜分析；全相位頻譜校正；水聲通信；多普勒

0 引言

快速傅里葉變換(Fast Fourier Transform, FFT)顯著減小了離散傅里葉變換的運(yùn)算量，被廣泛應(yīng)用于雷達(dá)、聲吶等信號(hào)處理領(lǐng)域。在工程實(shí)際中進(jìn)行FFT頻譜分析時(shí)，受到器件處理樣本長(zhǎng)度的限制，需要先將連續(xù)信號(hào)離散化、截?cái)嗖⒅芷谘油亍r(shí)域非整周期截?cái)嘣陬l域表現(xiàn)為能量的泄露，導(dǎo)致離散頻譜頻率、幅值和相位估計(jì)精度降低。

離散頻譜分析的誤差校正及高精度頻譜估計(jì)成為近幾十年內(nèi)現(xiàn)代信號(hào)處理的研究方向之一。Offelli C.等人研究了加窗函數(shù)對(duì)多頻率信號(hào)參數(shù)估計(jì)的影響[1]。Santamria I.等人對(duì)包括內(nèi)插FFT在內(nèi)的幾種頻率估計(jì)方法進(jìn)行了對(duì)比[2]。張強(qiáng)等人研究了雙頻率模型的頻譜校正形式及性能[3]。目前國(guó)內(nèi)外有關(guān)頻譜校正技術(shù)的研究重點(diǎn)落在單頻率諧波信號(hào)或者頻率間隔較大的多頻率信號(hào)的離散頻譜校正上。典型的方法有比值法[4]、能量重心法[5]和相位差法[6]。這幾種方法首先對(duì)截?cái)嗪笾芷谘油氐男盘?hào)進(jìn)行FFT分析，再對(duì)得到的離散頻譜進(jìn)行校正，其精度受頻譜泄露程度的影響較大，本文稱之為傳統(tǒng)頻譜校正方法。在圖像處理領(lǐng)域，常用重疊方式來改善由圖像分塊導(dǎo)致的邊界“方塊效應(yīng)”。Malvar 提出的雙正交重疊變換在圖像壓縮中得到廣泛應(yīng)用[7]。全相位數(shù)據(jù)預(yù)處理起源于王兆華對(duì)圖像處理中的重疊數(shù)字濾波的研究，主要用來解決信號(hào)非整周期截?cái)嗟膯栴}。預(yù)處理考慮包含某樣點(diǎn)的所有可能的數(shù)據(jù)截?cái)嗲闆r，即進(jìn)行最大程度的重疊處理，最終輸出這些處理結(jié)果的有機(jī)綜合。至今這種方法的應(yīng)用已經(jīng)由二維的圖像處理領(lǐng)域擴(kuò)展到了更普遍的一維信號(hào)處理領(lǐng)域，并形成了全相位快速傅里葉變換(all-phase FFT, apFFT)頻譜分析、全相位頻譜校正[8]和全相位濾波等理論及方法。apFFT具有非常優(yōu)良的抑制譜泄露性能以及相位不變性[9]，因而全相位頻譜校正受泄露的影響比傳統(tǒng)方法小，相位也不需要校正；在頻譜泄漏嚴(yán)重、譜線間干擾較大等情況下比傳統(tǒng)方法更具優(yōu)勢(shì)，具有廣泛應(yīng)用前景。

水聲通信中的多普勒頻移是由收發(fā)平臺(tái)之間的相對(duì)運(yùn)動(dòng)或信道中水的流動(dòng)引起的。相對(duì)速度與聲速之比可稱為多普勒頻移因子。由于水中聲速比空中電磁波波速約小了105倍，水聲移動(dòng)通信中多普勒頻移因子遠(yuǎn)大于無線電通信。而水聲通信中常用的信號(hào)處理方法對(duì)于多普勒導(dǎo)致的信號(hào)頻譜擴(kuò)展、頻率偏移等畸變適應(yīng)性較差，從而導(dǎo)致系統(tǒng)可靠性急劇下降[10]。工程實(shí)際中通常先估計(jì)多普勒頻移因子，再根據(jù)頻移因子進(jìn)行多普勒補(bǔ)償[11]。其中高精度估計(jì)頻移因子成為降低通信誤碼率的關(guān)鍵。而頻譜校正技術(shù)能夠高精度測(cè)頻，所以本文采用全相位頻譜校正技術(shù)來提高多普勒頻移因子估計(jì)的精度，進(jìn)而提高水聲通信系統(tǒng)的性能。

1 全相位頻譜校正原理及性能分析

1.1 全相位頻譜分析

全相位FFT(apFFT)頻譜分析是指對(duì)經(jīng)過全相位預(yù)處理后的數(shù)據(jù)再進(jìn)行FFT頻譜分析[12]。全相位數(shù)據(jù)預(yù)處理將長(zhǎng)度為的數(shù)據(jù)向量映射為長(zhǎng)度為的數(shù)據(jù)向量，且是中包含其中心點(diǎn)的長(zhǎng)度為的所有分段序列的疊加。設(shè)為單頻復(fù)指數(shù)信號(hào)，

(2)

apFFT頻譜為

1.2 全相位頻譜校正

全相位頻譜校正是指對(duì)apFFT譜分析結(jié)果進(jìn)行譜線頻率、幅度和相位補(bǔ)償修正，通常先估計(jì)出主譜線和真實(shí)譜線間的歸一化頻率誤差：

(5)

1.2.1 全相位能量重心法

全相位能量重心法利用對(duì)稱窗函數(shù)的離散譜能量重心在原點(diǎn)或原點(diǎn)附近這一特性對(duì)頻率進(jìn)行校正。設(shè)單頻率諧波信號(hào)FFT譜分析得到的離散頻譜的主譜線位置為，幅度譜值、功率譜值分別為、，apFFT譜對(duì)應(yīng)幅值為，且有

則依據(jù)離散功率譜線的能量重心位置是信號(hào)真實(shí)頻率的理論位置，可得該信號(hào)頻率校正公式：

(7)

由帕賽瓦爾定理可得幅度校正公式為

1.2.2 全相位比值法

全相位比值法先將apFFT振幅譜線歸一化，再求主譜線和最大旁譜線幅值的比值，記為：

(10)

將式(10)代入式(5)可得校正后頻率。幅度校正公式為

1.2.3 全相位時(shí)移相位差法

全相位時(shí)移相位差法基于apFFT譜分析的線性性質(zhì)和相位不變性，利用輸入離散信號(hào)的兩段序列apFFT譜分析得到的對(duì)應(yīng)主譜線的相位差進(jìn)行校正。取采樣數(shù)據(jù)中兩段存在延時(shí)的長(zhǎng)為的樣本分別進(jìn)行apFFT，取對(duì)應(yīng)主譜線的相位值分別為和，則由相位差可得頻率校正值為

幅度校正值可由式(11)求得。具體流程如圖1所示。圖1中，為主譜線上的頻偏估值。

圖1 全相位時(shí)移相位差法頻譜校正流程

Fig.1 Flowchart for spectrum correction of all-phase time-shift difference

相比傳統(tǒng)相位差法，全相位時(shí)移相位差法的相位估計(jì)與頻偏無關(guān)，直接取對(duì)應(yīng)相位譜值即可，計(jì)算量減??；另外，由于頻譜泄漏得到良好抑制，當(dāng)信號(hào)包含多頻成分時(shí)，估計(jì)精度受譜線間干擾影響較小。

1.2.4 FFT/apFFT綜合相位差法

FFT/apFFT綜合相位差法利用了兩種譜分析的主譜線信息[13]。首先對(duì)一段信號(hào)分別做FFT和apFFT分析，得到對(duì)應(yīng)主譜線的相位值，分別為和，其差值與頻率誤差成正比，比例系數(shù)為群延時(shí)，則頻率誤差為

將式(13)代入式(5)可得校正后頻率；取apFFT主譜線的相位值即為信號(hào)初相位；由式(6)可得幅度校正公式為

(14)

1.3 性能仿真分析

1.2節(jié)推導(dǎo)了幾種典型的全相位頻譜校正技術(shù)的計(jì)算公式，并與傳統(tǒng)頻譜校正技術(shù)進(jìn)行理論上的對(duì)比。下面將通過仿真實(shí)驗(yàn)對(duì)這兩類校正方法的性能進(jìn)行綜合比較分析。

表1 不同信噪比下幾種頻譜校正方法校正均方根誤差

2 基于全相位頻譜校正技術(shù)的水聲通信系統(tǒng)仿真研究

2.1 水聲通信系統(tǒng)設(shè)計(jì)

根據(jù)上述的理論推導(dǎo)與仿真驗(yàn)證，相比于傳統(tǒng)頻譜校正技術(shù)，基于全相位頻譜分析的全相位頻譜校正技術(shù)受到頻譜泄漏的影響更小，在提高頻率估計(jì)精度上更有優(yōu)勢(shì)。因此，考慮將全相位頻譜校正技術(shù)應(yīng)用于水下移動(dòng)通信系統(tǒng)中，利用其頻率估計(jì)精度高的性能對(duì)接收信號(hào)的多普勒進(jìn)行估計(jì)，從而采取相應(yīng)的補(bǔ)償措施，以提高通信系統(tǒng)可靠性。

與無線電通信相比，水聲通信系統(tǒng)的工作環(huán)境更為惡劣。發(fā)射機(jī)發(fā)射的聲波在復(fù)雜結(jié)構(gòu)的水聲通信信道中傳輸，主要產(chǎn)生能量的衰減及波形的畸變。聲信號(hào)傳播損失主要由幾何擴(kuò)展損失和介質(zhì)吸收損失構(gòu)成，即

(17)

式(17)中，為信號(hào)頻率，單位：kHz。可見，海水吸收損失近似與聲波頻率成正比。

聲信號(hào)能量的衰減主要影響水聲通信系統(tǒng)的通信距離。由于水聲信道具有多徑特性，信號(hào)經(jīng)信道傳輸在接收端會(huì)產(chǎn)生時(shí)間擴(kuò)展及頻率選擇性衰落。另一方面，當(dāng)收發(fā)平臺(tái)存在相對(duì)運(yùn)動(dòng)時(shí)，接收信號(hào)受多普勒效應(yīng)的影響會(huì)產(chǎn)生頻率的偏移及時(shí)域的伸縮。這些因素都會(huì)導(dǎo)致信號(hào)波形的畸變，從而影響通信系統(tǒng)的可靠性。

由以上分析，仿真條件及系統(tǒng)參數(shù)設(shè)置如下：系統(tǒng)通信距離為5 km，信道沖擊響應(yīng)如圖2所示，采用帶限高斯白噪聲。接收機(jī)相對(duì)發(fā)射機(jī)的運(yùn)動(dòng)速度為，且不超過40 kn，即在±20 m/s(相向運(yùn)動(dòng)時(shí)取正，相背運(yùn)動(dòng)時(shí)取負(fù))范圍內(nèi)。設(shè)計(jì)二進(jìn)制頻移鍵控水聲數(shù)字通信系統(tǒng)，數(shù)據(jù)碼元信息采用頻率為10 kHz和10.15 kHz的兩個(gè)單頻脈沖作為載波；用于多普勒頻移修正的導(dǎo)頻同步信號(hào)，采用頻率為6 kHz的CW脈沖；采樣頻率為51.2 kHz；同步脈沖及碼元信號(hào)脈寬分別為15 ms和20 ms，對(duì)應(yīng)的帶通濾波器的階數(shù)均為128，3 dB帶寬分別為600 Hz和100 Hz。

接收端處理流程如圖3所示。

主要步驟如下：

(1) 采用全相位頻譜校正技術(shù)對(duì)接收數(shù)據(jù)中的導(dǎo)頻信號(hào)進(jìn)行頻率估計(jì)；

(2) 由導(dǎo)頻信號(hào)的頻率估值及其發(fā)射頻率估算多普勒頻移因子；

(19)

(3) 根據(jù)多普勒頻移因子重構(gòu)本地載波，進(jìn)而對(duì)接收的數(shù)據(jù)信號(hào)進(jìn)行解調(diào)及判決。

2.2 水聲通信系統(tǒng)性能分析

由1.3節(jié)的仿真實(shí)驗(yàn)可知，全相位時(shí)移相位差法有較高的頻率估計(jì)精度，因此下面主要以該方法為例進(jìn)行系統(tǒng)性能分析，并以頻譜校正前的傳統(tǒng)FFT分析作為參考。設(shè)譜分析點(diǎn)數(shù)為=256，則傳統(tǒng)方法采用數(shù)據(jù)點(diǎn)數(shù)為，校正方法采用數(shù)據(jù)點(diǎn)數(shù)為(3-1)。

由圖4可見，當(dāng)SNR在0~10 dB范圍內(nèi)，校正前后的估計(jì)均方根誤差均隨SNR的增大而降低，且校正后降低速率更大；SNR大于10 dB時(shí)，校正前后誤差均趨于穩(wěn)定，校正前誤差趨于0.01，而校正后誤差趨于零。

表2 不同信噪比下采用各頻譜校正技術(shù)的系統(tǒng)誤碼率

通過上述分析可見，全相位頻譜校正技術(shù)對(duì)于降低水聲移動(dòng)通信系統(tǒng)誤碼率具有重要意義。實(shí)際應(yīng)用中，如果信噪比較低，且對(duì)誤碼率要求很高，可以采用全相位時(shí)移相位差法；FFT/apFFT綜合相位差法的運(yùn)算復(fù)雜度高于傳統(tǒng)相位差法，低于全相位時(shí)移相位差法，且估計(jì)性能接近全相位時(shí)移相位差法，適用于需要綜合考慮通信速率和誤碼率、且信噪比較高的條件。

3 總結(jié)

本文從理論和仿真的角度，對(duì)全相位頻譜校正技術(shù)進(jìn)行了綜合性能分析，并在此基礎(chǔ)上研究了該技術(shù)在水聲通信中的應(yīng)用。

仿真結(jié)果表明，與傳統(tǒng)頻譜校正技術(shù)相比，全相位頻譜校正技術(shù)受頻譜泄漏影響小，頻率估計(jì)精度高。

在水聲移動(dòng)通信系統(tǒng)中，采用全相位頻譜校正技術(shù)進(jìn)行多普勒頻移估計(jì)可以提高估計(jì)精度，從而降低系統(tǒng)誤碼率。全相位頻譜校正技術(shù)在改善水聲通信系統(tǒng)性能上具有重要的理論意義和應(yīng)用價(jià)值。

[1] Offelli C, Petri D. The Influence of Windowing on the Accuracy of Multifrequency Signal Parameter Estimation[J]. IEEE Transactions on Instrumentation and Measurement, 1992, 41(2): 256-261.

[2] Santamria I, Pantaleon C, Ibanez J. A Comparative Study of High-accuracy Frequency Estimation Methods[J]. Mechanical Systems and Signal Processing, 2000, 14(5): 819-834.

[3] 張強(qiáng), 張頻, 陳奎孚. 對(duì)兩個(gè)頻率相近成分作頻譜校正的非迭代形式研究[J]. 振動(dòng)與沖擊, 2012, 31(10): 24-28.

ZHANG Qiang, ZHANG Pin, CHEN Kuifu. Non- iterative spectrum correction for signals with closely spaced frequency components[J]. Journal of Vibration and Shock, 2012, 31(10): 24-28.

[4] 朱磊, 董亮, 孫振龍, 等. 離散頻譜梯形窗幅度比值校正法的頻率估計(jì)[J]. 電視技術(shù), 2011, 35(11): 102-104.

ZHU Lei, DONG Liang, SUN Zhenlong, et al. Discrete Spectral Amplitude Ratio Correction Method with Trapezoidal Window of Frequency Estimation[J]. Video Engineering, 2011, 35(11): 102-104.

[5] 丁康, 鄭春松, 楊志堅(jiān). 離散頻譜能量重心法頻率校正精度分析及改進(jìn)[J]. 機(jī)械工程學(xué)報(bào), 2010, 46(5): 43-48.

DING Kang, ZHENG Chunsong, YANG Zhijian. Frequency Estimation Accuracy Analysis and Improvement of Energy Barycenter Correction Method for Discrete Spectrum[J]. Journal of Mechanical Engineering, 2010, 46(5): 43-48.

[6] 胡文彪, 夏立, 向東陽, 等. 一種改進(jìn)的基于相位差法的頻譜校正方法[J]. 振動(dòng)與沖擊, 2012, 31(1): 162-166.

HU Wenbiao, XIA Li, XIANG Dongyang, et al. An improved frequency spectrum correction method based on phase difference correction method[J]. Journal of Vibration and Shock, 2012, 31(1): 162-166.

[7] Malvar H S. Biorthogonal and Nonuniform Lapped Transforms for Transform Coding with Reduced Blocking and Ringing Artifacts[J]. IEEE Transactions on Signal Processing, 1998, 46(4): 1043-1053.

[8] 王兆華, 黃翔東. 數(shù)字信號(hào)全相位譜分析與濾波技術(shù)[M]. 北京: 電子工業(yè)出版社, 2009.

WANG Zhaohua, HUANG Xiangdong. Digital Signal All-phase Spectrum Analysis and Filtering[M]. Beijing: Publishing House of Electronics Industry, 2009.

[9] 黃翔東, 王兆華, 羅蓬, 等. 全相位FFT密集譜識(shí)別與校正[J]. 電子學(xué)報(bào), 2011, 39(1): 172-177.

HUANG Xiangdong, WANG Zhaohua, LUO Peng, et al. Discrimination and Correction for Dense All-phase FFT Spectrums[J]. Acta Electronica Sinica, 2011, 39(1): 172-177.

[10] 許天增, 許鷺芬. 水聲數(shù)字通信[M]. 北京: 海洋出版社, 2010.

XU Tianzeng, XU Lufen. Underwater Acoustic Digital Communication[M]. Beijing: China Ocean Press, 2010.

[11] Kenneth A. Perrine, Karl F. Nieman, Terry L. Henderson, et al. Doppler Estimation and Correction for Shallow Underwater Acoustic Communications[C]// Proceedings of ASILOMAR 2010. Pacific Grove, CA, USA: 2010. 746 - 750.

[12] 賈方秀, 丁振良, 袁峰, 等. 基于全相位快速傅里葉變換譜分析的激光動(dòng)態(tài)目標(biāo)實(shí)時(shí)測(cè)距系統(tǒng)[J]. 光學(xué)學(xué)報(bào), 2010, 30(10): 2928- 2934.

JIA Fangxiu, DING Zhenliang, YUAN Feng, et al. Real-Time Laser Range Finding System for Moving Target Based on All -Phase Fourier Transform Spectrum Analysis[J]. Acta Optica Sinica, 2010, 30(10): 2928-2934.

[13] 黃翔東, 王兆華. 基于全相位頻譜分析的相位差頻譜校正法[J]. 電子與信息學(xué)報(bào), 2008, 30(2): 293-297.

HUANG Xiangdong, WANG Zhaohua. Phase Difference Correcting Spectrum Method Based on All-phase Spectrum Analysis[J]. Journal of Electronics & Information Technology, 2008, 30(2): 293-297.

Research on the application of all-phase spectrum correction technique in underwater acoustic communication

SUN Xiang-qian1, LI Qing2,3, FAN Zhan2,3

(1.Unit 91388, The People’s Liberation Army of China, Zhanjiang 524022, Guangdong, China;2. Acoustic Science and Technology Laboratory, Harbin Engineering University, Harbin 150001, Heilongjiang, China;3. College of Underwater Acoustic Engineering, Harbin Engineering University, Harbin 150001, Heilongjiang, China)

The Doppler effect is a main factor influencing the performance of underwater mobile communication, so it’s important to effectively estimate the frequency shift for the improvement of the system reliability. Resulting from non-inter-period sampling in time domain, the leakage error of discrete spectrum analysis leads to reducing accuracy of spectrum estimation. The satisfactory characteristics of all-phase spectrum analysis are lower spectrum leakage level and non-phase shift. Simulation experiments prove that all-phase spectrum correction technique has the advantage in precision over traditional methods. Furthermore, the application of all-phase spectrum correction in underwater acoustic communication is studied. The technique is adopted for Doppler estimation and compensation to reduce the error rate of system. Simulation results show that adopting all-phase spectrum correction for high frequency estimation accuracy improves the reliability of underwater mobile communication system.

all-phase spectrum analysis; all-phase spectrum correction; underwater acoustic communication; Doppler

TN929.3

A

1000-3630(2015)-02-0127-07

10.16300/j.cnki.1000-3630.2015.02.005

2014-03-06;

2014-06-07

國(guó)家自然科學(xué)基金資助項(xiàng)目(51279043, 61201411); 國(guó)家“863計(jì)劃”資助項(xiàng)目(2013AA09A503)

孫向前(1968－), 男, 河南安陽人, 高級(jí)工程師, 研究方向?yàn)樗暅y(cè)控總體技術(shù)。

李晴, E-mail: liqing_heu@163.com