巴特沃斯低通滤波器 matlab,利用matlab设计巴特沃斯低通滤波器
利用matlab设计巴特沃斯低通滤波器
三峡大学 课 程 设 计 报 告 专业班级 20091421 课 程 数字信号处理课程设计 学 号 2009142116 学生姓名 姜祥奔 指导教师 王露 2012年 5 月 平时成绩(20%) 报告成绩(40%) 答辩成绩(40%) 总成绩数字信号处理课程设计实验一:用双线性变换法和脉冲响应不变法设计IIR 数字滤波 器 采用双线性变换法设计一个巴特沃斯数字低通滤波器,要求:通带截止频 率 100Hz,阻带截止频率 200Hz,通带衰减指标 Rp 小于 2dB,阻带衰减 Rs 大 于 15dB,滤波器采样频率 Fs=500Hz。绘制频率响应曲线。 理论部分:(原理及设计过程) 第一步:论 的关系及数字域性能的公式表示。模拟频率与数字频率之间为 和 f 线性关系, 为抽样周期,故 T f T fT T s , 10 1 , 2 4 4 . 0 500 1 100 2 1 c c kHz f 对应于 修正后) 2 tan( 2 1 c T 8 . 0 500 1 200 2 1 st st kHz f 对应于 修正后) 2 tan( 2 1 st T 按衰减的定义 2 ) ( ) ( log 20 4 . 0 0 10 j j e H e H 15 ) ( ) ( log 20 8 . 0 0 10 j j e H e H 设 处频率响应幅度归一化为 1,即 ,则上两式变成 0 1 ) ( 0 j e H………………(1) 2 ) ( log 20 4 . 0 10 j e H ………………(2) 15 ) ( log 20 8 . 0 10 j e H 这就是数字滤波器的性能指标的表达式。 2 下面把数字低通滤波器的性能要求转变为“样本”模拟低通滤波器的性 能要求。由 ,按修正式 ,设没有混叠效应(即混 T ) ( ) ( T j H e H a j 叠效应设计完成后再进行校验) 则有……………(3) ), ( ) ( ) ( j H T j H e H a a j 利用(3)式,由(1) 、 (2)式可写出模拟低通滤波器的指标为…….(4) 2 ) 10 2 ( log 20 ) 4 . 0 ( log 20 2 10 10 j H T j H a a…(5) 15 ) 10 4 ( log 20 ) 8 . 0 ( log 20 2 10 10 j H T j H a a 3 计算“样本”模拟低通滤波器所需的阶数 N 及 3dB 截止频率 。巴特沃 C 思低通滤波器的幅度平方函数是N C a j H 2 2 ) ( 1 1 ) ( 以分贝形式表示上式,即……………(6) N C a j H 2 10 10 ) ( 1 log 10 ) ( log 20 把求出的性能指标关系(4)式、(5)式代入(6)式得 2 ) ( 1 log 10 2 1 10 N C 15 ) ( 1 log 10 2 2 10 N C 先用等号来满足指标,可得= N C 2 1 ) ( 1 2 . 0 10 = N C 2 2 ) ( 1 5 . 1 10解此两方程,得 N=1.3709, N 是滤波器阶次,必须取整数,为了满足或超过给 定指标,故应选取比求出的 N 大一点的整数,故取 N=2,代入通带条件(6- 118a)式,得 =1.038 ,即 C 3 10 N=2, =1.038 ……………(7) C 3 10 显然,利用(7)式的 N 和 ,阻带指标正好满足,而通带指标则低于 C 2dB,这正好对减少冲激响应不变法造成的频谱混叠效应是有利的。 程序部分: wp=100*2*pi; %数字滤波器的通带截止频率 ws=200*2*pi; %数字滤波器的阻带截止频率 Rp=2;As=15; % 输入滤波器的通阻带衰减指标 %转换为模拟滤 波器指标 Fs=500;T=1/Fs; wp1=wp*T; ws1=ws*T; Omgp=(2/T)*tan(wp1/2); %原型通带频率预修正 Omgs=(2/T)*tan(ws1/2); %原型通阻频率预修正 %模拟原型滤波器计算 [n, Omgc]=buttord(Omgp, Omgs,Rp,As, s ); %计算阶数 n和截止频率 [z0,p0,k0]=buttap(n); %归一化原型设计 [ba,aa]=zp2tf(z0,p0,k0); %把滤波器零极点转化为传递函数模型 [ba1,aa1]=lp2lp(ba,aa, Omgc); %变换为模拟低通滤波器系数 b,a %或者[ba1,aa1]=butter(n, Omgc,’s’); %用双线性变换法计算数字滤波器系数 [bd,ad]=bilinear(ba1,aa1,Fs) ; %用双线性变换法求数字滤波器系数 b,a %求数字系统的频率特性 [H,w]=freqz(bd,ad); plot(w*Fs/(2*pi),20*log10(abs(H))); ylabel( 增益/dB );xlabel( 频率(Hz) );title( 数字滤波器幅度响应 );grid on;数字滤波器的频率响应图: 实验二:用窗函数法设计FIR 数字滤波器 选择合适的窗函数设计一个 FIR 数字低通滤波器,要求:带通截止频率为 Wp=0.2π,阻带截止频率 Ws=0.3π,阻带衰减不小于 40dB,通带衰减不大于 3dB。描绘滤波器的幅频响应曲线。 理论部分:(原理及设计过程) 1.数字频率。 通带截止频率: 2 . 0 p w 阻带截止频率: 3 . 0 s w 阻带衰减: dB 40 2 2.求 。设 为理想线性相位滤波器 ) (n h d ) ( jw d e H
巴特沃斯低通滤波器 matlab,利用matlab设计巴特沃斯低通滤波器相关推荐
- 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序
环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. YID:856064261479453 ...
- 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强
环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. 现有:6960642614794538 ...
- 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数
环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. 现有:6960642614794538 ...
- 物联网控制原理与技术--基于Matlab/利用MATLAB进行频域分析(伯德图)的应用(超详细/设计/实验/作业/练习)
目录 课程名:物联网控制原理与技术 内容/作用:设计/实验/作业/练习 学习:利用MATLAB进行频域分析(伯德图) 一.前言 二.环境与设备 三.原理 四.内容 五.结果与分析 课程名:物联网控制原 ...
- 单闭环调速仿真matlab,利用Matlab仿真平台设计单闭环直流调速系统
内容简介: 毕业论文 利用Matlab仿真平台设计单闭环直流调速系统 共21页,3989字. 目 录 一.摘要 --------------------------2 二.总体方案设计 ------- ...
- 热传递 matlab,利用matlab程序解决热传导问题
利用matlab程序解决热传导问题 1哈佛大学能源与环境学院课程作业报告作业名称:传热学大作业--利用 matlab 程序解决热传导问题院系:能源与环境学院专业:建筑环境与设备工程学号:5201314 ...
- 神经网络人口预测matlab,利用MATLAB工具模拟实现人口预测实例分析.doc
摘要:人口预测对于国民经济计划的制定和社会战略目标的决策具有重要参考价值.本文主要利用MATLAB软件对中国1974-2015年的人口总数数据进行统计分析,建立BP神经网络预测模型.把它与一元非线性回 ...
- 二维静电场有限元matlab,利用MATLAB计算电磁场有关分布
<利用MATLAB计算电磁场有关分布>由会员分享,可在线阅读,更多相关<利用MATLAB计算电磁场有关分布(17页珍藏版)>请在人人文库网上搜索. 1.电磁场实验报告实验一 模 ...
- 投入产出实例matlab,利用MATLAB进行基础的投入产出表分析实验报告.doc
利用MATLAB进行基础的投入产出表分析实验报告 学 生 实 验 报 告 学 院: XX学院 课程名称:宏观经济数量分析方法与模型 专业班级: 姓 名: 学 号: 学生实验报告 学生姓名学号同组人实验 ...
- 多普勒谱matlab,利用MATLAB仿真多普勒效应.doc
利用MATLAB仿真多普勒效应 利用MATLAB仿真多普勒效应 某某某 摘 要:分析多普勒效应特性,建立数学模型,利用MATLAB软件对其进行仿真试验,进行定量分析,根据仿真试验结果绘制出听者接收到的 ...
最新文章
- 单片机和嵌入式哪个好?单片机会被嵌入式取代吗?
- morphologyEx 形态学
- 网络爬虫:URL去重策略之布隆过滤器(BloomFilter)的使用
- 大型网站架构 - 1.架构的演变过程
- 向量余弦值python_向量/矩阵的余弦值打印(元素明智的操作) 使用Python的线性代数
- Hibernate中的三种数据持久状态和缓存机制
- webservice使用EF生成的model序列化问题
- Android在WindowManagerService和ActivityManagerService中的Token
- Node编写API接口,ajax实现增删改查
- 回溯算法讲解--适用于leetcode绝大多数回溯题目
- 关于sts与eclipse版本对应的插件图
- 基于Web的MIS系统环境配置和相关含义(拷贝版)
- 求两个数的最大公约数和最小公倍数
- 基于SSM的设备租赁管理系统
- 因为一条SQL,我差点被祭天......,我太难了!
- vue返回上一页并不刷新
- android 短信数据库在哪个文件夹,安卓手机的短信备份后存到了哪个文件夹下
- ps如何快速消除黑眼圈或者眼袋
- Unity景深效果解析
- 数据可视化Matplotlib-中
热门文章
- MySQL数据库安装超级详细教程
- 深度学习模型迁移到VisualC++ demo
- 耗时6个月,整理了30款免费高评分软件,完爆付费
- java 反编译 报错_java反编译后再编译成class时提示缺少包
- sql注入之——sqlmap教程
- Kotlin Compose 标准组件
- 值得收藏的30道Python练手题(附详细答案)
- Python100经典练习题,附答案
- J2ME 手机模拟器插件 eclipseme.feature_1.7.9_site.zip
- 2022计算机Java二级考试四十五套题真题【收藏版】(一周裸考计划)