cheb1ord (Signal Processing Toolbox)

Chebyshev I型フィルタの次数選択

表示

[n,Wn] = cheb1ord(Wp,Ws,Rp,Rs)
[n,Wn] = cheb1ord(Wp,Ws,Rp,Rs,'s')

詳細

cheb1ordは、フィルタ設計仕様を満たすディジタルまたはアナログ Chebyshev I 型フィルタの最小次数を計算します。

ディジタル領域

[n,Wn] = cheb1ord(Wp,Ws,Rp,Rs) は、通過帯域でRp dB 以下の減衰、遮断帯域で Rs dB 以上の減衰をもつ Chebyshev 1型フィルタの最小次数 n を出力します。また、対応するカットオフ周波数ベクトル Wn も出力します。これらの出力引数 nとWnは、関数cheby1 で使用します。

つぎの表に従って、入力引数を選択して、遮断帯域と通過帯域を設定してください。

表 7-3: 遮断帯域と通過帯域のフィルタパラメータの記述
Wp
通過帯域コーナ周波数 Wp、すなわち、カットオフ周波数は、0と1の間の値で、スカラまたは2要素から構成されています。ここで、1は、正規化されたNyquist周波数、ラジアン/サンプルに対応します。

Ws
遮断コーナ周波数 Wsは、0と1の間の値で、スカラまたは2要素から構成されています。ここで、1は、正規化されたNyquist周波数に対応します。

Rp
dB 単位で表した通過帯域での許容される最大変化量。

Rs
dB 単位で表した遮断帯域で必要な減衰量の最大量。

**表 7-3: 遮断帯域と通過帯域のフィルタパラメータの記述**
`Wp`	通過帯域コーナ周波数 `Wp`、すなわち、カットオフ周波数は、0と1の間の値で、スカラまたは2要素から構成されています。ここで、1は、正規化されたNyquist周波数、ラジアン/サンプルに対応します。
`Ws`	遮断コーナ周波数 `Ws`は、0と1の間の値で、スカラまたは2要素から構成されています。ここで、1は、正規化されたNyquist周波数に対応します。
`Rp`	dB 単位で表した通過帯域での許容される最大変化量。
`Rs`	dB 単位で表した遮断帯域で必要な減衰量の最大量。

つぎの表を使って、種々のタイプのフィルタを設定します。

表 7-4: フィルタタイプに依る遮断帯域と通過帯域の設定
フィルタタイプ
遮断帯域と通過帯域の条件
遮断帯域
通過帯域

ローパス
Wp < Ws、共にスカラ
(Ws,1)
(0,Wp)

ハイパス
Wp > Ws、共にスカラ
(0,Ws)
(Wp,1)

バンドパス
Wsで指定される区間は、Wp (Ws(1) < Wp(1) < Wp(2) < Ws(2))で指定されるものです。
(0,Ws(1)) と (Ws(2),1)
(Wp(1),Wp(2))

バンドストップ
Wpで指定される区間は、Ws (Wp(1) < Ws(1) < Ws(2) < Wp(2))で指定されるものです。
(0,Wp(1)) と (Wp(2),1)
(Ws(1),Ws(2))

**表 7-4: フィルタタイプに依る遮断帯域と通過帯域の設定**
フィルタタイプ	遮断帯域と通過帯域の条件	遮断帯域	通過帯域
ローパス	`Wp` < `Ws`、共にスカラ	`(Ws,1)`	`(0,Wp)`
ハイパス	`Wp` > `Ws`、共にスカラ	`(0,Ws)`	`(Wp,1)`
バンドパス	`Ws`で指定される区間は、`Wp` (`Ws(1) < Wp(1) < Wp(2) < Ws(2)`)で指定されるものです。	`(0,Ws(1))` と `(Ws(2),1)`	`(Wp(1),Wp(2))`
バンドストップ	`Wp`で指定される区間は、`Ws` (`Wp(1) < Ws(1) < Ws(2) < Wp(2)`)で指定されるものです。	`(0,Wp(1))` と `(Wp(2),1)`	`(Ws(1),Ws(2))`

