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

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

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

figure(1)
freqz(b,a,512);
title('FIR高通濾波器');

f2=filter(b,a,x1);

figure(2)
subplot(2,1,1)
plot(x1)
title('FIR高通濾波器濾波前的時(shí)域波形');
subplot(2,1,2)
plot(f2);
title('FIR高通濾波器濾波后的時(shí)域波形');
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);    %濾波后的語音信號(hào)