崔付俊，劉世博，何立功，曹 霞，郭建中

(陜西師范大學(xué)物理學(xué)與信息技術(shù)學(xué)院，陜西西安 710062)

0 引 言

在超聲醫(yī)學(xué)中，生物組織對超聲效應(yīng)的受限程度使得超聲脈沖幅度不能太高。在傳統(tǒng)超聲診斷中，盡可能采用窄脈沖激勵(lì)寬帶超聲換能器來提高系統(tǒng)分辨率。為了提高信噪比和檢測深度，現(xiàn)有的實(shí)驗(yàn)都在最大限度地提高發(fā)射峰值聲功率。雖然峰值幅度達(dá)到了最大允許值，但平均聲功率還是很低。Matthew O'Donnell.等人發(fā)現(xiàn)[1]，超聲檢測峰值聲功率已經(jīng)很接近安全標(biāo)準(zhǔn)上限，可是平均聲功率卻還不到安全標(biāo)準(zhǔn)的1%。

編碼激勵(lì)和脈沖壓縮有效地解決了雷達(dá)系統(tǒng)中作用距離以及距離分辨力之間的矛盾，提高了平均發(fā)射功率[2]。1974年，Newhouse首次將編碼激勵(lì)引入到醫(yī)學(xué)超聲成像，并在此后的時(shí)間里做了大量的研究。文獻(xiàn)[3]表明，使用二維解析小波閾值對超聲脈沖反射波信號進(jìn)行濾噪，性能可以提高。傅娟等[4]提出相對幅度相等的包絡(luò)線性調(diào)頻編碼方法和預(yù)失真線性調(diào)頻編碼方法，不僅提高了軸向分辨力，而且使最大旁瓣幅度減小至?48 dB以下。劉波[5]通過仿真和實(shí)驗(yàn)對1～5周期載頻的13位巴克碼載頻調(diào)制信號激勵(lì)換能器，得到單位碼元載有3個(gè)周期載頻的激勵(lì)信號，通過換能器后，脈沖壓縮效果最佳。在超聲血流檢測與成像系統(tǒng)中，線性調(diào)頻信號能夠在不損失精度的前提下提高回波信號的信噪比，提高血流成像質(zhì)量[6]。Jennifer等[7]提出短促聲波激發(fā)的超聲導(dǎo)波可以減少采集時(shí)間，提高信噪比。魏東等[8]指出，改進(jìn)的非線性調(diào)頻脈沖壓縮方法可以有效提高脈沖壓縮效果，并且有利于改善空氣耦合超聲檢測的效果。

依據(jù)編碼調(diào)制及脈沖壓縮方法原理[9]，信號的距離分辨率取決于信號帶寬，單載頻脈沖信號的時(shí)寬T和帶寬B滿足TB=1，寬脈沖必然降低距離分辨率，系統(tǒng)的穿透力和分辨率存在矛盾[10]。在不增加峰值功率的前提下，編碼信號的TB積變大，增強(qiáng)了信號能量，提高了平均發(fā)射功率。對長脈沖回波信號進(jìn)行脈沖壓縮處理后，信號幅值大大增強(qiáng)，信噪比顯著提高。

本文考慮換能器頻帶特性，用適當(dāng)頻帶的線性調(diào)頻波激勵(lì)信號，研究換能器對激勵(lì)信號及其脈沖壓縮結(jié)果的影響。為了提高激勵(lì)換能器脈沖壓縮后的信噪比和分辨率，選用包絡(luò)調(diào)制線性調(diào)頻波激勵(lì)換能器。通過仿真和實(shí)驗(yàn)對比等幅及包絡(luò)調(diào)制兩種不同線性調(diào)頻波激勵(lì)換能器后的脈沖壓縮結(jié)果。

1 原 理

1.1 線性調(diào)頻信號

線性調(diào)頻信號是脈沖壓縮最常見的調(diào)制信號，其數(shù)學(xué)表達(dá)式為

