《用巴特莱特窗函数法设计数字FIR带通滤波器dsp课程设计》由会员分享，可在线阅读，更多相关《用巴特莱特窗函数法设计数字FIR带通滤波器dsp课程设计（21页珍藏版）》请在金锄头文库上搜索。
1、用巴特莱特窗函数法设计数字FIR带通滤波器dsp课程设计河击/本次摩课程设计课程设计名称：数字信号处理课程设计专业班级：通信学生姓名：学 号：指导教师：乔丽红课程设计时间：6程6-6.20电子信息工程专业课程设计任务书学生姓名专业班级学号题目用巴特莱特窗函数法设计数字FIR带通 滤波器课题 性质其他课题来源自拟指导 教师乔丽红同组姓名主要 内容用巴 FIR带通 400Hz, 5 550Hz,鬟 减 40dB MATLAE 波器传输信号 滤波器，： 器的输出 结论并给特莱特窗函数法设计一个数字 1滤波器，要求通带边界频率为 00Hz,阻带边界频率为350Hz, 费带最大衰减IdB,阻带最小衰 ，
2、抽样频率为2000Hz ,用 ；画出幅频特性，画出并分析滤 函数的零极点；、x(f) = x1 (t) + X,(r) = sin(2用I) + sin(2；r)经,过 其中力=450Hz, /;=600Hz,滤波 是什么？用Matlab验证你的 出 XQ),X2(f)，x(/),yQ)的图形。任务 要求1 .掌握用巴特莱特窗 字FIR带通滤波器的原理;2 .掌握用Kaiser方宿 滤波器长度的方法。3 .求出所设计磔器4 .用MATLAB画出所f函数法设计数 和设计方法。 必古计FIR数字:的Z变换。务频特性图并验证所设计的滤波器。1、程佩青著，数字信号处理教程,清华大学出版社，20012、
3、Sanjit K.Mitra著，孙洪，余翔宇 译，数字信号处理实验指导书参考文献（MATLAB版），电子工业出版社，2005年1月3、郭仕剑等，MATLAB 7.x数字信号处理，人民邮电出版社，2006年4、胡广书，数字信号处理理论算 法与实现，清华大学出版社，2003年指导教师签字:审查意见教研室主任签字：2014年06月12日说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页-需求分析和设计内容数字信号处理是把许多经典的理论体系作为自己的理论基础，同时又使自 己成为一系列新兴学科的理论基础。现如今随着电子设备工作频率范围的不断 扩大，电磁干扰也越来越严重，接
4、收机接收到的信号也越来越复杂。为了得到 所需要频率的信号，就需要对接收到的信号进行过滤，从而得到所需频率段的 信号，这就是滤波器的工作原理。对于传统的滤波器而言，如果滤波器的输入， 输出都是离散时间信号，则该灌波器的冲激响应也必然是离散的，这样的滤波 器定义为数字漉波器。它通过对采样数据信号进行数学运算来达到频域滤波的 目的.滤波器在功能上可分为四类，即低通(LP)、高通(HP)、带通(BP)、带阻 (BS)滤波器等，每种又有模拟滤波器(AF)和数字滤波器(DF)两种形式。对数字滤波器，从实现方法上，具有有限长冲激响应的数字滤波器被称为FIR 漉波器，具有无限长冲激响应的数字滤波器被称为HR流
5、波器。FIR数字滤波器的主要优点有：一、具有严格的线性相位特性；二、不存在 稳定性问题；三、可利用DFT来实现。这些优点使FIR数字滤波器得到了广泛 应用。窗函数法是一种设计FIR数字漉波器的基本方法，但它不是最佳设计方 法，在满足同样设计指标的情况下，用这种方法设计出的滤波器的阶数通常偏 大。在窗函数法的基础上，以所定义的逼近误差最小为准则来进行优化设计的 算法，由于其中的逼近误差可根据不同的设计要求进行定义，故此算法适应性 强，它既可用于设计选频型滤波器，又适用于非选频型滤波器的设计。常用的 窗函数有矩形窗函数、巴特莱特窗函数、三角窗函数、汉宁(Hann)窗函数、 海明(Hamming)窗
6、函数、布莱克曼(Blackman)窗函数、凯塞(Kaiser)窗函 数等。本设计通过MATLAB软件对FIR型滤波器进行理论上的实现，利用巴特莱 特窗函数设计数字FIR带通滤波器。FIR系统不像IIR系统那样易取得较好的通 带和阻带衰减特性，要取得较好的衰减特性,一般要求H (z)阶次要高,也即M要 大。FIR系统有自己突出的优点：系统总是稳定的;易实现线性相位;允许设计多 通带(或多阻带)滤波器，后两项都是IIR系统不易实现的。FIR数字滤波器的 设计方法有多种，如窗函数设计法、频率采样法和Chebyshev逼近法等。随着 Matlab软件尤其是Matlab的信号处理工作箱的不断完善,不仅数
7、字滤波器的计 算机辅助设计有了可能,而且还可以使设计达到最优化。本实验的数字滤波器的 MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定的输入信号x (n) 进行淀波，得到滤波后的输出信号y(n)。用巴特莱特窗函数法设计一个数字FIR带通滤波器，要求通带边界频率为 400Hz, 500Hz,阻带边界频率为350Hz, 550Hz,通带最大衰减IdB,阻带最小 衰减40dB,抽样频率为2000Hz,用MATLAB画出幅频特性，画出并分析滤波器 传输函数的零极点;信号 x(t) = Xj (r)+x2 (0 = sin(2|0 + sin(220 经过该滤波器，其中 /I=4
8、50Hz, f2 =600Hz,滤波器的输出y(t)是什么？用Matlab验证你的结论并 给出为Ox?。)/。),)，。)的图形。二设计原理及设计思路1.设计FIR数字滤波器的基本方法：jv-ljV-1FIR数字滤波器的系统函数无分母，为H(Z) = bi/=h(n)E：系统频 rU)/-0N-率响应可写成：(*) = Z(心5，令:(卬)、H(w)为幅度函数,夕(M，)称为相位函数。这与模和辐角的表示方法不同，H(w)为可为正可为负的实 数，这是为了表达上的方便。如某系统频率响应(e)二Sin4we-小，如果采用 模和幅角的表示方法，sin4w的变号相当于在相位上加上不(因-1二，)，从而
9、造成相位曲线的不连贯和表达不方便，用”(2好)则连贯而方便。窗函数法又 称傅里叶级数法，其设计是在时域进行的。函数一般是无限长且非因果的， 设计时需用一个合适的窗函数把它截成有限长的因果序列，使对应的频率响应 (的傅里叶变换)尽可能好地逼近理想频率响应。窗函数法的主要缺点是：一、 不容易设计预先给定截止频率的滤波器；二、满足同样设计指标的情况下所设 计出的滤波器的阶数通常偏大。一些固定窗函数的特性表名称主瓣宽度过度带宽最小阻带衰减矩形4/(2M+l)0. 92 兀/M20. 9dB巴特利 特4 乃/(2M+1)2. l/r/M25dB汉宁8/(2M+l)3.11/r/M43. 9dB海明8/(
10、2M+l)3. 32 -/M54. 5dB布莱克 曼12/(2M+1)5. 56 /M75. 3dB2. FIR数字滤波器设计的基本步骤如下:(1)确定技术指标在设计一个滤波器之前,必须首先根据工程实际的需要确定滤波器的技术 指标。在很多实际应用中，数字滤波器常被用来实现选频操作。因此，指标的形 式一般在频域中给出幅度和相位响应。幅度指标主要以2种方式给出。第一种 是绝对指标。他提供对幅度响应函数的要求，一般应用于FIR滤波器的设计。第 二种指标是相对指标。他以分贝值的形式给出要求。本文中滤波器的设计就以 线性相位FIR滤波器的设计为例。(2)逼近确定了技术指标后,就可以建立一个目标的数字滤波
11、器模型(通常采用理想 的数字滤波器模型)。之后,利用数字淀波器的设计方法(窗函数法、频率采样 法等)，设计出一个实际滤波器模型来逼近给定的目标。(3)性能分析和计算机仿真上两步的结果是得到以差分或系统函数或冲激响应描述的滤波器。根据这 个描述就可以分析其频率特性和相位特性，以验证设计结果是否满足指标要求; 或者利用计算机仿真实现设计的漉波器,再分析滤波结果来判断。三程序流程图利用kaiser方程求滤波器阶数四程序源代码clear all; clcFs=2000;fpl=400;fpu=500;fsl=350;fsu=550;wpl=2*pi*fpl/Fs; wpu=2*pi*fpu/Fs;%
12、通带截止频率 %wsl=2*pi*fsl/Fs;wsu=2*pi*fsu/Fs;% 阻带截止频率 B=wpl-wsl; %计算过渡带宽度阳=2*6“1(3.32*1-9+1;%计算所需卜(n)长度NO, ceil取大于等于WC=(wpl+wsl)/2/pi, (wpu+wsu)/2/pi ;%计算理想带通滤波器截止频率 ap=l;as=40;dp=l-10*(-ap/20);ds=10(-as/20);fenzi=-20*logl0(sqrt(dp*ds)-13;fenmu=14. 6*B/2/pi;N=ceil(fenzi/fenmu);n=0:N;hn二打30北0日16黄3+1);%调用立
13、1计算带通对1 频率/Hz);ylabel C 相位/rad);grid on;figure(2);zplane(hn, 1) ; ;%绘制零极点图titleC巴特莱特窗FIR数字带通滤波器系统函数的零极点图)；legend C 零点，极点)； grid on;n=0:1023;dt=l/10000;t=n*dt;f1=450;f2=600;x=sin(2*pi*fl*t)+sin(2*pi*f2*t);xl=sin(2*pi*fl*t);x2=sin(2*pi*f2*t);y=f liter (hn, 1, x);figure (3);subplot (3, 1, 1) ;plot (x);title C 输入信号x(t);xlabelC t/s，);ylabelC x(t) ,);grid on;axis(0, 50*pi, -2, 2);subplot(3, 1, 2);plot(xl);grid on;title (输入信号xl
