自適應(yīng)譜線增強(qiáng)器在生物雷達(dá)中的應(yīng)用

陳鐵軍 付瑞玲

(鄭州大學(xué)電氣工程學(xué)院，河南 鄭州 450001)

0 引言

生物雷達(dá)是地震、災(zāi)害救援的專用設(shè)備，它不僅能從攜帶有微弱人體生命體征(呼吸、體動等)的雷達(dá)回波信號中檢測出人體的生命參數(shù)，而且還可以穿透一定的介質(zhì)(磚墻或廢墟)，并間隔一定的距離探測人體的生命體征，發(fā)現(xiàn)存活、被困人員，使受困人員得到及時營救［1－3］。然而，雷達(dá)在接收生命信號的同時不可避免地會采集到諸如背景噪聲、地雜回波甚至災(zāi)害源等信號，這些信號容易與選擇的生命信號產(chǎn)生交疊。因此，對采集信號處理的效果直接決定了生命探測儀性能的優(yōu)劣。

為解決這一問題，基于變步長自適應(yīng)算法的各種信號處理的方法不斷涌現(xiàn)［4－7］。傳統(tǒng)的變步長自適應(yīng)算法均需要預(yù)先給定參考信號。但在救援環(huán)境復(fù)雜以及廢墟下存活人體的個體差異，在生物雷達(dá)信號處理中很難找到一個理想的信號作為參考信號。為解決這一問題，本文提出了基于誤差變化率的變步長最小均方誤差(least mean square，LMS)自適應(yīng)譜線增強(qiáng)算法。該算法不需要預(yù)先給定參考信號，而是以輸入信號本身的延遲作為參考信號［8－14］。

1 自適應(yīng)譜線增強(qiáng)器原理

自適應(yīng)譜線增強(qiáng)器最早是由Widrow等人于1975年提出的。目前，基于自適應(yīng)線性組合器的自適應(yīng)譜線增強(qiáng)器已廣泛應(yīng)用于頻譜估算、譜線估計(jì)以及窄帶檢測等領(lǐng)域。在窄帶信號加上寬帶信號的情況下，自適應(yīng)譜線增強(qiáng)方法無需獨(dú)立的參考信號就能將信號分離出來，從而提高了生命探測系統(tǒng)對微弱信號的檢測能力，其原理如圖1所示。

圖1 自適應(yīng)譜線增強(qiáng)器原理圖Fig.1 The principle of adaptive spectral enhancer

如果在輸入端加入的信號x(n)是1個窄帶信號sN(n)和1個寬帶噪聲sB(n)的混合，由于窄帶信號的自相關(guān)函數(shù)比寬帶噪聲自相關(guān)函數(shù)的時間相關(guān)半徑要短。因此，當(dāng)延遲時間Δ小于寬帶噪聲的時間相關(guān)半徑而大于窄帶信號的時間相關(guān)半徑時，寬帶噪聲sB(n)將與sB(n－Δ)不相關(guān)，而窄帶信號sN(n)與sN(n－Δ)仍然相關(guān)，因而自適應(yīng)濾波器的輸出將是窄帶信號sN(n)的最佳估計(jì)。sN(n)+sB(n)與sN(n)相減后得到的是sB(n)的最佳估計(jì)，從而能將窄帶信號sN(n)與寬帶噪聲sB(n)分離開來。這樣，自適應(yīng)濾波器的輸出y(n)就是所需要的有用信號。

固定步長的LMS算法迭代公式為:

式中:x(n)為輸入信號;d(n)為x(n)的延遲，作為參考信號;e(n)為誤差信號;W(n)為權(quán)重系數(shù);μ(n)為步長。

2 改進(jìn)的譜線增強(qiáng)算法

初始收斂速度、時變系統(tǒng)跟蹤能力以及穩(wěn)態(tài)失調(diào)是衡量自適應(yīng)濾波算法優(yōu)劣的三個最重要的技術(shù)指標(biāo)。在LMS算法中，最簡單的學(xué)習(xí)速率參數(shù)選擇是取μ(n)為常數(shù)，即:

