張寶榮1,2，金森林1

　?。?. 燕山大學(xué) 信息科學(xué)與工程學(xué)院, 河北 秦皇島 066004；2. 河北省特種光纖與光纖傳感重點實驗室, 河北 秦皇島 066004）

       摘要：提出了一種利用壓縮感知原理測量微弱信號的方法，測量信號由偽隨機(jī)序列調(diào)制，應(yīng)用改造的測量矩陣，在一次測量基礎(chǔ)上進(jìn)行二次測量，利用壓縮感知的恢復(fù)算法可以精確地確定信號在字典中的位置并且得到其幅度值。仿真實驗證明本文方法可以用于檢測信噪比高于-20 dB的微弱信號，將信號較完整地從噪聲干擾中恢復(fù)出來，信號幅度誤差很小。

　　關(guān)鍵詞：壓縮感知; 信號處理; 微弱信號檢測; 信噪比

　　中圖分類號：TN911.72文獻(xiàn)標(biāo)識碼：ADOI： 10.19358/j.issn.1674-7720.2017.08.024

　　引用格式：張寶榮，金森林.基于壓縮感知的微弱信號檢測方法［J］.微型機(jī)與應(yīng)用，2017,36（8）：76-78，82.

0引言

　　隨著現(xiàn)代信息技術(shù)的迅猛發(fā)展，人們不僅對信息量的要求劇增，而且對信息準(zhǔn)確度的要求也越來越高。在信號測量過程中，信號經(jīng)常會被湮沒在測量系統(tǒng)的噪聲和干擾中，如何準(zhǔn)確地測量這些微弱信號是一個亟需解決的問題。

　　圖1壓縮感知微弱信號測量系統(tǒng)由于信號、噪聲和干擾的特點不同，檢測方法也迥異。目前，微弱信號的檢測技術(shù)［12］主要有三種：(1) 相干檢測與頻域窄帶化技術(shù)，利用噪聲與信號符合同頻又同相的可能性很小的特性檢測微弱信號，鎖相放大器是這種方法的典型［3］；(2) 時域信號的積累平均法，基于信號的穩(wěn)定性和噪聲的隨機(jī)性，由于信號周期性地重復(fù)，噪聲卻不具有這個特性，因此只要重復(fù)的次數(shù)足夠多，就可以把噪聲中的微弱信號提取出來，而且重復(fù)的次數(shù)越多，提取微弱信號的能力越強(qiáng)［4］；(3) 新興的利用混沌理論和小波技術(shù)來檢測微弱信號［57］。

　　2006年，DONOHO L、CZNDS J和TAO T等人提出了壓縮感知(Compressive Sensing)理論［89］，這是一種全新的信號獲取和處理理論。壓縮感知理論認(rèn)為只要信號在某已知變換域具有稀疏性，就可以通過原信號在某投影域的投影近似無損地重構(gòu)原始信號，只要求投影域的基與已知變換域的基是不相干的。利用壓縮感知方法不僅可以對信號進(jìn)行采樣和壓縮，而且能夠去除微弱信號噪聲，更加精確地實現(xiàn)了原始信號的恢復(fù)和重構(gòu)。

　　本文在傳統(tǒng)微弱信號去噪方法的基礎(chǔ)上，基于壓縮感知的思想，將微弱信號與改進(jìn)的測量矩陣做相關(guān)，得到相關(guān)系數(shù)，然后通過OMP算法［10］將微弱信號重構(gòu)。

1應(yīng)用壓縮感知檢測微弱信號原理

　　為了測量系統(tǒng)的某些特性，把已知信號輸入到待測系統(tǒng)中，系統(tǒng)會影響信號的某些參數(shù)，測量輸出信號，比較輸出與輸入信號的差異，可以得到系統(tǒng)的特性。在信號通過系統(tǒng)時，系統(tǒng)內(nèi)的噪聲和系統(tǒng)外的干擾都會疊加到輸出信號中，降低了輸出信號的信噪比，增大了測量誤差。如何在低信噪比條件下準(zhǔn)確地測量信號是關(guān)鍵問題，因此設(shè)計了應(yīng)用壓縮感知原理測量微弱信號的測量系統(tǒng)，見圖1。

