ex0706.m

数字信号处理（机械工业出版社）的源码

% Chapter 7: Example 7.6 filter design
%	
Wp = 40/500; Ws = 150/500;
[n,Wn] = buttord(Wp,Ws,3,60)
[b,a] = butter(n,Wn);
freqz(b,a,512,1000); title('n=5 Butterworth Lowpass Filter')