资源描述：
《IIR滤波器,kalman滤波器综合设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、基于MATLAB的IIR数字带通滤波器的设计设计一:利用冲击响应不变法设计一个切比雪夫带通滤波器,设计的指标为:通带边缘频率w1=0.4pi,w2=0.8pi阻带边缘频率w3=0.3pi,w4=1.0pi通带波纹=0.5dB阻带衰减=20dB解:1.基本原理滤波器的传递函数的一般形式为:n当M<=N,N阶IIR系统;M>=N,H(z):N阶IIR系统+(M-N)阶的FIR系统,n以上两种表示等价，部分分式形式和零极点增益形式nIIR系统的逼近,就是找到滤波器的系数ak,bk，或者是系统的零极点和增益(z,p,k)。切比雪夫I型的幅度平方函数为:的特点如下:(1)当=0
2、时,N为偶数时,=,当N为奇数时，=0.(2)当=c时，即此时所有的幅度函数曲线都经过点,c即为切比雪夫滤波器的通带截止频率.(3)在通带内,在1之间等波纹地起伏。(4)在通带外,随着的增大,迅速单调的趋近于零。该滤波器在通带内具有等波纹起伏特性,在阻带内则单调下降且具有更大的衰减.相比于巴特沃斯滤波器,阶数N较小.2冲击响应不变法冲击响应不变法是使数字滤波器的单位冲击响应序列h(n)模仿模拟滤波器的单位冲击响应ha(t)。将模拟滤波器的单位冲击响应加以等间隔抽样,使h(n)正好等于ha(t)的抽样值，即满足:h(n)=ha(nT)既有:冲击响应不变法是将模拟滤波器的
3、s平面变换成数字滤波器的z平面。冲击响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好。一个线性相位的模拟滤波器通过脉冲响应不变法得到的仍然是一个线性相位的数字滤波器。但是，由于任何一个实际的模拟滤波器频率响应都不是严格限带的，变换后就会产生频率响应的混叠失真。这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。3．设计流程:(1)确定数字滤波器的技术指标。(2)将数字滤
4、波器的技术指标转换成相应的模拟滤波器的技术指标。(3)按照模拟低通滤波器的技术指标设过渡模低通滤波器。(4)用脉冲响应不变法，将模拟滤波转换成数字滤波器4.matlab实现主要步骤:(1)先由数字频率转化为模拟的频率,由于采用脉冲响应不变法,所以:Omega1=fs*w1;Omega2=fs*w2;Omega3=fs*w3;Omega4=fs*w4;采样频率fs=2000HZ由函数[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,'s');[z,p,k]=cheb1ap(N,Rp);求出模拟低通滤波器。(2)再通过[bt,at]=lp2bp(bb,aa,Wo,B
5、w)函数实现由模拟低通滤波器向带通滤波器的转换.其中:Wo=sqrt(Omega1*Omega2);Bw=Omega2-Omega1;进行频率变换后实现了带通模拟滤波器(3):采用脉冲响应不变法:利用函数:[bz,az]=impinvar(bt,at,fs)将上述带通模拟滤波器转换为数字带通滤波器.其中fs为采样频率(4):验证：预期输出信号:uzs=sin(2*pi*451*nn)加入噪声后的信号为:uz=0.5*cos(2*pi*21*nn)+uzs+2*sin(2*pi*1000*nn)观察经过该滤波器后的输出波形。其中Fn=100;dn=1/fs;n=0:Fn
6、-1;nn=n*dn;5.matlab实现结果:(1)根据滤波器的设计指标,首先实现一个模拟的低通滤波器,其曲线如图一所示:图一：模拟低通滤波器(2)将设计的模拟低通滤波器进行频率变换,变换成模拟的带通滤波器,其幅频,相频曲线如图2所示:图二:模拟带通滤波器(3)冲击响应不变法实现数字带通滤波器的带通曲线如图三所示:其幅频相频曲线如图四所示:图四:带通滤波器的幅频和相频曲线(1)为了对设计的数字带通滤波器的性能进行验证,本设计加入了噪声信号,如图五所示即为原来理想的信号波形uzs=sin(2*pi*451*nn)(2)如图六所示即为加入噪声信号后的波形曲线。噪声信号的
7、曲线方程为:uz=0.5*cos(2*pi*21*nn)+uzs+2*sin(2*pi*1000*nn)(7)如图七所示即为通过滤波器后输出的波形图五:理想波形图六:噪声信号图七:滤波后输出波形Matlab程序代码:w1=0.4*pi,w2=0.8*pi;%digitalsignalfreqRp=0.5;Rs=20;w3=0.3*pi,w4=1*pi;fs=2000;Omega1=fs*w1;Omega2=fs*w2;Omega3=fs*w3;Omega4=fs*w4;f1=Omega1/2/pi;f2=Omega2/2/pi;f3=Omega3/2
