切比雪夫1型滤波器讲解Word文档下载推荐.docx

1、低、带、高、带阻。（2）从实现方法上分:FIR、IIR（3）从设计方法上来分：Chebyshev（切比雪夫）,Butterworth（巴特沃斯）（4）从处理信号分：经典滤波器、现代滤波器2.2 模拟滤波器的设计指标设ha（j）是一个模拟滤波器的频率响应，则基于平方幅度响应J （） = Ha（j）的低通滤波器技术指标为： 0Ha （j） 1/A2,s其中为通带波动系数，p和s是通带和阻带边缘频率。A为阻带衰减系数从图知必须满足 其中参数和A是数字滤波器指标2.3切比雪夫1型滤波器2.3.1 切比雪夫滤波器介绍在巴特沃兹滤波器中，幅度响应在通带和阻带内都是单调的。因此，若滤波器的技术要求是用最大通

2、带和阻带的逼近误差来给出的话，那么，在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径是使逼近误差均匀地分布于通带或阻带内，或同时在通带和阻带内都均匀分布，这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质：其频率响应的幅度既可以在通带中是等波纹的，而在阻带中是单调的（称为I型切比雪夫滤波器），也可以在通带中是单调的，而在阻带中是等波纹的（称为II型切比雪夫滤波器）。I型切比雪夫滤波器的幅度平方函数是 其中1。而 是滤波器在截止频率的放大率 （注意: 常用的以幅度下降3分贝的频率点作为截

3、止频率的定义不适用于切比雪夫滤波器!） 是 阶切比雪夫多项式：从定义切比雪夫多项式可以直接得出由Tn（x）和Tn-1（x）求Tn+1（x）的递推公式。将三角恒等式得 Tn+1（x）=2（x）T（x）-Tn-1（x）当0x1时，cos是虚数，所以Tn（x）像双曲余弦一样单调地增加;Hn（w）对于0w/w01呈现出在1和1/（22）之间的波动；而对于w/w01单调地减小。需要用三个参量来确定该滤波器：，w0和N。在典型的设计中，用容许的通带波纹来确定，而用希望的通带截止频率来确定w0。然后选择合适的阶次n，以便阻带的技术要求得到满足。2.3.2切比雪夫1型滤波器的设计原理切比雪夫滤波器（又译车比雪

4、夫滤波器）是在通带或阻带上频率响应幅度等波纹波动的滤波器。在通带波动的为“I型切比雪夫滤波器”，在阻带波动的为“II型切比雪夫滤波器”。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。切比雪夫I型为例介绍其设计原理 幅度平方函数用H（j）表示 为小于1的正数，表示通带内幅度波动的程度，越大，波动幅度也越大。p称为通带截止频率。令=/p ，称为对p的归一化频率。 定义允许的通带内最大衰减p用下式表示定义允许的通带内最大衰减p用下式表示p=10lg maxHa（j ）2/ minH

5、a（j ）2 p 式中 maxHa（j ）2=1 minHa（j ）2=1/（1+2） 因此 p= 10lg（1+2） 2=10（0.1p）-1 这样，可以根据通带内最大衰减p，可求出参数。 阶数N影响过渡带的宽度，同时也影响通带内波动的疏密，因为N等于通带内的最大值和最小值的总个数。设阻带的起点频率为s，则有 2.3.3 切比雪夫低通滤波器的设计步骤（1） 确定低通滤波器的技术指标：边带频率p，它们满足（2）求滤波器阶数N和参数先求出1/k1，再求出阶数N，取符合条件的N的最小整数。（3）求归一化系统函数Ga（p）（4）将去归一化，得到实际的Ha（s） 函数说明：用MATLAB设计切比雪夫1

6、型低通滤波器（1）z,p,k=cheb1ap（N,Rs） 该格式用于计算N阶切比雪夫1型归一化模拟滤波器的零极点和增益因子，返回长度为N的列向量Z和P，分别给出N个零点和极点的位置。Rs是阻带最小衰减。（2）N,wso=cheb1ord（wp,ws,Rp,As）该格式用于计算切比雪夫1型数字滤波器的阶数N和阻带截止频率wso，调用参数分别为数字滤波器的通带频率和阻带边界频率的归一化值。 （3）N,wso=cheb1ord（wp,ws,Rp,As）该格式用于计算切比雪夫1型模拟滤波器的阶数N和阻带频率wso。wp.ws是实际模拟角频率。（4）B,A=cheby1（N,Rs,wso,ftype）该

7、格式用于计算N阶切比雪夫1型数字滤波器系统函数的分子和分母的多项式系数向量B和A。调用参数N和wso分别为切比雪夫1型数字滤波器的阶数和阻带截止频率的归一化值。（5）B,A=cheby1（N,Rp,wso,ftype,s）该格式用于计算N阶切比雪夫1型模拟滤波器系统函数的分子和分母多项式系数向量B和A。调用参数N和wso分别为切比雪夫1型模拟滤波器的阶数和阻带截止频率的归一化值。Impinvar功能：用“脉冲响应不变法”将模拟滤波器离散化。BZ,AZ =impinvar（B,A,Fs） Plot（X,Y） 功能：绘制线性二维图形grid on 对当前坐标系添加主要网格线。3 脉冲响应不变法脉冲

