数字滤波器第一讲.pptVIP

主讲人：李旭涛 《数字信号处理》课程 2010-2011学年度 秋季学期 ;主要内容;1. 线性移不变系统（LSI）;LSI系统的输入，冲激响应和输出之间的关系 ; 讨论离散时间系统，手段主要有两个： 序列的卷积——时域方法； Z变换分析法——频域方法。 二种方法统一于时频对应。 ;系统函数与传输函数 ; 进一步 单位脉冲响应为系统函数的反Z变换： 这就是分析LSI系统的频域方法。 ;离散时间系统的因果性和稳定性 ; 稳定系统意味着系统对每个有界输入，产生一个有界输出的系统。LSI系统是稳定系统的充要条件为： 系统稳定要求 的收敛域包含单位圆 ;离散时间系统的频率响应 ;2. 数字滤波器的基本概念 ;数字滤波器的特点 ;数字滤波器的实现形式 ;直接I型 ;级联型和并联型的基本思想是对系统函数进行因式分解（有理分解），分解为若干个子系统的并联和级联。 ;并联型：;值得注意的问题： 有限字长效应的存在； 级联结构具有累积误差大的缺点，但是便于准确实现滤波器的极零点，便于调整滤波器的频率响应的性能，所需要的存储器少； 并联结构的每一个子系统都是独立的，不受其他子系统的量化误差和舍入误差的影响，因此对于误差最不敏感，但是对于极零点的调整不易。;滤波器的Z域极零点分析 ;设系统稳定，将 代入，得传输函数 ;在Z平面上向量： 分别称为零点矢量和极点矢量，用极坐标表示为 ;对应的幅度谱： 相位谱为： 当频率 从0变化到 时，这些向量的终点B沿单位圆逆时针旋转一周。 如果滤波器的N个零点等间隔分布在单位圆上，设N=8，极零点分布如下图。当从0变化到时 ，每遇到一个零点，幅度为零。具有如下图所示幅度特性的滤波器称为梳状滤波器。 ;梳状滤波器的极零点分布及幅度特性 ;3. FIR数字滤波器 ;FIR数字滤波器的主要优点有： 1、易于实现线性相位； 2、由于不存在极点，系统是稳定的。 FIR数字滤波器的缺点： 需要较长的滤波器长度以满足短过渡带的要求，造成所需的计算量较大。 ;线性相位 ;因此， 的零点也是 的零点，反之??然。这说明线性相位的FIR滤波器其零点相对于单位圆对称分布，如图为线性相位FIR滤波器的零点分布 ;4. FIR数字滤波器的设计 ;理想滤波器;理想的低通滤波器;窗函数法设计FIR滤波器;设理想低通滤波器的幅频为1，相频为零：;;时域的乘积－频域的卷积;于是：;于是：;如何设计FIR滤波器？;频带指标;频带起伏指标;窗函数和这些技术指标之间的关系如何呢？也就是说具体的窗函数如何同这些指标相对应。;窗函数;Kaiser窗;Kaiser窗的设计参数为;常见的几种窗函数及其对应滤波器的频率响应;;设计步骤;不同的窗函数的设计在于其参数同指标要求之间的对应关系不同； 无外乎窗口的长度如何确定，窗口的形状参数如何同指标相对应； 工程设计中通常采用简化的经验设计公式，具体的内容可以查工程手册得到;频率采样法（时域的设计方法）;;clear all clf fc=15; %kHz% fs=40; %kHz% f = [-fs/2,-fc,-fc, 0, fc,fc,fs/2]; Df=[0, 0, 1, 0.5,1, 0, 0]; axis([-fs/2,fs/2,-0.2,1.2]) hold on plot(f,Df) xlabel(analog frequency f in kHz) title(desired AF frequency response D(f)) grid hold off pause;;;;clf N=21; M=(N-1)/2; delta_f=fs/N; fk=-fs/2*(1-1/N): delta_f : fs/2*(1-1/N); Dfk=((1+abs(fk)/fc)/2).*(sign(fk+fc)-sign(fk-fc))/2; axis([-fs/2,fs/2,-0.2,1.2]) hold on plot(fk,Dfk,*) xlabel(analog frequency f in kHz) title(samples of D(f)) grid f = [-fs/2,-fc,-fc, 0, fc,fc,fs/2]; plot(f,Df,r) hold off pause;;clf k=-M:M; Dk=Dfk; axis([-M,M,-0.2,1.2]) hold on bar(k,Dk,0) plot(k,Dk,*)