[学习总结]中南大学《数字信号处理》课程设计报告.doc
《[学习总结]中南大学《数字信号处理》课程设计报告.doc》由会员分享,可在线阅读,更多相关《[学习总结]中南大学《数字信号处理》课程设计报告.doc(36页珍藏版)》请在三一文库上搜索。
1、中南大学现代信号处理课程设计报告 学 院: 专业班级: 姓 名: 学 号: 指导老师: 设计时间: 前 言 现代信号处理课程设计与现代信号处理课程配套,是通信工程类专业的重要实践环节。数字信号处理是每一个电子信息科学工作者必须掌握的重要知识。它采用计算机仿真软件,以数值计算的方法对信号进行分析、变换、滤波、检测、估计与识别等加工处理,以达到提取信息便于使用的目的。通信工程专业的培养目标是具备通信技术的基本理论和应用技术,能从事电子、信息、通信等领域的工作。鉴于我校充分培养学生实践能力的办学宗旨,对本专业学生的培养要进行工程素质培养、拓宽专业口径、注重基础和发展潜力。特别是培养学生的创新能力,以
2、实现技术为主线多进行实验技能的培养。学生通过现代信号处理课程设计这一重要环节,可以将本专业的主干课程现代信号处理从理论学习到实践应用,对数字信号处理技术有较深的了解,进一步增强学生动手能力和适应实际工作的能力。 目 录前 言2第一章 设计目的及要求41.1 课程设计目的41.2 课程设计要求41.3 课程设计选做题目4第二章 设计题目及思想52.1 各题目的设计思想52.1.1 第一题52.1.2 第二题62.1.3 第三题62.1.4 第四题72.1.5 第五题72.1.6 第六题82.2 GUI界面设计思想9第三章 调试及结果分析123.1 GUI界面的结果与调试123.2 各题目的结果与
3、调试123.2.1 第一题123.2.2 第二题133.2.3 第三题143.2.4 第四题153.2.5 第五题163.2.6 第六题21第四章 课程设计总结244.1 课程设计所遇问题及解决方案244.2 课程设计心得体会25附录:26参考文献35第一章 设计目的及要求1.1 课程设计目的1全面复习课程所学理论知识,巩固所学知识重点和难点,将理论与实践很好地结合起来。2. 掌握信号分析与处理的基本方法与实现3提高综合运用所学知识独立分析和解决问题的能力;4熟练使用一种高级语言进行编程实现。1.2 课程设计要求一、使用MATLAB(或其它开发工具)编程实现上述内容,写出课程设计报告。滤波器设
4、计题目应尽量避免使用现成的工具箱函数。为便于分析与观察,设计中所有频谱显示中的频率参数均应对折叠频率归一化。二、课程设计报告的内容包括:(1)课程设计题目和题目设计要求;(2)设计思想和系统功能分析;(3)设计中关键部分的理论分析与计算,关键模块的设计思路;(4)测试数据、测试输出结果,及必要的理论分析和比较;(5)总结,包括设计过程中遇到的问题和解决方法,设计心得与体会等;(6)参考文献;(7)程序源代码清单。三、课程设计考核要求: 课程考核分三部分,一部分是上机率,占20%;第二部分是检查成绩,最后两次上机为检查时间,占50%;第三部分为课程设计报告,占30%。注意:(1)演示系统使用GU
5、I界面或混合编程实现集成打包发布,酌情加分。(2)若发现程序或课程设计报告雷同,一律不及格。1.3 课程设计选做题目选题1.音乐信号处理: 1)获取一段音乐或语音信号,设计单回声滤波器,实现信号的单回声产生。给出加入单回声前后的信号频谱。 2)设计多重回声滤波器,实现多重回声效果。给出加入多重回声后的信号频谱。 3)设计均衡器,使得得不同频率的混合音频信号,通过一个均衡器后,增强或削减某些频率区域。选题二. 对混有噪声的音乐信号进行滤波: 1)在一段音乐信号中混入两个频率的正弦型干扰信号,利用FFT计算该混合信号的频谱并确定干扰信号的频谱; 2)利用二阶带阻滤波器设计一个能滤出干扰信号的梳状滤
6、波器; 3)利用梳状滤波器滤除信号中的噪声,播放处理前后的信号,并比较处理前后的结果本次课程设计共有6道必做题目,并且我们并没有系统的学过MATLAB,所以我并没有做选做题。第二章 设计题目及思想2.1 各题目的设计思想2.1.1 第一题一、题目:1. 给定模拟信号:1)选择采样频率Fs = 5000Hz和合适的信号长度,采样得到序列 x1(n)。求并画出x1(n)及其序列傅里叶变换 |X1(ejw)|。2)选择采样频率Fs = 1000Hz和合适的信号长度,采样得到序列 x2(n)。求并画出x2(n)及其序列傅里叶变换 |X2(ejw)|。3)说明|X1(ejw)|与|X2(ejw)|间的区
7、别,为什么?二、思想:第一题主要的考查点是对连续时间信号的抽样,对抽样信号进行傅里叶变换以及不同采样频率对幅频特性曲线的影响等。主要的理论依据如下:(1) 对连续时间信号的抽样:设连续时间信号为x(t),则对其采样得到的离散时间信号为x(n)=x(nT),T为采样周期;(2) 离散时间信号的傅里叶变换:。 由于抽样后的离散时间信号题目没有给出,故应先对连续信号做抽样,后再对离散时间信号进行傅里叶变换。本次我直接采用输入公式的方式进行傅里叶变换。2.1.2 第二题一、题目:2. 已知两系统分别用下面差分方程描述:试分别写出它们的传输函数,并分别打印曲线。说明这两个系统的区别。二、思想:第二题的主
8、要考查点是由系统的方程得出系统的传输函数。而在MATLAB中有freqz函数可直接求出系统的传输函数,然后由abs函数以及angle函数可得出其幅频相频特性曲线。2.1.3 第三题一、题目:3. 已知已调信号,其中调制信号,载波为。1)选择合适的采样频率及信号长度,使用FFT分析该已调信号的频谱并画出其幅频特性和相频特性曲线图。2)对该已调信号进行解调,恢复原调制信号。二、思想:第三题的主要考查点是信号的调制与解调。这属于通信原理的主要原理知识。本题第一问需要使用FFT分析已调信号的频谱特性,在MATLAB中FFT可以直接作为函数调用,所以我选择直接调用FFT对其进行分析。第二问中要对已调信号
9、进行解调。其主要思想是解调已调信号需要先将已调信号乘以载波后选择适当的滤波器进行滤波,将其多余信号滤除。2.1.4 第四题一、题目:4. 已知三角波序列和反三角波序列: 用N=8点FFT分析序列x1(n)和x2(n)的幅频特性,观察两者的序列形状和频谱特性曲线有什么异同?绘出两序列及其幅频特性曲线。在x1(n)和x2(n)的末尾补零,用N=16点FFT分析这两个信号的幅频特性,观察幅频特性发生了什么变化?两情况的FFT频谱还有相同之处吗?这些变化说明了什么?二、思想:第四题的主要考查点是FFT快速傅里叶变换的用法。而这在MATLAB中,FFT可以直接作为函数调用,所以我选择用FFT快速傅里叶变
10、换进行分析三角波序列。不过首先要将三角波以及反三角波的取值范围和函数值计算出来,后才可直接对序列进行FFT变换并分析其频谱特性。2.1.5 第五题一、题目:5. 设有一信号,设计各种IIR数字滤波器以实现:1)低通滤波器,滤除的成分,保留成分2)高通滤波器,滤除的成分,保留成分3)带通滤波器,滤除的成分,保留成分4)带阻滤波器,滤除的成分,保留成分要求:1)求出各个滤波器的阶数,设计各滤波器。画出各滤波器的幅频和相频特性,计算滤波器的系统函数H(z)。 2) 画出滤波前后信号的时域、频域波形。二、思想:第五题的主要考查点是用IIR双线性变换法设计BW类型(低通、高通、带通以及帯阻)数字滤波器。
11、其主要理论依据如下:(1)正切变换实现频率压缩: (2)双线性变换: 或 (3)将所需数字滤波器的边界频率转换成相应类型的模拟滤波器的边界频率公式为:其设计主要按照课中所学滤波器设计的步骤进行分步设计。但由于其系统函数的公式无法直接输入,故采用buttord函数进行巴特沃斯滤波器的方法分析所给参数计算出其系统函数的系数矩阵。再借由系统函数的系数矩阵得出具体的系统函数。得出系统函数后便可得到相应的阶数以及幅频、相频特性。滤波前的信号频域波形由先对时域波形进行FFT变换后得出。对原信号的滤波则采用直接调用filter函数进行对原函数的滤波,而频域波形同样经由先FFT变换分析后得出。2.1.6 第六
12、题一、题目:6.1)用Hanning窗设计一线性相位带通数字滤波器,要求:N=15, 。观察它的实际3dB和20dB带宽。N=45,重复这一设计,观察幅频和相位特性的变化,注意长度N变化的影响;2)分别改用矩形窗和Blackman窗,设计(1)中的带通滤波器,观察并记录窗函数对滤波器幅频特性的影响,比较三种窗的特点;总结窗的不同长度和不同窗对滤波器的影响二、思想:第六题的主要考查点是几种主要的窗函数(矩形窗、汉宁窗、哈明窗等)对离散时间信号的截取,并用FFT对其频谱进行分析。其主要理论依据如下:(1)矩形窗(Rectangle Window) 其幅度函数为:(2)汉宁(Hanning)窗升余弦
13、窗 其频谱函数为: 其幅度函数为:(3)哈明(Hamming)窗改进的升余弦窗 其频谱函数为: 当时,其幅度函数可近似表示为:2.2 GUI界面设计思想本次课程设计由于结果中图形太多,若一次性完成所以题目的演示则需要开创许多个新的窗口进行图形建立,这样不方便验收及查看。所以最后决定采用图形用户(Graphical User Interface,简称 GUI)界面,下面对GUI界面进行简单的一些介绍。GUI 是 Graphical User Interface 的简称,即图形用户界面,通常人机交互图形化用户界面设计经常读做“goo-ee”,准确来说 GUI 就是屏幕产品的视觉体验和互动操作部分。
14、 GUI 是一种结合计算机科学、美学、心理学、行为学,及各商业领域需求分析的人机系统工程,强调人机环境三者作为一个系统进行总体设计。 在MATLAB主界面中输入“guide”,再敲击回车键,弹出GUI界面窗口,即可新建GUI界面,如下图所示:图2.5 GUI界面的新建窗口图在图2.5的窗口中点击“OK”,则可新建一个GUI界面,并进入其编辑窗口,如图2.6所示:图2.6 GUI界面的编辑窗口图 随后便可在此窗口内编辑演示界面,此次设计主要利用的是Static Text(静态文本)和Push Button(按钮)。在界面中输入相关的静态文本,并新建若干个按钮,再以此在每个按钮上点击右键,选择Vi
15、ew CallBacks调出Call Back函数,将相应设计题的源代码输入Call Back函数界面中对应按钮的段落后。如图2.7和2.8所示:图2.7 调用CallBack的操作流程图图2.8 CallBack代码编辑窗口图 如此便可以建立一个方便、美观又直观的演示界面,以便老师验收时操作简单、条理清晰、并且用时不长。第三章 调试及结果分析3.1 GUI界面的结果与调试 这次设计的演示界面采用的是GUI界面,这一点在前面提到过,在此不赘述了,只将界面的调试结果说明一下。演示界面如图3.1所示:图3.1 演示界面图 如图3.1所示,用户克根据需要选择题目,每个题中又有不同的选项,只要点击一下
16、对应的按钮,相关信号的处理结果就会出现在一个新的figure窗口中。3.2 各题目的结果与调试3.2.1 第一题 第一题的调试结果如图3.2所示,各图形的意义均有明确表明。图3.2 第一题的调试结果图 分析: 分别以所要求采样频率fs=5000Hz、1000Hz对进行等间隔采样,得到,(为采样周期),可知,采样信号的频谱函数是在原模拟信号频谱函数的周期延拓;若以频率为自变量,则以采样频率为延拓周期;对频带限于的模拟信号,根据采样定理,只有当时,采样后才不会发生频谱混叠失真。但是在Matlab上是无法计算连续函数,只有在当足够大的时候,我们才将频谱混叠忽略不计,从而可对采样序列进行傅里叶变换;最
17、后应用subplot()命令实现画图,并注意到要归一化。由以上实验得到的图形结果可以看到:,当采样频率越大的时候,采样信号的频谱越陡峭,而其失真情况也越来越小,根据采样定理“采样频率越接近信号频率,其失真情况就越小”可以得到上述结论。3.2.2 第二题 第二题的调试结果如图3.3所示,各图形的意义均有明确表明。图3.3 第二题的调试结果图分析:由于所求系统的系统方程已经由题目可知,则其零极点分布图也可由上图所知。故根据计算所得出其系统函数后做FFT变换所得幅频特性和相频特性均由上图可知。3.2.3 第三题 第三题的调试结果如图3.4所示,上图为第一问结果,下图为第二问结果。图3.4 第三题的调
18、试结果图分析:当已调信号进行不同点数的采样后再经不同间隔的FFT变换后,其展现的幅频特性以及相频特性类似于上图。因选取的采样点不同则结果不同。第二问中的滤波器的设计参数由题目中的调制信号和载波决定,由此设计出的低通滤波器的系统函数图如上。最后因滤波器的不严密而导致解调后的原信号有所失真,并不能完完全全的解调出原信号。3.2.4 第四题 第四题的调试结果如图3.5所示,图中所代表意义均有标明 图3.5 第四题的调试结果图分析:首先先根据题目中所给的三角波函数以及反三角波函数的函数式得出其波形。经由不同采样的FFT变换后有stem画出其图形如上。由图可看出三角波序列和反三角波序列在8点的FFT变换
19、后的图形基本一样,没有太大差别。而16点FFT后的图形则有不小差异,此说明对其FFT变换时,取的点数越多,图形差异越大。3.2.5 第五题 第五题的调试结果第一问中低通滤波器的图形如图3.6所示。图3.6 第五题的低通滤波器调试结果图高通滤波器的图形如图3.7所示。图3.7 第五题的高通滤波器调试结果图带通滤波器的图形如图3.8所示。图3.8 第五题的带通滤波器调试结果图带阻滤波器的图形如图3.9所示。图3.9 第五题的带阻滤波器调试结果图 分析:经由一下理论知识可依次求出各滤波器的参数及系统函数。(1)正切变换实现频率压缩: (2)双线性变换: 或 (3)将所需数字滤波器的边界频率转换成相应
20、类型的模拟滤波器的边界频率公式为:频域图形均由FFT变换后可得出。此次设计IIR数字滤波器是其参数选择并不精准,均由自己在可行范围中尝试得出最像的结果。但其滤波器的幅频特性和相频特性仍然具有失真现象。3.2.6 第六题 第六题的调试结果依次的图形如图4.0所示。图4.0 第六题的调试结果图 分析:此图由上往下依次为N=15及N=45的汉明窗、布莱克曼窗和矩形窗。由图形可以看出N越小时窗函数的主瓣范围越大,副瓣坡度越缓。滤波器滤波功能越好。第四章 课程设计总结4.1 课程设计所遇问题及解决方案 拿到设计题时,完全不知道从何处入手,一开始只会按照书上的解题思路一步步来算,一步步来解这些题目,解出来
21、后再设法将其输入编程中,从而的到想要的答案。可是每次如此,系统输出的结果总是和预想中的不一样,有些甚至不知道哪里出了错。不过后来,在翻阅了有关参考文献后,知道了题目中大部分内容都可以通过MATLAB强大的函数工具箱里的函数直接实现。开始设计后,总是用错函数,输错参数等,这类问题层出不穷,毕竟是第一次用MATLAB进行设计实践,很多问题还是要虚心向老师和同学请教的。在很多次的请教,修改后,程序基本成功地得到了运行,调试过程中又遇到了一些小问题,大多是由于粗心导致的,这些问题在调试过程中都得到了修正。这次课程设计过程中遇到的主要问题是频率归一化的问题,开始对归一化的概念不是很清楚。在问了班上能力较
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 学习总结 数字信号处理 学习 总结 中南 大学 数字信号 处理 课程设计 报告
链接地址:https://www.31doc.com/p-1973385.html