1、滤波器设计与实验1目录1 课题描述 .12 MATLAB 基础知识的介绍 .21.1 MATLAB 的介绍 .21.2 基本功能 .23 设计原理 .33.1 数字滤波器介绍 .33.2 IIR 数字滤波器设计原理 .53.3 巴特沃斯滤波器的介绍 .63.4 双线性变换法 .73.5 MATLAB 的应用 .114. 设计内容 114.1 用 MATLAB 编程实现 .154.2 设计结果分析 .175 总结 .206 参考文献 .21滤波器设计与实验21.课题描述数字滤波器,是指输入、输出均为数字信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数字滤波器分类的方法很多,可从以下几个方面进行分类:1根据功能分类,可以分为低通滤波器、高通滤波器、带通滤波器和带阻滤波器。2根据滤波器的单位脉冲响应分类,可以分为无限长脉冲响应数字滤波器和有限长脉冲响应数字滤波器,它们一般被称为IIR数字滤波器和FIR数字滤波器。3根据数字滤波器的实现方法分类,可以分为递归型数字滤波器RDF(Recursive Digital Filter)和非递归型数字滤波器NDF(N。
2、onrecursive Digital Filter2.MTTLAB 基 础 知 识 的 简 介2.1.1MATLAB 的 介 绍MATLAB 是 矩 阵 实 验 室 ( Matrix Laboratory) 的 简 称 , 是 美国 MathWorks 公 司 出 品 的 商 业 数 学 软 件 , 用 于 算 法 开 发 、 数 据 可视 化 、 数 据 分 析 以 及 数 值 计 算 的 高 级 技 术 计 算 语 言 和 交 互 式 环 境 ,主 要 包 括 MATLAB 和 Simulink 两 大 部 分 。2.2.基本功能MATLAB 是 由 美 国 mathworks 公 司 发 布 的 主 要 面 对 科 学 计 算 、可 视 化 以 及 交 互 式 程 序 设 计 的 高 科 技 计 算 环 境 。 它 将 数 值 分 析 、矩 阵 计 算 、 科 学 数 据 可 视 化 以 及 非 线 性 动 态 系 统 的 建 模 和 仿 真 等滤波器设计与实验3诸 多 强 大 功 能 集 成 在 一 个 易 于 使 用 的 视 窗 环 境 中 , 为 科 学 研 究 、工 程。
3、 设 计 以 及 必 须 进 行 有 效 数 值 计 算 的 众 多 科 学 领 域 提 供 了 一 种全 面 的 解 决 方 案 , 并 在 很 大 程 度 上 摆 脱 了 传 统 非 交 互 式 程 序 设 计语 言 ( 如 C、 Fortran) 的 编 辑 模 式 。MATLAB 和 Mathematica、 Maple 并 称 为 三 大 数 学 软 件 。 它 在数 学 类 科 技 应 用 软 件 中 在 数 值 计 算 方 面 首 屈 一 指 。 MATLAB 可 以进 行 矩 阵 运 算 、 绘 制 函 数 和 数 据 、 实 现 算 法 、 创 建 用 户 界 面 、 连 接 其 他 编 程 语 言 的 程 序 等 , 主 要 应 用 于 工 程 计 算 、 控 制 设 计 、 信号 处 理 与 通 讯 、 图 像 处 理 、 信 号 检 测 、 金 融 建 模 设 计 与 分 析 等 领域 。 MATLAB 的 基 本 数 据 单 位 是 矩 阵 , 它 的 指 令 表 达 式 与 数 学 、 工程 中 常 用 的 形 式 十 分 相 似 , 故 用 MATLAB 。
4、来 解 算 问 题 要 比 用C, FORTRAN 等 语 言 完 成 相 同 的 事 情 简 捷 得 多 , 并 且 MATLAB 也 吸收 了 像 Maple 等 软 件 的 优 点 ,使 MATLAB 成 为 一 个 强 大 的 数 学 软 件 。在 新 的 版 本 中 也 加 入 了 对 C, FORTRAN, C+ , JAVA 的 支 持 。 可以 直 接 调 用 ,用 户 也 可 以 将 自 己 编 写 的 实 用 程 序 导 入 到 MATLAB 函数 库 中 方 便 自 己 以 后 调 用 。3 设计原理31 数字滤波器介绍数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。数字滤波器和模拟滤波器有着相同的滤波概念,根据其频率响应特性可分为低通、高通、带通、带阻等类型,与模拟滤波器相比,数字滤波器除了具有数字信号滤波器设计与实验4处理的固有优点外,还。
5、有滤波精度高(与系统字长有关)、稳定性好(仅运行在 0 与 l 两个电平状态)、灵活性强等优点。时域离散系统的频域特性: ,其中 、分别是数字滤波器的输出序列和输入序列的频域特性(或称为频谱特性), 是数字滤波器的单位取样响应的频谱,又称为数字滤波器的频域响应。输入序列的频谱 经过滤波后,因此,只要按照输入信号频谱的特点和处理信号的目的, 适当选择 ,使得滤波后的 满足设计的要求,这就是数字滤波器的滤波原理。数字滤波器根据其冲激响应函数的时域特性,可分为两种,即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。IIR 数字滤波器的特征是,具有无限持续时间冲激响应,需要用递归模型来实现,其差分方程为:(1-1)系统函数为:(1-2)设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H(z),使其频率响应H(z)满足所希望得到的频域指标,即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。3.2 IIR数字滤波器设计原理滤波器设计与实验5频率交换冲激响应不变法双线性变换冲激响应不变法双线性变换IIR数字滤波器是一种离散时间系统,其系统函数为(1-。
6、3)假设MN,当MN时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数 和 ,它是数学上的一种逼近问题,即在规定意义上(通常采用最小均方误差准则)去逼近系统的特性。如果在S平面上去逼近,就得到模拟滤波器;如果在z平面上去逼近,就得到数字滤波器。设计高通、带通、带阻等数字滤波器通常可以归纳为如图所示的两种常用方法。方法 1: 首先设计一个模拟原型低通滤波器,然后通过频率变换成所需要的模拟高通、带通或带阻滤波器,最后再使用冲激不变法或双线性变换成相应的数字高通、带通或带阻滤波器。方法 2 :先设计一个模拟原型低通滤波器,然后采用冲激响应模拟原型低通数字高通、带通和带阻模拟高通、带通和带阻频率数字原型低通交换模拟原型低通模拟高通、带通和带阻滤波器设计与实验6不变法或双线性变换法将它转换成数字原型低通滤波器,最后通过频率变换把数字原型低通滤波器变换成所需要的数字高通、带通或带阻滤波器。方法一的缺点是,由于产生混叠是真,因此不能用冲激不变法来变换成高通或阻带滤波器,故一般采用第二种方法进行设计。本课程设计先构造一个切比雪夫模。
7、拟低通滤波器,然后将模拟低通滤波器转换成模拟带通滤波器,最后利用双线性变换将模拟带通滤波器转换成数字带通滤波器。3.3 巴特沃斯滤波器的介绍巴 特 沃 斯 滤 波 器 是 电 子 滤 波 器 的 一 种 。 巴 特 沃 斯 滤 波 器 的 特点 是 通 频 带 的 频 率 响 应 曲 线 最 平 滑 。 这 种 滤 波 器 最 先 由 英 国 工 程师 斯 替 芬 巴 特 沃 斯 ( Stephen Butterworth) 在 1930 年 发 表在 英 国 无 线 电 工 程 期 刊 的 一 篇 论 文 中 提 出 的 。巴 特 沃 斯 滤 波 器 的 特 点 是 通 频 带 内 的 频 率 响 应 曲 线 最大 限 度 平 坦 , 没 有 起 伏 , 而 在 阻 频 带 则 逐 渐 下 降 为 零 。 在 振 幅的 对 数 对 角 频 率 的 波 特 图 上 ,从 某 一 边 界 角 频 率 开 始 ,振 幅 随 着角 频 率 的 增 加 而 逐 步 减 少 ,趋 向 负 无 穷 大 。一 阶 巴 特 沃 斯 滤 波 器 的 衰 减 率 为 每 倍 频 6 分 贝 , 每 十 倍。
8、 频20 分 贝 。 二 阶 巴 特 沃 斯 滤 波 器 的 衰 减 率 为 每 倍 频 12 分 贝 、 三阶 巴 特 沃 斯 滤 波 器 的 衰 减 率 为 每 倍 频 18 分 贝 、 如 此 类 推 。 巴 特沃 斯 滤 波 器 的 振 幅 对 角 频 率 单 调 下 降 , 并 且 也 是 唯 一 的 无 论 阶 数 ,振 幅 对 角 频 率 曲 线 都 保 持 同 样 的 形 状 的 滤 波 器 。 只 不 过 滤 波 器滤波器设计与实验7阶 数 越 高 , 在 阻 频 带 振 幅 衰 减 速 度 越 快 。 其 他 滤 波 器 高 阶 的 振 幅对 角 频 率 图 和 低 级 数 的 振 幅 对 角 。图 二 巴 特 沃 斯 与 切 比 雪 夫 1, 切 比 雪 夫 2 和 椭 圆 滤 波 器 的滤 波 效 果 比 较3.4 双线性变换法目的:将模拟带通滤波器转换成数字高通滤波器为了克服冲激响应法可能产生的频率响应的混叠失真,这是因为从 S 平面到平面是多值的映射关系所造成的。为了克服这一缺点,可以采用非线性频率压缩方法,将整个频率轴上的频率范围压缩到-/T/T 之间,再用 z=esT转换到 Z 平面上。也就是说,第一步先将整个 S 平面压缩映射到 S1平面的-/T/T 一条横带里;第二步再通过标准变换关系 z=es1T将此横带变换到整个 Z 平面上去。这样就使 S 平面与 Z 平面建立了一一对应的单值关系,消除了多值变换性,也就消除了频谱混叠现象,映射关系如图 1-3o 1 1ZjImzRez/ Tj11/ TS1Sj o o滤波器设计与实验8图 1-3 双线性变换的映射关系为了将 S 平面的整个虚轴 j 压缩到 S1 平面 j1 轴上的-/T到 /T 段上,可以通过以下的正切变换实现 (1-13)式中,T 仍是采样间隔。当 1 由-/T 经过 0 变化到 /T 时, 由-经过 0 变化到+,也即映射了整个 j 轴。将式(1-9)写成(1-14)将此关系解析延拓到整个 S 平。
《巴特沃夫低通滤波器》由会员平***分享,可在线阅读,更多相关《巴特沃夫低通滤波器》请在金锄头文库上搜索。