1.2 脈沖壓縮

脈沖壓縮的原理是在發(fā)射端發(fā)送經(jīng)過調(diào)制的大時(shí)寬、大帶寬信號，在接收端通過相應(yīng)的匹配處理后，獲得與信號帶寬的倒數(shù)相對應(yīng)的窄脈沖信號，從而解決增強(qiáng)發(fā)射功率和提高距離分辨率之間的矛盾。

假設(shè)匹配濾波器的信號定義為s(t)，其頻譜函數(shù)是S(w)，匹配濾波器的頻域傳輸函數(shù)為輸入信號頻譜的共軛

其時(shí)域的沖激響應(yīng)函數(shù)為

它是輸入信號s(t)的鏡像s(?t)在時(shí)間上平移t0的結(jié)果。要使匹配濾波器物理可實(shí)現(xiàn)，平移時(shí)間必須大于等于信號的終止時(shí)間，s(t)必須在匹配濾波器輸出最大信噪比的時(shí)刻t0前消失(s(t)=0)，因此其最大信噪比輸出時(shí)刻t0必須在輸入信號全部結(jié)束后，物理意義是：若輸入信號沒有結(jié)束，就不能得到信號的全部輸入能量和最大信噪比，一般希望t0盡量小些，把t0取在輸入信號持續(xù)時(shí)間末尾。匹配濾波輸出信號是輸入信號的自相關(guān)函數(shù)，輸出最大信噪比只與輸入信號能量有關(guān)。

1.3 超聲檢測系統(tǒng)傳輸模型

電激勵(lì)信號e(t) 激勵(lì)換能器向被檢樣品發(fā)射超聲波，在傳播中遇到缺陷或界面產(chǎn)生反射回波，被換能器接收，轉(zhuǎn)換成電信號x(t)：

其中：樣品反射的傳輸特性為f(t)；設(shè)h1(t)和h2(t)分別為發(fā)射和接收換能器的沖激響應(yīng)函數(shù)，為簡化分析模型，收、發(fā)換能器的響應(yīng)等效為h(t)；n(t)為檢測隨機(jī)加性噪聲。把脈沖壓縮環(huán)節(jié)引入檢測模型，換能器接收的信號通過激勵(lì)信號e(t) 對應(yīng)的匹配濾波器e(t0?t),則脈沖壓縮結(jié)果為

式中：第二項(xiàng)是噪聲與匹配函數(shù)的互相關(guān)，明顯被抑制；第一項(xiàng)是激勵(lì)信號匹配濾波結(jié)果的自相關(guān)函數(shù)與換能器脈沖響應(yīng)的卷積。時(shí)域表現(xiàn)為壓縮后的脈寬被展寬，引起分辨率下降，頻域表現(xiàn)為激勵(lì)信號通過等效帶通濾波器產(chǎn)生能量衰減，引起脈沖壓縮的信噪比增益損失。

超聲檢測傳輸模型如圖1所示，為描述清楚以及我們實(shí)驗(yàn)的需要，這里把換能器發(fā)射和接收兩種響應(yīng)畫成兩個(gè)換能器表示。

圖1 超聲檢測傳輸模型Fig.1 The model of ultrasonic detection

2 實(shí)驗(yàn)研究

2.1 仿真實(shí)驗(yàn)研究

換能器沖激響應(yīng)模型可用高斯包絡(luò)的正弦信號來描述，e(t)經(jīng)h(t)傳輸后的激勵(lì)聲壓表達(dá)式為

其中：t是傳輸時(shí)間；θ是初始偏移相位；ω0是換能器中心角頻率；σ是與換能器相對帶寬BA有關(guān)的常數(shù)，σ=3.32/BA。

中心頻率為2.25 MHz、帶寬為1.2 MHz的換能器沖激響應(yīng)如圖2所示。換能器的頻譜響應(yīng)相當(dāng)于一個(gè)帶通濾波器，結(jié)果如圖3所示。

