基于kaiser的fir数字滤波器设计报告书.docVIP

课 程 设 计 课程设计名称：基于Kaiser窗的FIR数字低通 滤波器设计 专 业 班 级 ： 电信1003班 学 生 姓 名 ： 董成成 学 号 ： 201046830509 指 导 教 师 ： 李相国 课程设计时间： 2013年6月 电子信息工程 专业课程设计任务书 学生姓名 董成成 专业班级 电信1003班 学号 201046830509 题 目 基于Kaiser窗的FIR数字低通滤波器设计 课题性质 工程技术研究 课题来源 自拟课题 指导教师 李相国 同组姓名 主要内容 设计一个具有如下性能指标的FIR数字低通滤波器：=0.25，=0.4，=0.01，=0.001；窗函数为Kaiser窗；分析最后设计结果性能。 任务要求 温习窗函数法设计数字FIR滤波器的原理和设计方法； 求出最后所设计出的滤波器冲激响应系数； 绘制所设计滤波器的增益响应曲线。 参考文献 [1] 胡广书. 数字信号处理—理论、算法与实现[M]. 北京: 清华大学出版社, 1997. [2] R. Lyons. Understanding Digital Signal Processing [M]. 2nd ed. Prentice Hall PTR., 2004. [3] A.V.奥本海姆, R.W.谢弗 and J.R.巴克. 离散时间信号处理[M]. 第二版. 西安交通大学出版社, 2001. [4] S. K. Mitra. Digital Signal Processing: A Computer-Based Approach[M]. 3rd ed. McGraw-Hill, 2005. 审查意见 指导教师签字： 教研室主任签字： 年 月 日 说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页 1 需求分析 基于Kaiser窗的FIR低通滤波器是利用Kaiser函数截取无限冲击响应h(n)，最后利用函数freqz得到加窗后的滤波器的幅频响应和相频响应。 2 概要设计 开始读入 开始 读入窗口长度M 计算hd(n) 用Kaiser窗函数求w(n) 计算H(k)=frez[h(n)] 绘图子（函数）绘制H(k)幅度相位曲线 结束 图6-1 主程序框图 计算h(n)= hd(n) w(n) （1）代入参数值，计算、、、、M的值： omegap=0.25*pi;omegas=0.4*pi;deltap=0.01;deltas=0.001; omegac=(omegap+omegas)/2; deltaomega=omegas-omegap; As=-20*log10(deltas); beta=0.1102*(As-8.7); M=ceil((As-8)/(2.285*deltaomega)); （2）定义定义域，Kaiser函数求得，截取冲击响应h(n)得出滤波器时域值h,并求其频域特性： n=0:M; p=M/2; wk=kaiser(M+1,beta); h=(sin(omegac*(n-p))./(pi*(n-p+eps))).*wk; [H,w]=freqz(h,1); （3）定义定义域，画滤波特性 Nc=fix(omegac/(pi/length(w))); Hd=[ones(1,Nc),zeros(1,length(w)-Nc)]; N=length(h);L=N/2; b=2*(h(L:-1:1));k=[1:1:L];k=k-0.5;w1=[0:1:511]*pi/512; A2=cos(w1*k)*b; E=Hd-A2; subplot(3,1,1) stem(n,h,fill,MarkerSize,2) axis([0,N,-0.1,0.5]) subplot(3,1,2) plot(w,20*log10(abs(H))) axis([0 pi -100 10]) subplot(3,1,3) plot(w,E) axis([0 pi -0.001 0.001]) 3 运行环境 Window2003系统和Win7系统 4 开发工具和编程语言 Matlab软件、Matlab语言 5 详细设计 clear; ----【清除变量】---------------------------------- omegap=0