clc
close all
clear all
%给出技术指标:通带纹波ap=3dB,阻带衰减as=30dB,通带截止频率fp=150Hz,阻带截止频率fs=400Hz,抽样频率Fs=1000Hz
fp=150;fs=400;Fs=1000;
ap=3;as=30;
wp=2*pi*fp/Fs;
ws=2*pi*fs/Fs;
wap=2*Fs*tan(wp/2);was=2*Fs*tan(ws/2);
[n,wn]=buttord(wap,was,ap,as,'s');
%求出巴特沃思模拟滤波器的介数n及截止频率wn，注意不能省去s
[z,p,k]=buttap(n);
%创建巴特沃思低通滤波器原型
[num,den]=zp2tf(z,p,k);
%将系统的零级点增益模型转换为传递函数模型
[B,A]=lp2lp(num,den,wap)
%实现低通向低通的转变
[num1,den1]=bilinear(B,A,Fs)
%用双线性变换法把模拟滤波器转变成数字滤波器
[h,w]=freqz(num1,den1,500,Fs);
%求出频率响应
plot(w,abs(h));grid on;
%画出频率响应曲线
xlabel('频率/Hz');ylabel('幅值')
评论0