　　測量所用輸入信號是正交字典稀疏表示的信號x，設(shè)正交字典為Ψ，其稀疏表示為θ，即x=Ψθ。輸入信號x通過待測量的系統(tǒng)后的輸出信號x′=Ax+n，其中n表示測量過程中疊加的噪聲和干擾，A是待測系統(tǒng)的增益。壓縮感知的測量過程是利用測量矩陣Φ來獲取信號的信息，測量結(jié)果為y=Φx，信號的重構(gòu)依據(jù)求解最優(yōu)化問題minθ0s.t. y=ΦΨθ。如果選擇合適的輸入信號x和測量矩陣Φ，通過兩次測量的方法就可以將輸出信號中的噪聲和干擾去除，得到高信噪比的信號Ax。

　　1.1信號與正交字典的選取

　　本文選取了偽隨機(jī)序列中的m序列 ［11］作為測量原始信號。正交字典是一組正交基，任意信號都可以用正交字典來表示，若其表示中僅有幾項不為零, 則稱為信號的稀疏表示。根據(jù)m序列的性質(zhì)，m序列的循環(huán)移位仍然是m序列，再根據(jù)m序列的自相關(guān)和互相關(guān)特性，m序列及其N-1個循環(huán)移位序列可以組成一個正交字典。即：

　　Ψ=(m,m1,m2…,mN－1)T

　　其中mi,i=1,2,…,N－1表示m序列的循環(huán)移位序列，i是位移數(shù)。

　　1.2測量矩陣

　　測量矩陣Φ有多種選擇，只要滿足RIP條件［12］即可，其中隨機(jī)高斯矩陣是最常見的一種。測量矩陣Φ的行數(shù)m與信號的稀疏程度k有關(guān)，m應(yīng)該滿足m≥klog2(N/k)。

　　由于被檢測信號是微弱信號，信噪比低，當(dāng)測量矩陣為高斯矩陣時，不能夠保證高概率地從隨機(jī)投影獲得的少量測量值中將特征量準(zhǔn)確地提取出來。因此，測量矩陣Φ需要改造，將高斯矩陣中的任一行替換為原始信號，不失一般性，可以選取第一行。為了使測量矩陣Φ滿足RIP條件，替換信號kx，必須滿足k1，測量矩陣Φ成為以下結(jié)構(gòu)：

　　Φ=(kx,g2,g3…,gm)T

　　其中g(shù)i,i=2,3,…，m是隨機(jī)高斯矩陣的各行。

　　1.3信號處理方法

　　為了從噪聲中恢復(fù)信號，采用的信號處理方法是，首先應(yīng)用OMP算法從測量信號y1=Φ1x′中恢復(fù)含噪聲的信號x′的稀疏表示θ1，根據(jù)OPM算法的特性得到的θ1應(yīng)該是以下形式：

　　θ1=(kA,…,θi,…,θj,…0)T

　　其中kA為信號的稀疏表示系數(shù)，θi,θj表示噪聲不為零的稀疏表示系數(shù)。所以含噪聲的信號的第一次信號恢復(fù)結(jié)果為：x1=kAx+∑m=i,j,…mθmgm。求和項表示噪聲中與測量高斯矩陣最相似的那些向量。

　　為了除去噪聲，重新選取隨機(jī)高斯矩陣Φ，替換第一行后形成新的測量矩陣Φ2，對第一次恢復(fù)信號x1再進(jìn)行測量，y2=Φ2x1。對測量結(jié)果y2應(yīng)用OPM算法再次進(jìn)行恢復(fù)稀疏表示θ2，則：

　　θ2=(k2A,…,0)T

　　除去第一項其他項都為零，原因是新選取的高斯矩陣與原高斯矩陣相同的概率極小，即一次恢復(fù)信號x1中的噪聲既然與一次測量的高斯矩陣中的某些向量相似，則不可能與二次測量矩陣中的向量相似，所以稀疏表示的系數(shù)都是零。

　　恢復(fù)后的信號為x2=k2Ax，已經(jīng)除去了噪聲和干擾，k為已知，因此可以得到去噪的信號Ax。

　　2仿真與實驗

　　為了驗證上述理論分析，首先進(jìn)行了仿真，線性反饋移位寄存器是產(chǎn)生m序列的最簡便的方法，n位移位寄存器產(chǎn)生長度為2n-1的m序列。本文選取由10位移位寄存器產(chǎn)生的長度為1 023的m序列，反饋系數(shù)為0,0,1,0,0,0,0,0,0,1。序列如圖2所示。

