用matlab编程设计切比雪夫低通滤波器

1、目录1摘要.32设计原理.42.1切比雪夫滤波器介绍 .42.2滤波器的分类.52.3模拟滤波器的设计指标 .63切比雪夫I型滤波器.73.1切比雪夫滤波器的设计原理.73.2切比雪夫滤波器的设计步骤 .103.3用matlab编程设计切比雪夫低通滤波器 .114总结.185参考文献.18摘要随着信息和数字时代的到来，数字信号处理已成为当今一门 极其重要的学科和技术领域。在现代通信系统中，由于信号中经 常混有各种复杂成分， 因此很多信号的处理都是基于滤波器而进 行的。 所以，数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为根底的，这是因为模拟 滤波器的理论和设

2、计方方法都已开展的相当成熟，且有典型的模拟滤波器供我们选择。，如巴特沃思滤波器、切比雪夫滤波器 等。本次课程设计将运用MATLAB计一个基于切比雪夫低通滤 波器，并出所设计滤波器的幅度及幅度衰减特性。关键词：滤波切比雪夫模拟低通1切比雪夫滤波器介绍在巴特沃兹滤波器中，幅度响应在通带和阻带内都是单调的。因此，假设 滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话，那么，在靠近通带低频端和阻带截止频率以上的局部都会超出技术指标。一种比拟有效的途径是使逼近误差均匀地分布丁通带或阻带内， 或同时在通带和阻带内都 均匀分布，这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波 纹特性而不是单调

3、特性的逼近方法可以实现这一点。切比雪夫型滤波器就具 有这种性质：其频率响应的幅度既可以在通带中是等波纹的，而在阻带中是单调的(称为I型切比雪夫滤波器)，也可以在通带中是单调的，而在阻带 中是等波纹的(称为II型切比雪夫滤波器)，其中切比雪夫II型滤波器乂 称为逆切比雪夫滤波器。I型切比雪夫滤波器的幅度平方函数是1|Hc(j)|2= 1；2cN(“/“c)其中&是一个小丁1的正数，它与通带波纹有关，&越大，波纹也越大,式中为N阶切比雪夫多项式，定义为2 , 、从(2.2)式我们注意到，当 0X1时，cosx是虚数，所以CN(X)像双曲余弦一样单调地增加。参考(2.1),|Hc(j

4、O)|2对丁。5皿1呈现出在1和1/ (E2)之间的波动；而对 丁1单调地减小。需要用三个参量来确定该滤波器：p和 在典型的设计中，用容许的通带波纹来确定 而用希望的通带截止频率来c。然后选择适宜的阶次N,以便阻带的技术要求得到满足。(2.1)1 、CN(X)= cos(N cosX)(2.2)当N大丁或等丁1时，从定义切比雪夫多项式可以直接得出由CN(X)和CNX)求CN+(X)的递推公式。将三角包等式代入(2.2)式，得定义允许的通带最大衰减P用下式表示:2max Ha( jQ)ap= 10lg !，GXQPmin Ha(jQ)为了求切比雪夫滤波器在椭圆上极点的位置，我们首先要这样确定，在

5、大圆和小圆上以等角度 等间隔排列的那些点：这些点对丁虚轴呈对称分 布，并且没有一个点落在虚轴上；但当N为奇数时要有一个点落在实轴上， 而当N为偶数时，就都不会落在实轴上。切比雪夫滤波器的极点落在椭圆上，起纵坐标由相应的大圆上点的纵坐标来表示， 起横坐标由相应的小圆上点的 横坐标来表示。3切比雪夫I型滤波器3.1切比雪夫低通滤波器的设计原理切比雪夫滤波器的幅频特性具有等波纹特性。它有两种形式：振幅特性 在通带内是等波纹的，在阻带内是单调递减的切比雪夫I型滤波器，振幅特 性在阻带内是等波纹的，在通带内是单调递减的切比雪夫II型滤波器，如 图所示分别画出了滤波器的幅频特性和衰减函数。CN1(X)=2