選用中心頻率2.25 MHz、帶寬1.2 MHz的線性調(diào)頻信號s(t)，波形如圖4所示。將線性調(diào)頻信號s(t) 通過換能器后，即

圖2 換能器的單位脈沖響應(yīng)Fig.2 The unit impulse response of the transducer

圖3 換能器的單位脈沖響應(yīng)頻譜Fig.3 The spectrum of the unit impulse response

圖4 等幅線性調(diào)頻信號Fig.4 The equal amplitude chirp signal

g(t)波形如圖5所示，用其包絡(luò)y(t)調(diào)制等幅線性調(diào)頻信號s(t)，得r(t)=s(t)×y(t)。為了使s(t)和r(t)激勵(lì)換能器的能量相等，將r(t) 的幅值擴(kuò)大1.737倍，結(jié)果如圖6所示。用r(t)激勵(lì)換能器，即q(t)=r(t)*p(t)，波形如圖7所示。

對等幅線性調(diào)頻波s(t)、包絡(luò)調(diào)制線性調(diào)頻波r(t)及其激勵(lì)換能器后所得波形g(t)和q(t)進(jìn)行脈沖壓縮，其結(jié)果如圖8～11所示，對比這4幅圖中脈沖壓縮后信號的包絡(luò)脈寬，其仿真結(jié)果如圖12所示。因研究對象是幅度的趨勢，對縱坐標(biāo)均做了歸一化處理。

2.2 實(shí)驗(yàn)研究

圖5 等幅線性調(diào)頻信號激勵(lì)換能器后波形Fig.5 The equal amplitude chirp signal after transducer

圖6 包絡(luò)調(diào)制的調(diào)頻信號Fig.6 The envelope modulated chirp signal

圖7 包絡(luò)調(diào)制的調(diào)頻信號通過換能器后波形Fig.7 Envelope modulated chirp signal after transducer

圖8 等幅線性調(diào)頻信號未激勵(lì)換能器脈沖壓縮結(jié)果及包絡(luò)Fig.8 The pulse compression output of the signal shown in Fig.4

圖9 等幅線性調(diào)頻信號激勵(lì)換能器脈沖壓縮結(jié)果及包絡(luò)Fig.9 The pulse compression output of the signal shown in Fig.5

圖10 包絡(luò)調(diào)制線性調(diào)頻信號未激勵(lì)換能器脈沖壓縮結(jié)果及包絡(luò)Fig.10 The pulse compression output of the signal shown in Fig.6

圖11 包絡(luò)調(diào)制線性調(diào)頻信號激勵(lì)換能器脈沖壓縮結(jié)果及包絡(luò)Fig.11 The pulse compression output of the signal shown in Fig.7

實(shí)驗(yàn)裝置采用直徑為1.2 cm、中心頻率為2.25 MHz、?3 dB帶寬為1.2 MHz的超聲換能器發(fā)射信號，為了接收較為寬頻的信號，用在2.25 MHz附近響應(yīng)比較平坦的中心頻率為10 MHz的超聲換能器接收信號，其直徑為0.8 cm。將換能器平行擺放，在水槽中發(fā)射換能器的遠(yuǎn)場區(qū)域放置一根直徑為0.2 mm的銅靶線用于反射回波信號。通過任意波形發(fā)生器(RIGOL DG3000)輸出信號，將其作為換能器的激勵(lì)信號，它的脈沖周期為100 μs。用前置放大器(PANAMET 5800)將靶線反射回波信號進(jìn)行放大，然后由數(shù)字示波器(TELEDYNE LECROY HDO4000)采集信號，采樣頻率為 100 MHz，實(shí)驗(yàn)裝置如圖13所示。

激勵(lì)換能器的波形分別為s(t)和r(t)，如圖14、16所示。由數(shù)字示波器采樣反射波形，如圖15、17所示。

對圖14～17所示的波形進(jìn)行脈沖壓縮，其結(jié)果如圖18～21所示。對比脈沖壓縮后信號的包絡(luò)脈寬仿真結(jié)果如圖22所示。

