在超音波脈衝回波偵測中,匹配濾波器常用於從受雜訊干擾的接收訊號中偵測已知波形的回波。他的基本概念是將接收端訊號與發射端已知訊號的時間反轉共軛版本進行卷積;在數學上,此運算等價於將接收訊號與已知模板進行互相關分析。若發射訊號為 s(t),則對應的匹配濾波器脈衝響應可表示為:
:h(t)=s^{*}(T-t)
其中 s^{*}(t) 表示訊號的複數共軛,T 為觀測或取樣時間長度。若接收訊號為 r(t),匹配濾波器輸出為:
:y(t)=r(t)\ast h(t)
代入匹配濾波器的脈衝響應後,可得:
:y(t)=\int_{-\infty}^{\infty} r(\tau)s^{*}(T-t+\tau)\,d\tau
此式可視為接收訊號 r(t) 與已知發射波形 s(t) 的互相關運算。因此,當接收訊號中含有與發射訊號相似的回波成分時,匹配濾波器的輸出會在對應的時間位置產生較大的峰值,進而可用於估計回波到達時間與反射介面位置。
相似度量測與啁啾訊號
在醫用超音波與非破壞性檢測中,受測介質常包含散射、衰減與雜訊等因素,使微弱回波不易直接辨識。匹配濾波器可視為一種相似度量測方法,藉由比較接收訊號與已知發射波形的相似程度,將與目標波形相符的回波成分加以強化,而與目標波形不相關的雜訊則較不容易在相關輸出中形成明顯峰值。
若發射端使用線性調頻啁啾訊號,其訊號可表示為:
:s(t)=A\cos\left[2\pi\left(f_0t+\frac{k}{2}t^2\right)\right],\quad 0\leq t\leq T
其中 A 為振幅,f_0 為起始頻率,T 為訊號持續時間,k 為調頻率。若啁啾訊號的頻率由 f_0 掃描至 f_1,則:
:k=\frac{f_1-f_0}{T}
此類訊號的瞬時頻率會隨時間改變,可表示為:
:f(t)=f_0+kt
由於啁啾訊號能在較長時間內發射能量,常與匹配濾波器配合,用於超音波編碼激發與脈衝壓縮。
脈衝壓縮與回波定位
啁啾訊號雖然具有較長的持續時間,可增加發射能量與穿透能力,但直接使用長脈衝,會使軸向解析度下降。配合匹配濾波器後,接收端可將長時間的啁啾回波壓縮為較窄的峰值。此過程可用自相關函數描述:
:R_{ss}(\tau)=\int_{-\infty}^{\infty}s(t)s^{*}(t-\tau)\,dt
當 \tau=0 或回波與模板訊號對齊時,自相關函數會產生最大值。因此,匹配濾波器輸出的主瓣位置可用來判斷回波到達時間。若超音波在介質中的聲速為 c,回波往返時間為 \Delta t,則反射體深度約為:
:d=\frac{c\Delta t}{2}
式中的 2 代表超音波訊號需從探頭傳至反射體,再由反射體返回探頭。因此,匹配濾波器輸出峰值的位置可直接對應至反射體或組織界面的深度。
訊雜比提升機制
在加性白雜訊條件下,若已知欲偵測訊號的波形,匹配濾波器是能使特定取樣時刻輸出訊雜比最大化的線性濾波器。若接收訊號表示為:
:r(t)=\alpha s(t-t_0)+n(t)
其中 \alpha 為回波振幅,t_0 為回波延遲時間,n(t) 為加性雜訊。當匹配濾波器與回波訊號對齊時,訊號成分會相干累加,而隨機雜訊則不會以相同方式累加。對於功率譜密度為 N_0/2 的白雜訊,匹配濾波器在最佳取樣時刻的最大輸出訊雜比可表示為:
:\mathrm{SNR}_{\max}=\frac{2E_s}{N_0}
其中 E_s 為發射訊號能量:
:E_s=\int_{-\infty}^{\infty}|s(t)|^2\,dt
此式表示,在白雜訊環境中,匹配濾波器的偵測能力與訊號能量成正比。因此,使用較長的啁啾訊號可增加總能量,再透過匹配濾波器將能量壓縮集中於短時間區間,進而提升微弱回波的可偵測性。
在頻域中,匹配濾波器的轉移函數可表示為:
:H(f)=S^{*}(f)e^{-j2\pi fT}
其中 S(f) 為發射訊號 s(t) 的傅立葉轉換。此式表示匹配濾波器會根據已知發射訊號的頻譜進行加權與相位校正,使與發射波形相符的回波成分在輸出端形成較大的相關峰值。
超音波系統響應與波形設計
實際超音波系統中的探頭與前端電路並非理想元件,通常具有有限頻寬與帶通特性。因此,若發射訊號的頻譜與探頭頻率響應不匹配,可能會影響脈衝壓縮後的主瓣寬度與旁瓣大小。旁瓣過高時,可能遮蔽鄰近的微弱反射訊號,進而降低影像對比或回波判讀能力。
為了改善此問題,部分研究會設計非線性頻率調變訊號或使用加權、失配濾波器等方法,調整啁啾訊號的頻譜分布與脈衝壓縮特性。這類方法的目的通常是降低旁瓣、改善軸向解析度,或在特定系統頻寬限制下提高回波偵測品質。不過,旁瓣抑制、主瓣寬度與訊雜比之間通常存在取捨,因此實際設計需根據探頭頻寬、組織衰減、影像解析度與偵測目標進行調整。
參考資料
參見
*匹配濾波器
*啁啾
*醫學超音波檢查
评论 (0)