ついでに
b と　カットオフ周波数f0との関係は

w=2*PI*f0/サンプル周波数
r=tan(w/2)
b=(1-r)/(1+r)

です。 この式は 双１次変換 でアナログ1次フィルタとマッピングしたもので
高域では少しアナログと比べて性能がキツメですね。