切比雪夫1型滤波器

1、课题描述设计原理2.1滤波器的分类2.22.3模拟滤波器的设计指标切比雪夫1型滤波器目录 1 22.3.1切比雪夫1型滤波器的设计原理2.3.2切比雪夫1型滤波器的设计步骤脉冲响应不变法3.1脉冲响应不变法原理4设计内容*3*364.1设计步骤610.4.2用MATLAB编程实现64.3设计结果分析5总结106参考文献001 课题描述数字滤波器是数字信号处理的重要工具之一，它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序, 数字滤波器处理精度较高，体积小，稳定，重量轻，灵活，不存在阻抗匹配问 题，可以实现模拟滤波器无法实现的特殊功能。故本书课题使用MAT

2、LAB信号 处理箱和运算用切比雪夫法设计数字低通滤波器。利用脉冲响应不变法设计切比雪夫I数字低通滤波器，通带截止频率lOOhz，阻带截止频率150Hz，采样频率lOOOhz，通带最大衰减为0.5dB，阻带最小衰减为1OdB，画出幅频、相频响应曲线，并设计信号验证滤波器设 计的正确性。设计原理2.1.滤波器的分类从功能上分；低、带、高、带阻。从实现方法上分：FIR、IIR 从设计方法上来分：Chebyshev(切比雪夫)Butterworth (巴特沃斯)从处理信号分：经典滤波器、现代滤波器2.2模拟滤波器的设计指标设ha(j?)是一个模拟滤波器的频率响应，则基于平方幅度响应J (?)=Ha(j

3、 ?)的低通滤波器技术指标为：O< Ha (j?) I <1/A2, ? s<l ? I其中&为通带波动系数，P?和S?是通带和阻带边缘频率。 A为阻带衰减系数从图知必须满足其中参数£和A是数字滤波器指标2.3切比雪夫1型滤波器2.3.1切比雪夫滤波器介绍在巴特沃兹滤波器中，幅度响应在通带和阻带内都是单调的。因此，若滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话，那么，在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径是使逼近误差均匀地分布于通带或阻带内， 或同时在通带和阻带内都均匀分布，这样往往可以降低所要求的滤波器阶次。

4、通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质：其频率响应的幅度既可以在通带中是等波纹的，而在阻带中是单调的（称为I型切比雪夫滤波器），也可以在通带中是单调的，而在阻带中是等波纹的（称为II型切比雪夫滤波器）。丨型切比雪夫滤波器的幅度平方函数是T叱）匸严萌?其中I。而冋曲=E 是滤波器在截止频率如勺放大率(注意：常用的以幅度下降3分贝的频率点作为截止频率的定义不适用 于切比雪夫滤波器!) 几(彩)是和阶切比雪夫多项式:IL I 1 = cos (n arccos ) : 0 < ct； <3$cosh 1 n aiccosh 兰 J

5、; w > 地VSo 7从定义切比雪夫多项式可以直接得出由Tn (x)和Tn-1( x)求Tn+1(x)的递推公式。将三角恒等式得Tn+1 (x) =2 (x) T (x) -Tn-1 (x)当0VXV1 时，TnA2(x)在0和1之间变化；当x>1时,cos -是虚数,所以Tn(x)像双曲余弦一样单调地增加；IHn(w) I对于0<w/w0 <1呈现出 在1和1/ (2占2 )之间的波动；而对于 w/w01单调地减小。需要用三个参量来确定该滤波器：£, w0和N。在典型的设计中，用容许的通带波纹来确定£,而用希望的通带截止频率来确定 w0。然后选择

6、合适的阶次n,以便阻带的技术要求得到满足。2.3.2切比雪夫1型滤波器的设计原理切比雪夫滤波器(又译车比雪夫滤波器)是在通带或阻带上频率响应幅度等波纹 波动的滤波器。在通带波动的为“ I型切比雪夫滤波器”，在阻带波动的为“ II型切比雪夫滤波器”。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但 频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲 线之间的误差最小，但是在通频带内存在幅度波动。切比雪夫I型为例介绍其设计原理幅度平方函数用I H(j ) I表示&为小于1的正数，表示通带内幅度波动的程度，&越大，波动幅度也越大。？ P称为通带截止频率。令入=?

