数字信号处理上机实验指导书.docxVIP

数字信号处理实验指导书 《数字信号处理》上机实验指导书 陈纯锴 电子与信息工程学院 一、引言 “数字信号处理”是一门理论和实验密切结合的课程，为了深入地掌握课程内容，应当在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助学生深入地理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。所以，根据本课程的重点要求编写了四个实验。 第一章是全书的基础内容，抽样定理、时域离散系统的时域和频域分析以及系统对输入信号的响应是重要的基本内容。由于第一章大部分内容已经在前期《信号与系统》课程中学习完，所以可通过实验一帮助学生温习以上重要内容，加深学生对“数字信号处理是通过对输入信号的一种运算达到处理目的” 这一重要概念的理解。这样便可以使学生从《信号与系统》课程顺利的过渡到本课程的学习上来。 第二章、三章DFT、FFT是数字信号处理的重要数学工具，它有广泛的使用内容。采用实验二、实验三加深理解DFT的基本概念、基本性质。FFT是它的快速算法，必须学会使用。数字滤波器的基本理论和设计方法是数字信号处理技术的重要内容。学习这一部分时，应重点掌握IIR和FIR两种不同的数字滤波器的基本设计方法。IIR滤波器的单位冲激响应是无限长的，设计方法是先设计模拟滤波器，然后再通过S~Z平面转换，求出相应的数字滤波器的系统函数。这里的平面转换有两种方法，即冲激响应不变法和双线性变换法，后者没有频率混叠的缺点，且转换简单，是一种普遍应用的方法。FIR滤波器的单位冲激响应是有限长的，设计滤波器的目的即是求出符合要求的单位冲激响应。窗函数法是一种基本的，也是一种重要的设计方法。学习完第七章后可以进行实验四。 二、关于使用计算机语言 由于数字信号处理实验的主要目的是验证数字信号处理的有关理论，进一步理解巩固所学理论知识。所以，实现实验用的算法语言可以有许多种，但为了提高实验效率，要求学生用编程效率比C语言高好几倍的MATLAB语言。下面介绍MATLAB的主要特点。（有关MATLAB的启动、程序运行和有关信号处理工具箱函数等内容将放到最后附录中介绍。） MATLAB是一种交互式的以矩阵为基本数据结构的系统。在生成矩阵对象时，不要求明确的维数说明。所谓交互式，是指MATLAB的草稿纸编程环境。即用户每输入一条命令并按回车键，MATLAB系统便解释执行之，并显示执行结果。根据该结果，用户立即知道刚输入的命令的正确性，或利用中间结果进行其他处理等。 与C语言或FORTRON语言做科学数值计算的程序设计相比较，利用MATLAB可节省大量的编程时间。将其用于数字信号处理实验，则可大大提高实验效率，在有限的上机时间内，实验内容可增加几倍。例如，C语言FFT子程序有70多行，而用MATLAB只调用一个fft函数即可实现对序列进行FFT计算。 另外，MATLAB的工具箱及图形显示（打印）功能，可满足各层次人员直观、方便的进行分析、计算和设计工作，从而可大大节省时间。例如，序列的卷积、滤波，系统函数H(z)的幅频特性和相频特性等计算，均有现成的工具箱函数。而用其它算法语言完成这些计算的编程比较麻烦，且程序较长。 由于上述特点，在美国一些大学里，MATLAB已成为辅助教学的有益工具。MATLAB已成功地用于数字信号处理课程中的问题分析、实验、滤波器设计及计算机模拟。附录中所介绍的信号处理工具箱函数及绘图函数基本可满足本教材所要求的上机实验需要。对序列进行谱分析的MATLAB程序及运行结果见附录。 实验一 离散信号产生及频谱的绘制 一、实验目的: （1）熟悉Matlab环境。 （2）掌握 Matlab 中一些基本函数的建立方法 （2）通过编程绘制的幅度相位谱加深理解系统的特性 二、实验内容 1、编写程序，产生以下离散序列： （1）f(n)=δ(n) (-3n4) n1=-3;n2=4;n0=0; n=n1:n2; x=[n==n0]; stem(n,x,filled); axis([n1,n2,0,1.1*max(x)]); xlabel(时间(n));ylabel(幅度x(n)); title(单位脉冲序列); （2）f(n)=u(n) (-5n5) n1=-5;n2=5;n0=0; n=n1:n2; x=[n=n0]; stem(n,x,filled); axis([n1,n2,0,1.1*max(x)]); xlabel(时间(n));ylabel(幅度x(n)); title(单位阶跃序列); box （3） (0n16) n1=16;a=-0.1;w=1.6*pi; n=0:n1; x=exp((a+j*w)*n); subplot(2,2,1);plot(n,real(x)); title(复指数信号的实部); subplot(2,2,3);stem