13数字滤波器(1.docVIP

第1章 前言 1.1 课题研究的目的和意义 数字信号处理（Digital Signal Processing，简称DSP）是一门涉及多门学科并广泛应用于很多科学和工程领域的新兴学科，它是指利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理，以便提取有用的信息并进行有效的传输与应用，得到符合人们需要的信号形式。与模拟信号处理相比，数字信号处理具有精确、灵活、抗干扰能力强、可靠性高、体积小、易于大规模集成等优点。数字信号处理的实现在早期采用通用的处理器实现，但其速度低，价格昂贵。数字信号处理器（Digital Signal Processor）的出现则改变了这一状况，其内部带有专门用于数字信号处理的硬件，因此速度非常快，方便实时快速地实现各种数字信号处理算法。数字滤波器由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的 图3-1幅频和相频 用Simulink仿真窗函数设计已知滤波器的阶本例中，首先在中选择pass；在 (3-18) (3-19) 为了减小H（k）的通带边缘由于抽样点的变化而引起的起伏振荡，可以增加过渡点，加宽过渡带以减小通带的起伏。 每一个抽样值产生一个与sin()/sin()成正比，并位移（2k）/N的频率响应，而H（k）与内插函数的线性组合就是FIR滤波器的频率响应，增加一点过渡可以使阻带衰减提高到-44~54dB,二点过渡衰减-65~75,三点过渡衰减-85~95dB. 如果不能使过渡带太宽，同时要求增大阻带衰减，可以增加取样点数N，但这样会增加计算量、延时和误差。 频率取样型FIR滤波器设计步骤： (1)给定理想滤波器频率响应Hd(ejw)。 (2)根据过渡带宽和阻带衰减确定过渡点数和h(n)的长度N。 (3-20) (3)由IFFT计算IDFT得到： (3-21) 设计实例：率采样法设计一个带通滤波器，满足：低阻带边缘：w1s=0.2*;低通带边缘：w1p=0.35*;高通带边缘：w2p=0.65*;高阻带边缘：w2s=0.8*。设计过渡带中的频率样本值为t1和t2，取t1=0.109021,t2=0 设计程序如下： M=40; al=(M-1)/2; l=0:M-1; t1=0.109021; t2=0 Hrs=[zeros(1,5),t1,t2,ones(1,7),t2,t1,zeros(1,9),t1,t2,ones(1,7),t2,t1,zeros(1,4)]; k1=0:floor((M-1)/2); k2=floor((M-1)/2)+1:M-1; angh=[-al*(2*pi)/M*k1,al*(2*pi)/M*(M-k2)]; H=Hrs.*exp(j*angh); h=real(ifft(H,M)); freqz(h,1,512,1000) 实验得幅频相频特性如下图所示： 图3-5 幅频相频 3.3最优化方法设计FIR数字滤波器 MATLAB7.0信号处理工具箱提供了通用的工具函数firls和remez，采用不同的优化方法可以设计最优滤波器和多频带滤波器。 firls是fir1和fir2函数的扩展，其基本设计准则是利用最小二乘法使期望的频率响应和实际的频率响应之间的整体误差最小。Remez函数则利用Parks-McClellan算法设计线性相位FIR滤波器。该算法能使期望频率响应和实际频率响应之间的最大误差最小。Firls函数设计的滤波器在频率响应上出现等波纹，因此也叫等波纹滤波器。函数firls、remez调用的语法规则相同。 (1)firls。其调用格式为 Ⅰ. b=firls(n,f,a) Ⅱ. b=firls(n,f,a,w) Ⅲ. b=firls(n,f,a,’ftype’) Ⅳ. b=firls(n,f,a,w,’ftype’) 该函数用于设计n阶FIR滤波器，其幅频特性由f和a向量确定，f是频率点向量，其范围为0~1，是递增向量；a是指定频率点的幅度响应，与f必须同长度；b为函数返回的滤波器系数，长度为n+1,且具有偶对称关系： b(k)=-b(n+2-k), k=1,2,…,n+1 b=firls(n,f,a,w)则使用权系数w给误差加权。 b