带通滤波器的设计

《带通滤波器的设计》由会员分享，可在线阅读，更多相关《带通滤波器的设计（20页珍藏版）》请在装配图网上搜索。

1、带通滤波器设计作者：汤美玲陕西理工学院(物电学院)电子信息科学与技术专业2008级 陕西 汉中723000指导教师：蒋媛摘要：带通滤波器(bandpass filter)是从滤波器的特性上划分的，带通滤波器是指能通过某 一频率范围内的分量、但将其他范围的频率分量到极低水平的，与的概念相对。从实现的 网络结构或者从单位脉冲响应长度分类，可以分为无限长单位脉冲响应(IIR)滤波器和有 限长单位脉冲响应(FIR)滤波器。IIR数字滤波器的设计方法是利用模拟滤波器成熟的理 论及设计图表进行设计的，因而保留了一些经典模拟滤波器优良的幅度特性。但设计中只 考虑了幅度特性，没考虑相位特性，所设计的滤波器一般

2、是某种确定的非线性相位特性。 为了得到线性相位特性，对IIR滤波器必须另外增加相位相位校正网络，是滤波器设计变得 复杂，成本也高，又难以得到严格的线性相位特性。FIR滤波器在保证幅度特性满足技术要 求的同时，很容易做到有严格的线性相位特性。两者各有优点，择其而取之。后面的FIR 滤波器的设计中，为获得有限长单位取样响应，需要用窗函数截断无限长单位取样响应序 列。另外，在功率谱估计中也要遇到窗函数加权问题。由此可见，窗函数加权技术在数字 信号处理中的重要地位。关键词：带通滤波器，模拟，数字IIR,FIR,MATLAB软件Abstract： Bandpass filter (bandpass fi

3、lter) from the characteristics of the classification of the filter, belt filter is to point to by a frequency can within the scope of the frequency component, but will other range of frequency components to a very low level of attenuation filter, belt and the concept of elimination filter relative.

4、From the network structure or realize from the unit impulse response length classification, can be divided into an infinite long unit impulse response (IIR) filter and limited long unit impulse response (FIR filter. IIR the design of the digital filter method is to use the filter mature theory and s

5、imulation design charts for design, so keep some classic simulation filter excellent range characteristics. But design only considered the range characteristics, didnt consider phase characteristic, the design is a certain general filter nonlinear phase characteristic. In order to get the linear pha

6、se characteristic, for an additional filter must IIR phase phase correction network, is filter design complicated, the cost is high, and hard to get the strict linear phase characteristic. FIR filter in the guarantee range characteristics to meet technical requirements at the same time, very easy to

7、 do have the strict linear phase characteristic. Both have their advantages, pick the and of the take. The back of the FIR filters design, to acquire limited long unit sampling response, need to use the window function truncation infinite long unit sampling response sequence. In addition, in the pow

8、er spectrum estimation to meet a window function and weighted problem. This shows, window function weighted technology in the digital signal processing to the important position.Key words： Bandpass filter, simulation , digital , IIR , FIR , MATLAB software一.任务1. 基于IIR模拟带通滤波器的设计.2. 基于IIR数字带通滤波器的设计.3.

9、 基于窗函数的FIR带通滤波器的设计二要求1. 基本要求可显示任何汉字字符.可实现花样显示.2. 发挥部分不需要使用专门的字模软件提取固定汉字字模.可人性化设置.三. 说明时间要求：11月12日到11月24日.完成matlab设计程序、仿真，总结报告.四. 带通滤波器的设计原理、指标及方法步骤1. 带通滤波器的设计原理及窗函数法的MATLAB设计函数简介.带通滤波器的设计原理一个理想的滤波器应该有一个完全平坦的通带，例如在通带内没有增益或者衰减，并 且在通带之外所有频率都被完全衰减掉，另外，通带外的转换在极小的频率范围完成。实 际上，并不存在理想的带通滤波器。滤波器并不能够将期望频率范围外的所

10、有频率完全衰 减掉，尤其是在所要的通带外还有一个被衰减但是没有被隔离的范围。这通常称为滤波器 的滚降现象，并且使用每十的衰减幅度dB来表示。通常，滤波器的设计尽量保证滚降范围 越窄越好，这样滤波器的性能就与设计更加接近。然而，随着滚降范围越来越小，通带就 变得不再平坦一开始出现“波纹，。这种现象在通带的边缘处尤其明显，这种效应称为。 除了电子学和信号处理领域之外，带通滤波器应用的一个例子是在大气科学领域，很常见 的例子是使用带通滤波器过滤最近3到10天时间范围内的天气数据，这样在数据域中就只 保留了作为扰动的气旋。在频带较低的fl和较高的剪切频率f2之间是共振频率，这里滤波 器的增益最大，滤波

11、器的带宽就是f2和fl之间的差值。下图为理想模拟带通滤波器幅频特 性：图带通滤波器的幅频特性曲线.窗函数法的MATLAB设计函数简介(1) 窗函数法设计原理设希望逼近的滤波器频率响应函数为H*jw)，其单位脉冲响应是hd(n)。H (ejw) = h (n)e -jwnn=-s，/、1h (n) = J c H (ejw)ejwndc如果能够由已知的Hd (ejw)求出hd (n)，经过z变换可得到滤波器的系统函数。但通常以理想滤波器作为%(球)，其幅度特性逐段恒定，在边界频率处有不连续点，因而h ()是无限时宽的，且是非因果序列。为了构造一个长度为N的第一类线性相位FIR滤波器，只有将hd(

12、n)截取一段，并保证截取的一段关于n = (N-1)/2偶对称。设截取的一段用h(n)表示，即h(n) = h (n)R (n)式中，Rn(n)个矩形序列，长度为N，当a取值为(N-1)/2时，截取的一段顷) 关于n = (N -1)/2偶对称，保证所涉及的滤波器具有线性相位。我们实际所设计的滤波器的单位脉冲响应为h(n)，长度为N，其系统函数为H(z)，即H (z)=习 h(n)z -nn=0这样用一个有限长的序列h(n)去代替hn)，肯定会引起误差，表现在频域就是通常 所说的吉布斯效应。该效应引起过渡带加宽及通带和阻带内的波动，尤其使阻带的衰减减 小，从而不满足技术上的要求，这种吉布斯效应

