%Samp5_14
wp=[0.35 0.65]*pi;   %給出濾波器的通帶邊界頻率，用弧度表示
ws=[0.2 0.8]*pi;   %給出濾波器的阻帶邊界頻率，，用弧度表示
Rp=1;Rs=60;   %給出濾波器的通帶波紋和阻帶衰減
[n,wn]=cheb1ord(wp,ws,Rp,Rs,'s');  %根據濾波器的要求求濾波器的最小階數
[z,p,k]=cheb1ap(n,Rp);  %用求得的濾波器的最小階數設計Chebyshev I型濾波器
[b,a]=zp2tf(z,p,k); %將零點極點增益形式轉換為傳遞函數形式
Wo=sqrt(wn(1)*wn(2)); %計算中心點頻率
Bw=wn(2)-wn(1);  %計算頻帶寬度
[bt,at]=lp2bp(b,a,Wo,Bw);  %模擬原型低通濾波器轉換為帶通濾波器
[h,w]=freqs(bt,at,128);  %計算濾波器的復數頻率響應
plot(w/pi,20*log10(abs(h)));   %繪制濾波器的振幅響應
xlabel('角頻率/pi');
ylabel('|H(jw)|/dB')
grid on;hold on;
plot([0.2 0.2],ylim, 'r'); plot([0.8 0.8],ylim, 'r');%繪出阻帶界限以顯示達到的阻帶衰減
%plot([0.35 0.35],ylim); plot([0.65 0.65],ylim);%可繪出通帶界限以顯示達到的通帶衰減