ユーザのフィルタ仕様に通過帯域でのリップルの許容変化量と遮断帯域でのリップルの許容変化量がお互いに等しくない場合、表の中の設定に従ってローパスフィルタとハイパスフィルタを別々に設計し、2つのフィルタをカスケードに結合します。

アナログ領域

[n,Wn] = cheb1ord(Wp,Ws,Rp,Rs,'s') は、カットオフ周波数 Wn を満足するアナログ Chebyshev I型フィルタの最小次数 n を出力します。表 7-3 "遮断帯域と通過帯域のフィルタパラメータの記述"と同様に周波数 Wp と Ws を設定します。この場合、周波数はラジアン/秒でのみ、通過帯域または遮断帯域を無限大に設定することができます。

関数cheb1ordを使って、表 7-4 "フィルタタイプによる遮断帯域と通過帯域の仕様設定"に記述されているローパス、ハイパス、バンドパス、バンドストップフィルタを設定します。

例題

データサンプル周波数 1000 Hz で、0から40 Hz でリップルの最大変量を3 dB 以内の通過帯域とし、150 Hz からNyquist 周波数(500 Hz)で、最低でも60 dB 減衰するローパスフィルタを設計しましょう。

Wp = 40/500; Ws = 150/500;
Rp = 3; Rs = 60;
[n,Wn] = cheb1ord(Wp,Ws,Rp,Rs)
n = 
     4
Wn =
    0.0800
[b,a] = cheby1(n,Rp,Wn);
freqz(b,a,512,1000); 
title('n=4 Chebyshev Type I Lowpass Filter')

つぎに、60 Hz から 200 Hz でリップルの最大変量を3 dB 以内の通過帯域とし、両側 50 Hz の遷移部をもち、遮断帯域では、40 dB 以上の減衰をもつバンドパスフィルタを設計しましょう。

Wp = [60 200]/500; Ws = [50 250]/500;
Rp = 3; Rs = 40;
[n,Wn] = cheb1ord(Wp,Ws,Rp,Rs)
n =
     7
Wn =
    0.1200    0.4000
[b,a] = cheby1(n,Rp,Wn);
freqz(b,a,512,1000);
title('n=7 Chebyshev Type I Bandpass Filter')

アルゴリズム

cheb1ordは、参考文献[1]に記述されているChebyshevローパスフィルタの次数予測公式を使用します。この関数は、アナログとディジタルのいずれの場合も、アナログ領域で処理を行います。ディジタルの場合には、周波数パラメータをs領域に変換して、次数と固有周波数を推定し、その後でz領域に逆変換します。

cheb1ordは、まず、希望するフィルタの通過帯域周波数を1ラジアン/秒(ローパスフィルタおよびハイパスフィルタの場合)、または、-1および1ラジアン/秒(バンドパスフィルタおよびバンドストップフィルタの場合)に変換することにより、ローパスフィルタのプロトタイプを作成します。そして、遮断帯域仕様を満たすローパスフィルタに必要な最小次数を計算します。

参考

buttord
Butterworthフィルタの次数選択

cheby1
Chebyshev I型フィルタの設計(通過帯域リップル)

cheb2ord
Chebyshev II型フィルタの次数選択

ellipord
楕円フィルタの次数選択

kaiserord
KaiserウインドウをもつFIRフィルタ設計に対してのパラメータを推定

`buttord`	Butterworthフィルタの次数選択
`cheby1`	Chebyshev I型フィルタの設計(通過帯域リップル)
`cheb2ord`	Chebyshev II型フィルタの次数選択
`ellipord`	楕円フィルタの次数選択
`kaiserord`	KaiserウインドウをもつFIRフィルタ設計に対してのパラメータを推定

参考文献

[1] Rabiner, L.R., and B. Gold. Theory and Application of Digital Signal Processing. Englewood Cliffs, NJ: Prentice-Hall, 1975. Pg. 241.

cheb1ap cheb2ap