ex0707.m

来自「数字信号处理（机械工业出版社）的源码」· M 代码 · 共 9 行

% Chapter 7: Example 7.7 filter design
%	
Wp = [60 200]/500; Ws = [50 250]/500;
Rp = 3; Rs = 40;
[n,Wn] = buttord(Wp,Ws,Rp,Rs)
[b,a] = butter(n,Wn);
freqz(b,a,128,1000)
title('n=16 Butterworth Bandpass Filter')