7、/?p，称为对？ P的归一化频率。 定义允 许的通带内最大衰减a P用下式表示定义允许的通带内最大衰减 P a用下式表示ap=10lg max| Ha(j ? ) |八2/ minI Ha(j ?) I a2 I ? I <?p式中 max I Ha(j ? ) | 八2=1min I Ha(j ?) I 八2=1心+£A2)因此ap=10lg(1+ £2)£A2=10A(0.1卯)-1这样，可以根据通带内最大衰减a P，可求出参数£。阶数N影响过渡带的宽度,同时也影响通带内波动的疏密，因为N等于通带内的最大值和最小值的总个数。设阻带的起点频率为？

8、 s，则有2.3.3切比雪夫低通滤波器的设计步骤(1) 确定低通滤波器的技术指标：边带频率P,它们满足(2) 求滤波器阶数N和参数先求出1/k1 ,再求出阶数N，取符合条件的N的最小整数。(3) 求归一化系统函数 Ga(P)(4) 将去归一化，得到实际的Ha (s)函数说明： 用MATLAB设计切比雪夫1型低通滤波器(1) Z,p,k=cheb1ap(N,Rs)该格式用于计算N阶切比雪夫1型归一化模拟滤波器的零极点和增益因子，返回长度为N的列向量Z和P，分别给出N个零点和极点的位置。Rs是阻 带最小衰减。(2) N,wso=cheb1ord(wp,ws,Rp,As)该格式用于计算切比雪夫1型数

9、字滤波器的阶数N和阻带截止频率WSO , 调用参数分别为数字滤波器的通带频率和阻带边界频率的归一化值。(3) N,wso=cheb1ord(wp,ws,Rp,As)该格式用于计算切比雪夫1型模拟滤波器的阶数N和阻带频率wso °wp.ws是实际模拟角频率。(4) B,A=cheby1(N,Rs,wso, ' ftype ')该格式用于计算N阶切比雪夫1型数字滤波器系统函数的分子和分母的多 项式系数向量B和A。调用参数N和wso分别为切比雪夫1型数字滤波器的阶 数和阻带截止频率的归一化值。(5) B,A=cheby1(N,Rp,wso,' ftype '

10、， s')该格式用于计算N阶切比雪夫1型模拟滤波器系统函数的分子和分母多项 式系数向量B和A。调用参数N和wso分别为切比雪夫1型模拟滤波器的阶数 和阻带截止频率的归一化值。Impinvar功能：用“脉冲响应不变法”将模拟滤波器离散化。BZ,AZ =impinvar(BAFs)Plot(X,Y)功能：绘制线性二维图形 grid on功能：对当前坐标系添加主要网格线。3 脉冲响应不变法脉冲响应不变法原理脉冲响应不变法是实现模拟滤波器数字化的一种直观而常用的方法。它特别 适合于对滤波器的时域特性有一定要求的场合。脉冲响应不变法的核心原理是通过对连续函数ha(t)等间隔采样得到离散序列ha(

11、nT)。令h(n)= ha(nT)，T为采样间隔。它是一种时域上的转换方法。一个模拟滤波器的传递函数可以用有理分式表达式表示为:M./c a。ais . aM s H a®-: Nb。bs . bNS(3-1)通过反拉普拉斯变换我们就可以得到它的冲激相应:ha(t) L 1H a(s)(3-2)脉冲响应不变法就是要保证脉冲响应不变，即:h(n) ha( nT)(3-3)对上式的冲激相应序列h(n)作Z变换，就可以得到数字滤波器的传递函数:H(z) ZTha( nT)(3-4)一般来说，H a(s)的分母多项式阶次总是大于分子多项式的阶次。假定H a(s)的没有多重极点，则式(1)就可

