MATLAB辅助现代工程数字信号处理 第二版 教学课件 作者 李益华 第1 5章 第3章.pptVIP

第3章 数字滤波器的结构与分析 　　　　3.1 数字滤波器及其实现3.1.1 数字滤波器概述　　数字滤波器是数字信号处理的基础部分，与模拟滤波器相比，数字滤波器具有精度高、可靠性高、灵活性高、便于大规模集成和多维过滤等优点，已广泛应用于现代各类工程领域。　　数字滤波器是具有一定传输特性的数字信号处理装置，它的输入和输出都是离散数字信号，它借助于数字器件和一定的数值计算方法，对输入信号进行处理，改变输入信号，进而去掉信号中的无用成分而保留有用成分。如果在数字处理系统前、后分别加上A/D转换器和D/A转换器，就可以处理模拟信号。 　　数字滤波器的输入输出是一个时间序列。设H(z)为数字滤波器的系统函数，h(n)为其相应的脉冲序列，则在时域内有　　　　　　　　　　y(n)=x(n)*h(n)　　　　　　　　(3.1)　　在z域内有　　　　　　　　　　Y(z)=H(z)X(z)　　　　　　　　(3.2)式中，X(z)和Y(z)分别为输入x(n)和输出y(n)的Z变换。　　在频域内有　　　　　　　　　Y(jω)=H(jω)X(jω)　　　　　　(3.3)式中， H(jω)为数字滤波器的频率特性; X(jω)和Y(jω)分别为输入x(n)和输出y(n)的频谱。 3.1.2 数字滤波器的分类与实现　　数字滤波器有多种分类方法， 按功能，数字滤波器可分为低通滤波器、高通滤波器、带通滤波器和带阻滤波器。按脉冲响应的长度，数字滤波器可分为无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。　　IIR滤波器的差分方程为 　　FIR滤波器的差分方程为 图 3.1 数字滤波器的实现 　　一个输出序列是其过去N点的线性组合加上当前输入序列与过去N点输入序列的线性组合。输出y(n)除了与当前的输入x(n)有关外，同时还与过去的输入和过去的输出有关，系统是带有记忆的。　　对于上面的算式，可以化成不同的计算形式，如直接计算、分解为多个有理函数相加、分解为多个有理函数相乘等，不同的计算形式表现出不同的计算结构，而不同的计算结构可能会带来不同的效果，或者是实现简单、编程方便，或者是计算精度较高等。 　　由于数字信号是通过采样和转换得到的，而转换的位数是有限的，所以存在着量化误差。此外，计算机中数的表示也总是有限的，经此表示的滤波器的系数同样存在着量化误差，故在计算过程中因有限字长也会造成误差。　　量化误差主要有三种: A/D变换量化效应，系数的量化效应，数字运算的有限字长效应。 3.1.3 数字滤波器的运算结构　　数字滤波器无论采用硬件实现还是采用软件实现，首先都应确定数字滤波器的运算结构。运算结构可以用方框图表示，也可用信号流图来表示。为了简单起见，通常用信号流图来表示其运算结构。信号流图表示滤波器有加法、乘法及延迟三种基本运算单元。　　在信号流图中，只有输出支路的节点称为输入节点或源点; 只有输入支路的节点称为输出节点或阱点; 既有输入支路又有输出支路的节点叫做混合节点。通路是指从源点到阱点之间沿着箭头方向的连续的一串支路，通路的增益是该通路上各支路增益的乘积。回路是指从一个节点出发沿着支路箭头方向到达同一个节点的闭合通路，它象征着系统中的反馈回路。组成回路的所有支路增益的乘积通常称为回路增益。 　　在信号流图中，系统函数按梅逊(Mason)公式计算: 　　信号流图转置可以转变运算结构，并可以验证所计算的流图系统函数的正确与否。　　运算结构对滤波器的实现很重要，尤其对于一些定点运算的处理机，结构的不同将会影响系统的精度、误差、稳定性、经济性以及运算速度等许多重要的性能。对于IIR数字滤波器与FIR数字滤波器，它们在结构上各有自己不同的特点，下面分别加以讨论。 3.1.4 数字滤波器的MATLAB实现　　MATLAB信号工具箱提供了很多内部函数来设计数字滤波器，具体函数将在以后各章中逐一介绍，这里仅举几例来介绍数字滤波的MATLAB实现。　　【例 3.1】 用MATLAB的内部函数分别设计一个巴特沃斯和切比雪夫数字滤波器。　　MATLAB程序如下: 　　　　%MATLAB PROGRAM 3-1　　　　clc; 　　　　fp=2; 　 %2 kHz　　　　fs=10; 　　%10 kHz　　　　Rp=3; 　　%3 dB　　　　Rs=60; 　　%60 dB　　　　Wp=fp/fp; 　　　　Ws=fs/fp; 　　[n，Wn]=buttord(Wp，Ws，Rp，Rs，‘s′); %巴特沃斯滤波器阶数选择　　[z，p，k]=buttap(n); 　　　　　　　　　%巴特沃斯滤波器设计