　　當(dāng)被測信號通過檢測系統(tǒng)時，檢測系統(tǒng)中的噪聲會加到信號上，本文選用的是加性高斯噪聲，即含噪信號為x′=aamx+nam。共中aam是信號的幅值，nam為高斯噪聲的幅值。通過調(diào)整aam，nam的大小，可以得到不同信噪比的含噪信號如圖3所示。

　　根據(jù)上節(jié)的方法，采用的是OMP算法進(jìn)行信號的恢復(fù)，α=0.001，不同信噪比的信號恢復(fù)結(jié)果如圖4所示。在信噪比分別為-10 dB和-20 dB時，第一次測量所得到的特征量的位置由于噪聲的干擾，隨噪聲的增大而不準(zhǔn)確。而在進(jìn)行第二次測量后，噪聲特征量由于兩次測量的測量矩陣是隨機(jī)產(chǎn)生的，使得噪聲特征量的重合概率幾乎為0，但是信號特征量的位置依然不變，因此可以在第圖4二次測量后的信號恢復(fù)二次測量后準(zhǔn)確地找到信號特征量的位置，并獲取其幅值。

　　圖4中信號的幅值還存在一定的誤差，當(dāng)信噪比為-10 dB和-20 dB時，測量矩陣在不滿足RIP條件下和滿足RIP條件［13］下恢復(fù)信號的誤差如圖5所示。由圖5可知，在滿足RIP條件的情況下，對于信噪比為-20 dB的含噪聲信號，通過OMP算法，能夠準(zhǔn)確地找到特征量的位置，恢復(fù)信號的幅度誤差較小，而測量矩陣不滿足RIP條件的情況下，恢復(fù)信號誤差較大，當(dāng)信噪比再降低時，不能夠準(zhǔn)確地找到特征量的位置。

　　利用LabVIEW構(gòu)建了實驗系統(tǒng)，噪聲是電阻的熱噪聲，m序列由虛擬儀器產(chǎn)生，控制信號幅度，二者相加可以得到不同信噪比的被測信號，如圖6所示。

　　圖6含噪聲的m序列 (信噪比-10 dB)LabVIEW系統(tǒng)采集被測信號，應(yīng)用OMP算法進(jìn)行信號恢復(fù)，結(jié)果如圖7所示，與仿真結(jié)果是一致的。對不同信噪比的信號進(jìn)行了多次實驗，信號的恢復(fù)誤差見表1，隨著信噪比的降低，恢復(fù)信號的誤差增大，當(dāng)信噪比低于-25 dB后，恢復(fù)信號的誤差劇增，本文的方法已經(jīng)不適用。

3結(jié)論

　　本文簡要介紹了一種新的檢測微弱信號的方法，基于壓縮感知原理，根據(jù)噪聲和信號本身具有的特性，通過改造測量矩陣和在原有的測量基礎(chǔ)上再進(jìn)行一次測量，由重構(gòu)算法找到特征量的位置并測出其幅值的大小，進(jìn)而將原始信號恢復(fù)出來。通過仿真和實驗證實，在信噪比不低于-20 dB的情況下，重構(gòu)信號誤差較小，若在測量矩陣的選取和信號的處理方法上再做進(jìn)一步的改進(jìn)，有望在更低的信噪比條件下更精確地恢復(fù)信號。

　　參考文獻(xiàn)

　?。?］ 陳佳圭. 微弱信號檢測［M］. 北京：中央廣播電視大學(xué)出版社, 1987.

　　［2］ 夏均忠, 劉遠(yuǎn)宏, 冷永剛,等. 微弱信號檢測方法的現(xiàn)狀分析［J］. 噪聲與振動控制, 2011,31（3）:156-161.

　?。?］ 朱虹，林君，吳忠杰，等. 近紅外光譜儀中的數(shù)字鎖相技術(shù)研究［J］. 儀器儀表學(xué)報，2006，27(10): 1258-1261.

　　［4］ 王玉田，監(jiān)雄，王惠新，等. 基于Boxcar濾波器的熒光檢測二氧化硫信號的去噪研究［J］. 光譜學(xué)與光譜分析， 2012,32(12):3285-3289.

　?。?］ Xu Yanchun, Yang Chunling, Qu Xiaodong. A new approach to detecting weak signal in strong noise based on chaos system control［J］. Chinese Physics B， 2010, 19(3): 202-206.

　?。?］ 范劍，趙文禮，王萬強(qiáng). 基于Duffing振子的微弱周期信號混沌檢測性能研究［J］.物理學(xué)報， 2013,62(18):180502.