13、是由于将h刁(n)直接截断引起的，也称为 截断效应。通过构造窗函数(n)，用来减少截断效应，从而设计一个能满足技术要求的FIR 线性相位滤波器。以上即为窗函数法设计FIR滤波器的思想。总结为以下步骤： 通过傅里叶逆变换获得理想滤波器的单位脉冲响应h, (n)。从时域出发，截取有限长的一段冲击响应作为H(z)的系数，冲击响应长度N就是系 统函数H(z)的阶数。只要N足够长，截取的方法合理，总能满足频域的要求。一般这种 时域设计、频域检验的方法要反复几个回合才能成功。要设计一个线性相位的FIR数字滤 波器，首先要求理想频率响应H (ejw)。H (ejw)是的周期函数，周期为2兀，可以展 开成傅氏

14、级数：H * jw)=工 h (n)e-jwnn=一3其中h, (n)是与理想频响对应的理想单位抽样响应序列。但不能用来作为设计FIR DF 用的h(n)，因为h (n) 一般都是无限长、非因果的，物理上无法实现。为了设计出频响类 似于理想频响的滤波器，可以考虑用h(n)来近似h (n)。d 由性能指标确定窗函数(n)和窗口长度N。设x(n)是一个长序列，(n)是长度为N的窗函数，用(n)截断x(n)，得到N点序列x (n)，在频域上则有nXN 侦)=土 fCUJeE赢 求得实际滤波器的单位脉冲响应h(n)，h(n)即为所设计FIR滤波器系数向量。h(n) = h (n) & (n)一由此可见