6、XCN(X) -CN项X)(2.3)确定以切比雪夫I型为例介绍其设计原理，一一 、.一、“，F2.幅度平万函数用|H(jQ)表小11 + /CN(: )(2-3-1 ) p式中，&为小丁1的正数，表示通带内幅度波动的程度，8越大，波动 幅度也越大。Qp称为通带截止频率。令A = %，称为对Qp的归一化频 率。定义允许的通带内最大衰减p用下式表示pmax |Ha(jQ)icjcap=10lg !,Q|)式中2maxha(jQ)| =1(2-3-3)21min ha(jQ) =-2(2-3-4 )1；2因此:p =10lg(12)2H(j)(2-3-5 )I乩IlI M)|；2=10*1(

7、2-3-6 )(2-3-13 )这样，可以根据通带内最大衰减ap,可求出参数阶数N影响过渡带的宽度，同时也影响通带内波动的疏密,因为N等丁通带内的最大值和最小值的总个数。设阻带的起点频率为Os,那么有H(j-s)1 ;2CN2(2-3-7 )令令，s =，由s1,有p1CN-s=ch Narchs=Ha(jQs)(2-3-8 )可以解出arch1|2 1Ha(jd |arch ,s(2-3-9 )= QpchLrch13dB截止频率用 d 表示,12Hajs(2-3-10 )(2-3-11 )按照2-3-1 式，有-larch1N(2-3-12 )经过一系歹0推论Ga:fc)i W得归一化系统

8、函数为；2=10*1(2-3-6 )(2-3-13 )(3)求归一化系统函数Gap去归一化的系统函数为- 1N-ps =NF；2Np-ppi 43.2切比雪夫低通滤波器的设计步骤(1)确定低通滤波器的技术指标：边带频率Qp,通带最大衰减app p阻带最大衰减阻带截至频率 d,它们满足(2)求滤波器阶数N和参数(2-3-17 )(2-3-18 )K1archk/N =archs这样，先由(2-18 )式求出 K1*,代入(2-3-19),求出阶数N,最后取大丁或等丁N的最小整数o:-p = 10lgHa三j*2(2-3-15 ):s=10lg2Ha三j-s(2-3-16 )Ha(s) = Ga(

9、p)(2-3-14 )(2-3-19 )(2-3-20 )Ga P )=NcN以以 I . P - Pii 4(4)将去归一化，得到实际的HasHa(s)=G(p)SQ(2-3-21 )3.3用MATLAB计切比雪夫低通滤波器(1) Matlab的信号处理工具箱提供了频谱分析函数：fft函数、filter函数和freqz函数。fft函数filter函数功能：利用IIR滤波器和FIR滤波器对数据进行滤波。调用格式：y=filter(b,a,x)y,zf=filter(b,a,x)y=filter(b,a,x,zi)说明：filter采用数字滤波器对数据进行滤波，其实现采用移位直接皿型结构，因而适

10、用丁IIR和FIR滤波器。滤波器的系统函数为Hbpz b2z qz生bmz1 az &z a3zyz#即滤波器系数a=a0 a1 a2 .an,b=b0 b1 .bm,输入序歹0欠量为x。这里，标准形式为a0=1，如果输入欠量a时，a0 1, WJ MATLAB将自动进行归一化系数的操作；如果a0=0,那么给出出错信息。y=filter(b,a,x)利用给定系数欠量a和b对x中的数据进行滤波，结果放入y欠量中，y的长度取max(N,M)。y=filter(b,a,x,zi)可在zi中指定x的初始状态。y,zf=filter(b,a,x)除得到欠量y夕卜，还得到x的最终状态欠量zf(2-

11、3-20 )freqz函数功能：离散时间系统的频率响应格式：h,w=freqz(b,a,n)h,f=freqz(b,a,n,Fs) h=freqz(b,a,w)h=freqz(b,a,f,Fs)freqz(b,a,n)说明：freqz用于计算数字滤波器H(Z)的频率响应函数H(ej3 )。h,w=freqz(b,a,n)可得到数字滤波器的n点复频响应值，这n个点均 匀地分布在0,兀上，并将这n个频点的频率记录在w中，相应的频响值记 录在h中。要求n为大于零的整数，最好为2的整数次籍,以便采用FFT计算， 提高速度。缺省时n =512。h,f=freqz(b,a,n,Fs)用于对H(ej 3)在

