chebhigh2.m

来自「自己做的： （一）用双线性变换法设计并用实验系统实现一个三阶的契比雪夫Ⅰ型低通」· M 代码 · 共 14 行

function [Fs,N,ai,bi]=chebhigh2(fp,ap)
Fs=16000;T=1/Fs;fs=4000;
fp=input('请输入通带临界频率值fp=：');
ap=input('请输入通带衰减ap=:');
wp=fp/Fs*2*pi;
ws=fs/Fs*2*pi;
OmegaP=2*Fs*tan(wp/2);
OmegaS=2*Fs*tan(ws/2);
[N,Wn]=cheb2ord(OmegaP,OmegaS,ap,20,'s');
[b,a]=cheby2(N,ap,Wn,'high','s');
freqs(b,a);
[bi,ai]=bilinear(b,a,Fs);
x2=[Fs N ai bi];
disp(x2(1:4));