实验三窗函数的特性分析报告.doc

  1. 1、本文档共12页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验三窗函数的特性分析报告

本科学生实验报告 学号 *************** 姓名 *************** 学院 物电学院 专业、班级 *************** 实验课程名称 数字信号分析与处理 教师及职称 *************** 开课学期 2015 至 2016学年 下 学期 填报时间 2016 年 3 月 25 日 云南师范大学教务处编印 一、验设计方案 实验序号 实验三 实验名称 窗函数的特性分析 实验时间 2016/3/25 实验室 同析楼三栋313实验室 1.实验目的 分析各种窗函数的时域和频域特性,灵活应用窗函数分析信号频谱和设计FIR数字滤波器。 2. 实验原理、实验流程或装置示意图 在确定信号谱分析、随机信号功率谱估计以及FIR数字滤波器设计中,窗函数的选择对频谱分析和滤波器设计都起着重要的作用。在确定信号谱分析和随机信号功率谱估计中,截短无穷长的序列会造成频率泄露,影响频率普分析的精确度和质量。合理选取窗函数的类型,可以改善泄露现象。在FIR数字滤波器设计中,截短无穷长的系统单位脉冲序列会造成FIR滤波器的幅度特性产生波动,且出出现过渡带。 【例1.3.1】 写出分析长度N=51点矩形窗的时域波行和频谱的MATLAB程序。 [解] N=51; w=boxcar(N); W=fft(w,256); subplot(2,1,1); stem([0:N-1],w); subplot(2,1,2); plot([-128:127],abs(fftshift(W))); 运算结果如图1.3.1所示 图1.3.1 矩形窗的时域波形和频谱 3.实验设备及材料 计算机,MATLAB软件 4.实验方法步骤及注意事项 注意事项: 在使用MATLAB时应注意中英输入法的切换,在中文输入法输入程序时得到的程序是错误的; MATLAB中两个信号相乘表示为x.*u,中间有个‘.’,同样两个信号相除也是如此; 使用MATLAB编写程序时,应新建一个m文件,而不是直接在Comandante窗口下编写程序; 在使用MATLAB编程时,应该养成良好的编写习惯。 5.实验数据处理方法 图像法 6.参考文献 《信号分析与处理》 《MATLAB数值计算与方法》 二、报告 1.实验现象与结果 实验内容 1.分析并绘出常用窗函数的时域波形特性。 N=51; w=boxcar(N); subplot(3,1,1); stem([0:N-1],w); title(矩形窗的时域波形) w=hanning(N); subplot(3,1,2); stem([0:N-1],w); title(汉宁窗的时域波形) w=hamming(N); subplot(3,1,3); stem([0:N-1],w); title(汉明窗的时域波形) 实验运行结果如图 3.1 图3.1 矩形窗、汉宁窗及汉明窗的时域特性波形 N=51; w=blackman(N); subplot(3,1,1); stem([0:N-1],w); title(布莱克曼窗的时域波形); w=bartlett(N); subplot(3,1,2); stem([0:N-1],w); title(Bartlett窗的时域波形); beta=2*N; w=Kaiser(N,beta); subplot(3,1,3); stem([0:N-1],w); title(凯泽窗的时域波形); 实验结果如图3.2所示 图3.2 布莱克曼窗、Bartlett窗、凯泽窗时域特性波形 3. 研究凯泽窗(Kaiser)的参数选择对其时域和频域的影响。 (1)N=20; beta=4; w=Kaiser(N,beta); subplot(3,2,1); stem([0:N-1],w); title(第3题凯泽窗N=20时域波形); W=fft(w,256); subplot(3,2,2); plot([-128:127],abs(fftshift(W))); title(第3题凯泽窗N=20频域波形 ) N=60; w=Kaiser(N,beta); subplot(3,2,3); stem([0:N-1],w); title(第三题凯泽窗N=60波形); W=fft(w,256); subplot(3,2,4); plo

文档评论(0)

ee88870 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档