ローパス ベッセル フィルターの周波数応答

最大 $$10^4$$ ラジアン/秒の定数近似群遅延を持つ 5 次のアナログ ローパス ベッセル フィルターを設計します。freqs を使用してフィルターの振幅応答および位相応答をプロットします。

wc = 10000;
[b,a] = besself(5,wc);
freqs(b,a)

フィルターの群遅延応答を、アンラップされた位相応答の微分の負として計算します。カットオフ周波数までの範囲で近似的に定数であることを検証するために群遅延をプロットします。

[h,w] = freqs(b,a);
grpdel = -diff(unwrap(angle(h)))./diff(w);

clf
loglog(w(2:end),grpdel)
xlabel('Frequency (rad/s)')
ylabel('Group delay (s)')
xline(wc)
grid

参考

besself | freqs