15、，窗函数(n)不仅仅会影响原信号x(n)在时域上的波形，而且也会影响到频域内的形状。.（2） MATLAB 窗函数MATLAB软件工具箱主要提供了以下五种窗函数，如表4-1所示为五种窗函数的定义 式，表4-2为五种窗函数的性能。表4-1 MATLAB窗函数窗函数定义式矩形窗（Boxcar）1, 0 n N -1w(n) = R (n) = N0,其他三角窗(Triang)2kn +1,1 k nw(k) =Jn +12w(k)2(n - k +1) n +1 , k nn +12海明窗(Hamming)w(k) = 0.54 - 0.46cosf 2 兀一 0 n M-1N -1J汉宁窗(Ha

16、nning)w(k) = 0.5 1 - cos 2兀工0 n M-1I In +1JJ巴特利特窗（Bartlett）|Q,1 k n+1n12w(k) = n， i、i2(k D n +1 2 , k n n 12表4-2 MATLAB窗函数窗类型旁瓣峰值主瓣峰值最小阻带衰减矩形窗13dB4n/M21dB三角窗25dB8n/M25dB汉宁窗31dB8n/M44dB海明窗41dB8n/M53dB凯塞窗57dB12n/M74dB（3）窗函数法的MATLAB设计函数简介实际设计一般用MATLAB工具箱函数。可调用工具箱函数firl实现窗函数法设计步骤。 firl是用窗函数法设计线性相位FIR数字滤

17、波器的工具箱函数，实现线性相位FIR数字滤波 器的标准窗函数法设计。fir2为任意形状幅度特性的窗函数法设计函数，用fir2设计时，可以指定任意形状的H （ejw），它实质是一种频率采样法与窗函数法的综合设计函数。主要 d用于设计幅度特性形状特殊的滤波器。用help命令查阅其调用格式及调用参数的含义。2. 带通滤波器的设计指标逼近通带：3pi，%，通带最大衰减：apdB；逼近阻带：0s I , o，疽兀,阻带最小衰减：气dB。remezord调用参数:0,1,。，口p = 52,523.带通滤波器的设计方法及步骤.模拟带通滤波器的设计步骤00 。2 。2（1）通过映射关系式人=-七一，将希望设

18、计的带通滤波器指标转换为相应的W低通圆形滤波器Q（p）的指标。（2）设计相应的归一化低通系统函数Q（P）。S2 + O 2（3）用式Hp（s） = Q（p） lp& b s。将Q（p）转换成所希望设计得到的带通滤波器 W系统函数HBPs）。.数字带通滤波器的设计步骤（1）确定数字带通滤波器的技术指标：通带边界频率op、通带最大衰减a p、阻带截止频率，、阻带最小衰减a s。（2）将数字带通滤波器的技术指标转换成相应的模拟带通滤波器的技术指标。这里主要是边界频率o p和，的转换，ap和as指标不变。如果采用脉冲响应不变法，边界频率的 转换关系为oQ=ipT1 roQssT如果采用双线性变换法，边

19、界频率的转换关系为QPQs2o=tamT 22o=tanT 2 J（3）按照模拟带通滤波器的技术指标设计过渡模拟带通滤波器(4)用所选的转换方法，将模拟滤波器H履(s)转换成数字带通滤波器系统函数H(z)。.窗函数法设计FIR带通滤波器的设计步骤(1) 根据对阻带衰减及过渡带的指标要求，选择窗函数的类型，并估计窗口长度N。(2) 构造希望逼近的频率响应函数% (ejw)，即H (ejw) = H (w )e-jw(n-i)/2所谓的“标准窗函数法”，就是选择Hd (ejw)为线性相位理想带通滤波器。(3) 计算hd(n)。如果给出待求滤波器的频率响应函数为Hd (ejw)，那么单位脉冲响 应可

20、用下式求出：h (n) = 2“H (ejw)ejwndw(4) 加窗得到设计结果：h (n) = s。%(n -。d兀(n 一 a)五. 带通滤波器的设计1. IIR和FIR滤波器的比较从性能上说，IIR滤波器系统函数的极点可位于单位圆内的任何地方，因此零点和极点 相结合，可用较低的阶数获得较高的选择性，所用的存储单元少，计算量小，所以经济高 效。但这个高效是以相位的非线性为代价的。相反，FIR滤波器却可以得到严格的线性相位， 然而由于FIR滤波器系统函数的极点固定在原点，因而只能用较高的阶数达到高的选择性； 对于同样的滤波器幅频特性指标,FIR滤波器所要求的阶数一般比IIR滤波器高510倍

21、， 成本较高，信号延时也较大。下图为IIR数字滤波器的设计方法分类框图：图IIR数字滤波器的设计方法分类从结构上看，IIR滤波器必须采用递归结构，极点位置必须在单位圆内，否则系统将不 稳定。另外，由于运算过程中对序列的舍入处理，这种有限字长效应有时会引起寄生振荡。 相反，FIR滤波器主要采用非递归结构，不论在理论上还是实际的有限精度运算中都不存在 稳定性问题，运算误差引起的输出信号噪声功率也较小。FIR滤波器可采用FFT算法实现， 在相同阶数条件下，运算速度大大提高。从设计工具看，IIR滤波器可以借助成熟模拟滤波器设计成果，因此一般都有封闭形式 的设计公式可供准确计算，计算量小，对工具要求低。

22、FIR滤波器计算通带和阻带衰减等仍 无显示表达式，其边界频率也不易精确控制FIR滤波器设计只有计算程序可循，因此对计 算工具要求高。IIR滤波器虽然设计简单，但主要适用于设计具有片段常数特性的选频型滤波器，如低 通、高通、带通、带阻等，往往脱离不了几种典型模拟滤波器的频响特性的约束。而FIR 滤波器则要灵活的多，易于适应某些特殊的应用，两者各有所长，所以在实际应用时应该 全面考虑加以选择。2.基于IIR模拟带通滤波器的设计（1）利用MATLAB设计巴特沃斯带通滤波器：利用MATLAB函数buttap（语法Z,P,K =buttap（N），乙P,K分别表示滤波器的零点， 极点，增益）设计，验证极

23、点位置是否正确。技术指标为：通带截止频率Q。=14rad /s，通带的最大衰减七=3dB，阻带截止频率Q = 4 x 104 rad / s阻带的最小衰减气=35dB。MATLAB程序如下：a=10000;b=40000;Ap=3;As=35;%技术指标T1=sqrt(10A*P)-1);T2=sqrt(10A*S)-1);N=ceil(log10(T2/T1)/log10(b/a);%求阶数Z,P,K=buttap(N);%求零、极点，增益syms rad;%定义基本符号对象hs1=K/(i*rad/a-P(1)/(i*rad/a-P(2)/(i*rad/a-P(3);hs2=10*log1

24、0(abs(hs1)A2);ezplot(hs2,-60000,60000);%绘制符号表达式的二维曲线图巴特沃斯带通滤波器频率响应曲线(2) 利用MATLAB设计切比雪夫带通滤波器：现在以MATLAB的函数cheblap(语法：z,p,k =cheblap(n, rp),其中n为阶数，rp为 通带的幅度差，z,p,k分别代表零点、极点、增益)为工具设计切比雪夫I滤波器。设计技术 指标为：通带截止频率 Qc = 104 rad / s，通带的最大衰减 七=3dB，阻带截止频率 Q = 4 x 104 rad / s，阻带的最小衰减气=35dB。程序如下：a=10000;b=40000;Ap=3

25、;As=35;%技术指标T1=sqrt(10A*AP)-1);T2=sqrt(10A*AS)-1);N=ceil(acosh(T2/T1)/acosh(b/a);%求阶数Z,P,K=cheb1ap(N,AP);%求零、极点，增益syms rad;%定义基本符号对象hs1=K/(i*rad-P(1)/(i*rad-P(2)/(i*rad-P(3);hs2=10*log10(abs(hs1)A2);ezplot(hs2,-15,15);%绘制符号表达式的二维曲线图切比雪夫I型带通滤波器响应曲线(3) 利用MATLAB设计椭圆带通滤波器：用ellipap函数设计(语法：z,p,k =ellipap(

26、n, rp, rs),其中n为阶数，rp为通带的幅 度差，rs为阻带内的最小衰减，z,p,k分别代表零点，极点，增益)。技术指标为：通带截止 频率Q( Pd)= 1rad / s，通带的最大衰减Ap (pn) = 3dB，阻带截止频率Q (sd) = 4rad / s 阻带的最小衰减 A (sn) = 35db。程序如下：pd=1;sd=4;pn=1;sn=35;%技术指标rad0=sqrt(pd*sd);%求归一化基准频率k=pd/sd;%求选择因数nd=pd/rad0;%求通带截止频率归一化值nad=sd/rad0;%求阻带截止频率归一化值q0=(1-(1-k*k*(1+(1-k*k)A/

27、2;q=q0+2*q0A5+15*q0A9+15*q0A13;b=10A*sn);e=10A*pn);n=ceil(log10(16*(b-1)/(e-1)/log10(1/q);%求阶数z,p,k=ellipap(n,pn,sn);%求零、极点，增益syms rad;%定义基本符号对象hs1=k*(i*rad/rad0-z(1)*(i*rad/rad0-z(2)/(i*rad/rad0-p(1)/(i*rad/rad0-p(2)/(i*rad/rad0-p(3);hs2=10*log10(abs(hs1)A2);ezplot(hs2,-10,10);%绘制符号表达式的二维曲线grid on;

28、%画坐标网格线仿真图形如图所示：图 椭圆带通滤波器频率响应曲线从上面的IIR带通滤波器设计特性曲线图可以看出巴特沃斯滤波器的频率特性在通带 和阻带内部都是随着频率的变化而单调下降的，在通带内部能够完成设计出超过预定指标 的要求，但会造成滤波器的阶数匕较高；切比雪夫滤波器在通带范围内具有良好的等波纹 特性，能将指标的精度均匀分布在整个通带内，或者均匀分布在阻带内，但在通带范围外 的幅频特性是单调递减的；椭圆滤波器的系统函数在截止频率的附近含有有限个零点和极 点，能够更好的逼近理想低通滤波器的特性。在三种常用滤波器的设计方法中，选择哪种 形式的滤波器可以根据具体应用情况而定。通常椭圆滤波器的阶数低

29、，切比雪夫滤波器次 之，巴特沃斯滤波器最高，而参数量化的灵敏度与此相反，巴特沃斯滤波器最佳，切比雪 夫滤波器次之，椭圆滤波器对量化最为敏感。2.基于IIR数字带通滤波器的设计用MATLAB设计一个IIR数字带通滤波器：要求：通带中心频率为 =0.5兀，通带截止频率为 =0.4兀， =0.6兀，通 p0plp2带最大衰减为a p = 3dB ；阻带最小衰减为a广15dB，阻带截止频率2= 0.7兀。设计步骤：根据要求，采用切比雪夫I型数字滤波器涉及该滤波器； 设采样频率 Fs=2000Hz,则可计算出 fp0=500Hz, fp1=400Hz, fp2=600Hz, fs1=300 Hz, fs

30、2=700Hz； 给出一测试信号，测试该带通滤波器的功能。本例调用MATLAB函数直接设计数字带通滤波器。设计程序如下：%设计数字带通滤波器程序：function y=bandp(x,fp1,fp2,fs1,fs2,rp,rs,Fs)%带通滤波%使用注意事项：通带或阻带的截止频率与采样率的选取范围是不能超过采样率的一 半，%即，fp1,fp2,fs1,fs2 的值小于 Fs/2%x：需要带通滤波的序列%fp1：通带左边界%fp2：通带右边界%fs1：衰减截止左边界%fs2：衰减截止右边界%rp:边带区衰减DB数设置%rs:截止区衰减数DB设置%Fs:序列x的采样频率fp1-400;fp2-60

31、0;fs1-300;fs2-700;rp=3;rs=15;Fs=2000;wp1=*pi;wp2=*pi;ws1=*pi;ws2=*pi;wp=wp1 wp3;ws=ws1 ws2;%设计切比雪夫滤波器；n,wn=cheblord(ws/pi,rp,rs);bz1,az1=cheby1(n,rp,wp/pi);%查看滤波器的曲线h,w=freqz(bz1,az1,256,Fs);h=20*log10(abs(h);figure;plot(w,h);title(所设计滤波器的通带曲线 );grid on;y=filter(hz1,az1,x);end%带通滤波器测试程序Fs=2000;t=(1:

32、Fs)/Fs;ff1=100;ff2=400;ff3=700;x=sin(2*pi*ff1*t)+sin(2*pi*ff2*t)+sin(2*pi*ff3*t);figure;subplot(211);plot(t,x);subplot(212);hua_fft(y,Fs,1);function hua_fft(y,style,varargin)%当style=1，画幅值谱;当style=2，画功率谱；当style=其他的，那么画幅值谱和功率 谱%当style=1时，还可以多输入2个可选参数%可选输入参数是用来控制需要查看的频率段的%第一个是需要查看的频率段起点%第一个是需要查看的频率段的终点

33、%其他style不具备可选输入参数，如果输入发生位置错误Nfft=2Anextpow2(length(y); %找出大于y的个数的最大的2的指数值(自动运算最佳 FFT 步长 nfft)%nfft=1024;%人为设置FFT的步长nffty=y-mean(y);%去除直流分量y_ft=fft(y,nfft);%对y信号进行DFT，得到频率的幅值分量 y_p=y_ft.*conj(y_ft)/nfft;%conj(涵数是求y函数的共轭复数，实数的共轭复数是它本身。y_f=fs*(0:nfft/2-1)/nfft;if style=1if nargin=3plot(y_f,2*abs(y_ft(1

34、:nfft/2)/length(y);%matlab 的帮助里画 FFT 的方法%ylabel ( 幅值);xlabel ( 频率);title (信号幅值谱)；%plot(y_f,abs(y_ft(1:nfft/2);%论坛上画 FFT 的方法 elsef1=varargin1;fn=varargin2;ni=round(f1*nfft/fs+1);na=round(fn*nfft/fs+1);plot(y_f(ni:na),abs(y_ft(ni:na)*2/nfft);endelseif style=2plot(y_f,y_p(1:nfft/2);%ylabel( 功率谱密度);xlab

35、el(频率 );title( 信号功率谱); elseSubplot(211);plot(y_ft(1:nfft/2)/length(y);ylabel(幅值 );xlabel( 频率 );title ( 信号幅值谱);Subplot(212);plot(y_f，y_p(1:nfft/2);ylabel( 功率谱密度 );xlabel( 频率);title(信号功率谱);endend仿真结果如下图为该带通滤波器频谱图形，图为经过带通滤波后的信号时及频谱图形， 图原始输入信号时域及频谱波形。J既计it芨!的旬带Hi立图带通滤波器频谱图图经带通滤波后的信号时及频谱图图原始输入信号时域及频谱波形经图

36、示可看出该滤波器具有带通特性，它将通过它的信号在带通范围内的频率成分保 留，阻带范围内的频率成分截止，从而达到滤波作用。3. 基于窗函数的FIR带通滤波器的设计(1) 利用矩形窗进行设计程序如下：fs=20000;%设定采样频率fp1=4500;fp2=6500;%第一截止频率fs1=3000;fs2=7500;%第二截止频率As=40;%最小阻带衰减Ws1=(fp1+fs1)/fs;Ws2=(fp2+fs2)/fs;%截止频率归一化处理w=(fp1-fs1)/fs;%求归一化过渡带M=ceil(/*w)%计算所需滤波器的阶数 juxing=boxcar(M+1);%生成长度为M+1的矩形窗

37、boxb=fir1(M,Ws1,Ws2,juxing);%生成矩形窗设计的fir滤波器 freqz(boxb,1,fs,fs);%绘制幅频和相频响应曲线 运行结果：曲线在通带内为直线，效果较好，信号失真小。(2) 利用三角窗进行设计利用三角窗进行设计时，原理与矩形窗基本相同，只不过生成窗函数时采用triang()函 数生成三角窗，程序运行结果如下：图三角窗fir滤波器幅频和相频响应曲线(3) 利用凯塞窗进行设计程序如下：fs=20000;%设定米样频率fp1=4500;fp2=6500;%第一截止频率fs1=3000;fs2=7500;%第二截止频率n,wn,bta,ftype=kaisero

38、rd(fs1,fp1,fp2,fs2;0,1,0, ,fs)%求滤波器参数b=fir1(n,wn,ftype,kaiser(n+1,bta);%生成fir滤波器freqz(b,1,fs,fs);%绘制幅频和相频响应曲线程序运行结果：图凯塞窗fir滤波器幅频和相频响应曲线(4) 利用汉宁窗进行设计利用汉宁窗进行设计时，原理与矩形窗基本相同，只不过生成窗函数时采用hanning () 函数生成三角窗，程序运行结果如下：图汉宁窗fir滤波器幅频和相频响应曲线(5)利用海明窗进行设计hamming利用海明窗进行设计时，原理与矩形窗基本相同，只不过生成窗函数时采用()函数生成三角窗，程序运行结果如下：由

39、MATLAB模拟生成含有不同频率的数字信号，然后利用设计的滤波器对数字信号进 行滤波，为方便观察，模拟生成的信号只含有包含在阻带的两个频率(2000Hz，8000Hz) 一个包含于通带的频率(5000Hz)。程序如下：fs=20000;t=0:1/fs:2;x=sin(2*pi*2000*t)+sin(2*pi*5000*t).+sin(2*pi*8000*t);%生成混合信号xo=filter(b,2,xn);%用滤波器对信号进行滤波figure;nn=5000:5100;%取一段信号subplot(211);tt=nn/fs;plot(tt,x(nn);%绘制原始信号axis(,-4,4)

40、;ylabel(原始信号);xlabel(时间);subplot(212);plot(tt,xo(nn);%绘制滤波后的信号axis(,);ylabel(滤波后的信号);xlabel(时间)；图滤波器滤波性能测试波形从原始型号和滤波后的信号对比可以看出，在用设计的滤波器进行滤波后信号基本成 单一频率的正弦波，滤波结果令人满意。（6）滤波器的性能测试：FIR滤波器的一个突出优点便是它的稳定性，从信号与系统的理论可知，当一个系统的 Z域的传递函数的极点都在单位圆内时系统是稳定的，FIR滤波器的传递函数的分母为1， 即FIR滤波器一定是稳定的，MATLAB中也提供了专门用于绘制零极点图的函数:zpl

41、ane（）， 通过该函数可以轻松绘制系统的零极点图，以凯泽窗设计的FIR滤波器为例，绘制其零极 点图如下：图3-3凯塞窗FIR滤波器零极点图从该零极点图可以看出，FIR设计的滤波器是属于无极点的系统，因而系统一定是稳定 的。从上面对滤波器测试可以看出FIR滤波器的特点，第一个便是其相位曲线，在通带内， 相频曲线一定为直线，二是FIR滤波器由于阶数比较高，因而有一定的时延，三是FIR滤 波器一定是稳定的。六. 测试报告在设计带通滤波器时首先要计算出过渡带，然后查表得到不同窗函数所需要的阶数， 不同的窗函数所设计的滤波器的形状各有差异，尤其在主瓣宽度、旁瓣的形状以及主瓣与 旁瓣的高度差上有比较明显

42、得差别，实际应用中应根据实际情况，折衷处理，兼顾各项指 标。七. 总结本设计是对滤波器的一个综合性的设计，首先设计了模拟带通滤波器，有模拟带通 滤波器进而设计数字带通滤波器。并且比较性的选择使用了IIR和FIR，熟练地利用窗函数 法的MATLAB工具箱完成窗函数法设计FIR带通滤波器。通过本次设计，我理解了各种滤 波器的幅频特性，熟练掌握了 MATLAB软件下如何调用工具箱函数firl和fir2来设计程序 并完成仿真。加深了对数字信号处理（DSP）的完整复习和深入理解，同时增加了对其的学 习兴趣。在设计过程中，我得到了一次用专业知识、专业技能分析和解决问题、全面系统 的锻炼的机会，对于同一个设

43、计，自己经过思考，做出多种方案，并选择其中适用的，合 理的，方案进行设计，使我在常用编程设计思路技巧的掌握方面都能向前迈了一大步。同时在老师的悉心指导和严格要求下，极大地提高了设计与分析设计的能力，同时也 感受到理论与实践之间的差距，使我们对数字信号处理的应用有了更加深刻的认识参考文献1高西全.丁玉美.数字信号处理M西安电子科技大学出版社，2008.2陈怀琛.数字信号处理教程一MATLAB释疑与实现M北京电子工业出版社，2004.3高西全.丁玉美.数字信号处理一原理、实现及应用M北京电子工业出版社，2006.4刘顺兰.吴杰.数字信号处理M西安电子科技大学出版社，2003.5刘益成.孙祥娥.数字信号处理M北京电子工业出版社，2004.6楼顺天.李博菡.基于MATLAB的系统分析与设计现一信号处理M西安电子科技大学出版社，1998.