zbuttap1.m

来自「用MATLAB实现数字信号处理的程序」· M 代码 · 共 14 行

%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]);