MATLABでバタワースハイパスフィルタをつくる

4 次のバタワースハイパスフィルターをつくってみる。
1000 Hz でサンプリングされたデータに対し、カットオフ周波数 30Hzを指定。
その振幅応答と位相応答を表示。

[z, p, k] = butter(4, 30/1000,"high");
sos = zp2sos(z,p,k);
fvttool(sos,"Analysis", "freq");

MATLAB修行中。