fs=22050;
x1=wavread('音樂.wav');
wp=2*pi*1000/fs;
ws=2*pi*1200/fs;

wdelta=ws-wp;
N=ceil(8*pi/wdelta);       %取整
wn=(wp+ws)/2;

beta=3.39;
window=kaiser(N+1,beta); 
[b,a]=fir1(N,wn/pi,'low',window); 

figure(1)
freqz(b,a,512);
title('FIR低通濾波器');
f2=filter(b,a,x1);
figure(2)
subplot(2,1,1)
plot(x1)
title('FIR低通濾波器濾波前的時域波形');
subplot(2,1,2)
plot(f2);
title('FIR低通濾波器濾波后的時域波形');
F0=fft(f2,1024);
f=fs*(0:511)/1024;
figure(3)
y2=fft(x1,1024);
subplot(2,1,1);
plot(f,abs(y2(1:512)));
title('FIR低通濾波器濾波前的頻譜')
xlabel('頻率/Hz');
ylabel('幅值');
subplot(2,1,2)
F2=plot(f,abs(F0(1:512)));
title('FIR低通濾波器濾波后的頻譜')
xlabel('頻率/Hz');
ylabel('幅值'); 
sound(f2,44100);    %濾波后的語音信號