冲激响应对称冲激响应反对称在传输函数的直接型实现中,利用线性FIR.pptVIP

第八章 数字滤波器的结构 例：级联格型数字滤波器结构 将最后两个方程代入第一次建立的四个方程中可得: 联立最后两个方程得到 数字滤波器结构在方框图表示中延时的数量等于传输函数的阶数称这种结构是规范结构 否则，称为非规范结构 因此，有必要选择一种在有限字长条件下，具有最小量化效应的实现结构 一种方法是先确定多种等效结构，分析每一种结构的有限字长效应，最后选出具有最小量化效应的那一个 在某些情况下，也可以直接通过具有最小量化效应的构造生成 转置后的直接型结构更是容易的,如下 可以用类似的分解实现一个2型FIR传输函数 例如,一个长度为8的2型FIR传输函数,相关分解为 注意:长度为7的1型FIR结构需要4个乘法器,而用直接型来实现原始滤波器需要7个乘法器 长度为8的2型FIR滤波器需要4个乘法器,而用直接型来实现原始滤波器需要8个乘法器 这种节省也可以在具有反对称冲激响应的FIR滤波器中获得 我们涉及的因果IIR数字滤波器可以用形如z-1的实有理传输函数来描述,或者用常系数差分方程来表示 从差分方程的表达式可以看出,要计算第N个输出样本值,需要知道一些前面的样本值,换言之,需要一定形式的反馈 N阶IIR数字滤波器的传输函数是用2N+1个不同的系数描述的,通常需要2N+1个乘法器和2N个两输入加法器来实现 直接型IIR滤波器: 若乘法器的系数等于传输函数的系数,则这种IIR滤波器结构就称为直接型结构 将H1(z)和H2(z)进行级联,得到原始IIR传输函数H(z)的实现方式,如下图,这就是我们所说的直接I型结构 交换H1(z)和H2(z) 的位置，系统函数不变。框图变换得到下图： 直接II型 一个N阶IIR传输函数的实现可以被证实 通过不同的极零点多项式对可以得到不同的级联实现,如下图 由于有限字长效应的影响,每种级联实现的性能是不相同的 通常,可以通过因式分解把一个多项式分解 为一阶多项式和二阶多项式的积: 考虑三阶传输函数 例-直接II和级联型的实现 IIR传输函数可以通过部分分式展开以并联的形式来实现,给出的部分分式称为并联I型 假设为简单极点,传输函数H(z)可以表示为 IIR数字滤波器结构的并联实现 三阶IIR传输函数的两种基本并联实现如图所示 IIR数字滤波器结构的并联实现 例-一个部分分式展开 IIR数字滤波器结构的并联实现 相关并联I型实现如右图 IIR数字滤波器结构的并联实现 另一种是将H(z)直接分解为Z的多项式的比值的形式 MATLAB实现 级联型要求传输函数的因式分解,可以用M文件zp2sos sos = zp2sos(z,p,k)生成一个矩阵sos包含每个二阶等效传输函数H(Z)的系数,此传输函数H(Z)由它的零极点决定 MATLAB实现 它的第i行包括系数{pil}和{dil},它们分别是二次多项式第i项的分子分母 MATLAB实现 L 代表长度 总的传输函数的形式如下 MATLAB实现 注意: FIR传输函数可以看作是分子是带常系数的IIR传输函数，分母是一个FIR描述的的多项式的分式。 MATLAB实现 并联I和II型可以分别用函数residuez 和residue实现 Program 6_2 就是两个函数实现的 作业 习题 8.5(6.4), 8.10, 8.24(6.26), 8.26(6.28) M8.1, M8.2, M8.3 习题 8.5 分析图P8.4所示的数字滤波器结构并求出它的传输函数H(z)=A(z)/X(z). (a) 这是一个规范结构吗? (b) 若使H(z)在w=0处有单位增益,乘法器系数K的值应该是什么? (c)若使H(z)在w=∏处有单位增益,乘法器系数K的值应该是什么? (d) 这两个K值有区别吗?若没有,为什么? 习题 8.5 图 P8.4 习题 8.10 用如下形式实现FIR传输函数 H(z)=(1-0.6z-1)6=1-3.6z-1 + 5.4z-2 - 4.32z-3 + 1.944zz-4 - 0.4666z-5 + 0.0467z-6 两个不同的直接形式, 六个1阶部分的级联, 三个2阶部分的级联, 两个3阶部分的级联, 两个2阶部分和两个1阶部分的级联. 比较上面每个实现的计算复杂度 习题 8.24 因果LTI离散时间系统对输入x[n]=(0.2)nu[n ]的输出为y[n]=(0.4)nu[n]-0.3(0.4)n-1u[n-1] (a)确定这个系统的传输函数. (b) 确定描述这个系统的差分方程 (c)用不多于三个乘法器来给出这个系统的一个规范直接形式实现. (d)给出系统的并联I型实现. (e)以封闭的形式给出该系统的冲激响应. (f) 求系统在输入