设计一个低通巴特沃斯数字滤波器，采样频率 1000Hz，通带边界频率 100Hz，通带衰减为 1dB，阻带边界频率为 200Hz，阻带衰减为 60dB

设计一个低通巴特沃斯数字滤波器，采样频率 1000Hz ，通带边界频率

100Hz ，通带衰减为 1dB ，阻带边界频率为 200Hz ，阻带衰减为 60dB, 显

示出滤波器的阶数，画出该滤波器的幅频图和相频图，观察滤波器的幅

频曲线和相频曲线的特点。其它参数不变的情况下，通带衰减改为 3dB ，

分析滤波器的幅频曲线、相频曲线。

clear;close all;clc

%1.数字滤波器的技术指标

ap = 1 ;%通带最大衰减

as = 60;%阻带最小衰减

fp = 100;%通带截止频率

fs = 200;%阻带截止频率

Fs = 1000;%抽样间隔

T = 1/Fs;

%2.将数字指标转化成模拟滤波器技术指标（采用冲激响应不变法）

wap = 2*pi*fp;

was = 2*pi*fs;

%3.设计滤波器

[N,wac] = buttord(wap,was,ap,as,'s');%计算阶数N和3dB截止频率wac

[z,p,k] = buttap(N);%创建巴特沃斯低通滤波器 z零点p极点k增益

[Bap,Aap] = zp2tf(z,p,k);%由零极点和增益确定归一化Han(s)系数

[Bbs,Abs] = lp2lp(Bap,Aap,wac);%将s/wc 代替s，去归一化

[B,A] = impinvar(Bbs,Abs,Fs);%模拟域到数字域——冲激响应不变法

[H1,w] = freqz(B,A);%根据H(z)求频率响应

%绘制数字滤波器频响幅度谱

figure(1);

f = w*Fs/(2*pi);

subplot(2,1,1);

plot(f,20*log10(abs(H1)));%绘制幅度响应

title('冲激响应不变法——巴特沃斯BLPF(幅度)');

xlabel('频率/Hz');

ylabel('H1幅值/dB');

subplot(2,1,2);

plot(f,unwrap(angle(H1)));% 绘制相位响应

xlabel('频率/Hz');

ylabel('角度/Rad');

title('冲激响应不变法——巴特沃斯BLPF(相位)');