buttord数字滤波器-matlab

clear  ; clc; close all;

fs =16000;
fs2=fs/2;                       % 降采样后采样频率的一半
fp1= 3000;                     % 通带频率
fs1= 4000;                     % 阻带频率
wp1=fp1/fs2;                    % 归一化通带频率
ws1=fs1/fs2;                    % 归一化阻带频率
Ap=3; As=30;                    % 通带波纹和阻带衰减
[n,Wn]=buttord(wp1,ws1,Ap,As);  % 求滤波器原型阶数和带宽
[bn1,an1]=butter(n,Wn);         % 求数字滤波器系数
[H,f]=freqz(bn1,an1,1000,fs);   % 求数字滤波器幅频曲线

% y1=filter(bn1,an1,x);           % 对降采样后的数据进行滤波

figure(2)
plot(f,abs(H),'k');
grid; 
xlabel('频率/Hz'); ylabel('幅值')
title('巴特沃斯滤波器的幅值响应')
set(gcf,'color','w');