易　輝，侯孝民，馬　宏，吳　濤

(裝備學院，北京　101416)

?

一種基于譜分析的非帶限衛(wèi)星信號信噪比估計算法

易輝，侯孝民，馬宏，吳濤

(裝備學院，北京101416)

摘要：為了實現(xiàn)對非帶限衛(wèi)星信號信噪比的快速高精度估計，分析了傳統(tǒng)的基于譜分析的信噪比估計算法的對非帶限衛(wèi)星信號估計不準的原因，提出了一種新的基于譜分析的信噪比估計算法；該方法利用主瓣范圍內(nèi)信號功率與信號總功率之比和主瓣范圍內(nèi)噪聲功率與噪聲總功率之比求得信號功率和噪聲功率，從而得出信噪比；對BPSK和QPSK信號的仿真結(jié)果表明：新算法性能穩(wěn)定，可有效提高估計精度，在-10～10 dB的范圍內(nèi)，估計的偏差和均方根誤差基本都小于0.5 dB；提出的新的基于譜分析的信噪比估計算法可用于各種調(diào)制方式，計算復雜度小，可以滿足對非帶限衛(wèi)星信號信噪比估計的需求。

關鍵詞：信噪比估計；譜分析；非帶限；衛(wèi)星信號

信噪比(SNR，Signal to Noise Ratio)是衡量衛(wèi)星通信信號品質(zhì)的重要指標之一。對于某一種調(diào)制方式，信噪比與誤碼率有著直接對應的關系[1]，通過估計接收信號的信噪比，可以有效監(jiān)測接收信號的信號品質(zhì)，為信號解調(diào)提供重要參考依據(jù)。此外，在通信系統(tǒng)中，很多環(huán)節(jié)都需要信噪比的先驗知識優(yōu)化通信性能，如自適應調(diào)制系統(tǒng)中調(diào)制和編碼方式的選擇，移動通信網(wǎng)絡中的越區(qū)切換、碼分多址系統(tǒng)中的功率控制[2-3]等等。

根據(jù)是否需要輔助數(shù)據(jù)，信噪比估計算法可以分為數(shù)據(jù)輔助類算法(DA:Data-aided)和非數(shù)據(jù)輔助類算法(NDA:Non Data-aided)兩類。DA估計算法準確性較高，但是需要犧牲信道傳輸效率為代價。NDA方法在系統(tǒng)正常傳輸信息數(shù)據(jù)的同時進行信噪比估計，不影響信息傳輸效率，適用范圍較廣。

常用的信噪比估計算法有最大似然估計法(ML)[4-5]、分離符號矩算法(SSME)[6-7]、二階矩四階矩估計算法(M2M4)[8-9]、子空間估計算法[10-11]和基于功率譜的估計算法[3,12]。最大似然估計法精度高，但需要周期地發(fā)送訓練序列，運算復雜，是一種DA算法；分離符號矩算法只能對判決域信號進行處理，是一種NDA算法；二階矩四階矩估計算法是基于統(tǒng)計量的NDA算法，其估計隨著觀測數(shù)據(jù)的增加會越來越接近真實值；子空間估計算法在高信噪比的情況下性能較好，在低信噪比情況下性能下降，且其運算量較大，是一種NDA算法；基于譜分析的估計算法可以用FFT快速實現(xiàn)，運算量少，復雜度低，且不需要輔助數(shù)據(jù)，是一種NDA算法，在工程上應用廣泛。

衛(wèi)星通信通常呈現(xiàn)出信噪比低，數(shù)據(jù)量大的特點，要實現(xiàn)對接收衛(wèi)星信號的實時監(jiān)測和準確估計，基于譜分析的信噪比估計算法是一種可選方案。但是，對于非帶限的衛(wèi)星信號，傳統(tǒng)的基于譜分析的估計算法存在缺陷。文章通過分析傳統(tǒng)算法的不足，提出一種新算法，通過利用主瓣范圍內(nèi)信號功率與信號總功率之比和主瓣范圍內(nèi)噪聲功率與噪聲總功率之比來求得信號功率和噪聲功率，從而改善了基于譜分析的信噪比估計算法對非帶限衛(wèi)星信號信噪比估計的性能。

1傳統(tǒng)的基于譜分析的信噪比估計算法

