动态信号处理(第四、五、六章).ppt

目录 第四章 频谱分析 第五章 数字滤波 第六章 细化选带频谱分析 4.1 卷积与相关定理 卷积 相关分析 4.2 连续傅立叶变换 傅立叶级数 傅立叶变换 4.4 傅里叶变换的性质 线性性 时间尺度的变化 频率尺度的变化 对称性 时间位移 频率位移 偶函数 奇函数 复时间函数 小结 如果的自变量被移动一个常量，则作代换，傅立叶变换为： 于是有： 5.5 FIR非递归滤波器 1. FIR数字滤波器的设计方法 2.FRI低通数字滤波器的设计 3.FRI高通数字滤波器的设计 4.FRI带通数字滤波器的设计 5.FRI多通带带通数字滤波器的设计 6. 滤波输出序列的计算 5.6复解析带通滤波器 FIR非递归带通数字滤波器及其特点 FIR非递归复解析带通滤波器 基于复解析带通滤波器的复调制细化分析法 基于复解析带通滤波器的复调制细化谱分析原理和方法 基于复解析带通滤波器的复调制细化选带频谱分析算法 密集频率成分的离散频谱校正方法及其实现 4.FRI带通数字滤波器的设计 冲激响应函数为： 为带通上限频率， 为带通下限频率， 为1/2采样频率，2M+1为窗口的长度， M 称 为半阶数。 5.FRI多通带带通数字滤波器的设计 设有K个通带，第j个通带的上限频率为 ，下限频率为 ，其冲激响应函数为： 式中 为1/2采样频率，2M+1为窗口的长度，M称为半阶数。 6.滤波输出序列的计算 由于FIR滤波器冲击响应函数是偶对称（保证零相移），故只需计算滤波卷积因子的一半，滤波阶数为2M+1（M为滤波器的半阶数）。对一信号序列为x（n），n 0,1,2,…,N，其滤波输出序列为： 滤波输出序列比原始输入序列两头各少N点。滤波半阶数M越大，通带与阻带的过渡带越陡，滤波效果越好，但滤波卷积序列所需计算点数越大，计算量越大，综合考虑，一般取半阶数为50~300。 3.按照滤波器的实现方法分类 1 .递归滤波器。传统硬件分析仪器大都采用IIR无限冲击响应递归滤波器，即它的输出不仅取决于过去和现在的输入，而且也取决于过去的输出。其差分方程可表示为： ──滤波后的输出值序列； ──滤波前的输入值序列； ──为滤波系数序列。 此法优点是结构简单，每节滤波器的阶数少，可用递归方法提高滤波器的滤波特性，适用于硬件实现。但是，软件数字滤波采用这种方法的最大缺点是由于每节初始值是人为设定导致滤波输出序列的起始段失真 称为截头效应），其相位信息也失真。 Z变换得到频率响应函数为： 2 .非递归滤波器 数字滤波器大都采用FRI有限冲击响应非递归滤波器，即它的输出仅取决于过去和现在的输入，与过去的输出无关。其差分方程可表示为： ──滤波后的输出值序列； ──滤波前的输入值序列； ──为滤波系数序列。 Z变换得到频率响应函数为： 采用这种方法其缺点是每节滤波器的阶数多，不宜用硬件实现；优点是滤波一次完成，没有递归滤波的中间序列，易于软件实现且精度高，没有递归滤波的初始值需由人为赋给而产生的初始段滤波值不准确（截头效应）的缺点，其计算出的每一点都是对应原始序列的精确滤波值，且保留了相位信息；用于细化谱分析、离线重抽抗混滤波等运算时，只需计算重抽点或细化点滤波值即可，大大减少了运算量，提高了精度。 计算机软件滤波器一般常采用有限冲激响应（FIR）滤波器。 5.2滤波器的技术指标 实际上并不希望滤波器具有理想的幅频特性。因为具有理想幅频特性的滤波器在时域中的单位阶跃响应有过冲和振铃现象，而且是非因果性的，无法实现。而因果性是物理系统实现的必要条件，因此只能用一个因果系统去逼近它。另外还要考虑 逼近后 系统的难易和精度、速度的问题，所以实用中的数字滤波器总是在通带中允许一定的容限，阻带中允许另一个一定的容限，而且在通带与阻带之间有一个过渡带。 实用数字滤波器总是在通带中允许一定的容限，阻带中允许另一个一定的容限，而且在通带与阻带之间有一个过渡带。 这种情况我们在图6.2.1中用低通滤波器来说明。 下图中为 通带的容限， 为阻带的容限， 为通带的截止频率， 为阻带的下限频率，在 与 之间为过渡带。但在严格一些意义上，可以把0.8 （即精度不能满足通带的容限 要求）与 之间定义为过渡带。 a 滤波器线性幅频特性 b 滤波器对数幅频特性 低通滤波器的实际幅频特性 滤波器的主要技术指标有： 1 .截止频率 Hz ：一般定义为对数幅频特性上幅值为-3dB点的横坐标值。 2 .倍频程衰减率 -x/OCT ：定义为对数幅频特性横坐标值为 时的幅值。这个指标越大越好，一般硬件抗混滤波器大于36dB/OCT，现代频谱分析仪器的混滤波器大于100dB/OCT。 3 .通带内的波动率 滤波器精度 。有两种表示方法：