プチコン 非公式コミュニティ
トピック
こういち
◆.Id/aHiU36hu
2022/6/27 12:51
2022/6/27 12:51
質問
全極型IIRフィルタの周波数特性
ここ数ヶ月ほど悩んでる問題ですが、
リンク先の講義資料1の12ページあたりにあるパワースペクトルの求め方が分からないです。
パッと見フーリエ変換だなと思って雑にFFTを使ってみましたが、実行の度に出力が大きく変わってしまいました。
https://ocw.u-tokyo.ac.jp/lecture_231/
リンク先の講義資料1の12ページあたりにあるパワースペクトルの求め方が分からないです。
パッと見フーリエ変換だなと思って雑にFFTを使ってみましたが、実行の度に出力が大きく変わってしまいました。
https://ocw.u-tokyo.ac.jp/lecture_231/
コメント
こういち
2022/6/27 19:56
◆.Id/aHiU36hu
http://digitalfilter.com/smartdsp/dai2bu/chapiir/biquad02.html
http://www.ic.is.tohoku.ac.jp/~swk/lecture/yaruodsp/dfanalysis.html
どうやら、これを見る限りωは正規化された角周波数らしい。
値の範囲は0〜2πあるいは-π〜π。
恐らくサンプリング周波数が2πに対応する。
http://www.ic.is.tohoku.ac.jp/~swk/lecture/yaruodsp/dfanalysis.html
どうやら、これを見る限りωは正規化された角周波数らしい。
値の範囲は0〜2πあるいは-π〜π。
恐らくサンプリング周波数が2πに対応する。
こういち
2022/6/28 7:42
◆.Id/aHiU36hu
https://www.desmos.com/calculator?lang=ja
このサイトに
\sqrt{\frac{\left(1+\cos\left(2x\right)\right)^{2}+\left(\sin\left(2x\right)\right)^{2}}{\left(1-0.9\cos\left(x\right)+0.81\cos\left(2x\right)\right)^{2}+\left(-0.9\sin\left(x\right)+0.81\sin\left(2x\right)\right)^{2}}}
と打てばやらない夫のそれと同じものが出てくる。
分母は確かにフーリエ変換で間違いないと思うけど…
このサイトに
\sqrt{\frac{\left(1+\cos\left(2x\right)\right)^{2}+\left(\sin\left(2x\right)\right)^{2}}{\left(1-0.9\cos\left(x\right)+0.81\cos\left(2x\right)\right)^{2}+\left(-0.9\sin\left(x\right)+0.81\sin\left(2x\right)\right)^{2}}}
と打てばやらない夫のそれと同じものが出てくる。
分母は確かにフーリエ変換で間違いないと思うけど…
こういち
2022/6/28 22:40
◆.Id/aHiU36hu
フィルタの係数
2.116242857120330
-2.172584486211329
1.697035273949253
-0.650806197479927
0.173409469533881
-0.025655642361743
0.009731235000507
-0.5568342907644296
0.547592907644296
-0.335117867788814
0.128042559178572
0.0029580011934110
-0.014003558058372
0.014117614068190
-0.001219932531336
0.019441828241815
-0.123288785371843
0.302883788909530
-0.279135660772881
0.151630312909496
0.166907193868522
-0.441702933656371
0.360596216436326
-0.213756511891678
0.064778697282373
2.116242857120330
-2.172584486211329
1.697035273949253
-0.650806197479927
0.173409469533881
-0.025655642361743
0.009731235000507
-0.5568342907644296
0.547592907644296
-0.335117867788814
0.128042559178572
0.0029580011934110
-0.014003558058372
0.014117614068190
-0.001219932531336
0.019441828241815
-0.123288785371843
0.302883788909530
-0.279135660772881
0.151630312909496
0.166907193868522
-0.441702933656371
0.360596216436326
-0.213756511891678
0.064778697282373
こういち
2022/6/28 22:52
◆.Id/aHiU36hu
…！
窓関数が悪さしてる可能性あるかも？
ただ、窓関数を掛けないと、変な成分が出てしまう。
なら、FFTを使わずに定義通り計算したらどうなるだろうか…
窓関数が悪さしてる可能性あるかも？
ただ、窓関数を掛けないと、変な成分が出てしまう。
なら、FFTを使わずに定義通り計算したらどうなるだろうか…
こういち
2022/6/29 12:08
◆.Id/aHiU36hu
さて、窓関数が原因として、どうしようか…
C言語ならともかく、SmileBASICはコードを書かない方が速い言語。
素直に書くと数百倍ぐらい時間がかかりそう。
…IFFTってどういう定義だっけ。
C言語ならともかく、SmileBASICはコードを書かない方が速い言語。
素直に書くと数百倍ぐらい時間がかかりそう。
…IFFTってどういう定義だっけ。
こういち
2022/6/29 12:21
◆.Id/aHiU36hu
IFFTダメか。
…最終的にグラフを描画するのが目的なら、計算する点は400あれば十分。(つーか、これが限界)
あとはインターレースと併用すれば、気にならない程度にはなるんじゃないかな。
…最終的にグラフを描画するのが目的なら、計算する点は400あれば十分。(つーか、これが限界)
あとはインターレースと併用すれば、気にならない程度にはなるんじゃないかな。
こういち
2022/6/29 23:02
◆.Id/aHiU36hu
コメントを書く
- こちらは「プチコン3号」「プチコンBIG」など、プチコンシリーズに関する話題を扱ったコミュニティです
- プチコンシリーズにまったく関係ない書き込みはご遠慮下さい。削除の対象となります
- こちらにはその他のゲームや雑談のコミュニティはなく、作る予定もありません (ひとりで管理できないため)。ごめんなさい
- ユーザー登録なしで書き込みができます
- 秘密の合い言葉は成りすましの防止 (トリップ機能)、書き込みの編集時の本人認証に使用します
- 秘密の合い言葉に他人に推測されやすい言葉、他サービスと同じパスワードは入力しないでください。
- 書き込むと、投稿時に入力したお名前と秘密の暗号が記憶され、ログイン状態になります