傳統(tǒng)的基于譜分析的信噪比估計算法[2,3,12]對白噪聲信道下的帶限信號有較好的適用性。白噪聲信道下的接收信號中，噪聲功率譜均勻分布在整個觀測頻帶內(nèi)，而帶限信號的功率譜只在有限的頻帶范圍內(nèi)存在。利用這一特性，先對接收到的時域信號進行功率譜估計，根據(jù)功率譜估計出信號的頻帶范圍，將信號范圍外的功率譜當作噪聲譜，計算信號范圍外的噪聲譜的均值作為噪聲功率譜密度的估計值。然后用噪聲功率譜密度乘以觀測帶寬得到噪聲功率，用接收信號總功率減去噪聲功率得到信號功率，兩者相比就可以得到SNR估計值。

假設以 fS采樣率對接收信號進行采樣，得接收信號

(1)

其中，S(N)為發(fā)送信號，E(N)為方差為σ2的零均值白噪聲。對離散信號x(N)進行功率譜估計，得到接收信號的功率譜

(2)

其中，Px(f)表示接收信號的功率譜，PS(f)表示信號的功率譜，PE(f)表示噪聲功率譜。由于正負功率譜對稱，因此只需計算正半軸。

假設接收信號功率譜的正半軸如圖 1所示，虛線以下表示理想白噪聲的功率譜，虛線以上表示帶限信號的功率譜。經(jīng)過估計得到帶限信號頻帶范圍為[fL,fH]，根據(jù)信號頻帶范圍外的功率譜可以估計出噪聲功率譜密度：

(3)

從而可以估計出信噪比：

(4)

圖1　帶限信號功率譜圖

基于譜分析的估計算法的關鍵在于估計信號的頻帶范圍。文獻[2]中提出了一種利用平均功率譜密度單調(diào)性求頻帶范圍的方法，在頻帶邊緣下降陡峭時效果較好，在頻帶邊緣下降緩慢時會帶來誤差。對于帶限信號，求出有誤差的頻帶范圍后，將頻帶范圍適當放寬就可以較好地將信號限制在估計的頻帶范圍內(nèi)，從而估計帶外的噪聲功率譜密度。對于非帶限信號，為了準確估計帶外的噪聲功率譜密度，需要準確地估計信號的頻帶范圍，而非帶限信號的旁瓣相對主瓣較小，幅度下降也比較緩慢，再用上述方法估計頻帶范圍時會帶來較大誤差。此外，這種方法較為復雜，計算量較大。

對于非帶限信號，另一種簡單的做法是將信號的主瓣(無噪聲情況下的第一過零點)作為信號的頻帶范圍，把主瓣外的旁瓣信號當成噪聲處理，這是一種有偏估計。這種做法雖然簡單，但也會像文獻[2]方法那樣帶來誤差，甚至誤差會超過文獻[2]。這種做法在信噪比較低時，可以對接收信號的信噪比進行近似估計，但是當信噪比較高時，將會帶來較大的誤差。

通過分析可以發(fā)現(xiàn)，傳統(tǒng)算法之所以對非帶限信號估計性能不好，是因為在計算噪聲功率譜密度的時候，沒有準確估計信號頻帶范圍，將部分信號功率當作了噪聲功率進行計算，這就導致估計出來的噪聲功率偏大，而信號功率偏小，最終導致信噪比估計值比真實值偏小。

2新的基于譜分析的信噪比估計算法

針對傳統(tǒng)的基于譜分析的估計算法對非帶限信號估計不準的問題，提出一種新的基于譜分析的信噪比估計算法。假設信號S(N)為非帶限信號，理想的接收信號x(N)的功率譜如圖2所示，虛線以下表示理想白噪聲的功率譜，虛線以上表示信號的功率譜。

圖2　非帶限信號功率譜圖

同樣，對于接收信號先做功率譜估計，得到功率譜 Px(f)，由功率譜估計信號主瓣的頻帶范圍[fL,fH]。在確定信號主瓣頻帶范圍時，方法是通過對功率譜做差分得功率譜差分圖，如圖3所示。搜索功率譜差分圖的最大值，向負方向遍歷，求得第一零點的橫坐標作為 fL；搜索功率譜差分圖的最小值，向正方向遍歷，求得第一零點的橫坐標作為 fH；這種方法求出的是功率譜圖的第一凹點，也就是無噪聲情況下的第一過零點。

圖3　非帶限信號功率譜差分圖

在此基礎上，計算出信號主瓣頻帶范圍內(nèi)信號與噪聲的功率之和：

(5)

用接收信號的總功率減去信號主瓣頻帶范圍內(nèi)信號與噪聲的功率之和，得到主瓣頻帶范圍外信號與噪聲的功率之和：

