ついでに、16KHzでの減衰が大きすぎるなら (1-b)*(1+Z)/2(1-bZ) の実現式

out+=k*(in+old -2*out)
old = in


k = (1-b)/2
k = (1 -0.37)/2= 0.685 で 16KHz -3dB です。

>>790　より少しだけ重いですが、特性は 16KHzより下で減衰が減り　高い所でよく減衰してくれます。