24t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 9 行

clc;close all;clear all;
Rp=1;Rs=15;fs=2000;
wp1=0.2*pi;ws1=0.3*pi;
wp2=(2*fs)*tan(wp1/2);
ws2=(2*fs)*tan(ws1/2);
[n,wn]=cheb1ord(wp2,ws2,Rp,Rs,'s');
[b,a]=cheby1(n,Rp,wn,'s');
[bz,az]=bilinear(b,a,fs);
[db,mag,pha,grd,w]=freqz_m(bz,az);plot(w/pi,db);