IIR滤波器设计
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
Wp = 40/500; Ws = 150/500; %归一化频率计算
[n,Wn] = buttord(Wp,Ws,3,60); %计算butterworth滤波器的阶数n和截止频率wn
% Returns n = 5; Wn=0.0810;
[b,a] = butter(n,Wn); %设计butterworth滤波器,返回滤波器系数
freqz(b,a,512,1000); %画出幅频特性
title('n=5 Butterworth Lowpass Filter')
[n,Wn] = buttord(Wp,Ws,3,60); %计算butterworth滤波器的阶数n和截止频率wn
% Returns n = 5; Wn=0.0810;
[b,a] = butter(n,Wn); %设计butterworth滤波器,返回滤波器系数
freqz(b,a,512,1000); %画出幅频特性
title('n=5 Butterworth Lowpass Filter')
%%butterworth
[nb,wnb]=buttord(0.3,0.4,1,50);
[bb,ab]=butter(nb,wnb);
fvtool(bb,ab)
%%Elliptic
[nb,wnb]=ellipord(0.3,0.4,1,50);
[bb,ab]=ellip(nb,1,50,wnb);
fvtool(bb,ab)
[nb,wnb]=buttord(0.3,0.4,1,50);
[bb,ab]=butter(nb,wnb);
fvtool(bb,ab)
%%Elliptic
[nb,wnb]=ellipord(0.3,0.4,1,50);
[bb,ab]=ellip(nb,1,50,wnb);
fvtool(bb,ab)
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
Wp = [60 200]/500; Ws = [10 250]/500;
Rp = 3; Rs = 40;
[n,Wn] = buttord(Wp,Ws,Rp,Rs);
% Returns n =16; Wn =[0.1198 0.4005];
[b,a] = butter(n,Wn);
freqz(b,a,128,1000)
title('n=16 Butterworth Bandpass Filter')
Rp = 3; Rs = 40;
[n,Wn] = buttord(Wp,Ws,Rp,Rs);
% Returns n =16; Wn =[0.1198 0.4005];
[b,a] = butter(n,Wn);
freqz(b,a,128,1000)
title('n=16 Butterworth Bandpass Filter')
作者:wanglei_wan
本文版权归作者和博客园共有,欢迎转载,但未经作者同意必须保留此段声明,且在文章页面明显位置给出原文连接,否则保留追究法律责任的权利。