close all, clear all, clc fs = 300; fs2 = fs/2; B = [1 0 0 0 0 -1]/5; [H, w] = freqz(B, 1); plot(w/(2*pi)*fs, abs(H)), grid