第7章 无限脉冲响应滤波器的的设计 《数字信号处理》课件.ppt

若将公式的s换成jΩ=j2πf，即可计算其幅频特性|H(f)|。其MATLAB程序为： 程序转下页 (7.58) fp=3000;fs=6000;Ap=1;As=40;%模拟切比雪夫滤波器的技术指标 r=(10^(Ap/10)-1)^(0.5);%计算波动系数 程序接上页 程序转下页 N=acosh(sqrt(10^(As/10)-1)/r)/acosh(fs/fp);%计算切比雪夫滤波器的阶 N=ceil(N);%对计算得到的阶向上取整数 Wp=2*pi*fp;k=1:N;u=(2*k-1)/2/N*pi;v=asinh(1/r)/N;%准备计算模拟系统的极点 sk=-Wp*sin(u)*sinh(v)+j*Wp*cos(u)*cosh(v);%计算模拟系统的极点 B=Wp^N/r/2^(N-1);%计算模拟系统函数的分子系数 A=poly(sk)%从系统函数的极点到分母多项式的系数，poly与roots的作用相反 f=0:0.1:6000;W=2*pi*f;%给出离散模拟频率和离散模拟角频率 程序接上页 程序运行的结果如图7.9所示。它的通带最大衰减=-20lg(0.8912)≈1.00 (dB)，阻带最小衰减=-20lg(0.0054)≈45.35 (dB)，满足技术指标。 图7.9没有采用对数纵坐标，原因是对数纵坐标表现该幅频特性的波动不明显。 图7.9 H=freqs(B,A,W);%计算模拟滤波器的频谱 plot(W/2/pi,abs(H),r);grid;xlabel(f/Hz);ylabel(|H(f)|);axis([0,6e3,0,1.1]);%画幅频特性曲线 36节 7.2 间接设计数字滤波器 间接设计数字滤波器的方法是，借助模拟系统函数和数字系统函数的数学关系，也就是s和z之间的对应关系，对模拟滤波器进行数学变换。变换方法有两种。 7.2.1 脉冲响应不变法 脉冲响应不变法的原理和模拟信号变成数字信号的原理是一样的，将单位脉冲响应当作模拟信号，这个原理就好理解了。 为了从模拟系统函数Ha(s)获得数字系统函数H(z)，首先，将模拟系统函数写成 根据拉普拉斯变换的定义式(7.4)，这个Ha(s)的拉普拉斯反变换 按照h(n)=ha(nT)离散化这个模拟系统ha(t)，即采样ha(t)，就可得数字系统 (7.83) (7.84) (7.85) 然后，求这个离散系统h(n)的z变换，得到的就是数字系统函数 它就是从Ha(s)获得H(z)的脉冲响应不变法。 观察公式(7.86)，可得数字系统的极点zi和模拟系统的极点si的对应关系： 它是计算数字系统极点的公式。 公式(7.86)的幅度|H(ω)|与|H(Ω)|还不是完全相等，观察数字频谱与模拟频谱的关系就知道。 (7.86) (7.87) 让我们观察数字频谱与模拟频谱的关系 为了满足Ha(Ω)=T×H(ω)，公式(7.86)还要乘上T，即 复变量s和z存在对应关系： 它使得采样系统Hs(s)与数字系统H(z)等价，变量z与s之间的对应关系就像地图和自然地形的对应关系，科学上将这种数学转换方法称为映射(map)。 (7.81) (7.88) (7.91) 例题7.5 检测地球物理信号时，需要滤掉被测信号中的高频噪声。假设有用信号的频谱成分分布在频率f=0~500Hz的范围，请你设计一个4阶的巴特沃斯数字滤波器，让它完成这项任务。 解 截止频率通常是指半功率点截止频率，也就是3dB截止频率。这么理解的话，本题的设计可以分两步完成。 （1）设计模拟系统函数 已知有用信号的截止角频率Ωc=2π500弧度/秒，滤波器的阶N=4，根据公式 计算模拟系统函数Ha(s)的极点， (7.18) 并写出系统函数Ha(s)的部分分式表达式 将s=jΩ和极点代入该表达式，并用计算机和数学软件MATLAB辅助，可快速获得其幅频特性。相应的MATLAB程序为： (7.99) (7.100) 程序运行结果如图7.13所示。 N=4;Wc=2*pi*500;k=1:N;g=Wc^N;z=[];%系统的阶、截止角频率、极点序号、增益、零点 s=Wc*exp(j*pi/2*(1+(2*k-1)/N));%模拟滤波器的极点 [b,a]=zp2tf(z,s,g);%系统函数的零极增益型的系数变为传递函数型的系数 f=0:2:2000;H=freqs(b,a,2*pi*f);%计算模拟滤波器的频率f和频谱 plot(f,abs(H),r,linewidth,2),axis(