[y,fs,bits]=wavread('音樂.wav');
sound(y,fs);

figure(1)
plot(y)                   
title('聲音信號的時域波形')
xlabel('時間');
ylabel('幅度');
figure(2)
freqz(y)           %繪制信號采樣后的頻率響應圖
title('聲音經采樣后的頻率響應')


b=1;
while(b);
k=menu('選擇濾波器類:','1、FIR濾波器','2、IIR濾波器','3、退出   ');
if k==1
   i=menu('請選擇FIR濾波器類型:','低通','帶通','高通','退出');
    if i==1
           FIR_L;
       else if i==2
              FIR_B;
        else if i==3
                  FIR_H; 
              else if  i==4
                 close all;return ; 
              end
            end
        end
    end
else if k==2
        j=menu('請選擇IIR濾波器類型:','低通','帶通','高通','退出');
        if j==1
              IIR_L; 
        else if j==2
                 IIR_B;
            else if j==3
                      IIR_H; 
                  else if  i==4
                 close all;return ;   
                 end
               end
            end
        end
    end
if k==3
    close all;return ;
end
end
b=menu(' 還要繼續嗎? ','  是  ','  否  ');
if b==2
    b=0;
end
end