信息处理方向课程设计介绍.docxVIP

江南大学 物联网工程学院 PAGE  PAGE 1 信息处理方向课程设计 实验报告 班级: 姓名: 学号: 一、课程设计目的及要求 数字信号处理是一门理论性和实践性都很强的学科，通过课程设计可以加深理解掌握基本理论，培养学生分析问题和解决问题的综合能力，为将来走向工作岗位奠定坚实的基础，因此做好课程设计是学好本课程的重要教学辅助环节。 本指导书结合教材《数字信号处理教程》的内容，基于MATLAB程序语言提出课程设计的题目及要求，在做课程设计之前要求学生要尽快熟悉MATLAB语言，充分预习相关理论知识，独立编写程序，以便顺利完成课程设计。 二、课程设计任务 课程设计的过程是综合运用所学知识的过程。课程设计主要任务是围绕数字信号的频谱分析、特征提取和数字滤波器的设计来安排的。根据设计题目的具体要求，??用MATLAB语言完成题目所规定的任务及功能。设计任务包括：查阅专业资料、工具书或参考文献，了解设计课题的原理及算法、编写程序并在计算机上调试，最后写出完整、规范的课程设计报告书。课程设计地点在信息学院机房，一人一机，在教师统一安排下独立完成规定的设计任务。 三、课程设计题目 根据大纲要求提供以下四个课程设计题目供学生选择，根据实际情况也可做其它相关课题。 1. DFT在信号频谱分析中的应用 2. 数字滤波器的单位脉冲响应输出及其稳定性 3. 离散时间系统频域分析 4.数字滤波器的设计实现 四、实验内容及步骤 . 设计一 DFT在信号频谱分析中的应用 用MATLAB语言编写计算序列x(n)的N点DFT的m函数文件dft.m。并与MATLAB中的内部函数文件fft.m作比较。 dft.m程序： function Xk=dft(xn,N) if length(xn)N xn=[xn,zeros(1,N-length(xn))]; end n=0:N-1; for k=0:N-1 Xk(1,k+1)=sum(xn.*exp((-1)*j*n*k*(2*pi/N))); End 比较后发现与fft.m结果一样 2. 对离散确定信号 作如下谱分析： 截取使成为有限长序列N()，(长度N自己选)写程序计算出的N点DFT ,并画出相应的幅频图。 程序： n=0:11; xn=cos(0.48*pi*n)+cos(0.52*pi*n); Xk=fft(xn,12); subplot(2,1,1);stem(n,xn);grid;title(xn波形图); subplot(2,1,2);stem(n,abs(Xk));grid;title(幅频图); 将 (1)中补零加长至M点（长度M自己选），编写程序计算的M点DFT ,并画出相应的图。 程序： n=0:11; xn=cos(0.48*pi*n)+cos(0.52*pi*n); n1=0:29; xn1=[xn,zeros(1,18)]; Xk1=fft(xn1,30); subplot(2,1,1);stem(n,xn);grid; title(xn); subplot(2,1,2);stem(n1,abs(Xk1));grid;title(|Xk1|~k); 利用补零DFT计算 (1)中N点有限长序列频谱并画出相应的幅频图。 程序： n=0:11; xn=cos(0.48*pi*n)+cos(0.52*pi*n); n1=0:999; xn1=[xn,zeros(1,988)]; Xk1=fft(xn1,1000); plot(n1,abs(Xk1)); 3. 在计算机上用DFT对模拟信号进行谱分析时, 只能以有限大的采样频率对模拟信号的近似频谱。其误差主要来源于截断效应(频谱泄漏和谱间干扰)和频谱混叠失真。前者使谱分辨率降低，产生谱间干扰；后者使折叠频率附近的频谱产生较大失真。 实践证明，加大截取长度可提高频率分辨率；选择合适的窗函数可降低谱间干扰；而频谱混叠失真要通过提高采样频率和预滤波来改善。 按题目要求编写程序，验证截断效应及加窗的改善作用，参数选取如下： 采样频率fs=400Hz,T=1/fs； (2) 对作4096点DFT作为的近似频谱； (3) 取三种截取时间长度分别： Tp=0.04s, Tp=4*0.04s,Tp=8*0.04s;两种窗函数分别为矩形窗和Hamming窗。 实验结果记录： 程序： fs=800; T=1/fs; Tp=0.04; N=Tp*fs; N1=[N,4*N,8*N]; for m=1:3; n=1:N1(m); xn=cos(200*pi*n*