グラフ上マウス等でf0付近をドラッグすればf0を変更可
f0:=1000;w=2PI*f0/fs; w/(1-(1-w)Z)
f0:=1000;w=2PI*f0/fs; (1-Z)/(1-(1-w)Z)
L1
H1
f0:=1000;L1:=0.5;H1:=1.5;w=2PI*f0/fs; (w*Z*L1+(1-Z)*H1)/(1-(1-w)Z)
f0:=1000;w=2PI*f0/fs; (-(1-w)+Z)/(1-(1-w)Z)
f0:=1000;w=2PI*f0/fs; w/2*(1+Z)/(1-(1-w)Z)
fs/2でカットされるのでアナログよりも高域では急峻となります。
やはりf0>3KHzでアナログより高くなるので
w=2PI*f0/fs;をw=sin(2PI*f0/fs)とすれば若干改善しますがf0<fs/4までです
もっと良い補正方法は
f0:=1000;w=2PI*f0/fs; w =tan(w/ 2); b1=(1 - w) / (1 + w); a1=(1 - b1) / 2; a1*(1+Z)/(1-b1Z)
f0:=1000;w=2PI*f0/fs; (1-w/2)*(1-Z)/(1-(1-w)Z)
f0:=1000;w=2PI*f0/fs; w =tan(w/ 2); b1=(1 - w) / (1 + w); a1=(1 -b1)/2/w; a1*(1-Z)/(1-b1Z)
L0
H0
f0:=1000;L0:=0.5;H0:=1.5;w=2PI*f0/fs; w =tan(w/ 2); b1=(1 - w) / (1 + w); k1=(1 - b1) / 2; a0=k1*(L0+H0/w); a1=k1*(L0-H0/w); (a0+a1Z)/(1-b1Z)