%切比雪夫Ⅰ型模擬濾波器原型

n=input('N= ');           %由使用者輸入濾波器階數N
Rp=input('Rp= ');         %由使用者輸入濾波器通帶衰減指標
[z0,p0,k0]=cheb1ap(n,Rp); %計算n階模擬低通原型，得到左半平面零極點
b0=k0*real(poly(z0))     %求濾波器系數b0
a0=real(poly(p0))        %求濾波器系數a0
freqs(b0,a0);             %求系統的頻率特性
Pb=poly2str(b0,'s')       %給出b0決定的關于s多項式
Pa=poly2str(a0,'s')           %給出a0決定的關于s多項式