《数字信号处理课程设计.docVIP

《数字信号处理》 课程设计报告 任课教师：王晓君 指导教师：王晓君 学生学号：130701213 学生姓名：田继辉 所学专业：电信132 2016年1月8日 课程设计成绩评定表 学生姓名 田继辉 学号 130701213 成绩 专业班级 电信132 起止时间 2016.1.6至2016.1.9 设计题目 设计一 DFT在信号频谱分析中的应用 设计四 离散时间系统频域分析 指 导 教 师 评 语 指导教师： 年 月 日 目录 一、 设计题目 3 二、 设计目的 3 （一）DFT在信号频谱分析中的应用 3 （二）离散时间系统频域分析 3 三、 设计原理 3 （一）DFT在信号频谱分析中的应用 3 （二）离散时间系统频域分析 3 四、 实现方法 3 （一）DFT在信号频谱分析中的应用 3 （二）离散时间系统频域分析 4 五、 设计内容及要求 4 （一）DFT在信号频谱分析中的应用 4 1、编写计算序列x(n)的N点DFT的m函数文件 4 2、分析离散确定信号的幅频图 5 3、研究高密度谱与高分辨率频谱 8 （二）离散时间系统频域分析 11 1、编写计算N阶差分方程所描述系统频响函数的m函数文件 11 2、构造N阶数字低通滤波器的差分方程 11 3、改变系数转为数字高通、全通滤波器 12 六、 设计结果及改进建议 14 （一）DFT在信号频谱分析中的应用 14 （二）离散时间系统频域分析 14 七、 回答思考题 14 （一）DFT在信号频谱分析中的应用 14 （二）离散时间系统频域分析 15 八、 设计体会 15 九、 参考文献 16 设计题目 设计一 DFT在信号频谱分析中的应用 设计四 离散时间系统频域分析 设计目的 （一）DFT在信号频谱分析中的应用 1、熟悉DFT的性质。 2、加深理解信号频谱的概念及性质。 3、了解高密度谱与高分辨率频谱的区别。 （二）离散时间系统频域分析 1、学习离散时间系统频率特性的计算方法。 2、深刻理解离散时间系统频率特性与滤波特性的关系。 3、掌握离散时间系统的系统参数、系统零极点及系统频率特性间的关系。 设计原理 （一）DFT在信号频谱分析中的应用 1、学习用DFT和补零DFT的方法来计算信号的频谱。 2、用MATLAB语言编程来实现，在做课程设计前，充分预习课本DTFT、DFT及补零DFT的有关概念，熟悉MATLAB语言，独立编写程序。 （二）离散时间系统频域分析 1.用MATLAB语言编程分析数字滤波器的各种滤波特性及其相关的参数。 2.充分预习相关理论知识，理解低通、高通、全通滤波器的工作原理。 实现方法 MATLAB语言由美国的Clever Moler博士于1980年开发；集科学计算、图象处理；声音处理于一身，并提供丰富的Windows图形界面设计方法；被称为第四代编程语言，程序简洁、可读性很强而且调试十分容易；可以指令驱动模式和程序文件模式工作。 （一）DFT在信号频谱分析中的应用 因直接计算DFT变换的运算量很大，因此常用快速傅里叶算法FFT代替DFT变换。此题目调用到MATLAB内部的FFT函数。 （二）离散时间系统频域分析 首先构造N阶差分方程，使该差分方程为数字低通滤波器；进而改变差分方程的系数，使之为数字高通滤波器；重新改变差分方程的系数，使之为数字全通滤波器，从而计算N阶差分方程所描述系统频响函数H。 设计内容及要求 （一）DFT在信号频谱分析中的应用 1、编写计算序列x(n)的N点DFT的m函数文件 用MATLAB语言编写计算序列x(n)的N点DFT的m函数文件DFTmatlab.m，并与MATLAB中的内部函数文件fft.m作比较。 %DFTmatlab.m function[q]=DFTmatlab(xn,N) n=[0:1:N-1]; %n取0到N-1 k=[0:1:N-1]; WN=exp(-j*2*pi/N); %求WN nk=n*k; WNnk=WN.^nk; q=xn*WNnk %求出xn的DFT表达式 DFT与FFT计算时间比较的源程序： %TimeCompare.m %DFT与FFT计算时间比较 dft_time=zeros(1,100); %零矩阵：zeros(m,n),zeros(m) for n=1:100; t=clock; %clock:将当前时间作为日期向量返回 a=sin(n)+cos(n); DFTmatlab(a,n); dft_time(n