いま >>440をZ変換してあげてもいいんだけど、それじゃ勉強にならないだろうから
ttp://web.kyoto-inet.or.jp/org/create21/dsp06-2.htm

の13式から15式を作り
それを複素数の扱えるプログラム言語で書けばそのまま周波数特性が出ます