(6)

于是，可列出如下方程組：

(7)

其中，PS表示信號功率；PN表示噪聲功率；λ為信號主瓣頻帶范圍與總的觀測帶寬之比，也就是主瓣頻帶范圍內(nèi)噪聲功率與噪聲總功率的比值；η為無噪聲條件下信號的主瓣功率與信號總功率之比值，該比值只與信號的碼元波形有關，可根據(jù)不同碼元的功率譜計算。如，碼元波形為矩形脈沖的信號主瓣功率與信號總功率之比約為0.902 8。而在通信系統(tǒng)中，信號的碼元波形一般為已知條件。

解方程組，可得信噪比估計值：

(8)

新算法的本質(zhì)是利用信號主瓣范圍內(nèi)與主瓣范圍外的功率之比和噪聲主瓣范圍內(nèi)與主瓣范圍外功率之比來求得信號功率和噪聲功率。新算法通過方程組求解信號功率和噪聲功率，避免了因信號頻帶范圍估計不準確而帶來誤差的問題，在功率譜估計理想的情況下，可以準確求出信號和噪聲功率，不存在信號功率和噪聲功率劃分不清的問題，因而是一種無偏估計。即便噪聲不理想導致估計頻帶范圍存在一定的偏差，由于是在第一過零點附近，帶內(nèi)信號功率和帶外信號功率的比值也不會偏離真實值太多，而噪聲功率的比值并不會因頻帶范圍估計不準受到影響，因此新算法對頻帶范圍的估計誤差有一定的容忍性。

新算法也可用于處理帶限信號。用新算法處理帶限信號時，η取1即可，這種情況下，新算法與傳統(tǒng)算法是等效的。因此，傳統(tǒng)算法可以看作是新算法的一種特殊情況。

算法步驟歸納如下：

3仿真分析

本節(jié)主要針對BPSK和QPSK等常用非帶限衛(wèi)星信號進行仿真驗證，分析算法隨信噪比變化的趨勢以及適用的信噪比范圍。仿真中主要用估計均值和均方根誤差評價算法的估計性能。信噪比估計均值由下式計算得到：

(9)

(10)

在對新算法仿真的同時，仿真了傳統(tǒng)的基于譜分析的信噪比估計算法(取主瓣作為頻帶范圍)的性能作為對比。仿真條件：1.4 Mbps的碼元速率，碼元波形為矩形脈沖，載波為7 MHz，采樣速率56 MHz，數(shù)據(jù)長度為4 096個采樣點。進行500次蒙特卡洛仿真。

圖4和圖5為信噪比估計的均值圖。從圖中可以看出，正如前面所分析的那樣，傳統(tǒng)的基于譜分析的估計算法在信噪比較高時算法性能嚴重下降；而新的譜分析估計算法在-15～25 dB求得的均值都優(yōu)于傳統(tǒng)算法，一直在真值附近擺動，且在-12～13 dB估計值的均值偏離真值不超過0.5 dB。

圖4　BPSK信號信噪比估計均值

圖5　QPSK信號信噪比估計均值

圖6和圖 7為信噪比估計均方根誤差曲線圖。從中可以看出，當信噪比低于-10 dB時，傳統(tǒng)算法和新算法的均方根誤差都急速下降，這是因為當信噪比過低時，信號頻帶邊界的判斷會因噪聲影響而呈現(xiàn)出隨機性，導致均方根誤差過大。當信噪比大于0 dB時，傳統(tǒng)算法的均方根誤差開始增大，而新算法由于更好的分辨了噪聲能量和信號能量，在信噪比大于10 dB時，均方根誤差才開始呈現(xiàn)明顯增大。新的基于譜分析的信噪比估計算法在-10～10 dB均方根誤差都基本小于0.5 dB。

綜上，可以看出，新的基于譜分析的信噪比估計算法在估計均值和均方根誤差性能上都較傳統(tǒng)算法有較大的提升。無論在信噪比較低還是較高的情況下，算法的估計精度都得到了提高，在-10～10 dB能穩(wěn)定精確地估計信號的信噪比，估計精度可達0.5 dB。算法對于BPSK和QPSK兩種信號的估計性能基本一致。

圖6　BPSK信號信噪比估計均方根誤差曲線

圖7　QPSK信號信噪比估計均方根誤差曲線

4小結(jié)

