chb1.m

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

function [a,b]=chb1(Wp,Ws,Rp,As)
%根据技术指标设计切比雪夫滤低通波器
if Wp<=0
    error('通带截止频率Wp必须大于0')
end
if Ws<=0
    error('阻带截止频率Ws必须大于0')
end
if (Rp<=0)|(As<=0)
    error('通带最大衰减Rp和通阻最大衰减As必须都大于0')
end
ep=sqrt(10^(Rp/10)-1);
A=10^(As/20);
OmegaC=Wp;
OmegaR=Ws/Wp;
g=sqrt(A*A-1)/ep;
N=ceil((log10(g+sqrt(g*g-1)))/(log10(OmegaR+sqrt(OmegaR*OmegaR-1))));
[b,a]=u_chb1ap(N,Rp,OmegaC);