资源描述:
《matlab滤波器设计》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验四IIR滤波器的实现(IIR滤波器是无限脉冲响应数字滤波器的简称)一、授课目的熟悉用双线性变换法设计IIR数字滤波器和方法。二、授课内容(一)基本概念1、数字滤波器是指输入、输出均为数字信号,通过一定运算关系,改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数字滤波器与模拟滤波器概念相同,只是信号的形式和实现滤波的方法不同。一般数字滤波器从现实的网络结构或者从单位脉冲响应分类,可以分成无限脉冲脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。2、IIR滤波器设计的主要方法是先设计低通模拟滤波器,进行频率变换,将其转换为相应的高通、带通模拟滤波器,再将模拟滤
2、波器转换为相应的数字滤波器。3、N阶模拟滤波器系统函数的一般形式为4、N阶数字模拟滤波器系统函数的一般形式为(二)求滤波器频率响应的函数1、freqs函数:求模拟滤波器的频率响应freqs(b,a,w)计算由向量w(rad/s)指定的频率点上(或W直接为采样点数)频率响应,其中b和a分别为模拟滤波器系统函数H(s)的分子与分母。freqs函数将自动绘出幅频和相频曲线。实验4-1:系统传输函数为的模拟滤波器,绘出其幅频和相频谱。a=[1,0.4,1];b=[0.2,0.3,1];w=logspace(-1,1);%产生从10-1到101之间地50个等间距点,即50个频率点freqs
3、(b,a,w);运行得幅频特性和相频特性如下:2、freqz函数:求数字滤波器的频率响应freqz(b,a,w)计算由向量w指定的频率点上模拟滤波器H(z)的频率响应,其中b和a分别为数字滤波器系统函数H(z)的分子与分母。freqz函数自动将频率点均匀设置在频率范围[0,Π]上,并将自动绘出幅频和相频曲线。实验4-2:系统传输函数为的数字滤波器,绘出其幅频和相频谱。a=[1,0.4,1];b=[0.2,0.3,1];freqz(b,a,128);运行得幅频特性和相频特性如下:三、ButterWorth(巴特沃斯)模拟和数字滤波器的设计1、直接设计(1)buttord函数:But
4、terWorth滤波器阶数的选择(即ButterWorth滤波器的最小阶数的计算)ButterWorth模拟滤波器:调用格式:[n,Wn]=buttord(Wp,Ws,Rp,Rs)ButterWorth数字滤波器:调用格式:[n,Wn]=buttord(Wp,Ws,Rp,Rs,’s’),其中’s’为变元在给定滤波器性能的情况下(通带临界频率Wp、阻带临界频率Ws、通带由最大衰减Rp和阻带内最小衰减Rs),计算ButterWorth滤波器的阶数n和截止频率Wn。(2)butter函数:ButterWorth滤波器设计ButterWorth模拟滤波器:调用格式:[b,a]=butte
5、r(n,Wn)ButterWorth数字滤波器:调用格式:[b,a]=butter(n,Wn,’s’),其中’s’为变元根据滤波器的阶数n和截止频率Wn计算ButterWorth滤波器分子和分母系数(b为分子系数的矢量形式,a为分母系数的矢量形式)。实验4-3:采用直接设计方法设计一个ButterWorth数字滤波器满足以下参数:采样频率为1Hz,通带临界频率Wp=0.2Hz、通带内衰减小于1dB(Rp=1);阻带临界频率Ws=0.3Hz、阻带内衰减大于25dB(Rs=25)。[n,Wn]=buttord(0.2,0.3,1,25);[b,a]=butter(n,Wn);freq
6、z(b,a,128);运行得幅频特性和相频特性如下:2、脉冲响应不变法设计ButterWorth(巴特沃斯)数字滤波器(设计思路:先设计模拟滤波器,再模拟滤波器转换成数字滤波器)调用格式:[bz,az]=impinvar(b,a,Fs),在给定模拟滤波器参数b、a和采样频率Fs的前提下,把模拟滤波器的参数变为近似等价的数字滤波器参数,从而设计数字滤波器。实验4-4:采用脉冲响应不变法设计上述ButterWorth(巴特沃斯)数字滤波器[n,Wn]=buttord(0.2,0.3,1,25,'s');[b,a]=butter(n,Wn,'s');[bz,az]=impinvar(b
7、,a,1);freqz(bz,az,128);运行得幅频特性和相频特性如下:3、双线性不变法设计ButterWorth(巴特沃斯)数字滤波器(设计思路:先设计模拟滤波器,再模拟滤波器转换成数字滤波器)调用格式:[bz,az]=bilinear(b,a,Fs),在给定模拟滤波器参数b,a和采样频率Fs的前提下,把模拟滤波器的参数变为近似等价的数字滤波器参数,从而设计数字滤波器。实验4-5:采用双线性不变法设计上述ButterWorth(巴特沃斯)数字滤波器[n,Wn]=buttor