レート変換用のインターポーレーションフィルタ
実はディジタルフィルタの設計をちゃんと勉強したことないので、いろいろハマりました。SOX のマニュアルを参考に、Nuttall 窓をかけた FIR フィルタを作りました。これまた SOX を参考に、1024 次、カットオフ周波数 0.95 (2倍のインターポーレーションなので、実際は 0.475)のローパスフィルタを作り、そこに Nuttall 窓をかけました。
Octave のコマンドは、こんな感じー。
f = fopen("coef.txt", "wt") fprintf(f, "%g\n", nuttallwin(1025) .* fir1(1024, 0.475)) fclose(f)