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

取N=6。为求3dB截止频率Ωc,将Ωp和αp代入(6.2.17)式，得到Ωc=0.7032rad/s，显然此值满足通带技术要求，同时给阻带衰减留一定余量，这对防止频率混叠有一定好处。根据阶数N=6，查表6.2.1，得到归一化传输函数为为去归一化，将p=s/Ωc代入Ha(p)中，得到实际的传输函数Ha(s),第127页，共183页，星期日，2025年，2月5日④用脉冲响应不变法将Ha(s)转换成H(z)。首先将Ha(s)进行部分分式，并按照(6.3.11)式、(6.3.12)式，或者(6.3.13)式和(6.3.14)式，得到：第128页，共183页，星期日，2025年，2月5日图6.4.7例6.4.2图——用脉冲响应不变法设计的数字低通滤波器的幅度特性第129页，共183页，星期日，2025年，2月5日一般Ha(s)的极点si是一个复数，且以共轭成对的形式出现，在(6.3.1)式中将一对复数共轭极点放在一起，形成一个二阶基本节。如果模拟滤波器的二阶基本节的形式为极点为(6.3.11)可以推导出相应的数字滤波器二阶基本节(只有实数乘法)的形式为(6.3.12)第95页，共183页，星期日，2025年，2月5日如果模拟滤波器二阶基本节的形式为极点为(6.3.13)(6.3.14)第96页，共183页，星期日，2025年，2月5日例6.3.1已知模拟滤波器的传输函数Ha(s)为用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。解首先将Ha(s)写成部分分式：极点为那么H(z)的极点为第97页，共183页，星期日，2025年，2月5日按照(6.3.4)式，并经过整理，得到设T=1s时用H1(z)表示，T=0.1s时用H2(z)表示，则转换时，也可以直接按照(6.3.13),(6.3.14)式进行转换。首先将Ha(s)写成(6.3.13)式的形式，如极点s1,2=σ1±jΩ1,则第98页，共183页，星期日，2025年，2月5日再按照(6.3.14)式，H(z)为第99页，共183页，星期日，2025年，2月5日图6.3.3例6.3.1的幅度特性第100页，共183页，星期日，2025年，2月5日软件设计IIR数字滤波器一、MATLAB信号处理工具箱提供的butterworth模拟低通滤波器原型函数为：[z,p,k]=buttap(N)其中，N为滤波器的阶数；[z,p,k]为滤波器的零点、极点和增益。例1、绘制10阶butterworth模拟低通滤波器的平方幅频响应曲线。第101页，共183页，星期日，2025年，2月5日程序如下：n=0:0.01:2;N=10；[z,p,k]=buttap(N);[b,a]=zp2tf(z,p,k);%由零、极点求出系统函数中分子、分母的系数b和a[H,w]=freqs(b,a,n);%求出滤波器的幅频响应magH=(abs(H)).^2;plot(w,magH);axis([0201]);xlabel(w/wc);ylabel(H(jw));title(butterworth);第102页，共183页，星期日，2025年，2月5日第103页，共183页，星期日，2025年，2月5日二、MATLAB函数cheb1ap、cheb2ap调用格式为：[z,p,k]=cheb1ap(N,Rp)为切比雪夫I型滤波器，N为阶数，Rp为通带波纹；[z,p,k]为零点、极点和增益。第104页，共183页，星期日，2025年，2月5日例2：????绘制8阶切比雪夫I型模拟低通滤波器的平方幅频响应曲线。程序如下：n=0:0.01:2;N=8;Rp=1;[z,p,k]=cheb1ap(N,Rp);[b,a]=zp2tf(z,p,k);[H,w]=freqs(b,a,n);magH=(abs(H)).^2;subplot(211);plot(w,magH);axis([0201]);xlabel(w/wc);ylabel(H(jw));title(chebbyI);第105页，共183页，星期日，2025年，2月5日三、1、butter函数用于设计数字滤波器其调用格式