zbutter.m

来自「用MATLAB实现数字信号处理的程序」· M 代码 · 共 8 行

%Butter函数：设计巴特沃斯滤波器
N=5;                            %阶数的一半
Fs=1000;
Wn=[100 200]/(Fs/2);            %得到归一化角频率
[b,a]=butter(N,Wn,'stop');
[y,t]=impz(b,a,81);
stem(y);