信号分析与数据处理02幻灯片.pptVIP

数字信号处理 多 媒 体 教 学 课 件 前 言 序列Z变换的收敛域ROC 特点： 仅影响第k对零点，同样 仅影响第k对极点，子网络的零、极点也是整个系统的零、极点。便于调节滤波器的频率特性。 A 2. 并联型： 将H（Z）展成部分分式形式。 其中， 均为实数， 与 复共轭； 当MN时，不包含 项；M=N 时，该项为 G0。 一阶环节： 二阶环节： 传递环节： 其结构图如下： |H(ejω)| π -π ω Arg(H(ejω)) π -π ω 5.4 FIR数字滤波器的结构： 特点：*无反馈系统，系统只有一个极点 Z=0，为 M阶极点。ROC：|Z|0。 *无论差分方程的系数取任何有效的值，系统都是因果稳定的。 *冲击响应等于差分方程系数： h(n)=bn n=0,1,·····,M *设计时选定阶数M和系数bn使系统特性满足设计指标。 5.4.1 FIR数字滤波器计算直接实现流程图： z-1 z-1 z-1 b0 b1 b2 bM x(n) y(n) 1、直接实现方框图： h(0) h(1) h(2) h(N-2) h(N-1) 2、直接实现信号流图（横截型）： 5.4.2 FIR数字滤波器计算级联实现流程图： 系统函数可以表示成二阶因式的乘积： 得级联实现的信号流图： 特点：便于灵活调整零点位置。但设计时系数的计算量增大。 5.4.3 FIR数字滤波器的快速算法实现流程图： 因为FIR数字滤波器的冲击响应h(n)为有限长度，所以可以采用FFT 方法（快速算法）计算DFT： 当MN 时，可以采用分段卷积方法提高计算效率，减少计算量。 Arg(H(ejω)) π/2 -π/2 π -π ω |H(ejω)| π -π ω 数字信号处理多媒体教学系统 版权所有：yuning 2003。3 第2版 第 六 章IIR数字滤波器的理论与设计 数字信号处理 （1）软件实现。用于各种多媒体系统。 （2）单片信号处理微机（DSP）实现。如著名的TMS320系列信号处理机。具有固化的FFT程序，计算1024点浮点数据只需2mS；128阶FIR数字滤波器可处理10MHz的信号。用于实时多媒体系统，如可视电话，语言识别系统等。 4.3.1 FFT的实现方法： （3）阵列处理机实现。如信号处理阵列机；并行处理机等。用于大型图象信号处理，如气象卫星云图处理等。 4.3.2 FFT的应用 1、信号的频谱分析： 抗混叠 滤波器 抽样 A/D 频谱分析 FFT 窗口截断 x(t) x[n] X[k] 主要误差的处理：混叠误差（抗混叠滤波，提高采样频率）；泄漏误差（窗口函数的修正）；栅栏效应（采样点数N的选定，延长采样时间t=NT）。 A/D变换器的精度要求要根据实际情况，使频谱X(k)满足分析要求。 4.3 FFT的实现方法与应用 例：某FFT信号分析处理器采样频率为 5.12KHz。为了便于计算，取信号数据点数N是2的整数次方。问该处理器信号分析的频率范围是多少？若信号分析要求频谱的分辨率达10Hz，求信号采样持续时间至少为多长？信号数据点数应为多少？ |H(k)| fs/2 -fs/2 tp 2、快速卷积 FFT 计算（快速算法）： 一般应用中，有M N。这时，可以采用分段计算方法，进一步减少计算量。 分段计算，可以用重叠相加法和重叠保留法。 4.4 线性调频Z变换算法 DFT频谱是序列的Z变换在单位圆上的等间隔取样。但在一些应用中，需要计算某一段范围的较密集取样点的频谱；或非等间隔取样点的频谱；甚至可能要求频谱的取样点不在单位圆上，而在某一条螺旋线上。对于这样一些频谱计算要求，DFT计算无法满足。在这些情况下，采用线性调频Z变换算法（Chirp Z Trasnation CZT)是很有效的。 4.4.1 CZT算法的基本原理 CZT计算流程 4.4.2 CZT算法的快速算法： 4.4.2 CZT算法的快速算法计算步骤： 4.4.3 CZT算法的计算量： 4.5 实数序列的FFT高效算法和高效卷积： FFT计算为复数运算，所以输入序列 x(n)在运算时可以为复数数据。如果是实序列一般是把虚部置另。如果利用虚部数据可以提高计算效率。 所以，通过组合后计算一次N点FFT可以得出2个N点实数序列的DFT频谱。 4.5.1实数序列的FFT高效算法 如果希望计算2N点实数序列的FFT，也可以采用组合方法提高计算效率。 高效卷积可