%Samp5_9
[z,p,k]=buttap(6);   %設計階數為6的Butterworth模擬原型濾波器
[b,a]=zp2tf(z,p,k);  %將零點極點增益形式轉換為傳遞函數形式
[H,w]=freqs(b,a,0:0.01:2); %計算指定頻率點的復數頻率響應
subplot(2,1,1),plot(w,abs(H).^2); %繪制平方幅頻響應
xlabel('w/wc');ylabel('Butterworth |H(jw)|^2');
title('Btterworth低通原型濾波器(wc=1)')
w1=0.7;w2=1.5;   %帶阻濾波器的下限及上限頻率
w0=sqrt(w1*w2);bw=w2-w1;   %計算中心點頻率和阻帶寬度
[bt,at]=lp2bs(b,a,w0,bw);   %頻率轉換
[Ht,wt]=freqs(bt,at,0:0.01:2);  %計算帶阻濾波器指定點的復數頻率響應
subplot(2,1,2),plot(wt,abs(Ht).^2);%繪制平方幅頻響應
xlabel('w/wc');ylabel('Butterworth |H(jw)|^2');
title('Butterworth帶阻濾波器(wc=0.7~1.5) ')