xi=importdata('文件名');
number=變換點數;
fs=采樣頻率;
y=fft(xi,number);
再選定妳要提取的頻帶範圍（n：m）
yt=y（n：m）；
xt=ifft（yt）；
還有，可以直接用設計的FIR 或IIR濾波器在時域濾波