数字信号处理实验4（副本）.docVIP

贵州大学实验报告 学院：计算机科学与信息 专业： 通信工程 班级：通信091 姓名 何川 学号 0908060235 实验组 实验时间 2011--12--21 指导教师 刘洪 成绩 实验项目名称 用凯塞窗设计线性相位FIR滤波器 实验目的 通过本实验掌握窗函数法设计线性相位FIR低通滤波器的方法。 实验要求 集中授课的教学形式 实验原理 1、窗函数设计法原理参见《数字信号处理》教材。 2、凯塞窗是对于给定的阻带衰减，是主瓣具有最大能量意义下的最佳窗函数，因此具有最陡的过渡带。所以通常选用凯塞窗来进行FIR滤波器的设计。 凯塞窗函数： 其中为形状参数，是第二类修正零阶贝塞尔函数，其幂函数展开式为： 给定滤波器的过渡带宽度和阻带衰减，则滤波器的长度和形状参数可由下列经验公式给出： 其中。 实验仪器 计算机一台； Matlab软件 实验步骤 1、编写程序。计算理想低通滤波器单位取样响应的M函数文件dlp.m。 其中， 2、编写计算凯塞窗函数的M函数文件kwind.m。 计算贝塞尔函数的参考程序besl.m如下： 3、编写程序，调用dlp.m和kwind.m文件，设计系统FIR滤波器。 实验内容 利用凯塞窗来设计线性相位的FIR滤波器。 实验数据 1、function h=dlp(p,s,N) ws=(p+s)/2; a=(N-1)/2; for n=0:N-1 na=n-a; if na==0 h(n)=ws/pi; else h(n+1)=sin(ws*na)/na/pi; end end end 2、function s=besl(x) eps=10^(-9); n=1; s=1; D=1; while D(eps*s) T=x/(2*n); n=n+1; D=D*T^2; s=s+D; end function w=kwind(B,N) I=besl(B); for n=0:N-1 x=B*(1-(1-(2*n/(N-1))^2))^(1/2); w(n+1)=besl(x)/I; end end 实验总结 通过本次实验，使我初步掌握了窗函数法设计线性相位FIR低通滤波器的方法，为以后更近一步的学习打下了坚实的基础！ 指导教师意见 签名： 年 月 日 注：各学院可根据教学需要对以上栏木进行增减。表格内容可根据内容扩充。