butwthband.m

来自「用matlab编写的数字信号处理程序」· M 代码 · 共 18 行

clc,clear all,close all;
Wp=[0.2 0.3];                   %通带两个边界角频率
Ws=[0.1 0.4];                   %阻带两个边界角频率
Rp=-1;                          %通带允许起伏
Rs=-25;                         %阻带最小衰减  
[N,Wn]=buttord(Wp,Ws,Rp,Rs);    %用buttord确定低通原型参数
[B,A]=butter(N,Wn);             %用butter生成带通滤波器
[H,w]=freqz(B,A,2*N);
subplot(2,1,1);
plot(w/pi,abs(H));
xlabel('角频率（单位pi）');
ylabel('滤波器幅度');
title('幅频特性');
subplot(2,1,2);
plot(w/pi,angle(H));
xlabel('角频率（单位pi）');
ylabel('对应相位');
title('相频特性');