胡速謀,屠健康
(1.中國(guó)人民解放軍91033 部隊(duì),山東 青島 266000;2.中國(guó)人民解放軍92325 部隊(duì),山西 大同 037008)
甚低頻是工作頻率為3~30 kHz,波長(zhǎng)范圍在10~100 km 的無(wú)線電磁波。甚低頻信道有空氣信道和海洋信道兩種信道傳播環(huán)境。甚低頻通信具有信號(hào)傳播穩(wěn)定、衰減小、穿透海水能力較強(qiáng)的特點(diǎn),因此在對(duì)潛通信和遠(yuǎn)洋通信領(lǐng)域有著廣泛的應(yīng)用。
最小移頻鍵控(Minimum Shift Keying,MSK)信號(hào)是當(dāng)前世界發(fā)達(dá)國(guó)家海軍甚低頻通信的主要調(diào)制方式,該調(diào)制方式具有功率譜集中、頻帶利用率高、包絡(luò)恒定、帶外輻射干擾小等優(yōu)點(diǎn)。但信號(hào)經(jīng)過(guò)地面及淺層的土壤、大氣層、電離層以及海水等多種介質(zhì)的遠(yuǎn)距離傳輸后,再到達(dá)水下甚低頻接收端的過(guò)程中,受到多徑效應(yīng)以及多普勒頻移等因素影響,不可避免地產(chǎn)生失真和畸變,從而引起碼間串?dāng)_。碼間串?dāng)_會(huì)導(dǎo)致信道中攜帶的信息在解調(diào)時(shí)誤碼較高,嚴(yán)重影響甚低頻通信的質(zhì)量;而采用合適的均衡技術(shù)就能有效減少或者清除信號(hào)的碼間串?dāng)_影響,提高甚低頻系統(tǒng)的性能。
單載波頻域均衡(Single Carrier-Frequency Domain Equalization,SC-FDE)系統(tǒng)是近幾年發(fā)展較快的,能夠有效對(duì)抗多徑效應(yīng)的方法。單載波頻域均衡技術(shù)借鑒多載波傳輸系統(tǒng)正交頻分復(fù)用(Orthogonal Frequency Division Multiplexing,OFDM)的基本模式,在接收端通過(guò)離散傅里葉變換快速算法(Fast Fourier Transform,F(xiàn)FT)與離散傅里葉反變換快速算 法(Inverse Fast Fourier Transform,IFFT)進(jìn) 行數(shù)字信號(hào)處理,實(shí)現(xiàn)頻域均衡,可有效恢復(fù)甚低頻數(shù)字信號(hào)。具體地,對(duì)于MSK 調(diào)制信號(hào),由于其屬于二進(jìn)制的全響應(yīng)連續(xù)相位調(diào)制信號(hào),若采用Viterbi 算法進(jìn)行時(shí)域的最大似然檢測(cè),算法復(fù)雜程度會(huì)很高;而SC-FDE 系統(tǒng)提供了一種低復(fù)雜度的頻域?qū)勾a間串?dāng)_的有效方案,系統(tǒng)實(shí)現(xiàn)簡(jiǎn)單、成本大大降低,因此研究利用SC-FDE 線性均衡技術(shù)來(lái)減少或者消除信道的碼間串?dāng)_具有重大的應(yīng)用價(jià)值[1]。
SC-FDE 系統(tǒng)是近幾年興起的有效對(duì)抗多徑干擾的方法。該系統(tǒng)有效借鑒了多載波傳輸系統(tǒng)正交頻分復(fù)用(Orthogonal Frequency Division Multiplexing,OFDM)技術(shù),并在系統(tǒng)接收端運(yùn)用了FFT 和IFFT 算法,提高了系統(tǒng)的運(yùn)算效率;同時(shí),結(jié)合了傳統(tǒng)的單載波傳輸技術(shù)的優(yōu)點(diǎn),將單載波系統(tǒng)接收端的時(shí)域均衡濾波器變換成頻域均衡濾波器進(jìn)行數(shù)字信號(hào)處理,大大降低了均衡的復(fù)雜程度,提高了系統(tǒng)的抗多徑干擾性能[2]。SC-FDE 技術(shù)主要有以下兩個(gè)方面特點(diǎn):
(1)與OFDM 相比,SC-FDE 技術(shù)克服了峰均功率比(Peak to Average Power Ratio,PAPR)和對(duì)相位噪聲的敏感性特點(diǎn),大大降低了功率放大器等模擬器件的成本[3];
(2)與單載波系統(tǒng)相比,系統(tǒng)實(shí)現(xiàn)簡(jiǎn)單,均衡器的復(fù)雜度大大降低,且抗多徑干擾效果與多載波傳輸系統(tǒng)相當(dāng)[4]。
SC-FDE 系統(tǒng)完整的系統(tǒng)結(jié)構(gòu)原理[5]如圖1 所示。SC-FDE 系統(tǒng)與OFDM 系統(tǒng)相似,主要區(qū)別在于SC-FDE 系統(tǒng)中FFT 與IFFT 的位置全部在接收端。
圖1 SC-FDE 系統(tǒng)原理
在發(fā)射端,信源產(chǎn)生的二進(jìn)制比特流經(jīng)過(guò)映射分塊操作得到長(zhǎng)度為N的符號(hào)序列塊:x(n)={x0(n),x1(n),…,xN-1(n)}。插入導(dǎo)頻后,將每個(gè)塊的最后長(zhǎng)度為G的數(shù)據(jù)符號(hào)循環(huán)前綴UW 復(fù)制到塊首,形成長(zhǎng)度為X=N+G的發(fā)射符號(hào)序列數(shù)據(jù)塊,然后將無(wú)線信道、多徑衰落信道和加性高斯白噪聲(Additive White Gaussian Noise,AWGN)信道的高斯白噪聲干擾一起送入系統(tǒng)接收端[6]。
在接收端,先將接收到的信號(hào)分成長(zhǎng)度為X數(shù)據(jù)塊;然后刪除每個(gè)數(shù)據(jù)塊中長(zhǎng)度為G的數(shù)據(jù)符號(hào)循環(huán)前綴UW;再進(jìn)行N點(diǎn)的FFT 計(jì)算,并將信道變換到頻域,在頻域內(nèi)進(jìn)行均衡處理;最后,通過(guò)N點(diǎn)的IFFT 變換到時(shí)域進(jìn)行判決,從而恢復(fù)信源信號(hào)[7]。
SC-FDE 均衡技術(shù)可用于不同的調(diào)制技術(shù),如正交相移鍵控(Quadrature Phase Shift Keying,QPSK),π/4-DQPSK,高階的正交振幅調(diào)制(Quadrature Amplitude Modulation,QAM)或相移鍵控(Phase-Shift Keying,PSK)、MSK 等。用于分塊傳輸?shù)腗SK 調(diào)制的SC-FDE 系統(tǒng)模型有其特有的傳輸方式,需要注意的是映射和逆映射部分。為了方便研究,通常設(shè)想甚低頻系統(tǒng)發(fā)射端和接收端精確同步狀態(tài),并在甚低頻MSK 信號(hào)準(zhǔn)靜態(tài)信道中,隨機(jī)抽取不同間隔的抽頭作為信道使用模型,使得每個(gè)抽頭的增益獨(dú)立同分布,且信道沖激響應(yīng)持續(xù)時(shí)間為碼元間隔的整數(shù)倍,信道的沖激響應(yīng)被理想估計(jì)[8]。
發(fā)射端信源數(shù)據(jù)的第i塊經(jīng)過(guò)映射得到長(zhǎng)度為N的調(diào)制符號(hào)塊xn,可表示為xn=[xcp,x1,x2,…,xN],其中循環(huán)前綴CP 為xcp=[xN-G+1,xN-G+2,…,xN],G為CP的長(zhǎng)度。構(gòu)成發(fā)送符號(hào)序列Xn=xn+xcp。發(fā)送端中總的信道沖激響應(yīng)可表示為hn=[h1,h2,…,hLch],Lch 為最大多徑數(shù),當(dāng)其小于等于循環(huán)前綴CP 的長(zhǎng)度時(shí),循環(huán)前綴CP 可以消除多徑信道引起的塊間干擾。在接收端,經(jīng)過(guò)去除循環(huán)前綴CP 并進(jìn)行匹配濾波之后,接收信號(hào)rn可以表示為:
式中:?表示循環(huán)卷積運(yùn)算;vn(n=0,1,…,N-1)為信道加性高斯白噪聲,服從均值是0,方差是N0的復(fù)高斯分布。對(duì)式(1)中的rn進(jìn)行N點(diǎn)FFT 變換后,可以得到:
式中:·運(yùn)算表示兩個(gè)向量中對(duì)應(yīng)元素相乘;Hn是信道頻域響應(yīng);Vn是加性高斯白噪聲vn的傅里葉變換。
然后,對(duì)頻域信號(hào)Rn均衡處理,設(shè)頻域均衡濾波器的系數(shù)是Wn,經(jīng)過(guò)頻域均衡之后,得到:
再然后對(duì)Zn進(jìn)行N點(diǎn)的IFFT 后得到:
最后經(jīng)過(guò)判決和逆映射恢復(fù)為信源信息。實(shí)際中,不同的均衡算法,對(duì)應(yīng)的均衡系數(shù)是不一樣的。
MSK 是移頻鍵控(Frequency-Shift Keying,F(xiàn)SK)的一種改進(jìn)型,它是二進(jìn)制全響應(yīng)的連續(xù)相位調(diào)制(Continue Phase Modulation,CPM)信號(hào),最小移頻鍵控(MSK)中的“最小”指的是獲得正交的調(diào)制信號(hào)的最小的調(diào)制指數(shù)h=1/2。MSK 信號(hào)在一個(gè)碼元期間恰好相差1/2 周,即相差π。對(duì)MSK 連續(xù)相位調(diào)制信號(hào)來(lái)說(shuō),加循環(huán)前綴時(shí),為了保持其相位的連續(xù)性,本文采用最簡(jiǎn)單且常用的是相位歸零法[9]。
MSK 調(diào)制指數(shù)h=1/2,有4 種可能的累積相位狀態(tài)如圖2 所示。
圖2 MSK 相位累積狀態(tài)
在MSK 調(diào)制信號(hào)數(shù)據(jù)塊{xn},n=0,1,…,N-1 中,xn在(1,-1)之間。假設(shè)初始相位為0,CP 長(zhǎng)度是G,那么必須在n=N-G的時(shí)刻將相位歸零。根據(jù)相位歸零法,因?yàn)镸SK 信號(hào)是二進(jìn)制全響應(yīng)信號(hào),記憶長(zhǎng)度L=1,因此只需要1 個(gè)長(zhǎng)度的尾符號(hào)使相位回到0 狀態(tài)。對(duì)于頻域均衡來(lái)說(shuō),N點(diǎn)FFT 變換需要進(jìn)行N2運(yùn)算,共需要進(jìn)行Nlog2N次的運(yùn)算,運(yùn)算復(fù)雜度大大降低;所以選擇N為2 的冪次,即選擇的數(shù)據(jù)符號(hào)為偶數(shù)。這樣N-L(L=1)為奇數(shù),在N-1 時(shí)刻累積相位到達(dá)1 或3 狀態(tài),只需通過(guò)1 個(gè)尾符號(hào)xl可以讓累積相位回到0 狀態(tài)[10]。
圖3 是MSK 信號(hào)加了循環(huán)前綴的完整的幀結(jié)構(gòu)。假設(shè)MSK 相位網(wǎng)格路徑開(kāi)始并結(jié)束于相位為0 的狀態(tài),n=0 時(shí)相位狀態(tài)為0。在n=N-G-1 時(shí)刻后,需要一個(gè)尾符號(hào)xl1,使在n=N-G時(shí)刻相位狀態(tài)回到0;在最后的n=N-2 時(shí)刻后,需要一個(gè)尾符號(hào)xl2,使在n=N-1 時(shí)刻相位狀態(tài)仍回到0;然后再將后面G個(gè)符號(hào)復(fù)制到前面作為循環(huán)前綴,此時(shí)構(gòu)成的MSK 信號(hào)相位是連續(xù)的。
圖3 MSK 信號(hào)加循環(huán)前綴完整幀
由于MSK 信號(hào)的連續(xù)相位變化特性,本文不能直接對(duì)接收到的采樣點(diǎn)進(jìn)行FFT 運(yùn)算,這會(huì)給MSK 單載波均衡帶來(lái)很大困難。CPM 信號(hào)運(yùn)用到單載波頻域均衡的離散線性算法包括Gram-Schmidt正交分解法和Laurent 分解法。對(duì)于MSK 信號(hào)來(lái)說(shuō),在AWGN 信道下,Gram-Schmidt 正交分解法的復(fù)雜程度略高于Laurent 分解法[11],因此本文采用運(yùn)算復(fù)雜程度相對(duì)較低的Laurent 分解法。
Laurent 分解代表了將MSK 信號(hào)所形成的復(fù)基帶信號(hào)按相應(yīng)的條件分解為2L-1 個(gè)脈沖幅度調(diào)制(Pulse Amplitude Modulation,PAM)脈沖信號(hào)的線性組合(記憶長(zhǎng)度為L(zhǎng))。MSK 是二進(jìn)制全響應(yīng)的CPM 信號(hào),其調(diào)制指數(shù)h=1/2,且頻率成型函數(shù)g(t)是持續(xù)時(shí)間為T(mén) 的矩形脈沖,脈沖信號(hào)的記憶長(zhǎng)度L=1,在這種情況下,MSK 信號(hào)復(fù)基帶可表示為只有一個(gè)Laurent 波形。
MSK 信號(hào)的Laurent 分解可以表示為:
復(fù)數(shù)據(jù)序列為:
式(7)中的yn可以用遞推形式表示為:
MSK 信號(hào)經(jīng)過(guò)Laurent 分解線性展開(kāi)之后就可以應(yīng)用FFT 算法,進(jìn)而可以完成頻域均衡。
由于甚低頻信道的特點(diǎn),使用時(shí)域均衡會(huì)有很大的復(fù)雜性,不利于系統(tǒng)實(shí)現(xiàn)。本文研究的SCFDE 線性均衡技術(shù)與OFDM 系統(tǒng)和單載波傳輸系統(tǒng)相比,實(shí)現(xiàn)簡(jiǎn)單,成本低大大降低,可有效對(duì)抗甚低頻MSK 信號(hào)中的碼間干擾;但由于MSK 調(diào)制信號(hào)是全響應(yīng)的連續(xù)相位信號(hào),需要插入尾符號(hào),并進(jìn)行Laurent 分解,進(jìn)而應(yīng)用SC-FDE 技術(shù),從而達(dá)到良好的均衡效果[12]。