main3.m

按照给定指标采用模拟低通滤波器（巴特沃思低通）去逼近理想模拟低通滤波器

clear;
close all;

fp=input('通带截止频率fp(Hz)=');
fs=input('阻带截止频率fs(Hz)=');
Fs=input('采样频率Fs(Hz)=');
Ap=input('通带最大衰减Ap(dB)=');
As=input('阻带最小衰减As(dB)=');

Wp=2*pi*fp;  %rad/sec
Ws=2*pi*fs;

Fs=2*Fs;
%数字要求转换成模拟要求
Omegap=Fs*tan(Wp/Fs);
Omegas=Fs*tan(Ws/Fs);

[b,a,z,p,k]=afd_butt(Omegap,Omegas,Ap,As);

%进行双线性Z变换
[b1,a1]=mybilinear(z,p,k,Fs);
[db,mag,pha,grd,w]=freqz_m(b1,a1);%数字滤波器响应子程序
figure
subplot(211);plot(w/pi,db);title('所设计的数字滤波器幅度响应(dB)');grid on; axis([0,1,-100,10]);

%使用toolbox得到的数字滤波器
[N, Wc] = buttord(0.2,0.6,3,20);
[b2, a2] = butter(N, Wc)
[db,mag,pha,grd,w]=freqz_m(b2,a2);
subplot(212);plot(w/pi,db);title('使用toolbox得到的数字滤波器幅度响应(dB)');grid on; axis([0,1,-100,10]);