レート変換用のインターポーレーションフィルタ

実はディジタルフィルタの設計をちゃんと勉強したことないので、いろいろハマりました。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)