matlab模拟巴特沃斯滤波器设计,巴特沃斯滤波器matlab实现
描述
巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦,没有起伏,而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。
一阶巴特沃斯滤波器的衰减率为每倍频6分贝,每十倍频20分贝。二阶巴特沃斯滤波器的衰减率为每倍频12分贝、三阶巴特沃斯滤波器的衰减率为每倍频18分贝、如此类推。巴特沃斯滤波器的振幅对角频率单调下降,并且也是唯一的无论阶数,振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高,在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。
设计步骤
如设计一个数字低通滤波器,其技术指标为:
通带临界频率fp ,通带内衰减小于rp;
阻带临界频率fs,阻带内衰减大于αs;采样频率为FS
1、将指标变为角频率 wp=fp*2*pi;ws= fs*2*pi;
2、将数字滤波器的频率指标{Wk}由wk=(2/T)tan(Wk/2)转换为模拟滤波器的频率指标{wk},由于是用双线性不变法设计,故先采取预畸变。
3、将高通指标转换为低通指标,进而设计高通的s域模型
4、归一化处理
由以上三式计算出N,查表可得模拟低通滤波器的阶数,从而由下式确定模拟高通滤波器的参数。
仿真程序的设计与调试
数字域指标变换成模拟域指标
其程序为:
fp = 400 fs= 300;
Rp = 1; Rs = 20;
wp =fp*2*pi;
ws =fs*2*pi;
FS=1000;T=1/FS;
程序执行结果为:wp=2.5133e+003 ws=1.8850e+003 与实际计算结果相符。
数字域频率进行预畸变
其程序为:
wp2=2*tan(Wp/2)/T;
ws2=2*tan(Ws/2)/T;
经过预畸变,可以发现频率变为: wp2= 6.1554e+003
ws2= 2.7528e+003
模拟滤波器的设计
其程序为
%设计模拟滤波器
[N,Wn] = buttord(wp2,ws2,Rp,Rs,‘s’)
武汉理工大学《数字信号处理》报告
9
[z,p,k]=buttap(N); %创建Buttord低通滤波器原型
[Bap,Aap]=zp2tf(z,p,k); %由零极点转换为传递函数的形式
figure(1) freqs(Bap,Aap); %模拟低通滤波器的频率响应
title(‘模拟滤波器(低通原型)的频率响应’)
[Bbs,Abs]=lp2hp(Bap,Aap,Wn); %模拟低通变高通
figure(2)
freqs(Bbs,Abs);
title(‘模拟滤波器的频率响应’)
程序执行后可以发现其频率响应为: N=4,其波形如下图
模拟滤波器的频率响应
由上图分析可得:其符合高通的一般特征,与预期的效果一样。 而在此条件下,Butterworth滤波器低通原型的波形如下图。
模拟滤波器(低通原型)的频率响应
在设计的过程中,涉及一个频率变换的问题,即将模拟低通原型变为高通,其函数及用法如下:
[b,a]=lp2hp(Bap,Aap,Wn);
功能:把模拟滤波器原型转换成截至频率为 Wn 的高通滤波器。 其中,Bap,Aap分别为低通传递函数的分子向量和分母向量;
b,a分别为高通传递函数的分子向量和分母向量。
模拟滤波器变成数字滤波器
其程序为:
[Bbz,Abz]=bilinear(Bbs,Abs,FS); %用双线性变换法设计数字滤波器 freqz(Bbz,Abz,512,FS);
程序运行的结果为:如下图
数字滤波器的频率响应
由于使用的是双线性不变法设计的,其相位为非线性。此处主要是基于要获得严格的频率响应,以及较准确地控制截止频率的位置,故画出了详细的幅频响应。(如下图)
详细的幅频响应
分析该图可知其在0.6(即300Hz)处的衰减为40dB,而在0.8(即400Hz)处的衰减极小,应小于1dB。由此可见,此设计符合要求设计的参数。
而在调试的过程中发现:通带衰减越小,滤波器的性能越好 阻带衰减越大,滤波器的性能越好 其曲线也越陡峭,选择性越好,当然所用的滤波器阶数也越高。
当阻带衰减变为40dB(之前为20dB),通带不变时,其波形如下图。对比上图可知,其在阻带临界频率处衰减变为了40dB,曲线变陡峭了。
详细的幅频响应(阻带衰减为40dB)
当通带变为5dB时,阻带不变时,其波形如下图。对比图3-3可知,其在通带处的衰减变为了5dB,曲线平滑了一些。
详细的幅频响应(通带衰减为5dB)
理论计算数字滤波器的仿真
wp=0.8*pi;
ws=0.6*pi;
OmegaP=2*1000*tan(wp/2);
OmegaS=2*1000*tan(ws/2);
lamdas=OmegaP/OmegaS;
N=0.5*log10((10.^(20/10)-1)/(10.^(1/10)-1))/log10(lamdas);
%笔算的结果为N=3.6947;故取N=4 %
此处为计算高通的传递函数 Wn= 4.8890e+003 az=[0 0 0 0 1];
bz=[1 2.613 3.414,2.613,1]; [Bbs,Abs]=lp2hp(az,bz,Wn) %用双线性不变法处理
[Bbz,Abz]=bilinear(Bbs,Abs,1000); 其运行结果为:N=3.6947;图形如下图
理论计算的滤波器的幅频响应
综上所述,本滤波器以四阶即实现了预期的设计目标:采样频率为1000Hz,通带临界频率fp =400Hz,通带内衰减小于1dB(αp=1);阻带临界频率fs=300Hz,阻带内衰减大于20dB(αs=25),其在通带内的性能更好。
打开APP精彩内容
点击阅读全文
matlab模拟巴特沃斯滤波器设计,巴特沃斯滤波器matlab实现相关推荐
- [Matlab]FIR滤波器设计:(FIR滤波器的结构)
[Matlab]FIR滤波器设计:(FIR滤波器的结构) FIR(Finite Impulse Response)滤波器:有限长单位冲激响应滤波器,又称为非递归型滤波器,是一种在数字信号领域应用非常广 ...
- matlab滤波器设计工具箱带阻滤波器,用matlab信号处理工具箱进行fir滤波器设计的三种方法...
用matlab信号处理工具箱进行fir滤波器设计的三种方法 摘 要 介绍了利用 MATLAB 信号处理工具箱进行 FIR 滤波器设计的三种方法:程序设计法. FDATool 设计法和 SPTool 设 ...
- 音频均衡器 matlab code,急求高手 设计声音均衡器 滤波器
本帖最后由 lzhwzs 于 2010-12-6 14:02 编辑 注意: 必须上交MATLAB编码 这次的任务包括设计一个简单的双频带的声音均衡器,然后在Matlab中实现. 你已经被提供了一个Ma ...
- matlab中滤波器设计与实现,基于MATLAB的IIR滤波器设计与实现
描述 由于电力系统中非线性电子元件的大量使用,使得谐波污染问题日益严重,当谐波含量超过一定限度时就可能对电网和用户造成极大的危害,且增加线路损耗,降低线路传输能力,干扰通信信号等.因此,应该积极寻找一 ...
- matlab iir滤波器设计 实验报告,基于MATLAB的IIR滤波器设计与实现
由于电力系统中非线性电子元件的大量使用,使得谐波污染问题日益严重,当谐波含量超过一定限度时就可能对电网和用户造成极大的危害,且增加线路损耗,降低线路传输能力,干扰通信信号等.因此,应该积极寻找一种治理 ...
- 基于matlab的模糊控制器的设计与仿真,基于MATLAB的模糊控制洗衣机的设与仿真详解.doc...
基于MATLAB的模糊控制洗衣机的设计与 仿真 卫瑶瑶,王胜红 (南京农业大学工学院,210031) 摘要:根据模糊控制的原理对传统洗衣机进行改造,设计了模糊控制系统.通过MATLAB仿真,采用取最大 ...
- 滤波器设计软件_滤波器设计——电路仿真软件的滤波器参数提取(下)
本文章仅代表个人观点,如有错误缺漏,欢迎指正. 接上一节的内容. 按照上一节的方法,用仿真得到的频率和耦合曲线,修正电磁仿真软件的尺寸,经过5次迭代,就可以满足回波和带外抑制的要求.迭代过程如下(继续 ...
- matlab 模拟风,三维空间相关风场的计算机模拟及Matlab程序实现
第 14 卷第 2 期2008 年 6 月 空 间 结 构 SPATIAL STRUCTURES Vol. 14 No. 2 Jun. 2008 收稿日期:200610208. 作者简介:刘文洋(19 ...
- 【光学】基于matlab模拟涡旋光与球面波的干涉含Matlab源码
1 简介 基于涡旋光与球面波的干涉原理,提出一种物体微位移的光学测量方法.改进马赫泽德干涉光路,其中一束光照射至空间光调制器产生涡旋光束作为参考光,另一束光经透镜变为球面波后照射至物体上,两束光干涉后 ...
最新文章
- 图像的稀疏表示——ScSPM和LLC的总结
- 领域驱动设计,让程序员心中有码(七)
- aspx 获取html值,aspx解析之后的html代码获取
- c语言mysql 学生信息管理系统_学生信息管理系统学生时代小作品源码(C语言版)...
- 数据库视图作用?什么时候用视图?
- 项目推进计划表_项目为啥延期之一:计划问题
- 中心极限定理_中心极限定理和Python图解
- Hello!CDSN!
- 怎样修改用户的计算机配置文件,计算机本地用户配置文件如何迁移至域账户
- html5怎么给标题居中,html5标题居中 Html5如何使div里面文字在水平垂直居中对齐...
- POJ 1129 Channel Allocation(四色定理)
- Windows文件系统过滤驱动开发教程-楚狂人
- BookCollectionSoftware - 图书管理(知识点糅合::类,对象,抽象类,接口,封装,继承,多态,顺序表)- java - 细节狂魔
- 2017年最后一篇文章:微信小程序游戏跳一跳刷榜原理解析!
- 【文智背后的奥秘】系列篇——文本聚类系统
- rust 飞天指令_rust腐蚀游戏指令 游戏命令大全
- 斐讯dns服务器未响应,斐讯k2无线网络初始密码是?
- H.266/VVC相关技术学习笔记21:帧间预测中五种Merge模式的熵编码方式
- Python实战技巧(9)Python发送邮件
- 深度估计 DenseDepth 笔记
热门文章
- java api font_java – PDFBox API:如何更改字体以处理Acro...
- LaTex 之 爱思唯尔(ELSEVIER)模板使用技巧
- RGB,CMYK,HSB各种颜色表示的转换
- 运行中出现键盘backspace键不起作用的情况
- linux中LVM是什么类型,Linux_Linux中的LVM详解,什么是LVM?什么是LVM?LVM(Logi - phpStudy...
- 分享几个FTP帐号-来自www.lgm.hk
- 如何成为一个linux开发人员
- 写自己的故事II + 2019年年度总结
- ZESH | 2023春夏系列 尽享静谧诗意的东方之美
- 《狂暴幻想》3.18上线链游玩家|Q版养成、挂机冒险