最後に実数だけとるのを作った。フィルタとして設計するhttps://nbviewer.jupyter.org/github/cho45/pseudo-audio-signal/blob/master/docs/02-real-fir.ipynbフィルタを左右対称にすればほぼ実数フィルタになるはずだが、低域の減衰が急峻すぎるせいかこれでも虚数が出てきてしまう。しかたないのでこのまま虚数を捨てて実数フィルタとして評価した。うまくいくけど、かなり重くなってしまう。特性から誤差が発生して気持ちがわるい。次数を増やせばそれだけ近似はできるが、どんどん重くなっていく。フィルタで補正する↑補正なし↓補正https://nbviewer.jupyter.org/github/cho45/pseudo-audio-signal/blob/master/docs/03-iir-fir.ipynbアナログフィルタ(周波数連続という意味)からデジタルフィルタ(離散という意味)をつくる方法の一つに双一次変換がある。簡単にできるので、これを試す。今回の場合、アナログの伝達関数はわかっているので、この規格に書いてある伝達関数を双一次変換してIIRフィルタの係数を得る。以下を変換してデジタルフィルタを得たい。微妙に括弧がついてるのかわからないが、結局これは次数ごとに書きなおすと見なれた以下の形になる。受けとる形になっている。