IIRFilter初始化用於ColorDopplerwallfilter的研究.docVIP

IIR Filter 初始化 用於Color Doppler wall filter 的研究 王聖維 吳積霖 前言 在彩色Doppler中，我們通常需要一個效能很好的濾波器，用來濾除頻率較低，因組織運動關係而產生的低Doppler 頻率位移。如果我們使用IIR濾波器來濾波，則在初始化時的暫態，將會使得此濾波器的輸出失真。 對於使用於彩色Doppler的wall Filter，我們通常希望此濾波器，能有下列的優良特點： (1)希望我們可以依照被測物的特性來調整自己希望的cutoff frequency。 (2)希望在stopband，能夠將信號有效率的濾除(suppression)，如此一來， 我們方能準確由頻通帶(passband)的信號來估算被測物的速度。 (3)希望暫帶(transition band)附近的響應是非常窄的，如此一來，才能確 保較低速的信號不會隨雜訊被濾掉，保持儀器對低速物體的靈敏度。 (4)此濾波器應該具備低頻譜失真(Minimal spectrum distortion)的特性， 以確保測量結果的正確。 而在信號處理的考量上來看，FIR濾波器原則上較易實現，而且提供了cutoff frequency的可程式化，但是，卻有一個很大的缺點，就是對於相等階數(equivalent order)的濾波器而言，FIR濾波器的暫通帶(transition band （非暫態）)太長。 相反的，IIR濾波器的暫通帶雖然是較FIR還短，也具有cuttoff frequency可程式化的特性，但是其暫態響應(transition response)卻十分長。 此外，FIR在硬體的實現上，需要十分長的長度來implement FIR，尤其當輸入filter的資料有著十分短的資料長度，或者是由一段一段短的資料節段所構成時，為了達到sharp cutoff，我們需要十分長段的資料來初始化濾波器。 相反的，因為IIR濾波器的長度通常較小，因此，很適用於輸入資料並不長的情況，唯一的缺點是IIR濾波器的暫態響應太長，幸好目前已有幾種方法，用來初始化IIR的初始值，使其暫態響應並不會長到我們無法接受。 因此，在此專題中，我們嘗試著以IIR filter的設計方式，來比較幾種不同的IIR filter之初始化，對於其filter cutoff的造成影響。以及這幾種filter的初始化法，對於組織回音信號的排除是否良好，對於血液流速測量的精確性，以及濾波器長度對於效能表現是否有好的效率比，以及穩定度上作一比較。 IIR Filter 初始化的幾種方式 對IIR filter而言，為了希望整個filter能夠及早進入穩態，脫離暫態，我們需要初始化的步驟。因為，IIR filter在方塊圖上是有回饋的，初始化的值將會一保留著他的影響直到filter到達穩態。也因此，慎選IIR filter 初始化的值是十分重要的。 對於IIR，我們有幾種不同初始化方法： (1).ZERO initialization: 這個方法是最直觀的狀況，也就是將輸入IIR的初始值都設為零。在此我們以兩階IIR為例。設一二階濾波器系統方程式表示如下： y[n]=DK[n-1]+b0x[n] where D=[b1 b2 –a1 –a2], 其中，b0,b1,b2,a1,a2為濾波器之係數。 由此，則此濾波器的內部狀態可以用下列向量來表示 K[n]=[ x[n] x[n-1] y[n] y[n-1] ]T 由此，此濾波器的的初始狀態可以下式來表示，即在n=0時，對K[n-1]算值。 K[-1]=[ x[-1] x[-2] y[-1] y[-2] ]T 而我們針對IIR Filter 作初始化的目標，就是要找出一個向量使他的暫態響應最小化。對Zero Initialization而言，即是將此向量全部設為零。 K[-1]=[ 0 0 0 0 ]T 由論文上所言，以zero initialization方式來初始化，所以對於我們要作color Doppler的wall filter而言，似乎不適用;不僅暫通帶很長，而且截止帶的信號的衰減只有8 db左右。並不是十分理想的初始化方法。 (2).STEP initialization Step initialization的作法是假設濾波器輸入的第一個取樣值，早在時間為無窮遠以前便存在，也就是以輸入的第一個取樣值為初始向量。在想像上就好像是一個DC之值來作初