数字信号处理 第五章(数字滤波器结构)

参考教程:​​​​​​https://www.bilibili.com/video/BV1qTCDYpEFr?spm_id_from=333.788.videopod.sections&vd_source=8f8a7bd7765d52551c498d7eaed8acd5

一、数字滤波器概述

1、数字滤波器的概念

(1)数字滤波器(DF)作用是对输入信号进行滤波,换句话说,DF是由差分方程描述的一类特殊的离散时间系统。

(2)数字滤波器把输入序列通过一定的运算变换成输出序列,不同的运算处理方法决定了滤波器实现结构的不同。

2、数字滤波器的工作原理

3、数字滤波器结构的表示方法

(1)数字滤波器表示所需的元素:

①表示方法:可选方框图表示法、信号流图表示法。

②运算类型:相加、乘以常数、延时。

③基本运算单元:加法器、乘法器、单位延时。

(2)可通过流图或方框图看出系统的运算步骤和运算结构,为方便起见,下面统一用流图分析数字滤波器的结构。

(3)举例:

4、研究数字滤波器结构的意义

(1)滤波器的基本特性(如有限长冲激响应FIR与无限长冲激响应IIR)决定了结构上有不同的特点。

(2)不同结构所需的存储单元及乘法次数不同,前者影响复杂性,后者影响运算速度。

(3)有限精度(有限字长)实现情况下,不同运算结构的误差及稳定性不同。

(4)好的滤波器结构应该易于控制滤波器性能,适合于模块化实现,便于时分复用。

二、IIR数字滤波器结构

1、直接I型及直接II型结构

(1)IIR DF的特点:

(2)直接I型IIR数字滤波器:

(3)直接II型IIR数字滤波器:

        直接I型IIR数字滤波器使用了较多延时单元,编程时需要分配内存去存储它们,这无疑是增加了空间复杂度,对此可考虑将延时单元合并,这样就由直接I型IIR数字滤波器转化为直接II型IIR数字滤波器

        直接II型IIR数字滤波器的结构为两个网络级联,第一个横向的M节延时网络实现零点,第二个反馈的N节延时网络实现极点,共需N(N≥M)级延时单元,所需延时单元数最少,又称为典范型

        不过,直接II型IIR数字滤波器一样具有直接型实现的一般缺点

(4)举例:

①例1:

②例2:

2、级联型结构

(1)问题的引出与解决:

(2)级联型IIR数字滤波器系统函数的因式分解:

(3)级联型IIR数字滤波器结构的特点:

①每个二阶节系数单独控制一对零点或一对极点,有利于控制频率响应,滤波器调整方便。

②硬件实现中,可采用一个二阶环分时复用;软件实现中,可以调用同一子函数。

③对有限字长的敏感程度比直接型低。需要注意的是,不同的配对方式具有不同的有限字长敏感程度。

3、并联型结构

(1)问题的引出与解决:

(2)系统函数的部分分式展开:

(3)并联型IIR数字滤波器结构的特点:

①可以单独调整极点位置,但不能像级联那样直接控制零点,故不能用于对零点位置精度要求高的滤波器,例如陷波器和窄带带阻滤波器。

②各基本网络是并联的,产生的运算误差互不影响,不像直接型和级联型有误差累积,故并联型结构的运算误差较小。

③硬件实现时,各子系统可以并行实现,运算速度快。

(4)举例:

4、转置型结构

(1)如果将一个实系数LSI系统的流图中的所有支路方向翻转,支路增益不变,将输入输出的位置交换,则形成原网络的转置型结构,且系统的传输函数不变。

(2)举例:

三、FIR数字滤波器结构

1、直接型结构

(1)FIR DF的特点:

(2)根据FIR数字滤波器的系统函数H(z)可以推导出系统的差分方程为

(3)根据系统的差分方程可直接作出信号流图,这就是FIR数字滤波器直接型结构,它简单直观,运算速度快,但不能直接控制零点。

2、级联型结构

(1)当需要控制滤波器的传输零点时,可将FIR数字滤波器的系统函数分解成二阶实系数因子的形成,即

        该式可由多个二阶节级联实现,如下所示