論文介紹了基于譜分析的信噪比估計算法的原理，分析了算法存在的缺陷及其對非帶限信號的不適用性，通過分析產(chǎn)生導致這種現(xiàn)象的原因，提出了一種新算法。新算法通過細致區(qū)分功率譜估計圖中的噪聲功率和信號功率，改善了基于譜分析的信噪比估計算法的性能。且算法適用于各種不同的調(diào)制方式，計算量小，滿足了對非帶限衛(wèi)星信號信噪比估計的需求。仿真結(jié)果驗證了新算法的正確性。

參考文獻：

[1]樊昌信，曹麗娜.通信原理[M].國防工業(yè)出版社,2010:212-213.

[2]彭耿，黃知濤，陸鳳波，等.中頻通信信號信噪比的快速盲估計[J].電子與信息學報，2010(1):102-106.

[3]張金成，彭華，趙國慶.信噪比估計算法研究[J].信息工程大學學報，2011(5):535-543.

[4]KERR R B.On Signal and Noise Level Estimation in a Coherent PCM Channel[J].IEEE Transactions on Aerospace and Electronic Systems.1966,AES-2(4):450-454.

[5]孟慶萍，周新力，田偉.基于先驗信息輔助的8PSK信號信噪比估計[J].北京郵電大學學報，2012(4):94-98.

[6]SIMON M K,MILEANT A.SNR estimation for the baseband assembly[R].The Telecommunications and Data Acquisition Progress Report 42-85,January-March 1986.1986:118-126.

[7]SIMON M,DOLINAR S.Improving signal-to-noise ratio estimation for autonomous receivers[J].The Interplanetary Network Progress Report，2004,42:159.

[8]BENEDICT T,SOONG T.The joint estimation of signal and noise from the sum envelope[J].IEEE Transactions on Information Theory，1967,13(3):447-454.

[9]MATZNER R.An SNR estimation algorithm for complex baseband signals using higher order statistics[J].Facta Universitatis (Nis)，1993,6(1):41-52.

[10]ANDERSIN M,MANDAYAM N B,YATES R D.Subspace based estimation of the signal to interference ratio for TDMA cellular systems[J].Wireless Networks，1998,4(3):241-247.

[11]詹亞鋒，曹志剛，馬正新.無線數(shù)字通信的盲信噪比估計[J].清華大學學報(自然科學版)，2003,43(7):957-960.

[12]陳璐，畢大平，徐梁昊.時域重疊多輸入/輸出雷達信號分離算法[J].探測與控制學報，2015(5)：51-56.

[13]胡冰舟，張蓉，雷維嘉，等.一種基于平均周期圖的頻域信噪比估計算法[J].電訊技術，2014(10):1385-1390.

(責任編輯楊繼森)

本文引用格式：易輝，侯孝民，馬宏，等.一種基于譜分析的非帶限衛(wèi)星信號信噪比估計算法[J].兵器裝備工程學報，2016(5):104-107.

Citation format:YI Hui,HOU Xiao-min,MA Hong,et al.An SNR Estimator Based on Spectrum Analysis for Non-Band-Limited Satellite Signals[J].Journal of Ordnance Equipment Engineering,2016(5):104-107.

An SNR Estimator Based on Spectrum Analysis for Non-Band-Limited Satellite Signals

YI Hui,HOU Xiao-min,MA Hong,WU Tao

(Equipment Academy,Beijing 101416,China)

Abstract:In order to estimate the SNR (Signal to Noise Ratio)of non-band-limited satellite signals in high precision,a new SNR estimator based on spectrum analysis was proposed after an investigation of the traditional SNR estimator based on spectrum analysis.The new estimator utilized the ratio of signal power in main lobe to total signal power and the ratio of noise power in main lobe to total noise power to calculate signal power and noise power,and then got the SNR.Simulations for BPSK (Binary Phase Shift Keying)and QPSK(Quadrature Phase Shift Keying)signals indicate that the new estimator is stable and can improve the estimation precision effectively.The bias and root mean square error of the new estimator are almost less than 0.5 dB when the true SNR range is from -10 dB to 10 dB.The new estimator can be used to different modulations.It has low complexity and can satisfy the SNR estimation demand of non-band-limited satellite signals.

Key words:SNR estimation; spectrum analysis; non-band-limited; satellite signal

doi:【信息科學與控制工程】10.11809/scbgxb2016.05.025

收稿日期：2015-10-25；修回日期：2015-11-29

作者簡介：易輝(1991—)，男，碩士研究生，主要從事測控與通信技術研究。

中圖分類號：TN571

文獻標識碼：A

文章編號：2096-2304(2016)05-0104-05