数字信号处理王震宇张培珍编第.ppt

上两式可以化成不同的计算形式(其中M≤N)，如直接计算、分解为多个有理函数相加、分解为多个有理函数相乘、交换运算次序等等。不同的计算形式和顺序也就表现出不同的计算结构。 但从运算上看，只需要加法、单位延迟、乘常数三种运算，因此数字滤波结构中有三个基本运算单元，即加法器，单位延时器，乘常系数乘法器。 系统函数|H(z)|在|z|0处收敛，极点全部在z=0处(即FIR一定为稳定系统)，结构上主要是非递归结构，没有输出到输入反馈，但频率采样结构也包含有反馈的递归部分。 FIR滤波器直接型结构特点如下。 (1) 实现简单，但结构相对复杂，需要N-1个延时器和N个常系数乘法器。 (2) 系数量化会受到有限字长效应的影响，从而产生较大误差。 num=[-4 6 5 6 -4]; %输入分子系数向量 den=[1 0 0 0 0]; %输入分母系数向量 [z,p,k]=tf2zp(num,den); %求出各子系统的零极点 [sos,A]=zp2sos(z,p,k); %求出各二阶节乘系数 disp(零点:);disp(z); disp(极点:);disp(p); disp(增益系数:);disp(k); disp(二阶节:);disp(real(sos)); 根据第7章的讨论，用频率采样表达系统函数的内插公式为 上式既包含极点，也包含零点，所以这时滤波器具有递归结构。频率采样型结构,是两部分级联而成，即 其中H1(z) 为全零点滤波器，即其零点位于单位圆的等间隔点上 (1) 由两个网络级联构成，第一个横向结构M节延时网络实现零点，第二个有反馈的N节延时网络实现极点，物理概念清楚、简单、易于理解。 (2) 需要的延迟单元太多，共需M+N个延迟器，当M=N时为2N个。 (3) 系数ai和bi对滤波器性能的控制不直接，一个ai和bi的改变会影响系统的零点或极点分布，对极、零点的控制难。 (4) 容易出现不稳定或产生较大误差。 按照差分方程画出如图8.19所示的直接Ⅰ型和直接Ⅱ型结构。 利用程序将H(z)的分子和分母进行因式分解。 a=[1 -3 0 -4 9] b=[5 -12 -2 4 2] [z,p,k]=tf2zp(a,b); [sos,A]=zp2sos(z,p,k) 级联结构如图所示 程序运行结果 留数： -0.0586 -0.5233 -1.8590 - 2.1155i -1.8590 + 2.1155i 极点： 2.3987 0.7851 -0.3919 - 0.2425i -0.3919 + 0.2425i 常数：4.5000 程序运行结果： 分子多项式的系数向量： -3.7181 -2.4831 0 分母多项式的系数向量： 1.0000 0.7839 0.2124 将复根合并成系数二阶节，所以并联型结构系统函数为 重点习题 2（2）；6；7；9，10 前面讨论了有限长序列系统函数H(z)在单位圆上作N等分采样，这个采样值也就是h(n)的离散傅里叶变换值H(k)，即 4. 频率采样型 8 4. 频率采样型 8 4. 频率采样型 8 N=7和N=8时全零点滤波器的零点分布图 滤波器H1(z)又被称为梳状滤波器，其频率响应和幅频特性表达式分别为 (8.12) 4. 频率采样型H1(z) 8 另一个滤波器的系统函数为 是由N个单极点的一阶滤波器并联构成，极点正好与梳状滤波器的一个零点(i=k)相抵消，从而使频率ω=2πk/N上的频率响应等于H(k)。 4. 频率采样型H2(z) 8 4. 频率采样型H2(z) 8 所有极点都是在单位圆上,由于系数量化的影响，极点就不能被梳状滤波器的零点所抵消，系统稳定性变差。 为了克服稳定性变差的缺点，可以首先进行修正，将所有谐振器的极点从单位圆向内收缩，使它处在一个靠近单位圆但半径稍小的圆上r≤1，同时，梳状滤波器的零点也移到半径为r的圆上，也即将频率采样由单位圆移到修正半径r的圆上 。 4. 频率采样型结构存在问题 8 FIR数字滤波器有哪四种基本网络结构 ？每种结构的特点。 （1）直接型 （2）级联型 （3）线性相位型 （4）频率采样型 知 识