简单低通滤波器设计及matlab仿真.docx
《简单低通滤波器设计及matlab仿真.docx》由会员分享,可在线阅读,更多相关《简单低通滤波器设计及matlab仿真.docx(9页珍藏版)》请在冰豆网上搜索。
简单低通滤波器设计及matlab仿真
东北大学
研究生考试试卷
考试科目:
课程编号:
阅卷人:
考试日期:
姓名:
XI学号:
注意事项
1.考前研究生将上述项目填写清楚.
2.字迹要清楚,保持卷面清洁.
3.交卷时请将本试卷和题签一起上交.
4.课程考试后二周内授课教师完成评卷工作,公共课成绩单与试卷交研究生院培养办公室
专业课成绩单与试卷交各学院,各学院把成绩单交研究生院培养办公室.
东北大学研究生院培养办公室
数字滤波器设计
技术指标:
fp=
频率:
Fs二200Hz
目标:
1、根据性能指标设计一个巴特沃斯低通模拟滤波器。
2、通过双线性变换将该模拟滤波器转变为数字滤波器。
原理:
一、模拟滤波器设计
每一个滤波器的频率范围将直接取决于应用目的,因此必然是千差万别。
为了使设计规范化,需要将滤波器的频率参数作归一化处理。
设所给的实际频率为(或,),归一化后的频率为,对低通模拟滤波器令=>,则…。
令归一化复数变量为卩,‘则一o所以巴特沃思模拟低通滤波器的设计可按以下三个步骤来进行。
(1)将实际频率规一化
(2)求Qc和N
<2mp/10
NIgl1o°V/r°11,gS
这样Qc和N可求。
根据滤波器设计要求=3dB,则C=1,这样巴特沃思滤波器的设计就只剩一个
参数N,这时
G(j)211
(3)确定强因为‘根据上面公式有
…zN2N
N2N
1
(1)"P*0”
pkexp(j2kN12)k
k2N=4
,_「2,•••,2N
这样可得
Gk(p)
G,:
P(PPk)(PPNIk)p22poo.(2kN1)1
2N
求得“•后,用••代替变量P,即得实际需要得-o
二、双线性变换法
双线性变换法是将•平面压缩变换到某一中介S1平面的一条横带里,再通过标准变换关系一将此带变换到整个Z平面上去,这样就使S平面与Z平面之间建立一一对应的单值关系,消除了多值变换性。
为了将S平面的,轴压缩到S1平面的’轴上的•到一段上,可以通过以下的正切变换来实现:
21
tan(iT)
1i
这样当1由•「经0变化到・4寸,由经过0变化到,也映射到了整个,轴。
将这个关系延拓到整个•平面和s1平面,则可以得到
stan(s1)sit
•w
再将S1平面通过标准变换关系映射到Z平面,即令•-…得到
同样对Z求解,得到
2
zt28s
T
这样的变换叫做双线性变换。
为了验证这种映射具有s平面的虚轴映射到z平面单位圆上的特性,考虑・,,ze得
21e■
Tie,
■tan
T2
除了使s平面的虚轴映射到单位圆上之外,s平面的左半部分映射到单位圆的内部,s平面的右半部分映射到单位圆的外部。
设计过程一、计算滤波器阶数N和截止频率Oc
根据公式:
S
lgk
代入数据,计算可得:
ps二0.5
K
ps二0.99885
N=3.3151
所以取N二4,即滤波器为四阶滤波器
计算可得
计算3dB截止频率,根据公
Oc=112.6096
二、计算系统函数
已知四阶归一化低通巴特沃斯模拟滤波器系统函数为:
“a・432
3s42.6131s33.4142s22.6131s1
将S用S/Qc替代,求出系统函数:
Ha(s)二b/(s+a3s+a2s+ais+ao)其
中:
b=2.5063*10^(11)
33=1.8489*10"(5)
32=1.7092*10"(6)
ai=9.2560*10^(8)
ao二2.5063*10^(11
)
三、程序实现模拟滤波器(MATLAB)
cIear;cIoseaII
%滤波器指标
%计算阶数N和3dB截至频率fc
%设计低通巴特沃斯模拟滤波器
fp=100;fs二200;Rp二3;As二20
t
[N,fc]二buttord(fp,fs,Rp,
As,'s')[B,A]=butter(N,fc,'s');
pIot(f,20*Iog10(abs(hf)),'r')grid;xIabeI('频率Hz');ylabel('幅度dB')titleC模拟低通滤波器’);axis([0,250,-25,5])Iine([0,250],[-3,-3]);line([100,100],[-25,5]);Iine([0,250],[-20,-20]);
line([200,200],[-25,5]);响应曲线如下
四、通过双线性法变换将模拟滤波器转变为数字滤波器
首先根据公式
H(z)H(s)21Z.
'T1「将H(s)转换为H(z)计算可得:
H(z)分子系数为:
Mz二0.07860.11790.07860.0197
0.0197
-1.72710.08711.8682-0.9138
H(z)分母系数为:
Nz=23Y
0•叩獅0786Z~'+O・行79Z"+0.0786Z‘+0・0197Z"4
1-1.7271r1+0.0871Z_2+1.8682Z_3-O.9138Z-4
五、程序实现数字滤波器(MATLAB)
T=0.005;M二2.5063*10^(11);
N=[1,1.8489*10^(5),1.7092*10^(6),9.2560*10^(8),2.5063*10
"(11)][Mz,Nz]^bilinear(M,NJ/T);%对模拟滤波器双线性变换
Mz,Nz
[hi,w1]=freqz(M乙Nz);%数字滤波器的幅频响应
figure
pIot(w1/pi,20*Iog10(abs(hi)),'r');grid;
xlabelCCO/n');ylabel('幅度(dB)');titleC数字低通滤波器
');axis([0,1.1,-160,20])
响应曲线如下
六、设计级联型滤波器
根据
0.0197+0.0786Z_1+0.1179Z2+0.0786r3+0.0197V=-1
-2-3-4
1-1.7271Z_1+0.0871Z为<8682Z_3-0.9138Z-4
编写matlab程序求级联型滤波器系数。
程序如下:
a二[0.0197,0.0786,0.1179,0.0786,0.0197];
b二[1.0000,-1.7271,0.0871,1.8682,-0.9138];
[sos,g]二tf2sos(a,b)%求级联型结构系数
结果如下:
g=0.0197
所以可由级联型结构系数写出系统函数:
H(z)=0.0197
-1-21+1.6348Z-
'+0.6942Z
级联型结构图如下图所示:
0.0197
Y(n)
-1
-0.3905
0.6025
Z-1
2.355
1.4405
-1
Z-1
2.11761.6348
-1.5166
Z-1
0.6942
结果分析
1、模拟滤波器性能分析
观察响应曲线,在通带边界频率100Hz处,幅度最大衰减为3dB;在阻带边界频率200Hz处,幅度衰减达到20dBo因此模拟滤波器设计符合要求。
2、双线性变换所得数字低通滤波器分析
双线性变换优点:
是频率坐标变换是线性的,即T,如果不考虑频率混叠现象,用这种方法设计的数字滤波器会很好的重现原模拟滤波器的频率特性。
另一个优点是数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应,时域逼近性好。
双线性变换法的缺点:
数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变。
例如,一个模拟微分器,它的幅度与频率是线性关系,但通过双线性变换后,就不可能得到数字微分器。
另外,一个线性相位的模拟滤波器经双线性变换后,滤波器就不再有线性相位特性。
虽然双线性变换有这样的缺点,但它目前仍是使用得最普遍、最有成效的一种设计工具。
这是因为大多数滤波器都具有分段常数的频响特性,如低通、高通、带通和带阻等,它们在通带内要求逼近一个衰减为零的常数特性,在阻带部分要求逼近一个衰减为X的常数特性,这种特性的滤波器通过双线性变换后,虽然频率发生了非线性变化,但其幅频特性仍保持分段常数的特性。
双线性变换比脉冲响应法的设计计算更直接和简单。
由于S与Z之间的简单代数关系,所以从模拟传递函数可直接通过代数置换得到数字滤波器的传递函数。
置换过程:
H⑵■弘(s)2W
频响:
这些都比脉冲响应不变法的部分分式分解便捷得多