12、0,Fs/2上等间隔采样n点，采样点频率及相应频响值分别记录在f和h中。由用户指定FS(以HZ为单位)值。h=freqz(b,a,w)用于对H(ej 3 )在0,2兀上进行采样， 采样频率点由 矢量w指定。h=freqz(b,a,f,Fs)用于对H(ej 3)在0,FS上采样，采样频率点由矢 量f指定。freqz(b,a,n)用于在当前图形窗口中绘制幅频和相频特性曲线。fft函数函数功能：对信号进行傅里叶变换。格式：fft(X)fft(X,N)fft(X,DIM)或fft(X,N,DIM)说明:fft(X)是对输入信号X的离散傅里叶变换。fft(X,N)是N点的傅里叶变换，如果X少于N点那么补

13、0凑齐位数，长于N点那么截断。如果x是个矩阵，列的长度将会以同样的方式调整，fft会对每列进行傅里叶 变换，并返回一个相同维数的矩阵。fft(X,DIM)或fft(X,N,DIM)是离散傅里叶变换在DIM尺度上的应用。DIM可适应于任意维度的fft运算。(2)确定低通滤波器的技术指标：N-椭圆滤波器最小阶数;Wn为椭圆滤波器的带宽;Wp-W圆滤波器通带截止角频率;Ws-ffi圆滤波器阻带起始角频率；Rp-通带波纹(dB);Rs-阻带最小衰减(dB);设计要求是低通滤波器，需要屏蔽的是15Hz和30Hz的波形，所以可令fp=5Hz，设fs=8Hz, Rp40dB(3)连续信号的产生及采样：cle

14、arf1=5;f2=15;f3=30;N=100;fs=100;n=(0:N-1); %采样点数为Nx1=sin(2*pi*f1*n/fs);x2=sin(2*pi*f2*n/fs);x3=sin(2*pi*f3*n/fs); %模拟信号转化为数字信号x=x1+x2+x3; %信号叠加plot(n,x); %作出时间-幅值图像xlabel(时间(s);ylabel(幅值);grid;连续信号仿真图时间(s)(4)低通滤波器的设计Wp=8;Ws=12;Rp=2;Rs=40;%S置指标参数N,wp1=cheb1ord(Wp/(fs/2),Ws/(fs/2),Rp,Rs);%利用cheblord函

15、数求滤波器最小阶数b,a=cheby1(N,Rp,wp1); %计算切比雪夫I型模拟低通滤波器系统函 数系数H,f=freqz(b,a,512,100); %求幅值H,频率fplot(f,abs(H);切比雪夫I型模拟低通滤波器仿真图17nr11rrT11r0.90.6/ u/ V V0,7.0.6-D.5-一0 4-0.3-0.2-0.1 -Vn1J|111H10510152。253。35404550(5)对滤波后的信号进行分析和变换sf=filter(b,a,x);%对原信号进行滤波plot(n,sf);xlabel(时间(s);ylabel(幅值);axis(0 1-1 1);X=fft

16、(x,512); %滤波前信号的傅里叶变换SF=fft(sf,512); %滤波后信号的傅里叶变换f=(0:255)/256*(fs/2);plot(f,abs(X(1:256) SF(1:256); % xlabel(频率(Hz);ylabel(傅立叶变换图);滤波前后信号图像比照grid;legend(before,after)信号通过椭圆低通滤波器的仿真图注：图中蓝色曲线代表滤波前的幅频曲线，绿线代表滤波后的幅频曲线3.4设计结果分析(1)求得阶数N=通带边界频率为wp = 0.1600 ,说明阶数越大，模 拟低通滤波器越理想，且接近矩形。(2)从图上可以看出：切比雪夫I型滤波器是在通带频率响应幅度等波 纹的滤波器，通带内具有等波纹起伏特性，阻带内单调下降且衰减更大:幅 度在低频处通过，随着频率的增加，逐渐衰减至零；综上所看：本设计符合 要求。(3)切比雪夫滤波器的特点是在通带内，具有相等的波纹。截频衰减陡度比同阶数巴特沃斯特性更陡度比同阶数程时的衰减就超过6NdB在阶数N一定时，波纹越大，截频衰减陡度越陡。相位响应也是非线性，但较之比巴 特沃斯为差。总结经过了将近一个周的课程设计， 我学会的如何使用matlab来进行滤波器的 设计，通过这次课程设计我不仅学到了许多关于matlab的知识，认识到了matlab的强大功能，更