圖12 脈沖壓縮包絡(luò)脈寬對比仿真結(jié)果Fig.12 The envelope waveforms of pulse compression outputs

圖13 實(shí)驗(yàn)裝置示意圖Fig.13 Schematic diagram of experimental apparatus

圖14 等幅線性調(diào)頻信號Fig.14 The equal amplitude chirp signal

2.3 數(shù)據(jù)分析

圖15 等幅線性調(diào)頻信號激勵(lì)換能器后波形Fig.15 The equal amplitude chirp signal after transducer

圖16 包絡(luò)調(diào)制的調(diào)頻信號Fig.16 The envelope modulated chirp signal

圖17 包絡(luò)調(diào)制的線性調(diào)頻信號通過換能器后波形Fig.17 Envelope modulated chirp signal after transducer

圖18 等幅線性調(diào)頻信號未激勵(lì)換能器脈沖壓縮結(jié)果及包絡(luò)Fig.18 The pulse compression output of the signal shown in Fig.14

圖19 等幅線性調(diào)頻信號激勵(lì)換能器脈沖壓縮結(jié)果及包絡(luò)Fig.19 The pulse compression output of the signal shown in Fig.15

圖20 包絡(luò)調(diào)制線性調(diào)頻信號未激勵(lì)換能器脈沖壓縮結(jié)果及包絡(luò)Fig.20 The pulse compression output of the signal shown in Fig.16

圖21 包絡(luò)調(diào)制線性調(diào)頻信號激勵(lì)換能器脈沖壓縮結(jié)果及包絡(luò)Fig.21 The pulse compression output of the signal shown in Fig.17

從仿真和實(shí)驗(yàn)角度進(jìn)行回波脈沖壓縮性能的比對研究，由圖8～12、18～22得出各信號脈沖壓縮后歸一化主瓣峰值及?3 dB主瓣寬度，如表1所示。

計(jì)算距離旁瓣水平通常有兩種方法，一種叫峰值旁瓣水平(Peak Sidelobe level, PSL)，用來衡量距離旁瓣的最高幅度；另一種叫平均旁瓣水平(Integrated Sidelobe level, ISL)，用來衡量距離旁瓣的總體能量。本文采用PSL來計(jì)算距離旁瓣水平，其定義如表達(dá)式(8)所示：

由圖12、22和表達(dá)式(8)可計(jì)算出各信號脈沖壓縮的峰值旁瓣水平，結(jié)果如表2所示。從表2可以看出，包絡(luò)調(diào)制線性調(diào)頻信號在中心頻段能量較低，低頻段和高頻段能量較高，而等幅線性調(diào)頻信號所有頻段能量分布都相等。激勵(lì)換能器相當(dāng)于乘以換能器單位脈沖頻譜(如圖3所示)，在低頻段和高頻段的能量衰減較大，所以，包絡(luò)調(diào)制線性調(diào)頻信號激勵(lì)換能器的能量損失較大。

表1 換能器對各激勵(lì)信號脈沖壓縮結(jié)果的性能影響分析Table 1 The impacts of the transducer on the pulse compression performances for different excitation signals

表2 換能器對各信號脈沖壓縮結(jié)果旁瓣峰值水平分析Table 2 The impacts of the transducer on the peak sidelobe levels in pulse compression results for different excitation signals

圖22 脈沖壓縮包絡(luò)脈寬對比仿真結(jié)果Fig.22 The envelope waveforms of pulse compression outputs

3 結(jié) 論

(1) 本文用等幅線性調(diào)頻信號和包絡(luò)調(diào)制線性調(diào)頻信號激勵(lì)換能器，分析激勵(lì)前后脈沖壓縮的主瓣峰值和主瓣寬度，探究換能器對脈沖壓縮效果的影響，得出激勵(lì)后主瓣峰值降低，主瓣寬度增加，即激勵(lì)換能器后脈沖壓縮的信噪比和分辨率都有所降低。

