% 2/9/2010 Dr. Durant EE3220 W9D2 - linear phase FIR filter design
format compact
format short g
fs = 22050;
f = [0 5000 7000 11025] / (fs/2); % normalized frequency bands
a = [1 1 0 0]; % passband and stopband gain
for N = [100 10] % compare short with long length effects
b = firpm(N,f,a); % Remez Exchange / Parks-McClellan algorithm
% B=FIRPM(N,F,A) % length, frequency bands, amplitude; fs/2 corresponds to 1
figure, stem(b)
figure, freqz(b)
figure, zplane(b) % input coefficients, row understood to be in B/A form, not Z/P form
end