Overview of Digital Filters (INCT)
1.11
2005/10/17 07:12:22 GMT-5
2008/01/14 03:17:37.390 US/Central
Yoji
Yamada
yyama@ishikawa-nct.ac.jp
Hironori
Takaryo
s043105@ishikawa-nct.ac.jp
Yoji
Yamada
yyama@ishikawa-nct.ac.jp
Hironori
Takaryo
s043105@ishikawa-nct.ac.jp
線形シフト不変システムとインパルス応答
ディジタル信号処理を実行する処理ブロックのモデルとして、図1に示すブロック図を考える。
システムの入出力関係 図1で、
x ( n )
および
y ( n )
をそれぞれ入力信号および出力信号とする。
ここで、
x ( n )
から
y ( n )
を生じる処理ブロックを、システムと呼ぶことにする。
システムが線形性、シフト不変性[1]などの条件を満たすとき、このシステムを
線形シフト不変システムという。
また、システムを実時間処理で実現するには、
因果性[1]が要求される。本テーマで扱うシステムはすべて因果性を満たす
ものとする。 線形シフト不変システムにディジタル単位インパルス
δ ( n )
δ ( n ) = {
1
0
( n = 0 )
( n ≠ 0 )
を入力したときの出力
h ( n )
を、インパルス応答と定義する。
インパルス応答は、システムの性質を示す重要な値である。
畳み込み演算とディジタルフィルタ
あるシステムが線形性とシフト不変性を満たす場合、そのシステムのインパル
ス応答
h ( n )
が与えられれば、任意の入力
x ( n )
に対する出力
y ( n )
は、
式(2)により計算することができる。
y ( n ) =
∑
k = − ∞
∞
x ( k ) h ( n − k )
=
∑
k = − ∞
∞
h ( k ) x ( n − k )
さらに、ある入力信号に対して畳み込み演算を実行することにより、入力信号
からある特定の周波数成分だけを抜き出すなどのフィルタ効果があることが知
られている。加えてこのフィルタ効果は、インパルス応答を変化させることに
より様々に変化させることができる。
このことから、畳み込み演算を実行するシステムは、ディジタルフィルタ
ともよばれている。インパルス応答が有限時間だけ持続するシステム
をFIR(Finite Impulse Response)システム(またはFIRフィルタ)、
無限時間持続するフィルタのことをIIR(Infinite Impulse Response)システム
(またはIIRフィルタ)と分類している。本テーマでは、FIRフィルタ
を扱うこととする。
インパルス応答と周波数特性
以上のようなフィルタ効果を与えるシステムの特性を、そのシステムの周波数
特性として表現することができる。インパルス応答
h ( n )
を有するシステム
の周波数特性
H (
e
j ω
)
は、
ω
を正規化角周波数として、
次式で定義される。
H (
e
j ω
) =
∑
n = − ∞
∞
h ( n )
e
− j ω n
システムがFIRフィルタの場合、インパルス応答
h ( n )
は有限時間しか持
続しない。また、因果性を満たすシステムの場合、因果性の定義より、
h ( n ) = 0 ( n < 0 )
であるので、因果性を満足しインパルス応答長が
N
で
あるFIRシステムの周波数特性は次式となる。
H (
e
j ω
) =
∑
n = 0
N − 1
h ( n )
e
− j ω n
N
のことをフィルタ長ということもある。
周波数特性の計算
周波数特性には、その大きさを表す振幅特性と、入出力信号間の位相の変化を
表す位相特性とがある。また、位相特性を別の形で表現した特性として、
群遅延特性がある。
(4)式より、振幅特性
H (
e
j ω
)
は、以下のように実部と虚部に分けて表すことができる。
H (
e
j ω
) =
∑
n = 0
N − 1
h ( n ) cos ω n + j {
−
∑
n = 1
N − 1
h ( n ) sin ω n
}
ここで、右辺第一項は実部、右辺第二項は虚部である。よって振幅特性
H (
e
j ω
)
の絶対値は、(5)式の絶対値をとって以下のように計算できる。
|
H (
e
j ω
)
| =
R
e
2
+
I
m
2
式(6)で
R
e
および
I
m
は、それぞれ(5)式における実部と虚部を表すものとする。振幅特性の大きさは、しばしば[dB]で表示される。この場合の振幅(ampitude) [dB]は、常用対数をとって次式で表される。
20
log
10
|
H (
e
j ω
)
|
(7)式で、正規化角周波数
ω
を細かく変化させて
ω
を横軸、振幅特性を縦軸としてグラフに示せば、振幅特性の
グラフを得ることができる。また、周波数を正規化周波数
f =
ω
2 π
で表現することも多い。正規化周波数
f
は、周波数
F
[Hz]をサンプリング周波数
F
s
[Hz]で割った値である。 位相特性
θ ( ω )
は、振幅特性
H (
e
j ω
)
を極座標表示したときの偏角である。
数値計算により
θ ( ω )
の値を求める場合には、次式を用いることができる（ただし、その処理系での逆正接関数 －たとえば atan( )関数－ の主値の範囲について理解しておくことが必要であろう）。
θ ( ω ) =
tan
− 1
I
m
R
e
[ r a d ]
さらに、群遅延特性
τ
の定義は、位相特性の傾きとして次式で与えられる。
τ = −
d θ ( ω )
d ω
[ S a m p l e ]
以上をまとめると、インパルス応答
h ( n )
の値が既知ならば、周波数特性(振幅、
位相、群遅延特性)を計算できることになる。あるインパルス応答について、実際に周波数特性を計算した例
を図2に示す。図2(a)がインパルス応答で、図2(b)から(d)がそれぞれ振幅、
位相、群遅延特性である。 周波数特性の計算例
伝達関数についての補足
式(2)をz変換すると、次式となる。
Y ( z ) = H ( z ) X ( z )
上式の
H ( z )
を伝達関数と呼ぶ。伝達関数
H ( z )
はインパルス応答
h ( n )
の
z
変換
であり、
ディジタルフィルタの入出力関係を
z
領域で表した式である。式(3)と式(11)
を比較すると、畳み込み演算が
z
領域では積で表現できることがわかる。 ここで、上式に
z =
e
j ω
を代入すると次式となる。
Y (
e
j ω
) = H (
e
j ω
) X (
e
j ω
)
X (
e
j ω
)
および
Y (
e
j ω
)
は、それぞれ入力信号
x ( n )
および出力信号
y ( n )
の離散時間フーリエ変換(周波数スペクトル)である。式(11)は、畳み込み演算
を周波数領域で表現した式に相当する。