第7章 无限脉冲响应数字滤波器的的设计.ppt

* * 设IIR滤波器由K个二阶网络级联而成，系统函数用H(z)表示 6-86 式中，A是常数；ai、bi、ci、di是待求的系数；而在(6-84)式中共有(4k+1)个待定的系数，求它们的原则是使E最小。 按照(6-83)式，E是(4k+1)个未知数的函数，用下式表示 上式θ表示4k个系数组成的系数向量。 * * 令 那么 为选择A使E最小，令 6-88 6-87 * * 这里只考虑幅度误差，不考虑A的符号，将Ag作为常数代入(6-85)式中。然后将E(θ,A)对4K个系数分别求偏导，令其等于零，共有4K个方程，可以解4K个未知数。 设θk是θ的第k个分量(ak或bk或ck或dk) 因为 式中 表示对Hi函数取共轭。所以 6-89 6-90 * * 将上式具体写成对ak、bk、ck、dk的偏导，得到 式中，k=1,2,…,K；i=1,2,…,N。同理求得 6-91 6-92 6-93 6-94 * * 由于系统函数是一个有理函数，极、零点均以共轭成对的形式存在，对于极点z1，一定有下面关系 上式表明，极点z1和它的共轭极点 ，均以其倒数 和 代替后，幅度特性的形状不变化，仅是幅度的增益变化了|z1|2。一般极点这样搬移以后，需要继续进行前面的优化迭代算法。 * * 例6-17 设计低通数字滤波器，其幅度特性如图6-44(a)所示。截止频率ωs=0.1πrad。 (a)要求的幅度特性 (b)k=1、k=2时的幅度特性 图6-44 例6-17图 * * 解：考虑到通带和过渡带的重要，在ω=0.2π区间，每隔0.01π取一点ωi值，在0.2π~π区间每隔0.1π取一点ωi值，并增加一点过渡带，在ω=0.1π处|Hd(ejω)|=0.5。 N=29，取k=1，系统函数为 * * 待求的参数是A，a1，b1，c1，d1。设初始值θ=(0,0,0,-0.25)T，经过90次迭代，求得E=1.2611，系统函数零、极点位置为 零点 0.678 344 30±j0.734 744 18 极点 0.756 777 93±j1.321 391 6为使滤波器因果稳定，将极点按其倒数搬入单位圆内，再进行62次优化迭化，求得结果为 零点 0.821 911 63±j0.569 615 01 极点 0.891 763 90±j0.191 810 84 Ag=0.117 339 78，E=0.567 31 * * 其幅度特性如图6-44(b)实线所示。k=2时幅度特性如图6-44(b)虚线所示。该图表明k=2比k=1幅度特性改善了，且幅度平方误差E也小了。 6-95 误差函数用下式表示 * * 在时域直接设计IIR数字滤波器 1．帕德逼近法 利用时域办法来设计一个IIR滤波器，使它的单位抽样响应逼近所要求的期望滤波器单位抽样响应。 设滤波器是因果性的(物理上可实现)，其系统函数为 6-96 帕德逼近法是找M+N+1个系数ai,bi，以使得在0k≤p-1的范围内，使h(k)逼近一个所希望的单位抽样响应hd(k) * * 已证明，可以找到一组ai，bi，使得以下的加权均方误差最小 6-97 式中ω(k)是误差序列的正加权函数，因为h(k)是滤波器参量ai，bi的非线性函数，故只能用迭代法来得到ε的最小值。当P=N+M+1时，解(N+M+1)个线性方程，可以得到使ε为极小的滤波器参量。 设a0=1，(6-94)式可写成 代入P=M+N+1，则有 * * 令等式两边z的同次幂的系数相等，即可得M+N+1个方程 … 由于iM时，bi=0，因而一般表达式可分为两种表示(一般情况下MN) 0≤k≤M M k≤M+N=p-1 6-98 6-99 * * 因为希望h(k)充分逼近hd(k)，所以用hd(k)代替h(k)，即令h(k)=hd(k)(k=0,1,…, M+N)，这样来求解(6-96)式及(6-97)式，得到N个ai；及M +1个bi；即⑴令hd(k)=h(k)(k= M+1, M+2,…, M+N)，利用(6-97)式求得N个系数ai(i=1,2,…, N)；⑵求得ai后，令hd(k)=h(k) (k=0,1,…, M)，利用(6-96)式求得M +1个系数bi(i=1,2,…, M)。由于是无限长脉冲响应系统，一定满足M≤N，故方程组可以求解，各ai，bi得出后，也就确定了系统函数H(z)。 * * 可以看出