数字信号处理-丁玉美_高西全_编著-第4章.ppt

第4章　时域离散系统的网络结构及 数字信号处理的实现 4.3.3 IIR中的级联结构和并联结构　　IIR基本结构有直接型、 级联型和并联型。 一般低阶的用直接型， 高阶的用级联型或并联型。 　　在设计级联型结构时， 需要将分子式和分母式进行因式分解， 阶数高时可借助于MATLAB语言用计算机解决。 设计并联结构时要进行部分分式展开。 　　部分分式展开要求分子多项式的阶数低于分母多项式的阶数， 否则是一个假分式（分子多项式的阶数不低于分母多项式的阶数）， 要将其化为整数和真分式之和， 然后再对真分式进行部分分式展开。 部分分式的各系数通过待定系数法解决。 部分分式的一般表达式为 上面两式说明系统H1(z)和H2(z)分别在x(n)=Aδ(n)的激励下可以分别产生正弦波和余弦波。H1(z)和H2(z)的极点为p1,2=e±jω0, 这正是在单位圆上的两个极点， 极点的相角为ω0。 这样，H1(z)和H2(z)可以分别称为正弦波和余弦波发生器， 画出H1(z)实现结构图如图4.4.3所示， 共需要两个乘法器、 两个加法器和两个移位器。 运行时要用x(n)=Aδ(n)作激励。也可以令图中的v(n)起始条件为v(0)=A， v(－1)=0，v(－2)=0， 代替输入信号x(n)=Aδ(n)。在实际应用中有时需要两个正交相位正弦波， 可以将H1(z)和H2(z)进行组合， 同时产生正弦波和余弦波， 实现结构如图4.4.4所示。 图4.4.6 题2解图 题3解图（一） 题4解图 题6图 　　12. 已知FIR滤波器系统函数在单位圆上16个等间隔采样点为: 　　H(0)=12， H(3)~H(13)=0　　 H(1)=－3－j　　，　　H(14)=1－j　　 H(2)=1+j，　　　　　 H(15)=－3+j　　 试画出它的频率采样结构, 取修正半径r =0.9, 要求用实数乘法器。　　解： 题14解图 　　16. 画出题15图中系统的转置结构， 并验证两者具有相同的系统函数。 　　解： 按照题15图， 将支路方向翻转， 维持支路增益不变， 并交换输入输出的位置， 则形成对应的转置结构， 画出题15图系统的转置结构如题16解图所示。 　　 将题16解图和题15图对照， 它们的直通通路和反馈回路情况完全一样， 写出它们的系统函数完全一样， 这里用Masson公式最能说明问题。 题17图 题19*解图 　　20. 已知FIR滤波器的系统函数为：　　 （1） H(z)=1+0.8z－1+0.65z－2　　 （2） H(z)=1－0.6z－1+0.825z－2－0.9z－3试分别画出它们的直接型结构和格型结构， 并求出格型结构的有关参数。　　解： 已知FIR滤波器的系统函数， 设计相应的格型结构需要用到的公式如下： 　　　　　　　ak=h(k) 题21解图 题22解图（二） 题22解图（三） 式中， 分母分子多项式各有一个零点z=1， 相互抵消， 因此该系统仍然稳定， 属于FIR系统。 由系统函数得到单位脉冲响应为 　　　h(n)=5δ(n)+5δ(n－1)+5δ(n－2)+3δ(n－3) 　　　　　+3δ(n－4)+3δ(n－5) H(k)=DFT［h(n)］ k=0, 1, 2, …, 5 　　按照上式画出频率采样修正结构如题19*解图所示。 图中系数 　　　a0k=2Re［H(k)］, a1k=－2Re［rH(k)W6－k］ 求系数程序ex519.m如下：  　　%程序ex519.m 　　hn=［5， 5， 5， 3， 3， 3］； r=0.95；  　　Hk=fft(hn， 6)；  　　for k=1： 3，  　　 hk(k)=Hk(k)； Wk(k)=exp(－j*2*pi*(k－1)/6)；  　　end 　　H0=Hk(1)； H3=Hk(4)； r0k=2*real(hk)； 　　　r1k=－2*real(r*hk.*Wk) 　　程序运行结果：  　　H(0) = 24 　　H(3) = 2 　　r0k = 48 4 0 　　r1k = －45.6000 3.8000 0 　　得到 　　　　　α01=48, α02=4, α11=－45.2, α12=38 进一步的说明： 此题h(n)的长度为6， 由单位圆上采样6点得到频率采样结构， 满足频率采样定理。