式中:λmax為輸入信號自相關(guān)矩陣的最大特征值。

這種方法會引起收斂與穩(wěn)定性能的矛盾，即大的學(xué)習(xí)速率能夠提高濾波器的收斂速率，但穩(wěn)態(tài)性能就會隨之降低;反之，為了提高穩(wěn)態(tài)性能而采用小的學(xué)習(xí)速率時，跟蹤速度和收斂就會變慢。因此，學(xué)習(xí)速率的選擇應(yīng)該兼顧穩(wěn)態(tài)性能與收斂速率。一種簡單而有效的方法就是在不同的迭代過程中使用不同的學(xué)習(xí)速率參數(shù)，即采用時變的學(xué)習(xí)速率。

文獻(xiàn)［8］～［11］給出了幾種變步長自適應(yīng)算法，基本解決了固定步長存在的矛盾。以上所有算法均需要理想信號作為參考信號，從而得到越來越小直至為0的誤差信號來控制μ(n)的變化［8－11］。但是自適應(yīng)譜線增強(qiáng)系統(tǒng)的理想誤差并不趨于0，而是近似于寬帶噪聲，所以即便系統(tǒng)已經(jīng)完全跟蹤上時，步長因子μ(n)并不為0，仍將繼續(xù)迭代。

本文提出了一種基于誤差變化率的變步長自適應(yīng)LMS譜線增強(qiáng)算法。該算法采用誤差變化率來控制步長因子μ(n)的變化，以保證當(dāng)系統(tǒng)達(dá)到穩(wěn)態(tài)時停止迭代。步長更新公式如下:

式中:k為步長控制因子;U(n)為步長向量;W(n)=［ω1(n)ω2(n)…ωM(n)］為濾波器的權(quán)系數(shù)，其中M為濾波器的階數(shù);α為控制函數(shù)形狀的常數(shù);m為調(diào)整因子;b(n)為誤差變化量的絕對值。

3 算法收斂性能分析

步長μ(n)與各參數(shù)的關(guān)系曲線如圖2所示。

圖2 步長μ(n)與各參數(shù)的關(guān)系曲線Fig.2 The relative curves between step size μ(n)and various parameters

從圖2(a)可以看出，當(dāng)m選擇過小時，曲線比較尖銳，b(n)在接近0處變化時(算法已達(dá)到或?qū)⒁_(dá)到穩(wěn)態(tài))，μ(n)變化太大，不具有緩慢變化的特性，易造成較大的穩(wěn)態(tài)噪聲;當(dāng)m選擇過大時，曲線底部形態(tài)比較平緩，會出現(xiàn)b(n)還未接近0處變化時μ(n)已成為0的情況，這樣又可能造成較大的穩(wěn)態(tài)噪聲。從圖2(b)可以看出，k從整體上影響算法的收斂速度。當(dāng)k選擇過大時，步長調(diào)整過早進(jìn)入緩慢變化區(qū)域;當(dāng)k選擇過小時會增大穩(wěn)態(tài)誤差。從圖2(c)可以看出，α控制函數(shù)的形狀。當(dāng)α選擇過大時，誤差變化率接近0仍有較大步長，穩(wěn)態(tài)誤差增大;當(dāng)α選擇過小時，步長較小且變化緩慢，收斂速度降低。因此，參數(shù)m、k和α應(yīng)根據(jù)具體的系統(tǒng)環(huán)境與要求進(jìn)行選擇。

4 算法應(yīng)用

因?yàn)樯锢走_(dá)探測到的人的呼吸信號比較規(guī)律且接近正弦波［12］，所以本文采用正弦信號來模擬被探測人員的呼吸信號。通過對所采集數(shù)據(jù)的分析，雜波可以近似用高斯雜波表示。根據(jù)正常人的呼吸頻率區(qū)間，設(shè)被探測人員的呼吸信號為被高斯白噪聲污染的正弦信號。該信號的信噪比為5 dB、頻率為0.4Hz的正弦信號。為了取得較小的穩(wěn)態(tài)失調(diào)量和較短的收斂時間，各參數(shù)設(shè)置如下:α=600、k=2、m=4、Δ=20，同時，設(shè)置濾波器階數(shù)為100。

