firbandpass.m

matlab的滤波器

fb1=1200;
fb2=3000;
fc1=1000;
fc2=3200;
fs=18000; 
wb1=2*pi*fb1/fs;
wc1=2*pi*fc1/fs;
wb2=2*pi*fb2/fs;
wc2=2*pi*fc2/fs;
deltaw=(wc2-wb2)/(2*pi);
N0=ceil((100-7.95)/(14.36*deltaw))+1;
N=N0+mod(N0+1,2)
beta=0.1102*(100-8.7);
[b,a]=fir1(N-1,[wb1/pi,wb2/pi],kaiser(N,beta));
[x1,fs]=wavread('music');
x2=fftfilt([b,a],x1);
sound(x2,44100);
f=fft(x2,length(x1));
printf(x2,abs(f),abs(length(x1)),b,a);