数字信号处理IIRFIR数字滤波器设计教程解析.docVIP

一、课程设计的目的与要求 1.掌握利用双线性变换法设计IIR数字滤波器的原理和具体方法。 2.深入理解利用双线性变换法设计IIR数字滤波器的优缺点及使用范围。 3.加深对窗函数设计FIR数字滤波器的基本原理的理解。 4.学习用MATLAB语言的窗函数法编写设计FIR数字滤波器的程序。 5.了解MATLAB有关窗函数法设计常用的子函数。 二、设计正文δ1 = 1-0.89125 = 0.10875 δ2 = 0.17783 2.用Kaiser窗法设计一个高通滤波器，满足技术指标∣H(ejω)∣≤0.021,∣ω∣≤0.35π 0.979≤∣H(ejω)∣≤1.021, 0.5π≤ω≤π。 用绝对指标要求来表示为：ωp = 0.5π，ωs = 0.35π，δ1 =δ2 =δ=0.021。 2.2设计原理 2.2.1．IIR数字滤波器： IIR数字滤波器的设计一般是利用目前已经很成熟的模拟滤波器的设计方法来进行设计，通常采用模拟滤波器原型有butterworth函数、chebyshev函数、bessel函数、椭圆滤波器函数等。IIR数字滤波器的设计步骤： ① 按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标； ② 根据模拟滤波器技术指标设计为响应的模拟低通滤波器； ③ 跟据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器； ④ 如果要设计的滤波器是高通、带通或带阻滤波器，则首先把它们的技术指标转化为模拟低通滤波器的技术指标，设计为数字低通滤波器，最后通过频率转换的方法来得到所要的滤波器。 在MATLAB中，经典法设计数字滤波器主要采用以下步骤： 图1.1 IIR数字滤波器设计步骤 2.2.2．FIR数字滤波器： FIR数字滤波器的系统函数无分母，为??，系统频率响应可写成：，令=，H(w)为幅度函数，称为相位函数。这与模和辐角的表示方法不同，H(w)为可为正可为负的实数，这是为了表达上的方便。如某系统频率响应 =sin4w,如果采用模和幅角的表示方法，sin4w的变号相当于在相位上加上（因-1= ），从而造成相位曲线的不连贯和表达不方便，用则连贯而方便。 窗函数法又称傅里叶级数法，其设计是在时域进行的。?一般是无限长且非因果的，设计时需用一个合适的窗函数把截成有限长的因果序列，使对应的频率响应（的傅里叶变换）尽可能好地逼近理想频率响应。窗函数法的主要缺点是：一、不容易设计预先给定截止频率的滤波器；二、满足同样设计指标的情况下所设计出的滤波器的阶数通常偏大。%03.1.6 clear,close all; %技术指标 derta1=1-0.89125; derta2=0.17738; wp=0.2*pi; ws=0.3*pi; %相对技术指标 rp=-20*log10(0.89125/1); rs=-20*log10(0.17738/1); T=1; Fs=1/T; %双线性变换确定模拟滤波器的边界频率 Op=(2/T)*tan(wp/2); Os=(2/T)*tan(ws/2); %确定阶数N和截止频率Wn [N,Wn]=buttord(Op,Os,rp,rs,s); %确定模拟滤波器 [num,den]=butter(N,Wn,s); %双线性变换获得数字滤波器的参数 [Num,Den]=bilinear(num,den,Fs); %设定显示范围，间隔 wmax=0.5*pi; w=[0:1:500]*wmax/500; %显示频率响应曲线 freqz(Num,Den,w); fprintf(\n直接型系数\n); Num Den [sos,g]=tf2sos(Num,Den); fprintf(\n级联型系数【NumG|DenG】\n); sos g %验证 %检验信号 f=sin(0.1*pi*[0:100])+sin(0.5*pi*[0:100]); figure(2); plot(f),title(输入波形); figure(3); plot(filter(Num,Den,f)),title(输出波形); 2.4.2.FIR数字滤波器： %使用kaise窗设计FIR高通滤波器 %03。1。6 clear,close all; %技术指标 ws=0.35*pi,wp=0.5*pi,derta=0.021; dw=wp-ws; A=-20*log10(derta); %确定beta和窗长M if(A50) beta=0.1103*(A-8.7), elseif(A=21A=50) beta=0.5842*(A-21)^0.4+0.07886*(A-21), else beta=0, end M=ceil((A-8)/(2.285*dw))+1; M %由使用者判断并输入