%Buttap函數：設計巴特沃斯濾波器
Wc=10000;
Ws=40000;
Ap=3;
As=35;
Np=sqrt(10^(0.1*Ap)-1);
Ns=sqrt(10^(0.1*As)-1);
n=ceil(log10(Ns/Np)/log10(Ws/Wc));
[z,p,k]=buttap(n);
syms rad;
Hs1=k/(i*rad/Wc-p(1))/(i*rad/Wc-p(2))/(i*rad/Wc-p(3));
Hs2=10*log10((abs(Hs1))^2);
ezplot(Hs2,[-60000,60000]);