butt.m

来自「这个是我数字信号处理这门课的课程设计,文件包括用matlab编写的代码和设计报告」· M 代码 · 共 20 行

function [a,b]=butt(Wp,Ws,Rp,As,T)
%根据技术指标设计数字巴特沃斯滤低通波器
if Wp<=0
    error('通带截止频率Wp必须大于0')
end
if Ws<=0
    error('阻带截止频率Ws必须大于0')
end
if (Rp<=0)|(As<=0)
    error('通带最大衰减Rp和通阻最大衰减As必须都大于0')
end
Wp=Wp/T;
Ws=Ws/T;
ksp=((10^(Rp/10)-1)/(10^(As/10)-1))^0.5;
asp=Ws/Wp;
N=ceil(-(log10(ksp))/(log10(asp)));
Omegac=Wp/((10^(Rp/10)-1)^(1/(2*N)));
[b,a]=u_buttap(N,Omegac);
[b,a]=imp_invr(b,a,T);