%MAT63.M TWO PASSBANDS-FIR WITH 63 COEFFICIENTS USING MATLAB Fs=10 kHz
nu=[0 0.1 0.12 0.18 0.2 0.3 0.32 0.38 0.4 1];  %normalized frequencies
mag=[0 0 1 1 0 0 1 1 0 0];     %magnitude at normalized frequencies 
c=remez(62,nu,mag);            %invoke remez algorithm for 63 coeff
bp63=c';                       % coeff values transposed
save mat2bp.cof bp63 -ascii;   %save in ASCII file with coefficients
[h,w]=freqz(c,1,256);          %frequency response with 256 points
plot(5000*nu,mag,w/pi,abs(h))  %plot ideal magnitude response