(2)级联型FIR数字滤波器结构的特点:

①级联型FIR数字滤波器结构所需的系数比直接型多,因而所需乘法运算比直接型多。

②级联型FIR数字滤波器结构的易于系统零点的控制,因而常在需要控制传输零点时使用。

3、频域采样型结构

(1)频域采样型滤波器结构的基本原理:

(2)频域采样型FIR数字滤波器结构的特点:

①可直接根据系统的频率响应的采样值构造滤波器,易于设计实现。

②谐振器柜中的每个一阶网络的系数均为复数,计算相对复杂。

③适用于窄带滤波器,即仅有少数H(k)不为0的滤波器。

④由于系数响应字长的有限,系统容易变得不稳定。

(3)频域采样型FIR数字滤波器的优化:

(4)频域采样型结构实现的说明:

①当谐振器柜中剩下少数几个所需要的谐振器时,该结构可以比直接型结构少用乘法器,但使用的存储器会多一些。

②可以共同使用多个并列的滤波器实现,在实现时,可以共用梳状滤波器及谐振柜,将各谐振器的输出适当加权组合就能组成各种所需的滤波器。

4、快速卷积型结构

(1)快速卷积型滤波器结构的基本原理:

(2)快速卷积型滤波器的结构框图:

5、线性相位型结构

(1)线性相位系统与广义线性相位系统:

①线性相位系统:

②广义线性相位系统:

③线性相位系统和广义线性相位系统统称为常数群延迟系统。

(2)线性相位系统单位脉冲响应的特点:

(3)基于线性相位系统单位脉冲响应的特点,可以进一步节省运算过程中使用的乘法器:

