% IIR_bs_2.m
% Calculating IIR filter coefficients
% Numerator coefficients
B0 = 31135;
B1 = -29652 * 2;
B2 = 31135 ;
% Denominator coefficients
A0 = 32767; % Always 1.0000
A1 = -29652 * 2;
A2 = 29503;
b = [B0 B1 B2];
a = [A0 A1 A2];
Fs = 48000 % Sampling frequency
[H,f] = freqz ( b, a, 512, Fs)
magnitude = 20 * log10(abs(H));
plot ( f, magnitude), grid on
xlabel ('Frequency (Hz)');
ylabel ('Magnitude Response (dB)');
title ('2000 Hz to 2800 Hz Bandstop Filter Using BLT. 48000 Hz Sampling');
set (gca, 'xlim', [0, 12000])
set (gca, 'XTick', [0, 2000, 4000, 6000, 8000, 10000, 12000])
set (gca, 'XTickLabel', [0 2000 4000 6000 8000 10000, 12000])
set (gca, 'ylim', [-42, 6])
set (gca, 'YTick', [-42, -36, -30, -24, -18, -12, -6, 0, 6])
set (gca, 'YTickLabel', [-42, -36,-30, -24, -18, -12, -6, 0, 6])
%zplane (b,a)