IR滤波器,kalman滤波器综合设计I.doc

基于MATLAB的IIR数字带通滤波器的设计 设计一: 利用冲击响应不变法设计一个切比雪夫带通滤波器,设计的指标为: 通带边缘频率 w1 = 0.4pi,w2 = 0.8pi 阻带边缘频率 w3 = 0.3pi,w4 = 1.0pi 通带波纹 = 0.5dB 阻带衰减 = 20 dB 解: 1.基本原理 滤波器的传递函数的一般形式为: 当M=N, N阶IIR系统; M=N, H(z): N阶IIR系统+(M-N)阶的FIR系统, 以上两种表示等价，部分分式形式和零极点增益形式 IIR系统的逼近, 就是找到滤波器的系数ak, bk，或者是系统的零极点和增益(z,p,k)。 切比雪夫I型的幅度平方函数为: 的特点如下: 当=0时,N为偶数时,=,当N为奇数时，=0. 当=c时， 即此时所有的幅度函数曲线都经过点, c即为切比雪夫滤波器的通带截止频率. 在通带内,在1之间等波纹地起伏。 在通带外,随着的增大, 迅速单调的趋近于零。 该滤波器在通带内具有等波纹起伏特性,在阻带内则单调下降且具有更大的衰减.相比于巴特沃斯滤波器,阶数N较小. 2 冲击响应不变法 冲击响应不变法是使数字滤波器的单位冲击响应序列h(n)模仿模拟滤波器的单位冲击响应ha(t)。将模拟滤波器的单位冲击响应加以等间隔抽样,使h(n)正好等于ha(t)的抽样值，即满足: h(n)= ha(nT) 既有: 冲击响应不变法是将模拟滤波器的s平面变换成数字滤波器的z平面。冲击响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好。一个线性相位的模拟滤波器通过脉冲响应不变法得到的仍然是一个线性相位的数字滤波器。 但是，由于任何一个实际的模拟滤波器频率响应都不是严格限带的，变换后就会产生频率响应的混叠失真。这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。 3．设计流程: (1)确定数字滤波器的技术指标。 (2)将数字滤波器的技术指标转换成相应的模拟滤波器的技术指标。 (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,Bw)函数实现由模拟低通滤波器向带通滤波器的转换. 其中: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-1; nn= n*dn; 5.matlab实现结果: (1)根据滤波器的设计指标,首先实现一个模拟的低通滤波器,其曲线如图一所示: 图一：模拟低通滤波器 (2) 将设计的模拟低通滤波器进行频率变换,变换成模拟的带通滤波器,其幅频,相频曲线如图2所示: 图二:模拟带通滤波器 (3)冲击响应不变法实现数字带通滤波器的带通曲线如图三所示: 其幅频相频曲线如图四所示: 图四:带通滤波器的幅频和相频曲线 为了对设计的数字带通滤波器的性能进行验证,本设计加入了