濾波器有四種表示形式:差分方程、脈衝函式、傳輸函式和頻率響應。其中前兩個是定義在實數域中,後兩個定義在複數域中。
首先看一下離散傅利葉變換的定義。
我們由濾波器的差分形式退出濾波器的頻率響應形式。假設有濾波器的差分方程為: σn
k=0a
ky[n
−k]=
σmk=
0bkx
[n−k
] 對每一項進行離散傅利葉變換,得到: a0
y(ω)
+a1e
−jωy
(ω)+
a2e−
j2ωy
(ω)+
⋯+an
e−jn
ωy(ω
)=b0
x(ω)
+b1e
−jωx
(ω)+
b2e−
j2ωx
(ω)+
⋯+bm
e−jm
x(ω)
得到: h(
ω)=y
(ω)x
(ω)=
b0+b
1e−j
ω+b2
e−j2
ω+⋯+
bme−
jma0
+a1e
−jω+
a2e−
j2ω+
⋯+an
e−jn
ω 稱h
(ω) 為濾波器的頻率響應。h(
ω)是複數,可以用極座標形式表示h(
ω)=|
h(ω)
|ejθ
(ω) ,
|h(ω
) 表示濾波器在數字頻率
ω 處的增益,θ(
ω)是相位差。
正弦輸入:
現在把余弦(或正弦)訊號(頻域上的)寫成簡寫形式,方便在頻域中計算。比如輸入x[
n]=a
cos(
nω0+
θx) 的x(
ω)可以記作a|
θx。假設在頻率ω0
下,頻率響應的增益為
h ,相位差為
θ,那麼頻率響應可以簡寫成h|
θ ,那麼dft的相乘現在很容易計算: y(
ω)=h
(ω)x
(ω)=
(h|θ
)(a|
θx)=
ha|θ
+θx
而且很容易再轉換到實數域上: y[
n]=h
acos(n
ω0+θ
+θx)
幅度響應和相位響應
幅度響應是增益|h
(ω)|
和ω之間的關係圖,相位響應是相位差θ(
ω)和|om
ega 之間的關係圖。頻率響應可以用幅度響應和相位響應表示。
注意,所有的幅度響應都是週期為2π
的偶函式;所有的相位響應都是週期為2π
的奇函式。所以兩個響應圖的自變數範圍只需要在[0
,π] 上就可以表示所有情況。
濾波器的形狀由幅度響應給出,幅度響應說明了濾波器的型別,如低通、高通、帶通、帶阻。
但是給定濾波器的頻率響應h(
ω),我們是沒有辦法畫出它的幅度響應曲線的。那怎麼判斷濾波器的型別呢?可以通過極零點位置確定濾波器的形狀。假設濾波器的傳輸函式為: h(
z)=k
(z−z
0)(z
−p0)
(z−p
1)其中
k 是增益,zj
是零點,pj
是極點,那麼濾波器的形狀為: |h
(ω)|
=k|e
jω−z
0||e
jω−p
0||e
jω−p
1|=k
×(半徑
為1的上
半圓上ω
對應的點
到z0的
距離)(
半徑為1
的上半圓
上ω對應
的點到p
0的距離
)(半徑
為1的上
半圓上ω
對應的點
到p1的
距離)
當極點和零點比較靠近單位圓時,會導致幅值大小的劇烈變化,標明濾波器的選擇性比較好。
matlab 數字訊號處理 FIR濾波器
matlab 數字訊號處理 fir濾波器 當我們處理一段訊號時,如圖1所示 往往需要通過消除雜訊,放大有效訊號的部分的處理,此時可以通過時域和頻域的方式來濾除雜訊。fir與iir的區別在於,fir具有嚴格的線性相位特性,但運算時間過長。而iir濾波器則可以大大的計算執行時間。頻域上消除雜訊成分則是通...
數字訊號處理 數字濾波器
輸入訊號中有用的頻率成分和洗完濾除的成分個占有不同的頻帶,通過濾波器選頻實現濾波目的 如高通濾波器,帶通濾波器,低通濾波器等等 訊號和干擾的頻譜會相互重疊,需要根據隨機訊號的統計特性,在某種準則下最大限度地抑制干擾,恢復訊號,達到濾波目的.寫出系統幅頻特性函式 h e h z big 令 omega...
論高通濾波器,帶阻濾波器,陷波濾波器
首先,對一副影象進行如下二維傅利葉變換。我們將u 0和v 0帶上式,我們可以得到如下式子。根據上式,可以到f 0,0 的值是非常大的。這裡,我們將 f 0,0 稱為直流分量,直流分量比其他的成分要大好幾個數量級。所以,這也就是傅利葉譜為什麼需要使用對數變換才能看清楚的原因。這裡,對於高通濾波器而言,...