全新IIR无动画.pptVIP

主要内容 2)按频率特性分类 理想滤波器的频率响应 数字滤波器的分类 IIR滤波器 ：无限长单位脉冲响应滤波器 (Infinite Impulse Response) 系统的单位脉冲响应h(n)延伸到无穷长 数字滤波器的分类 FIR滤波器 ：有限长单位脉冲响应滤波器(Finite Impulse Response) 系统的单位脉冲响应h(n)是一个有限长序列 数字滤波器的分类 数字滤波器的分类 数字滤波器的技术要求 1.设计数字滤波器的任务 寻求一个因果稳定的线性时不变系统，使其系统函数 H(z)具有指定的频率特性。 2.分析(通常为选频滤波器) 数字滤波器的性能要求（低通为例） 数字滤波器的性能要求（低通为例） 设计IIR滤波器的方法 IIR数字滤波器的系统函数可表示为有理分式 两种设计方法： 模拟滤波器的设计 模拟低通滤波器的设计思路 巴特沃思低通滤波器 巴特沃斯模拟低通滤波器的设计 用脉冲响应不变法设计IIR数字低通滤波器 脉冲响应不变法优点 用双线性变换法设计IIR数字低通滤波器 脉冲响应不变法的转换性能分析 (1) H(z)的因果稳定性 如果Ha(s)因果稳定，转换后得到的H(z)仍是因果稳定 (2) 频率响应特性 根据时域采样理论，数字滤波器频率响应函数与模拟滤波器频率响应函数的关系为： 数字滤波器频率响应是模拟滤波器频率响应的周期延拓函数 如果模拟滤波器ha(t)具有带限特性，而且T 满足时域采样定理，则数字滤波器频率响应完全模仿了模拟滤波器频率响应 脉冲响应不变法的转换性能分析 1）频率坐标变换是线性的，即ω=ΩT ，故这种方法会很好地重现原模拟滤波器的频率特性； 2) 脉冲响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好。 如果模拟滤波器ha(t)不是带限于-π/T~π/T之间，则会在π/T奇数倍附近产生频谱混叠，从而映射到z平面上，在ω=± π 附近产生频谱混叠。因此脉冲响应不变法不适合用于设计高通滤波器和带阻滤波器。 脉冲响应不变法存在的问题 脉冲响应不变法从Ha(s)转换成H(z)的设计步骤 (1) 将模拟滤波器Ha(s)写成部分分式形式，求出极点si (2) 求H(z)的极点 (3) 脉冲响应不变法直接从Ha(s)转换成H(z)的推导思路 【例4】已知模拟滤波器的系统函数 ，用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z) H(z)的极点为 脉冲响应不变法的MATLAB实现 [Bz,Az]=impinvar(B,A,Fs); Fs=1;T=1/Fs; wp=0.2*pi/T;ws=0.35*pi/T;rp=1;rs=10; [N1,wc1]=buttord(wp,ws,rp,rs,s); [B1,A1]=butter(N1,wc1,s); %绘制低通滤波器的损耗函数曲线 fk=0:Fs/2/512:Fs/2; wk=2*pi*fk; Hk=freqs(B1,A1,wk); %求模拟滤波器的频率响应 subplot(2,2,1); plot(fk,20*log10(abs(Hk))); grid on xlabel(f(Hz));ylabel(幅度(dB)) axis([0,0.5,-50,5]);title((a) 模拟滤波器（T=1s） ); [Bz1,Az1]=impinvar(B1,A1,Fs); Hz=freqz(Bz1,Az1,wk); subplot(2,2,2); plot(wk/pi,20*log10(abs(Hz)));grid on; %xlabel(\omega/\pi);ylabel(幅度(dB)) xlabel(w/π);ylabel(幅度(dB)); axis([0,1,-50,5]);title((b) 数字滤波器（T=1s） ); 【例5】用脉冲响应不变法设计低通巴特沃斯数字滤波器，其指标参数为 B1=[0 0 0 0 0.4872] A1=[1 2.1832 2.3832 1.524 0.4872] Bz1=[0 0.0456 0.1027 0.0154 0] Az1=[1 -1.9184 1.6546 -0.6853 0.1127] 一、双线性变换法的基本思想 非线性压缩到 ，再用 转换到z平面 双线性变换法的推导思路 令 则 令