(2) 包絡(luò)調(diào)制線性調(diào)頻信號激勵(lì)換能器脈沖壓縮的信噪比和分辨率高于等幅線性調(diào)頻信號。但包絡(luò)調(diào)制線性調(diào)頻信號激勵(lì)換能器后能量損失較大。

(3) 等幅線性調(diào)頻信號和包絡(luò)調(diào)制線性調(diào)頻信號激勵(lì)換能器脈沖壓縮峰值旁瓣水平都有所提高，說明通過換能器后旁瓣被衰減大于主瓣衰減。

參考文獻(xiàn)

[1] Matthew O'Donnell. Coded excitation system for improving the penetration of real-time phased-array imaging systems[J]. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 1992, 39(3): 341-351.

[2] 陸錦輝, 是湘全, 丁慶海, 等. 隨機(jī)二相碼脈沖壓縮雷達(dá)信號分析[J]. 電子學(xué)報(bào), 1996, 24(6): 125-127.LU Jinhui, SHI Xiangquan, Ding Qinghai, et al. Random two phase code pulse compression radar signal analysis[J]. Journal of electronic, 1996, 24(6): 125-127.

[3] Mohammad R. Hoseini. Denoising ultuasonic pulse-echosignal using two-dimendional analytic wavelet thresh olding. Department of Mechanical Engineer[C]// University of Alberta, Canada. April,IEEE, 2012, 45(3): 255-267.

[4] 傅娟, 韋崗, 黃慶華. 基于幅度加權(quán)的預(yù)失真線性調(diào)頻超聲編碼激勵(lì)[J]. 電子與信息學(xué)報(bào), 2013, 35(2): 494-498.FU Juan, WEI gang, HUANG Qinghua. In the weighted amplitude predistortion linear frequency modulated ultrasonic coding incentive[J]. Journal of Electronics & Information Technology,2013, 35(2): 494-498.

[5] 劉波. 基于脈沖壓縮與二相編碼激勵(lì)的超聲檢測提高信噪比研究[J]. 電子學(xué)報(bào), 2011, 39(11): 2701-2707.LIU Bo. Based on pulse compression and two phase encoding incentive Gui ultrasonic testing to improve SNR of inquiry[J]. Acta Electronica Sinica, 39(11): 2701-2707.

[6] 鄭培創(chuàng). 超聲血流檢測與成像技術(shù)研究[J]. 廣州: 華南理工大學(xué),2012, 39(10): 270-275.ZHENG Peichuang. Research on ultrasound blood flow detection and imaging technology[J]. Guangzhou: South China University of Technology, 2012, 39(10): 270-275.

[7] Jennifer E. Michaels. Chirp excitation of ultrasonic guided waves Georgia Institute of Techology[C]. Atlanta, IEEE, 2013, 53(1):265-270.

[8] 魏東, 周正干. 改進(jìn)的非線性調(diào)頻脈沖壓縮方法在空氣耦合超聲檢測中的應(yīng)用[J]. 機(jī)械工程學(xué)報(bào), 2012, 48(16): 8-13.WEI Dong, ZHOU Zhenggan. Improved nonlinear frequency modulan pulse compression method in the application of air coupling ultrsonic testing[J]. Journal of Mechanical Engineering,2012, 48(16): 8-13.

[9] Misaridis T, Jensen JA. Use of modulated excitatinsignals in medical ultrasound. part I: Basic concepts and expected benef i ts[J].IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 2005, 52(2): 177-191.

[10] 顧紅, 蘇衛(wèi)民, 付耀先, 等.一種提高二相編碼雷達(dá)檢測目標(biāo)動態(tài)范圍的方法[J]. 電子學(xué)報(bào),2002, 30(12): 1752-1754.GU Hong, SU Weimin, FU Yaoxian, et al. A method of enhancing the dynamic range of binary phase coded radars[J]. Acta Electronica Sinica, 2002, 30(12): 1752-1754.