- ベストアンサー
バンドパスIIRフィルタの低域におけるエラー
hot-ICの回答
- hot-IC
- ベストアンサー率87% (7/8)
次数が小さいからではないでしょうか? Labviewでは次数はいくつになっていますか? あと、Matlab centralなどを見ると、Matlab用のオクターブ分析のコードがいくつか あります。 http://www.mathworks.com/matlabcentral/fileexchange/69-octave 試してみたところOctaveでも一部動作するようです。 また、作者の方の論文も参考になりそうです。 http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.57.5728 JISは国際規格と互換性があるはずですが、細かいところはご自身で確認されたほうが良いと思います。
関連するQ&A
- IIR型の低域通過フィルタについて
次のような問題を出されたのですが、解き方がさっぱり分かりません。 どなたか解答例をお願いします。 遮断周波数がπ/3となるIIR型の低域通過フィルタを設計したい。 (a)アナログバタワースフィルタから、双1次z変換によってデジタルフィルタの設計をしたい。このとき、アナログフィルタにおける遮断周波数はどうなるか。 (b)上記の遮断周波数で設計した3次バタワースフィルタの伝達関数は添付画像の式のようになる。対応するデジタルフィルタの伝達関数H(z)を示せ。 (c)|H(z)|の周波数特性を求めよ。
- 締切済み
- 科学
- IIR型の低域通過フィルタについて
次のような問題を出されたのですが、中間試験に備えて解き方をマスターしたいと思っています。 どなたか解答例をお願いします。 遮断周波数がπ/3となるIIR型の低域通過フィルタを設計したい。 (a)アナログバタワースフィルタから、双1次z変換によってデジタルフィルタの設計をしたい。このとき、アナログフィルタにおける遮断周波数はどうなるか。 (b)上記の遮断周波数で設計した3次バタワースフィルタの伝達関数は添付画像の式のようになる。対応するデジタルフィルタの伝達関数H(z)を示せ。 (c)|H(z)|の周波数特性を求めよ。 どなたか、回答よろしくお願いします。
- 締切済み
- 科学
- IIRフィルタの応答を高速にする方法
入力信号をIIRフィルタに通すと、一般に、出力が落ち着くまで時間がかかります。 この時間を短縮する方法はありますでしょうか。 例えば、一定振幅 A Hz + B HzのSin波を、中心周波数がA Hzで帯域幅が狭く次数の高いIIR BPFに通すと、出力信号の振幅は、しばらく0が続き、途中からグワっと大きくなり、やがて一定のA Hzに落ち着きます。 帯域幅や次数を変えないで、信号を入れた直後から、振幅一定のA Hzの出力を得たい(応答を早める)のですが、可能でしょうか。 アドバイス、お願いいたします。
- ベストアンサー
- 物理学
- FIRフィルタの遅延量補正とIIRフィルタの安定性について教えて頂きた
FIRフィルタの遅延量補正とIIRフィルタの安定性について教えて頂きたいのですが。 双一次z変換を用いていくつかのIIRフィルタを作る事が出来ました。Scilbが持つIIR関数と比較して同じ結果になっています。IIRフィルタがフィルタ係数によっては不安定=発信したりする、というのはScilab等の結果と同じになる→フィルタ係数も同じ→不安定にならない。と思ってよいのでしょうか?。例えばですが、あるサンプリング周波数の波形データをフィルタする為に、同じサンプリング周波数で作ったImplus波形をIIRフィルタに放り込んで周波数応答を確認したら、実際にフィルタする波形データの最初から最後までその周波数応答でフィルタされると私は理解しているのですが。 私のフィルタを使う用途の場合は、フィルタの計算速度は特に制限は無く、速いに越した事はありませんが得られた結果が正確な事が大切なので、安定性や係数誤差の事を考えるとIIRフィルタではなくFIRフィルタでも構わないのですが、FIRフィルタの場合だとタップ数が多くなるので波形の遅延量が無視出来ません。フィルタを通した波形は次定数fastのレベル波形にして100Hz位で出すとは言え、IIRよりは数秒かコンマ数秒遅延した結果になると思います。そのような場合、遅延量を補正するとしたらどうしたら良いのでしょうか?。個人的にはタップ数が(IIRの場合は次数が)遅延量と思っているのですが。 宜しくお願い致します。
- ベストアンサー
- フリーウェア・フリーソフト
- IIRフィルタについて
IIRフィルタで疑問があり、過去のページを調べてもわからなかった ので新たに質問を作りました。 今検討しているのは下記HPの図1のIIRフィルタの1次版(a1,b1,b0のみ) で伝達関数(Z領域)、サンプリング周波数、各係数が分かっています。 (逆に言うとこれしか分かってません) ここで上記の資産を流用して新たにサンプリング周波数 だけを変更して使いまわしたいのですが係数をどのように 変更したらよいのか分からずに困っています。 (アナログの伝達特性などはわかりません) どなたか分かる方、教えていただけないでしょうか。 よろしくお願いします。 http://adsp2191.hp.infoseek.co.jp/misc/009_difficult_sp.shtml もしくは http://www.ics.es.yamanashi.ac.jp/~hanawa/asp/ppt/asp-6.ppt(39枚目スライドの右図)
- ベストアンサー
- 物理学
- デジタルフィルタ(IIR)
無限インパルス応答(IIRフィルタ)について、サルでもわかるようにメカニズム・カットオフ周波数について教えて下さい。又は、サルでもわかるように解説された書籍・HPなど御存知でしたら教えて頂けると幸いです。
- 締切済み
- その他(学問・教育)
- ナイキスト周波数に関して
ナイキスト周波数に関する質問です。 本に 「ナイキスト周波数は、最少標本化周波数である。」 ネットに 「ナイキスト周波数は、サンプリング周波数の1/2である。」 と書かれていました。 どちらも、成り立つのであるなら、 ナイキスト周波数が、原信号の帯域幅の2倍の周波数であるので、 標本化をする際、サンプリング周波数は帯域幅の4倍からしか取れなという事になるのですか? そんなはずは、ないと思うのですがどのように考えたら良いのでしょうか? よろしくお願いします。
- ベストアンサー
- 情報工学
お礼
回答ありがとうございます。 labviewでの次数は同じく4にしています。Scilab、Octave共に次数を上げるとよりナイキスト周波数に近い帯域もエラーになってしまいます。 ご紹介頂いたファイルの中身を見たら、低いOct帯域でdecimateしている記述があるので、私と同じ問題でそうしているのか、計算時間の短縮の為なのかは判りませんが、私もdecimateしてみます。ちなみにMatlabのFilter Design Toolboxには中心周波数と次数だけで結果が出てくる関数があるようなので、Labviewもそうですが流石に高いだけあると思います。しかしオクターブバンドフィルタは基本的にJISで定める計算式で中心周波数を出し、同じく計算式で上限下限遮断周波数を指定したバンドパスフィルタなので、苦労するとは思いませんでした。 どうもありがとうございました。