中南大学《数字信号处理》课程设计报告

1中南大学《现代信号处理》课程设计报告学院:专业班级：姓名：学号：指导老师：设计时间：2前言《现代信号处理》课程设计与《现代信号处理》课程配套，是通信工程类专业的重要实践环节。数字信号处理是每一个电子信息科学工作者必须掌握的重要知识。它采用计算机仿真软件，以数值计算的方法对信号进行分析、变换、滤波、检测、估计与识别等加工处理，以达到提取信息便于使用的目的。通信工程专业的培养目标是具备通信技术的基本理论和应用技术，能从事电子、信息、通信等领域的工作。鉴于我校充分培养学生实践能力的办学宗旨，对本专业学生的培养要进行工程素质培养、拓宽专业口径、注重基础和发展潜力。特别是培养学生的创新能力，以实现技术为主线多进行实验技能的培养。学生通过《现代信号处理》课程设计这一重要环节，可以将本专业的主干课程《现代信号处理》从理论学习到实践应用，对数字信号处理技术有较深的了解，进一步增强学生动手能力和适应实际工作的能力。3目录前言............................................................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.2GUI界面设计思想............................................9第三章调试及结果分析............................................123.1GUI界面的结果与调试.......................................123.2各题目的结果与调试.........................................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参考文献..........................................................354第一章设计目的及要求1.1课程设计目的1．全面复习课程所学理论知识，巩固所学知识重点和难点，将理论与实践很好地结合起来。2.掌握信号分析与处理的基本方法与实现3．提高综合运用所学知识独立分析和解决问题的能力；4．熟练使用一种高级语言进行编程实现。1.2课程设计要求一、使用MATLAB（或其它开发工具）编程实现上述内容，写出课程设计报告。滤波器设计题目应尽量避免使用现成的工具箱函数。为便于分析与观察，设计中所有频谱显示中的频率参数均应对折叠频率归一化。二、课程设计报告的内容包括：(1)课程设计题目和题目设计要求；(2)设计思想和系统功能分析；(3)设计中关键部分的理论分析与计算，关键模块的设计思路；(4)测试数据、测试输出结果，及必要的理论分析和比较；(5)总结，包括设计过程中遇到的问题和解决方法，设计心得与体会等；(6)参考文献；(7)程序源代码清单。三、课程设计考核要求：课程考核分三部分，一部分是上机率，占20%；第二部分是检查成绩，最后两次上机为检查时间，占50%；第三部分为课程设计报告，占30%。注意：（1）演示系统使用GUI界面或混合编程实现集成打包发布，酌情加分。（2）若发现程序或课程设计报告雷同，一律不及格。1.3课程设计选做题目选题1.音乐信号处理：1）获取一段音乐或语音信号，设计单回声滤波器，实现信号的单回声产5生。给出加入单回声前后的信号频谱。2）设计多重回声滤波器，实现多重回声效果。给出加入多重回声后的信号频谱。3）设计均衡器，使得得不同频率的混合音频信号，通过一个均衡器后，增强或削减某些频率区域。选题二.对混有噪声的音乐信号进行滤波：1）在一段音乐信号中混入两个频率的正弦型干扰信号，利用FFT计算该混合信号的频谱并确定干扰信号的频谱；2）利用二阶带阻滤波器设计一个能滤出干扰信号的梳状滤波器；3）利用梳状滤波器滤除信号中的噪声，播放处理前后的信号，并比较处理前后的结果本次课程设计共有6道必做题目，并且我们并没有系统的学过MATLAB，所以我并没有做选做题。第二章设计题目及思想2.1各题目的设计思想2.1.1第一题一、题目：1.给定模拟信号：ettxa1000)(1）选择采样频率Fs=5000Hz和合适的信号长度，采样得到序列x1(n)。求并画出x1(n)及其序列傅里叶变换|X1(ejw)|。2）选择采样频率Fs=1000Hz和合适的信号长度，采样得到序列x2(n)。求并画出x2(n)及其序列傅里叶变换|X2(ejw)|。3）说明|X1(ejw)|与|X2(ejw)|间的区别，为什么？二、思想：第一题主要的考查点是对连续时间信号的抽样，对抽样信号进行傅里叶变换以及不同采样频率对幅频特性曲线的影响等。6主要的理论依据如下：（1）对连续时间信号的抽样：设连续时间信号为x(t)，则对其采样得到的离散时间信号为x(n)=x(nT)，T为采样周期；（2）离散时间信号的傅里叶变换：()[()]()jjnnXeFTxnxne。由于抽样后的离散时间信号题目没有给出，故应先对连续信号做抽样，后再对离散时间信号进行傅里叶变换。本次我直接采用输入公式的方式进行傅里叶变换。2.1.2第二题一、题目：2.已知两系统分别用下面差分方程描述：)1()()(1nxnxny)1()()(2nxnxny试分别写出它们的传输函数，并分别打印weHjw~)(曲线。说明这两个系统的区别。二、思想：第二题的主要考查点是由系统的方程得出系统的传输函数。而在MATLAB中有freqz函数可直接求出系统的传输函数，然后由abs函数以及angle函数可得出其幅频相频特性曲线。2.1.3第三题一、题目：3.已知已调信号)9cos()cos()(tttx，其中调制信号)cos(t，载波为)9cos(t。1）选择合适的采样频率及信号长度，使用FFT分析该已调信号的频谱并画出其幅频特性和相频特性曲线图。2）对该已调信号进行解调，恢复原调制信号。二、思想：第三题的主要考查点是信号的调制与解调。这属于通信原理的主要原理知识。本题第一问需要使用FFT分析已调信号的频谱特性，在MATLAB中FFT可以直接作为函数调用，所以我选择直接调用FFT对其进行分析。第二问中要7对已调信号进行解调。其主要思想是解调已调信号需要先将已调信号乘以载波后选择适当的滤波器进行滤波，将其多余信号滤除。2.1.4第四题一、题目：4.已知三角波序列1()xn和反三角波序列2()xn：11,03()8,470,nnxnnn其它24,03()3,470,nnxnnn其它用N=8点FFT分析序列x1(n)和x2(n)的幅频特性，观察两者的序列形状和频谱特性曲线有什么异同？绘出两序列及其幅频特性曲线。在x1(n)和x2(n)的末尾补零，用N=16点FFT分析这两个信号的幅频特性，观察幅频特性发生了什么变化？两情况的FFT频谱还有相同之处吗？这些变化说明了什么？二、思想：第四题的主要考查点是FFT快速傅里叶变换的用法。而这在MATLAB中，FFT可以直接作为函数调用，所以我选择用FFT快速傅里叶变换进行分析三角波序列。不过首先要将三角波以及反三角波的取值范围和函数值计算出来，后才可直接对序列进行FFT变换并分析其频谱特性。2.1.5第五题一、题目：5.设有一信号2()=1+cos(n)+cos(n)43xn，设计各种IIR数字滤波器以实现：1）低通滤波器，滤除2cos(n)3的成分，保留成分1+cos(n)42）高通滤波器，滤除1+cos(n)4的成分，保留成分2cos(n)33）带通滤波器，滤除21+cos(n)3的成分，保留成分cos(n)44）带阻滤波器，滤除cos(n)4的成分，保留成分21+cos(n)3要求：1)求出各个滤波器的阶数，设计各滤波器。画出各滤波器的幅频和相8频特性，计算滤波器的系统函数H（z）。2)画出滤波前后信号的时域、频域波形。二、思想：第五题的主要考查点是用IIR双线性变换法设计BW类型（低通、高通、带通以及帯阻）数字滤波器。其主要理论依据如下：（1）正切变换实现频率压缩：12tan(/2)/TT（2）双线性变换：11(2/)[(1)(1)]sTzz或(2/)/(2/)zTsTs（3）将所需数字滤波器的边界频率转换成相应类型的模拟滤波器的边界频率公式为：(2/)tan(/2)T其设计主要按照课中所学滤波器设计的步骤进行分步设计。但由于其系统函数的公式无法直接输入，故采用buttord函数进行巴特沃斯滤波器的方法分析所给参数计算出其系统函数的系数矩阵。再借由系统函数的系数矩阵得出具体的系统函数。得出系统函数后便可得到相应的阶数以及幅频、相频特性。滤波前的信号频域波形由先对时域波形进行FFT变换后得出。对原信号的滤波则采用直接调用filter函数进行对原函数的滤波，而频域波形同样经由先FFT变换分析后得出。2.1.6第六题一、题目：6.1)用Hanning窗设计一线性相位带通数字滤波器，要求：N=15,。观察它的实际3dB和20dB带宽。N=45，重复这一设计，观察幅频和相位特性的变化，注意