12、分解为:Ha(S)(3-5)式中S，A均为复数，Si是H a(s)的极点。其拉氏变换为NSFha(t)Ae u(n)i 1一般来说，Ha(S)的分母多项式阶次总是大于分子多项式的阶次。数字滤波器的传递函数H(z)经过合并简化，成为一般形式的有理分式传递函数1Co Gzdo Hz 1ClZ LN.d N z(3-6)在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中，我们已经知道1H(Z)|zesT-I K2Ha(s jyK)(3-7)按照z= eST的关系,每一个s平面上宽度为2 /T的水平条带将重迭映射到z平面上。因此脉冲响应不变法将 s平面映射到z平面，不是一个简单的对应的关系。对

13、于高采样频率(T小)的情况，数字滤波器在频域可能有极高的增益。为此我们采用H(z)(3-8)Nh(n) Tha( nT)TAeSITnu( n)i 1(3-9)在脉冲响应不变法设计中，模拟频率与数字频率之间的转换关系是线性的(*t)。同时，它可以保持脉冲响应不变，h(t) = ha(nT)。因此，这一方法往往用于低通时域数字滤波器设计及相应的模拟系统数字仿真设计。4实验内容4.1设计步骤(1)按一定规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标;根据转换后的技术指标使用滤波器阶数函数，确定滤波器的最小阶数N和截止频率Wc;(3)利用最小阶数N产生模拟切比雪夫低通滤波原型。(4)利用冲

14、激响应不变法或双线性不变法把模拟切比雪夫滤波器转换成数字切比雪夫低通滤波器。4.2试验程序低通滤波器程序fp=100;fs=150;Rp=0.5;Rs=10;Fs=1000;Wp=2* pi *fp;Ws=2* pi*fs;N,wc=cheb1ord(Wp,Ws,R p,Rs,'s');%计算阶数N和截止频率wcwcz,p ,k=cheb1a p(N,R p);%计算系统函数零极点和增益Ba p,Aa p=z p2tf(z, p,k);%将系统函数的零极点转化为系统函数一般形式的系数b,a=l p2lp (Ba p,Aap ,wc);%将模拟低通滤波器原型，转换为模拟低通滤波器

15、b z, az=i mpinv ar(b,a,Fs);%采用脉冲响应不变法，将模拟低通滤波器，转换为数字低通滤波器figure(1) sub plot(2,1,1), plot(W/pi,20*log10(abs(H) xIabelC 幅频/Hz');ylabel('幅度/dB');grid on;H,W=freqz(b z,az,Fs);%计算频率H和频率响应Wsub plot(2,1,2), plot(W/pi,180/pi*u nwrap(an gle(H)xlabel('相频/Hz');ylabel('相位/o');grid on

16、;验证滤波器的正确性程序figure(2);f1=50;n=150;m=0:( n-1);t=m/Fs;%采样点数sub plot(2,3,1)x=si n(2* pi*f1*t);plot(t,x);title('f1输入信号');f2=200;x=si n(2*pi *f2*t);sub plot(2,3,2)plot(t,x);title('f2输入信号');sub plot(2,3,3) x=si n( 2* pi *f1*t)+si n(2* pi*f2*t);plot(t,x);title('叠加后输入信号');y=filtfilt(

17、bz,az,x);sub plot(2,3,4) plot(t,y);title('滤波之后的信号') sub plot(2,3,5) plot(t,sin(2*pi*f1*t);title(' 希望输出的信号');4.3实验结果分析求得阶数N=3 ,通带边界频率 wc=628.3125 ,频率在 0-100hz (0-0.2)完全通过，在100-150hz (0.2 -0.3 )通过幅度递减，150hz (0.3 - 后被截止。验证信号f1=50hz,f2=200hz，观察验证后的波形，可以看出f1通过，而f2被截止，数字低通滤波器的阻带截止频率为150hz，所以设计的数字低通滤波器符合技术指标5实验总结与体会通过对数字信号处理的课程设计， 熟悉了 MATLAB的运行环境，了解更多有关于Matlab软件的知识，初步掌握了 MATLAB语言在数字信号处理中一些基本库函数的调用和编写基本程序等应用；通过对数字低通滤波器的设计让我熟悉了滤波器设计的一般原理，对滤波器有了一个感性的认识；学会了数字低通滤波 器设计的一般步骤；加深了对脉冲响应不变法的理解和认识。在设计的