daitong.m

来自「用MATLAB编程实现巴特沃兹低通、带通滤波器」· M 代码 · 共 12 行

fs=1000;  %采样频率
Wp=[200 300]/(fs/2);
Ws=[190 310]/(fs/2);  %通带截止频和阻截频
Rp=1;Rs=15;                          %通带衰减和阻带衰减  
[N, Wn] = BUTTORD(Wp, Ws, Rp, Rs);   %求出滤波器的阶数和固有频率
[B,A] = butter(N,Wn,'bandpass');     %返回滤波器的传输函数
[H,W]=freqz(B,A,1024);      %H为离散系统响应，W则包含范围内N个频率等分点
plot(W*fs/(2*pi),abs(H));
title('巴特沃兹带通滤波器');
xlabel('频率');
ylabel('幅度');