%用fdatool設計濾波器，生成一個.fda文件，不知道怎么導入數據，查找資料說的是用simulink的濾波器模塊可以導入.fda生成的濾波器
%所以只用了程序來實現  
c=48;                        %定義濾波器階數
fc1=8400;fc2=13200;fs=48000;
w1=2*pi*fc1/fs;
w2=2*pi*fc2/fs;               %參數轉換，將模擬濾波器的技術指標轉換為數字濾波器的技術指標
window=hamming(c+1);        %使用hamming窗函數  根據不同的衰減要求選擇窗函數
h=fir1(c,[w1/pi w2/pi],window);  %使用標準響應的加窗設計函數fir1
freqz(h,1,512);                %數字濾波器頻率響應
y=wavread('yuanshi.wav');
ynew=filter(h,1,y);        %濾波處理

t=(0:length(ynew)-1)/fs;figure;   
subplot(2,1,1);plot(t,y);xlabel('time(s)');

subplot(2,1,2);plot(t,ynew);xlabel('Time(s)');