因果濾波器

信号处理裡的因果濾波器是线性时不变的因果系统。「因果」(causal)一詞是指濾波器輸出只和過去和現在的輸入有關。若濾波器的輸出也和未來的輸入有關,則為非因果(non-causal),若濾波器的輸出只和未來輸入有關,則為反因果(anti-causal)。可實現(可以實時運算)的系統和濾波器必須是因果因統,因為實時運算的系統無法靠未來的輸入來產生目前的輸出。以效果上,最可以代表時間t輸入的輪出,會在出現稍晚於時間t的輸出。数字滤波器裡的常見設計實務就是利用將非因果的脈衝響應變慢或是時間移位,以創建可實現的濾波器。若需要縮短脈衝響應的話,一般會將脈衝響應乘以窗函数來達成。

像最大相位濾波器就是反因果濾波器的一個例子,本身可以是有界輸入有界輸出穩定性,而反因果濾波器的倒數也是穩定的非因果濾波器。

例子
以下定義是輸入信號s(x),的移動平均,此處省略係數:

:f(x) = \int_{x-1}^{x+1} s(\tau)\, d\tau\ = \int_{-1}^{+1} s(x + \tau) \,d\tau\,

其中x可以表示空間坐標,就像在影像處理中的應用。但若x表示時間(t)\,,則這樣定義的移動平均是非因果的,因為f(t)和未來的輸入有關,就像s(t+1)\,。可實現的版本如下:

:f(t-1) = \int_{-2}^{0} s(t + \tau)\, d\tau = \int_{0}^{+2} s(t - \tau) \, d\tau\,

是上述濾波器延時後的版本。

任何一個線性濾波器(就像此處的移動平均)其特徵可以用其冲激响应h(t)表示,而其輸出為以下的卷积

:
f(t) = (h*s)(t) = \int_{-\infty}^{\infty} h(\tau) s(t - \tau)\, d\tau. \,

若用這種表示方式,因果代表

:
f(t) = \int_{0}^{\infty} h(\tau) s(t - \tau)\, d\tau

若上述的等式要成立,一般來說會需要h(t) = 0,針對所有t 
g(t) = {h(t) + h^{*}(-t) \over 2}

是非因果的。換句話說,g(t)是埃爾米特函數,因此傅立葉變換G(ω)是實值函數。可以得到以下的關係

:
h(t) = 2\, \Theta(t) \cdot g(t)\,

其中Θ(t)是单位阶跃函数。

因此h(t)和g(t)的傅立葉變換有以下的關係:
:
H(\omega) = \left(\delta(\omega) - {i \over \pi \omega}\right) * G(\omega) =
G(\omega) - i\cdot \widehat G(\omega) \,

其中\widehat G(\omega)\,是頻域下的希爾伯特轉換(不是時域的)。\widehat G(\omega)\,的正負號會依傅立葉變換的定義而不同。

將上式進行希爾伯特轉換,可得到以下H和其希爾伯特轉換的關係:
:
\widehat H(\omega) = i H(\omega)

參考資料
*
*

评论 (0)

  • 还没有评论,来抢沙发吧。