%巴特沃斯模擬濾波器原型

n=input('N= ');     %由使用者輸入濾波器階數(shù)N
%計(jì)算n階模擬低通原型，得到左半平面零極點(diǎn)
[z0,p0,k0]=buttap(n); 
b0=k0*real(poly(z0))     %求濾波器系數(shù)b0
a0=real(poly(p0))        %求濾波器系數(shù)a0
freqs(b0,a0);            %顯示系統(tǒng)的頻率特性
Pb=poly2str(b0,'s')      %給出b0決定的關(guān)于s多項(xiàng)式
Pa=poly2str(a0,'s')      %給出a0決定的關(guān)于s多項(xiàng)式