Q&A
まだ上記の関数を試したりはできていません。
使用する前に、各引数の意味や、使い方を知りたいと思い質問しました。
それぞれの私の認識があっているかなどご指摘いただきたいと思っています。よろしくお願いします。
実現したいこと
Pythonのライブラリである”buttord(wp, ws, gpass, gstop)”の
各引数の、意味や定義方法について教えてください。
N, Wn = signal.buttord(wp, ws, gpass, gstop)
現状の認識は以下になっています。それぞれ認識があっているかご教示ください。
・wp, ws：カットしたい周波数の開始と終了をナイキスト周波数で割ることで正規化
・gpass：カットオフ周波数より、大きさを1/2にしたい場合は-3dBとなるため、
「3」に設定することがメジャー
・gstop：全く分かっていません
・N：フィルタの次数
・Wn：臨界周波数 (※この意味もよく分かっていません)
[使用環境]
・Python3.8
・windos10 64bit
[参考]
・「PythonのSciPyでローパスフィルタをかける！」https://watlab-blog.com/2019/04/30/scipy-lowpass/
ここは質問に対するコメントを書くところです。
質問に追加したいことがあるのであれば、質問を編集して追記したほうがいいですよ。
承知しました。ご指摘ありがとうございます。
回答1件
0
ベストアンサー
僕も詳しいわけではありませんし、質問者さんがどれくらい知識があるのかもわからないので、ざっくりとした説明です。
バターワースフィルタを定義するには、 次数 (order)と 減衰開始のキーポイント(3db減数するポイント)が必要です。
また、フィルタが欲しい場合、xxHzまで通して yyHz以上は通したくないという要求があるはずです。この情報から、バターワースフィルタのパラメータを算出してくれるのが、buttord(wp, ws, gpass, gstop) 関数です。
- wp : 減衰開始周波数 (ナイキスト周波数で正規化したもの)
- gpass : 減衰開始とする減衰率
3db下ったところから減衰開始とすることが多いので、ここは通常3dbとします。 - wp : この周波数意向は通したくないの周波数 (ナイキスト周波数で正規化したもの)
- gstop: 通さない周波数はこの減衰量(db)以下にしたいという値
バターワースフィルタの特性は 以下のページの説明を見るとわかるでしょう。
https://cc.cqpub.co.jp/system/contents/2067/
一番最初の図だと、 wp=1kHz, ws=10kHz, gpass=3db, gstop=70db です。
これで計算して出てきた値は、
- N: 次数
- Wn: 3db下ったところの周波数。 gpassに3を与えていれば、wpと同じ値
で、これを scipy.signal.butter()関数に与えると、所望のフィルタが得られます。
投稿2023/02/22 07:01
総合スコア14249
あなたの回答
tips
プレビュー