8、响应不变法原理脉冲响应不变法是实现模拟滤波器数字化的一种直观而常用的方法。它特别适合于对滤波器的时域特性有一定要求的场合。脉冲响应不变法的核心原理是通过对连续函数ha（t）等间隔采样得到离散序列ha（nT）。令h（n）= ha（nT） ，T为采样间隔。它是一种时域上的转换方法。一个模拟滤波器的传递函数可以用有理分式表达式表示为： （3-1）通过反拉普拉斯变换我们就可以得到它的冲激相应： （3-2）脉冲响应不变法就是要保证脉冲响应不变，即： （3-3）对上式的冲激相应序列作变换，就可以得到数字滤波器的传递函数： （3-4）一般来说，的分母多项式阶次总是大于分子多项式的阶次。假定的没有多重极点，则

9、式（1）就可分解为： （3-5）式中，均为复数，是的极点。其拉氏变换为数字滤波器的传递函数H（z）经过合并简化，成为一般形式的有理分式传递函数 （3-6）在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中，我们已经知道 （3-7）按照z=的关系，每一个s平面上宽度为/T的水平条带将重迭映射到z平面上。因此脉冲响应不变法将s平面映射到z平面，不是一个简单的一一对应的关系。对于高采样频率（T小）的情况，数字滤波器在频域可能有极高的增益。为此我们采用 （3-8） （3-9）在脉冲响应不变法设计中，模拟频率与数字频率之间的转换关系是线性的（）。同时，它可以保持脉冲响应不变， =。因此，这一方法

10、往往用于低通时域数字滤波器设计及相应的模拟系统数字仿真设计。4实验内容4.1 设计步骤（1）按一定规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标；（2）根据转换后的技术指标使用滤波器阶数函数，确定滤波器的最小阶数N和截止频率Wc；（3）利用最小阶数N产生模拟切比雪夫低通滤波原型。（4）利用冲激响应不变法或双线性不变法把模拟切比雪夫滤波器转换成数字切比雪夫低通滤波器。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,Rp,Rs,s）; %计算阶数N和截止频

11、率wcNwcz,p,k=cheb1ap（N,Rp）; %计算系统函数零极点和增益Bap,Aap=zp2tf（z,p,k）; %将系统函数的零极点转化为系统函数一般形式的系数b,a=lp2lp（Bap,Aap,wc）; %将模拟低通滤波器原型，转换为模拟低通滤波器bz,az=impinvar（b,a,Fs）; % 采用脉冲响应不变法，将模拟低通滤波器，转换为数字低通滤波器figure（1）H,W=freqz（bz,az,Fs）; %计算频率H和频率响应Wsubplot（2,1,1）,plot（W/pi,20*log10（abs（H）xlabel（幅频/Hzylabel（幅度/dBgrid on;

12、subplot（2,1,2）,plot（W/pi,180/pi*unwrap（angle（H）相频/Hz相位/o验证滤波器的正确性程序figure（2）;f1=50;n=150;m=0:（n-1）;t=m/Fs;%采样点数subplot（2,3,1）x=sin（2*pi*f1*t）;plot（t,x）;title（f1输入信号f2=200;x=sin（2*pi*f2*t）;subplot（2,3,2）f2输入信号subplot（2,3,3）x=sin（2*pi*f1*t）+sin（2*pi*f2*t）;叠加后输入信号y=filtfilt（bz,az,x）;subplot（2,3,4）plot（

13、t,y）;滤波之后的信号subplot（2,3,5）plot（t,sin（2*pi*f1*t）;希望输出的信号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软件的知识，初步掌

14、握了MATLAB语言在数字信号处理中一些基本库函数的调用和编写基本程序等应用；通过对数字低通滤波器的设计让我熟悉了滤波器设计的一般原理，对滤波器有了一个感性的认识；学会了数字低通滤波器设计的一般步骤；加深了对脉冲响应不变法的理解和认识。在设计的过程中，我也认识到了自己所学知识的不足。以前上课都是学一些基本的东西，自以为知识掌握的很熟练，实际只是了解了最表面的东西，好多理论只是不能很好的应用于实践，现在运用学到得的东西做出有实际应用价值的东西,对所学知识点进一步的理解，并进行系统化。这也让我再次认识到知识是无尽的，只有不断的充实自己、完善自己的知识理论体系，才能够更好的胜任自己以后的工作。设计过程中知识的不足也让我更加坚定了终身学习的决心。6 参考文献数字信号处理（第三版） 高西全 丁玉美MATLAB基础及应用（第三版） 于润伟 朱晓慧 MATLAB辅助现代工程数字信号处理（第二版） 李益华