第l章 数字信号处理引言 1.1 引言 1.2 数字信号处理起源 1.3 信号域 1.4 信号分类 1.5 DSP:一个学科 第2章 采样原理 2.1 引言 2.2 香农采样原理 2.3 信号重构 2.4 香农插值 2.5 采样方法 2.6 多通道采样 2.7 MATLAB音频选项 第3章 混叠 3.1 引言 3.2 混叠 3.3 圆判据 3.4 IF采样 第4章 数据转换和量化 4.1 域的转换 4.2 ADC分类 4.3 ADC增强技术 4.4 DSP数据表示方法 4.5 量化误差 4.6 MAC单元 4.7 MATLAB支持工具 第5章 z变换 5.1 引言 5.2 z变换 5.3 原始信号 5.4 线性系统的z变换 5.5 z变换特性 5.6 MATLAB z变换设计工具 5.7 系统稳定性 5.8 逆z变换 5.9 赫维赛德展开法 5.10 逆z变换MATLAB设计工具 第6章 有限冲激响应滤波器 6.1 引言 6.2 FIR滤波器 6.3 理想低通FIR滤波器 6.4 FIR滤波器设计 6.5 稳定性 6.6 线性相位 6.7 群延迟 6.8 FIR滤波器零点位置 6.9 零相位FIR滤波器 6.10 最小相位滤波器 第7章 窗函数设计法 7.1 有限冲激响应综述 7.2 基于窗函数的FIR滤波器设计 7.3 确定性设计 7.4 数据窗 7.5 基于MATLAB窗函数的FIR滤波器设计 7.6 Kaiser窗函数 7.7 截尾型傅里叶变换设计方法 7.8 频率采样设计法 第8章 最小均方设计方法 8.1 有限冲激响应综述 8.2 最小二乘法 8.3 最小二乘FIR滤波器设计 8.4 MATIAB最小均方设计 8.5 MATLAB设计对比 8.6 PRONY方法 第9章 等波纹设计方法 9.1 等波纹准则 9.2 雷米兹交换算法 9.3 加权等波纹FIR滤波器设计 9.4 希尔伯特等波纹FIR滤波器 9.5 等波纹滤波器阶次估计 9.6 MATLAB等波纹FIR滤波器实现 9.7 LpFIR滤波器设计 9.8 基于Lp范数的MATLAB滤波器设计 第10章 FIR滤波器特例 10.1 引言 10.2 滑动平均FIR滤波器 10.3 梳状FIR滤波器 10.4 L波段FIR滤波器 10.5 镜像FIR滤波器 10.6 补码FIR滤波器 10.7 频率抽样滤波器组 10.8 卷积平滑FIR滤波器 10.9 非线性相位FIR滤波器 10.10 Farrow FIR滤波器 第11章 FIR的实现 11.1 概述 11.2 直接型FIR滤波器 11.3 转置结构 11.4 对称FIR滤波器结构 11.5 格型FIR滤波器结构 11.6 分布式算法 11.7 正则符号数 11.8 简化加法器图 11.9 FIR有限字长效应 11.10 计算误差 11.11 缩放 11.12 多重MAC结构 第12章 经典滤波器设计 12.1 引言 12.2 经典模拟滤波器 12.3 模拟原型滤波器 12.4 巴特沃斯原型滤波器 12.5 切比雪夫原型滤波器 12.6 椭圆原型滤波器 12.7 原型滤波器到最终形式的转换 12.8 其他IIR滤波器形式 12.9 PRONY(PADE)法 12.10 尤尔—沃尔 第13章 无限冲激响应滤波器设计 13.1 引言 13.2 冲激响应不变法 13.3 冲激响应不变滤波器设计 13.4 双线性z变换法 13.5 翘曲 13.6 MATLAB IIR滤波器设计 13.7 冲激响应不变与双线性z变换IIR对比 13.8 最优化 第14章 状态变量滤波器模型 14.1 状态空间系统 14.2 状态变量 14.3 模拟仿真 14.4 MATLAB仿真 14.5 状态变量模型 14.6 基变换 14.7 MATLAB状态空间 14.8 转置系统 14.9 MATLAB状态空间算法结构 第15章 数字滤波器结构 15.1 滤波器结构 15.2 直Ⅰ、Ⅱ型结构 15.3 直Ⅰ、Ⅱ型IIR滤波器的MATLAB相关函数 15.4 直Ⅰ、Ⅱ型结构的MATLAB实现 15.5 级联型结构 15.6 一阶、二阶子滤波器 15.7 一阶、二阶子滤波器的MATLAB实现 15.8 并联型结构 15.9 级联/并联型结构的MATLAB实现 15.10 梯型/格型IIR滤波器 第16章 定点效应 16.1 背景 16.2 定点系统 16.3 溢出(饱和)效应 16.4 算法误差 16.5 系数敏感度 16.6 二阶子滤波器 16.7 标准IIR滤波器 16.8 缩放 16.9 极限环振荡 第17章 IIR结构分析 17.1 溢出防范 17.2 Lp范数界 17.3 L2溢出预防 17.4 L2范数测定 17.5 L2范数的附加说明 17.6 L∞范数界 17.7 L1范数界 17.8 噪声功率增益 17.9 基于状态空间的噪声分析 17.10 相似变换 第18章 多采样率系统简介 18.1 背景 18.2 抽取 18.3 插值 18.4 采样率转换 18.5 多相表示法 18.6 子带滤波器 18.7 MATLAB 第19章 多采样率滤波器 19.1 引言 19.2 离散傅里叶变换(DFI)滤波器组 19.3 L波段滤波器 19.4 正交镜像滤波器 19.5 多相表达式 19.6 掩频滤波器 19.7 级联积分梳状滤波器(CIC) 附录一 MATLAB 附录二 词汇表 附录三 中英文对照 参考文献
11种经典软件滤波的原理和实现 1、限幅滤波法(又称程序判断滤波法) A、方法: 根据经验判断,确定两次采样允许的最大偏差值(设为A) 每次检测到新值时判断: 如果本次值与上次值之差<=A,则本次值有效 如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值 B、优点: 能有效克服因偶然因素引起的脉冲干扰 C、缺点 无法抑制那种周期性的干扰 平滑度差 2、中位值滤波法 A、方法: 连续采样N次(N取奇数) 把N次采样值按大小排列 取中间值为本次有效值 B、优点: 能有效克服因偶然因素引起的波动干扰 对温度、液位的变化缓慢的被测参数有良好的滤波效果 C、缺点: 对流量、速度等快速变化的参数不宜 3、算术平均滤波法 A、方法: 连续取N个采样值进行算术平均运算 N值较大时:信号平滑度较高,但灵敏度较低 N值较小时:信号平滑度较低,但灵敏度较高 N值的选取:一般流量,N=12;压力:N=4 B、优点: 适用于对一般具有随机干扰的信号进行滤波 这样信号的特点是有一个平均值,信号在某一数值范围附近上下波动 C、缺点: 对于测量速度较慢或要求数据计算速度较快的实时控制不适用 比较浪费RAM 4、递推平均滤波法(又称滑动平均滤波法) A、方法: 把连续取N个采样值看成一个队列 队列的长度固定为N 每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据.(先进先出原则) 把队列中的N个数据进行算术平均运算,就可获得新的滤波结果 N值的选取:流量,N=12;压力:N=4;液面,N=4~12;温度,N=1~4 B、优点: 对周期性干扰有良好的抑制作用,平滑度高 适用于高频振荡的系统 C、缺点: 灵敏度低 对偶然出现的脉冲性干扰的抑制作用较差 不易消除由于脉冲干扰所引起的采样值偏差 不适用于脉冲干扰比较严重的场合 比较浪费RAM 5、中位值平均滤波法(又称防脉冲干扰平均滤波法) A、方法: 相当于“中位值滤波法”+“算术平均滤波法” 连续采样N个数据,去掉一个最大值和一个最小值 然后计算N-2个数据的算术平均值 N值的选取:3~14 B、优点: 融合了两种滤波法的优点 对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差 C、缺点: 测量速度较慢,和算术平均滤波法一样 比较浪费RAM 6、限幅平均滤波法 A、方法: 相当于“限幅滤波法”+“递推平均滤波法” 每次采样到的新数据先进行限幅处理, 再送入队列进行递推平均滤波处理 B、优点: 融合了两种滤波法的优点 对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差 C、缺点: 比较浪费RAM 7、一阶滞后滤波法 A、方法: 取a=0~1 本次滤波结果=(1-a)*本次采样值+a*上次滤波结果 B、优点: 对周期性干扰具有良好的抑制作用 适用于波动频率较高的场合 C、缺点: 相位滞后,灵敏度低 滞后程度取决于a值大小 不能消除滤波频率高于采样频率的1/2的干扰信号 8、加权递推平均滤波法 A、方法: 是对递推平均滤波法的改进,即不同时刻的数据加以不同的权 通常是,越接近现时刻的数据,权取得越大。 给予新采样值的权系数越大,则灵敏度越高,但信号平滑度越低 B、优点: 适用于有较大纯滞后时间常数的对象 和采样周期较短的系统 C、缺点: 对于纯滞后时间常数较小,采样周期较长,变化缓慢的信号 不能迅速反应系统当前所受干扰的严重程度,滤波效果差 9、消抖滤波法 A、方法: 设置一个滤波计数器 将每次采样值与当前有效值比较: 如果采样值=当前有效值,则计数器清零 如果采样值<>当前有效值,则计数器+1,并判断计数器是否>=上限N(溢出) 如果计数器溢出,则将本次值替换当前有效值,并清计数器 B、优点: 对于变化缓慢的被测参数有较好的滤波效果, 可避免在临界值附近控制器的反复开/关跳动或显示器上数值抖动 C、缺点: 对于快速变化的参数不宜 如果在计数器溢出的那一次采样到的值恰好是干扰值,则会将干扰值当作有效值导入系统 10、限幅消抖滤波法 A、方法: 相当于“限幅滤波法”+“消抖滤波法” 先限幅,后消抖 B、优点: 继承了“限幅”和“消抖”的优点 改进了“消抖滤波法”中的某些缺陷,避免将干扰值导入系统 C、缺点: 对于快速变化的参数不宜 第11种方法:IIR 数字滤波器
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Zevalin爱灰灰

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值