cheb1hpf.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 15 行

function [b,a]=cheb1hpf(wp,ws,Rp,As)
wplp=0.2*pi;
alpha=-(cos((wplp+wp)/2))/(cos((wplp-wp)/2));
wslp=-angle(-(exp(-j*ws)+alpha)/(1+alpha*exp(-j*ws)));
T=1;Fs=1/T;
OmegaP=(2/T)*tan(wplp/2);
OmegaS=(2/T)*tan(wslp/2);
% Design Analog Chebyshev Prototype Lowpass Filter:
[cs,ds]=afd_chb1(OmegaS,OmegaP,Rp,As);
% Perform Bilinear transformation to obtain digital lowpass
[blp,alp]=bilinear(cs,ds,Fs);
% Transform digital lowpass into highpass filter
Nz=-[alpha,1];Dz=[1,alpha];
[b,a]=zmapping(blp,alp,Nz,Dz);