模拟低通巴特沃斯滤波器.docx
- 文档编号:4380660
- 上传时间:2022-12-01
- 格式:DOCX
- 页数:9
- 大小:220.53KB
模拟低通巴特沃斯滤波器.docx
《模拟低通巴特沃斯滤波器.docx》由会员分享,可在线阅读,更多相关《模拟低通巴特沃斯滤波器.docx(9页珍藏版)》请在冰豆网上搜索。
模拟低通巴特沃斯滤波器
《数字信号处理》课程设计报告
设计课题滤波器设计与实现
专业班级
姓名
学号
报告日期2012年11月
《数字信号处理》
课程设计任务书
题目
滤波器设计与实现
学生姓名
学号
专业班级
设
计
内
容
与
要
求
一、设计内容:
设计一个模拟低通巴特沃斯滤波器,技术指标:
通带截止频率
,通带最大衰减
;阻带起始频率
,阻带最小衰减
,画出其幅度谱和相位谱。
二、设计要求
1设计报告一律按照规定的格式,使用A4纸,格式、封面统一给出模版。
2报告内容
(1)设计题目及要求
(2)设计原理(包括滤波器工作原理、涉及到的matlab函数的说明)
(3)设计内容(设计思路,设计流程、仿真结果)
(4)设计总结(收获和体会)
(5)参考文献
(6)程序清单
起止时间
2012年12月3日至2011年12月11日
指导教师签名
2011年12月2日
系(教研室)主任签名
年月日
学生签名
年月日
目录
1课题描述
2设计原理
2.1滤波器的分类
2.2模拟滤波器的设计指标
2.3巴特沃斯低通模拟滤波器
2.3.1巴特沃斯低通模拟滤波器的设计原理
2.3.2巴特沃斯低通模拟滤波器的设计步骤
3设计内容
3.1用MATLAB编程实现
3.2设计结果分析
4总结
5参考文献
1课题描述
模拟滤波器的理论和设计方法已发展的相当成熟,且有多种典型的模拟滤波器供我们选择。
这些滤波器都有严格的设计公式,现成的曲线和图表供设计人员使用,而且所设计的系统函数都满足电路实现条件。
本课题使用MATLAB信号处理箱和BW(巴特沃斯)设计低通模拟滤波器。
2设计原理
2.1滤波器的分类
典型模拟滤波器有巴特沃斯滤波器,切比雪夫滤波器,椭圆滤波器,贝塞尔滤波器等。
选频型模拟滤波器按幅频特性可分成低通,高通,带通和带阻滤波器。
2.2模拟滤波器的设计指标
分别为:
通带截止频率=1000rad,阻带截止频率=5000rad,通带最大衰减=1dB,阻带最小衰减=30dB
N:
巴特沃斯低通滤波器的幅度平方函数|Ha(Jω)|用下式表示:
|Ha(jΩ)|=
公式中,N称为滤波器的阶数。
Wc:
3dB截止频率
B,A分别表示:
N阶巴特沃斯数字滤波器系统函数分子和分母多项式的系数向量
2.3巴特沃斯低通模拟滤波器
2.3.1巴特沃斯低通模拟滤波器的设计原理
巴特沃斯低通滤波器的幅度平方函数|Ha(jΩ)|2用下式表示:
|Ha(jΩ)|2=
公式中,N称为滤波器的阶数。
在Ω=0时,|Ha(jΩ)|=1;
Ω=Ωc时,|Ha(jΩ)|=1/
,Ωc是3dB截止频率。
在Ω=Ωc附近,随Ω加大,幅度迅速下降。
幅度特性与Ω和N的关,N越大,通带越平坦,过渡带越宽,过渡带与阻带幅度下降的速度越快,总的频响特性与理想低通滤波器的误差越小。
以S代替jΩ,将幅度平方函数|Ha(jΩ)|2写成s的函数:
Ha(S)Ha(-S)=
复变量S=α+jΩ,此公式表明幅度平方函数有2N个极点,极点Sk用下公式表示:
Sk=
公式中K=0,1,2…… ,2N-1。
2N个极点等间隔分布在半径为Ωc的圆上,间隔是
/Nrad。
为形成因果稳定的滤波器,2N个极点中只取S平面左半平面的N个极点构成Ha(S),而右半平面的N个极点构成Ha(-S)。
Ha(S)的表达式为 Ha(S)=
/
(S-Sk)
设N=3,极点有六个,它们分别是
S0=
S1=-ΩcS2=
S3=
S4=Ωc S5=
取S平面的左半平面的极点S0,S1,S2组成系统函数Ha(S),即Ha(S)=
由于不同的技术指标对应的边界频率和滤波器幅频特性不同,为使设计公式和图表统一将频率归一化。
巴特沃斯滤波器采用对3dB截止频率Ωc归一化,归一化后的系统函数为
令
=
,λ称为归一化频率,
称为归一化复变量,这样的巴特沃斯滤波器的归一化低通原型系统函数为
公式中,
称为归一化极点,用下公式表示,
k=0,1,…;,N-1
显然,Sk=Ωc*Pk
这样,只要根据技术指标求出阶数N,按照
公式求出N个极点,再按
得到归一化低通原型系统函数Ga(
),如果给定Ωc,再去归一化,即将
,代入
中,便可得到期望设计的系统函数Ha(S)。
所使用的MATLAB函数:
1.巴特沃斯滤波器阶数选择函数
(1)[N,wc]=buttord(wp,ws,αp,αs)
作用:
计算巴特沃斯数字滤波器的阶数N和3dB截止频率wc,wc为数字频率,单位rad。
说明:
调用参数wp,ws分别为数字滤波器的通带、阻带截止频率的归一化值,要求:
0≤wp≤1,0≤ws≤1。
αp,αs分别为通带最大衰减和组带最小衰减(dB)。
当ws≤wp时,为高通滤波器;当wp和ws为二元矢量时,为带通或带阻滤波器,这时wc也是二元向量。
(2)[N,Ωc]=buttord(Ωp,Ωs,αp,αs,‘s’)
作用:
计算巴特沃斯模拟滤波器的阶数N和3dB截止频率Ωc。
说明:
Ωp,Ωs,Ωc均为实际模拟角频率。
模拟频率f:
每秒经历多少个周期,单位Hz,即1/s,信号的真实频率,可用于模拟信号和数字信号;
模拟角频率Ω:
每秒经历多少弧度,单位rad/s,通常只于模拟信号;
数字频率w:
每个采样点间隔之间的弧度,单位rad,通常只用于数字信号。
关系:
Ω=2pi*f;w=Ω*T=2pi*f/F。
(F=1/Ts为采样频率,Ts为采样间隔)
2.完整巴特沃斯滤波器设计函数
(1)格式:
[b,a]=butter(N,wc,‘ftype’)
作用:
计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量b、a。
说明:
调用参数N和wc分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值,一般是调用buttord格式
(1)计算N和wc。
系数b、a是按照z-1的升幂排列。
(2)格式:
[B,A]=butter(N,Ωc,‘ftype’,‘s’)
作用:
计算巴特沃斯模拟滤波器系统函数的分子、分母多项式系数向量。
说明:
调用参数N和Ωc分别为巴特沃斯模拟滤波器的阶数和3dB截止频率(实际角频率),可调用buttord
(2)格式计算N和Ωc。
系数B、A按s的正降幂排列。
tfype为滤波器的类型:
◇ftype=high时,高通;Ωc只有1个值。
◇ftype=stop时,带阻;Ωc=[Ωcl,Ωcu],分别为带阻滤波器的通带3dB下截止频率和上截止频率。
◇ftype缺省时:
若Ωc只有1个值,则默认为低通;若Ωc有2个值,则默认为带通;其通带频率区间Ωcl<Ω<Ωcu。
3.求离散系统频响特性的函数freqz()
格式:
[H,w]=freqz(b,a,N)
说明:
b和a分别为离散系统的系统函数分子、分母多项式的系数向量,返回量H则包含了离散系统频响在0~pi范围内N个频率等分点的值(其中N为正整数),w则包含了范围内N个频率等分点。
调用默认的N时,其值是512。
可以先调用freqz()函数计算系统的频率响应,然后利用abs()和angle()函数及plot()函数,绘制出系统的频响曲线。
2.3.2巴特沃斯低通模拟滤波器的设计步骤
(1)确定滤波器的设计指标:
;
(2)运用函数
计算巴特沃斯低通滤波器的阶数N和归一化3db截止频率
;
(3)运用函数
求得低通滤波器的系统函数的分子、分母多项式形式;
(4)作图显示滤波器的幅频特性和相位特性。
3设计内容
3.1用MATLAB编程实现
t=1;
rp=1000;%带通截止频率
fp=fp*t;
rs=5000;%阻通截止频率
fp=fs*t;
rp=1;
rs=30;
wp=2*pi*fp;%利用Nyquist频率进行归一化
ws=2*pi*fs;
[n,wc]=buttord(wp,ws,rp,rs,’s’);%求数字滤波器的最小阶数和截止频率
[B,A]=butter(n,wc,'s');%设计高通数字滤波器系数b,a
fk=0:
14000/512:
14000;wk=2*pi*fk;
Hk=freqs(B,A,wk);
Subplot(2,1,1);plot(20*log10(abs(Hk)));gridon;%求系统频响特性,W为数字角频率,单位rad
xlabel('频率/Hz');ylabel('幅值');
Subplot(2,1,2);plot(180/pi*unwrap(angle(H)));grid;%绘出频率响应曲线
xlabel('频率/Hz');ylabel('相位');
3.2设计结果分析
Ωc是3dB截止频率。
在Ω=Ωc附近,随Ω加大,幅度迅速下降。
幅度特性与Ω和N的关系是幅度下降的速度与阶数N有关,N越大,通带越平坦,过渡带越窄,过渡带与阻带幅度下降的速度越快,总是频响特性与理想低通滤波器的误差越小。
综合可得:
本实验符合结论
4总结
本次课程实验中,我从刚开始的知识模糊终于在自己学习后对程序的编译有了一定的了解,不再完全都不了解。
所以说这次的实验对我的意义还是很大的。
首先要确定需要的函数,然后根据所给条件对数据进行处理后得到函数Ha(s),所以说只有自己也慢慢了解才能让它更好的使用在以后的生活中。
5参考文献
邹理和著:
《数字滤波器》,国防工业出版社,北京,1979。
A.AntoniouDigitalFilter:
AnalysisandDesign,McGraw-HillCo.,NewYork,1979.
T.S.Huang,TwoDimensional DigitalSignalProcessing,Springer-Verlag,Berlin,1981.
高西全丁玉美编著:
《数字信号处理》,西安电子科技大学出版社,2008.8
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 模拟 低通巴特沃斯 滤波器
![提示](https://static.bdocx.com/images/bang_tan.gif)