exa4_15.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 12 行

% exa4-15_butterband.m , for example 4-15
% using butterworth filter to design bandpass DF

clear all;
N= 10; 
Wn = [100 200]/500;
[b,a] = butter(N,Wn,'bandpass');
freqz(b,a,128,1000)
figure(2)
[y,t] = impz(b,a,101);
stem(t,y)