引き続き発散現象調査。
対策として、フィルタ係数のクリップ処理を入れてしまうのはどうかと考えている。
なるべく、というか厳密にフィルタの出力を変えないように係数をクリップしたいのだが、どうすればいいのか。その議論ってあんまりされてないように見えるんやがな。(FIRフィルタ係数の同値類ってなんだ?)
色々見ると、(あんまり記述無いけど)
にあるように、全係数を定数で割ると案の定ゲインが定数倍だけ小さくなる。しかし周波数特性はそのまま。
周波数特性は保たれるからええか、というノリで定数倍の方策を取る。幅は18bit([-4,4]の範囲)とした。16, 17bitだと顕著に性能が悪い。しっかしオーバーフローが怖いかも。
テスト追加に入る。