模拟滤波器及IIR数字滤波器的设计.docx
《模拟滤波器及IIR数字滤波器的设计.docx》由会员分享,可在线阅读,更多相关《模拟滤波器及IIR数字滤波器的设计.docx(17页珍藏版)》请在冰豆网上搜索。
模拟滤波器及IIR数字滤波器的设计
实验三模拟滤波器及IIR数字滤波器的设计
一、模拟滤波器的设计
1.设计一个巴特沃斯模拟低通滤波器,以满足:
通带截止频率fp=5Hz,通带最大衰减
:
=2dB,阻带截止频率fs=12Hz,阻带最小衰减:
&=30dB。
要求绘出滤波器的
幅频特性曲线。
(幅度用分贝值表示)
理论分析:
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s')
其中,参数Wp和Ws分别是通带边界频率和阻带边界频率,Wp和Ws的单位是rad/s。
Rp和Rs分别为通带最大衰减和阻带最小衰减(dB)。
返回的参数N和Wn分别为滤波器的
阶数和3dB截止频率。
对于带通和带阻滤波器,Wp和Ws都是二维向量,向量的第一个元
素对应低端的边界频率,第二个元素对应高端的边界频率。
[B,A]=butter(N,Wn,'s')
其中,N和Wn分别为滤波器的阶数和3dB截止频率。
利用此函数可以获得低通和带
通滤波器系统函数的分子多项式(B)和分母多项式(A)的系数。
H=freqs(B,A,w)
其中,B和A分别表示滤波器系统函数的分子多项式和分母多项式的系数。
该函数返
回矢量w指定的那些频率点上的频率响应,w的单位是rad/s。
不带输出变量的freqs函数,
将绘制出幅频和相频曲线。
源程序:
wp=2*pi*5;
ws=2*pi*12;
rp=2;
rs=30;
[N,Wn]=buttord(wp,ws,rp,rs,'s');
[B,A]=butter(N,Wn,'s');
w=0:
300;
h=freqs(B,A,w);
H=20*log10(abs(h));
plot(w,H);
title('巴特沃斯低通滤波器的幅频特性');
xlabel('频率/Hz');
ylabel('幅度/db');
实验结果:
-70
-80
:
-p=3dB,阻带截止频率fs=10Hz,阻带最小衰减〉s=15dB。
要求绘出滤波器的
幅频特性曲线。
(幅度用分贝值表示)
理论分析:
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s')
其中,参数Wp和Ws分别是通带边界频率和阻带边界频率,Wp和Ws的单位是rad/s。
Rp和Rs分别为通带最大衰减和阻带最小衰减(dB)。
返回的参数N和Wn分别为滤波器的
阶数和3dB截止频率。
对于带通和带阻滤波器,Wp和Ws都是二维向量,向量的第一个元
素对应低端的边界频率,第二个元素对应高端的边界频率。
[B,A]=butter(N,Wn,'high','s')
可以获得高通滤波器系统函数的分子多项式(B)和分母多项式(A)的系数。
H=freqs(B,A,w)
其中,B和A分别表示滤波器系统函数的分子多项式和分母多项式的系数。
该函数返
回矢量w指定的那些频率点上的频率响应,w的单位是rad/s。
不带输出变量的freqs函数,
将绘制出幅频和相频曲线。
源程序:
wp=2*pi*20;
ws=2*pi*10;
rp=3;
rs=15;
[N,Wn]=buttord(wp,ws,rp,rs,'s');
[B,A]=butter(N,Wn,'high','s');
w=0:
400;
h=freqs(B,A,w)
H=20*log10(abs(h));
plot(w,H);
title('巴特沃斯高通滤波器的幅频特性');
xlabel('频率/Hz');
ylabel('幅度/db');
巴特沃斯高通滤波器的幅频特性
实验结果:
0
-20
-40
b-60
度
幅-80
-100
-120
-140
050100150200250300350400
频率/Hz
3.设计一个巴特沃斯模拟带通滤波器,以满足:
通带范围为10Hz〜25Hz,阻带截止频率
分别为5Hz、30Hz,通带最大衰减为3dB,阻带最小衰减为30dB。
要求绘出滤波器的
幅频特性曲线。
(幅度用分贝值表示)
理论分析:
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s')
其中,参数Wp和Ws分别是通带边界频率和阻带边界频率,Wp和Ws的单位是rad/s。
Rp和Rs分别为通带最大衰减和阻带最小衰减(dB)。
返回的参数N和Wn分别为滤波器的
阶数和3dB截止频率。
对于带通和带阻滤波器,Wp和Ws都是二维向量,向量的第一个元素对应低端的边界频率,第二个元素对应高端的边界频率。
[B,A]=butter(N,Wn,'s')
其中,N和Wn分别为滤波器的阶数和3dB截止频率。
利用此函数可以获得低通和带
通滤波器系统函数的分子多项式(B)和分母多项式(A)的系数。
H=freqs(B,A,w)
其中,B和A分别表示滤波器系统函数的分子多项式和分母多项式的系数。
该函数返回矢量w指定的那些频率点上的频率响应,w的单位是rad/s。
不带输出变量的freqs函数,
将绘制出幅频和相频曲线。
源程序:
wp=[2*pi*102*pi*25];
ws=[2*pi*52*pi*30];
rp=3;
rs=30;
[N,Wn]=buttord(wp,ws,rp,rs,'s');
[B,A]=butter(N,Wn,'s');
w=0:
1000;
h=freqs(B,A,w);
H=20*log10(abs(h));
plot(w,H);
title('巴特沃斯带通滤波器的幅频特性');
xlabel('频率/Hz');
ylabel('幅度/db');
实验结果:
4.设计一个巴特沃斯模拟带阻滤波器,以满足:
通带截止频率分别为10HZ、35HZ,阻带
截止频率分别为15HZ、30HZ,通带最大衰减为3dB,阻带最小衰减为30dB。
要求绘出滤波器的幅频特性曲线。
(幅度用分贝值表示)
理论分析:
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s')
其中,参数Wp和Ws分别是通带边界频率和阻带边界频率,Wp和Ws的单位是rad/s。
Rp和Rs分别为通带最大衰减和阻带最小衰减(dB)。
返回的参数N和Wn分别为滤波器的
阶数和3dB截止频率。
对于带通和带阻滤波器,Wp和Ws都是二维向量,向量的第一个元素对应低端的边界频率,第二个元素对应高端的边界频率。
[B,A]=butter(N,Wn,'stop','s')
可以获得带阻滤波器系统函数的分子多项式(B)和分母多项式(A)的系数。
H=freqs(B,A,w)
其中,B和A分别表示滤波器系统函数的分子多项式和分母多项式的系数。
该函数返回矢量w指定的那些频率点上的频率响应,w的单位是rad/s。
不带输出变量的freqs函数,
将绘制出幅频和相频曲线。
源程序:
wp=[2*pi*102*pi*35];
ws=[2*pi*152*pi*30];
rp=3;
rs=30;
[N,Wn]=buttord(wp,ws,rp,rs,'s');[B,A]=butter(N,Wn,'stop','s');
w=0:
400;
h=freqs(B,A,w);
H=20*log10(abs(h));
plot(w,H);
title('巴特沃斯带阻滤波器的幅频特性');
xlabel('频率/Hz');
ylabel('幅度/db');
实验结果:
0F「一
-50-
b-loo-「
度
幅-150--
J
f
-200--
\\
I;
-250--
II
-300[1[ccc1:
050100150200250300350400
频率/Hz
二、用脉冲响应不变法和双线性变换法设计IIR数字滤波器
1.要求分别用脉冲响应不变法和双线性变换法设计一个数字低通滤波器,以满足:
通带截
止频率为0.2二,阻带截止频率为0.3二,通带最大衰减为1dB,阻带最小衰减为15dB,
采样间隔设为1s。
理论分析:
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s')
其中,参数Wp和Ws分别是通带边界频率和阻带边界频率,Wp和Ws的单位是rad/s。
Rp和Rs分别为通带最大衰减和阻带最小衰减(dB)。
返回的参数N和Wn分别为滤波器的
阶数和3dB截止频率。
对于带通和带阻滤波器,Wp和Ws都是二维向量,向量的第一个元
素对应低端的边界频率,第二个元素对应高端的边界频率。
[B,A]=butter(N,Wn,'s')
其中,N和Wn分别为滤波器的阶数和3dB截止频率。
利用此函数可以获得低通和带
通滤波器系统函数的分子多项式(B)和分母多项式(A)的系数。
[bz,az]=impinvar(b,a,Fs)
可以实现用脉冲响应不变法将模拟滤波器转换为数字滤波器。
其中b和a分别是模拟滤
波器的系统函数H(s)的分子多项式和分母多项式的系数,Fs是脉冲响应不变法中的采样频
率,单位为Hz,如果Fs没有说明,其缺省值为1Hz。
运算的结果bz和az分别表示数字滤
0.91
波器的系统函数H(z)的分子多项式和分母多项式的系数。
[bz,az]=bilinear(b,a,Fs)
可以实现用双线性变换法将模拟滤波器转换为数字滤波器。
参数含义同上。
利用freqz函数计算数字滤波器的频率响应
源程序:
wp=0.2*pi;
ws=0.3*pi;
rp=1;
rs=15;
Fs=1;
[N,Wn]=buttord(wp,ws,rp,rs,'s');
[B,A]=butter(N,Wn,'s');
[bz,az]=impinvar(B,A,Fs);
[h,w]=freqz(bz,az);plot(w/pi,20*log10(abs(h)));
title('数字低通滤波器的幅频特性’);
xlabel('频率/Hz');
ylabel('幅度/db');
实验结果:
0
-10
-20
-30
b
度-40
幅
-50
-60
-70
-80
00.10.20.30.40.50.60.70.8
频率/Hz
源程序:
wp=0.2*pi;
ws=0.3*pi;
rp=1;
rs=15;
Fs=1;
[N,Wn]=buttord(wp,ws,rp,rs,'s');
[B,A]=butter(N,Wn,'s');
[bz,az]=bilinear(B,A,Fs);
[h,w]=freqz(bz,az);
plot(w/pi,20*log10(abs(h)));
title('数字低通滤波器的幅频特性’);
xlabel('频率/Hz');
ylabel('幅度/db');
实验结果:
b-150
度
幅-200
-350
00.10.20.30.40.50.60.70.80.91
频率/Hz
波器的下列性能指标:
通带截止频率p=2二2000rad/s,通带最大衰减:
3dB,
阻带截止频率"s=2二3000rad/s,阻带最小衰减>s=15dB,设采样频率
Fs=10000Hz。
假设该数字低通滤波器有一个输入信号
x(t)二sin2二0.5cos2二f2t,其中,f^1000Hz,f^4000Hz。
试将滤波器的
输出信号与输入信号进行比较。
理论分析:
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s')
其中,参数Wp和Ws分别是通带边界频率和阻带边界频率,Wp和Ws的单位是rad/s。
Rp和Rs分别为通带最大衰减和阻带最小衰减(dB)。
返回的参数N和Wn分别为滤波器的
阶数和3dB截止频率。
对于带通和带阻滤波器,Wp和Ws都是二维向量,向量的第一个元素对应低端的边界频率,第二个元素对应高端的边界频率。
[B,A]=butter(N,Wn,'s')
其中,N和Wn分别为滤波器的阶数和3dB截止频率。
利用此函数可以获得低通和带通滤波器系统函数的分子多项式(B)和分母多项式(A)的系数。
[bz,az]=impinvar(b,a,Fs)
可以实现用脉冲响应不变法将模拟滤波器转换为数字滤波器。
其中b和a分别是模拟滤
波器的系统函数H(s)的分子多项式和分母多项式的系数,Fs是脉冲响应不变法中的采样频
率,单位为Hz,如果Fs没有说明,其缺省值为1Hz。
运算的结果bz和az分别表示数字滤
波器的系统函数H(z)的分子多项式和分母多项式的系数。
利用filter函数计算数字滤波器的输出
源程序:
wp=2*pi*2000;
ws=2*pi*3000;
rp=3;
rs=15;
Fs=10000;
t=0:
0.0001:
0.1;
[N,Wn]=buttord(wp,ws,rp,rs,'s');
[B,A]=butter(N,Wn,'s');
[bz,az]=impinvar(B,A,Fs);
x=sin(2*pi*1000.*t)+0.5*cos(2*pi*4000.*t);
y=filter(bz,az,x);
subplot(2,1,1);
plot(t,x);
subplot(2,1,2);
plot(t,y);
实验结果:
1
I1
L
E
1I
1
■
n
-I
t]iI
n
I!
J\
!
1
.I
Jij
■卄
\\1,^Mij]*.M
ifb
ki\}Ilif
m'I,f护
ii
Uf3/
¥灯
yf.[ifl卩v|
|v'ti
N0
F|Ji
1
1
1
!
-
':
I
...
rir
『
4
rd
rr
1i
00.0010.0020.0030.0040.0050.0060.0070.0080.0090.01
\
I
GLL1
I
、\/
\f
J
\AAAA
\八
-1\
f
\\\\i
1
1
\f1,I\
\1-
/ip
\
\/1|Ii1t/
1
1
1
1\1\i
1i/\1\f
1
\f
1-
I\1\
i/
1l
\\
1J
11
IJ
I
I!
-
I\
1i
'r
Itf
(J\]■/\}
VV『■<
I1
"1\I
rr
1.5
1
0.5
0
-0.5
-1
1.5
1
0.5
0
-0.5
-1
00.0010.0020.0030.0040.0050.0060.0070.0080.0090.01