第5章时域离散系统的网络结构剖析.ppt

　　18. 对于题18图中的系统，要求： 　（1） 确定它的系统函数 　（2） 如果系统参数为 　（a） b0=b2=1, b1=2, a1=1.5, a2=－0.9 　（b） b0=b2=1, b1=2, a1=1, a2=－2 画出系统的零极点分布图，并检验系统的稳定性。 题18图 　　19*. 假设滤波器的系统函数为  在单位圆上采样六点，选择r＝0.95，试画出它的频率采样结构，并在计算机上，用DFT求出频率采样结构中的有关系数。 　　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 试分别画出它们的直接型结构和格形结构，并求出格型结构的有关参数。 　　21. 假设FIR格型网络结构的参数k1=－0.08, k2=0.217, k3=1.0, k4=0.5，求系统的系统函数并画出FIR直接型结构。 　　22. 假设系统的系统函数为   要求: 　　（1） 画出系统的直接型结构以及描述系统的差分方程; （2） 画出相应的格形结构，并求出它的系数; （3） 判断系统是否是最小相位。 图5.7.3 例5.7.1图 　　略去由全零点格型网络结构转换到FIR直接型网络结构的公式，如需要了解该内容，请参考文献[19]。 　　实际上，调用MATLAB函数实现直接型网络结构与格型网络结构之间的相互转换非常容易。tf2latc实现直接型到格型结构变换，latc2tf 实现格型到直接结型结构变换。 　　K=tf2latc(hn): 求FIR格型结构的系数向量K=[k1, k2, …, kN], hn为FIR滤波器的单位脉冲响应向量，并关于hn(1)=h(0)归一化。应当注意，当FIR系统函数在单位圆上有零极点时，可能发生转换错误。 　　hn=latc2tf(K) 将FIR格型结构转换为FIR直接型结构。K为FIR格型结构的系数向量K，hn为FIR滤波器的单位脉冲响应向量，即FIR直接型结构系数向量。显然，该函数可以用于求格型结构的系统函数的系数。 　　例 5.7.1的求解程序如下： 　　hn=[1, -0.9, 0.64, -0.576]; 　　K=tf2latc(hn) 运行结果： 　　K=[-0.6728 0.1820 -0.5760] 与上面的递推结果相同。 5.7.2 全极点格型网络结构 　　全极点IIR系统的系统函数用下式表示:  　　　　　　　　　　　 (5.7.14) 　　　　　　　　　　　　　　　　　　　　(5.7.15) 　 式中, A(z)是FIR系统，因此全极点IIR系统H(z)是FIR系统A(z)的逆系统。下面先介绍如何将H(z)变成A(z)。 　　假设系统的输入和输出分别用x(n)、y(n)表示，由 (5.7.14) 式得到全极点IIR滤波器的差分方程为 如果将x(n)、y(n)的作用相互交换，差分方程则变成下式:   　　　　　　　　　　　　　　　　　　　　　(5.7.17)  观察上式，它描述的是具有系统函数H(z)=A(z)的FIR系统，而(5.7.16) 式描述的是H(z)=1/A(z)的IIR系统。按照(5.7.16) 式描述的全极点直接型结构如图5.7.4所示。 (5.7.16) 则 图5.7.4 全极点IIR系统的直接型结构 　　基于上面的事实，我们将FIR格型结构通过交换公式中的输入输出作用，形成它的逆系统，即全极点格型IIR系统。重新定义输入输出  再将FIR格型结构的基本公式(5.7.1)、(5.7.2)重写如下： （5.7.18） （5.7.19） def def 　　由于重新定义了输入输出，将el(n)按降序运算，rl(n)不变，即  （5.7.20） （5.7.21） （5.7.22） （5.7.23） 　　按照上面四个方程画出它的结构如图5.7.5所示。为了说明这是一个全极点IIR系统，令N=1, 得到方程为 （5.7.24） （5.7.25) （5.7.26） （5.7.27） 图5.7.5 全极点IIR格型结构 当x(n)和y(n)分别作为输入和输出时，（5.7.27）式就是一个全极点的差分方程，由（5.7.24）～（5.7.27）式描述的结构就是一阶的单极点格型网络，如图5.7.6(a)所示。如果N＝2，可得到下面方程组: （5.7.28） （5.7.29） （5.7.30）