wp=0.2*pi;
    ws=0.3*pi;
    Rp=1;
    Rs=15;
    Ts=0.02*pi;
    Fs=1/Ts;
%選擇濾波器的最小階數
    [N,Wn]=buttord(wp,ws,Rp,Rs,'s');
    %創建butterworth模擬濾波器
    [Z,P,K]=buttap(N);
    %把濾波器零極點模型轉化為傳遞函數模型
    [Bap,Aap]=zp2tf(Z,P,K);
    %把模擬濾波器原型轉換成截至頻率為Wn的低通濾波器
    [b,a]=lp2lp(Bap,Aap,Wn);
    %用雙線性變換法實現模擬濾波器到數字濾波器的轉換
    [bz,az]=bilinear(b,a,Fs);
    %繪制頻率響應曲線
    [H,W]=freqz(bz,az);
    plot(W*Fs/(2*pi),abs(H));
    grid
    xlabel('頻率／Hz')
    ylabel('頻率響應幅度')