文獻(xiàn)［11］和本文算法提取的正弦信號波形如圖3所示，兩種算法的收斂曲線比較如圖4所示。

由圖3、圖4可知，本文算法具有良好的跟蹤性能、較小的穩(wěn)態(tài)誤差以及良好的收斂性能。

5 結(jié)束語

本文分析了傳統(tǒng)的LMS算法在生物雷達(dá)信號處理過程中的應(yīng)用，根據(jù)固定步長算法以及文獻(xiàn)中的變步長算法的不足，提出了一種新的變步長自適應(yīng)LMS譜線增強(qiáng)算法。通過仿真試驗(yàn)可以看出，當(dāng)系統(tǒng)發(fā)生突變時，該算法也能快速收斂，具有很好的魯棒性，達(dá)到了很好的效果;在低信噪比的條件下，該算法也具有很好的跟蹤性能。

［1］Linjc K，Wollschlaeger P.Microwave apexcardiography［J］.IEEE Transaction on Microwave Theory Technology，1979，27(6):618 －620.

［2］Misra D K.Scattering of electromagnetic waves by human body and its applications［D］.East Lansing:Michigan State University，1984:292－297.

［3］倪安勝.傷員探測中生命體識別技術(shù)研究［D］.西安:第四軍醫(yī)大學(xué)，2003:13－15.

［4］付玉峻，商文忠，張晗.一種自適應(yīng)方向性技術(shù)在生命探測儀中的應(yīng)用［J］.武警學(xué)院學(xué)報(bào)，2008(24):34 －35.

［5］張楊，王建琪，荊西京.基于自適應(yīng)濾波的生物雷達(dá)干擾抑制方法［J］.儀器儀表學(xué)報(bào)，2009，30(4):807 －811.

［6］夏淑芳.自適應(yīng)濾波技術(shù)在雷達(dá)式生命探測儀中的應(yīng)用研究［J］.科技風(fēng)，2009(13):260.

［7］張陸游，張永順，潘克戰(zhàn).變步長自適應(yīng)算法在雷達(dá)旁瓣對消系統(tǒng)中的應(yīng)用［J］.空軍工程大學(xué)學(xué)報(bào):自然科學(xué)版，2008(4):25－28.

［8］Yasukawa H，Shimada S，F(xiàn)urukrawa I，et al.Acoustic echocanceller with high speech quality［C］//ICASSP，1987:2125 －2128.

［9］葉華，吳伯修.變步長自適應(yīng)濾波算法的研究［J］.電子學(xué)報(bào)，1990，18(4):63 －69.

［10］罩景繁，歐陽景正.—種新的變步長自適應(yīng)濾波算法［J］.數(shù)據(jù)采集與處理，1997，12(3):171 －194.

［11］高鷹，謝勝利.一種變步長LMS自適應(yīng)濾波算法及分析［J］.電子學(xué)報(bào)，2001，29(7):1094 －1097.

［12］李杰，王健琪，荊西京，等.非接觸生命參數(shù)檢測系統(tǒng)動/靜目標(biāo)的識別技術(shù)［J］.第四軍醫(yī)大學(xué)學(xué)報(bào):自然科學(xué)版，2005，26(20):1915－1918.

［13］黃春艷，張?jiān)迄i，黃紅艷.基于灰色關(guān)聯(lián)度的圖像混合噪聲的自適應(yīng)濾波算法［J］.微電子學(xué)與計(jì)算機(jī)，2010，27(2):126 －128.

［14］畢國堂，王曉輝，周燕，等.基于粗集的圖像混合噪聲濾波算法［J］.計(jì)算機(jī)工程與設(shè)計(jì)，2009，30(21):4898－4900.