[数字信号处理实验指导书.docVIP

数字信号处理实验指导书 前言 随着电子计算机、大规模集成电路（LSI）和超大规模集成电路（VLSI），以及微处理器技术的迅猛发展，数字信号处理无论在理论上还是在工程应用中，都是目前发展最快的学科之一，并且日趋完善和成熟。 《数字信号处理》作为高等院校电气信息学科学生的必修课，不仅是电类专业和信息类专业的一门基础课，而且也是工科学生受益较多的一门课程。该课程主要介绍离散信号与离散系统分析的基本理论,离散付立叶变换、快速付立叶变换和数字滤波器的设计、数字信号处理的实现等内容。通过该课程的学习,要求学生掌握数字信号处理的基本理论、基本方法和基本技能,为今后学习和应用数字信号处理技术奠定基础。 由于《数字信号处理》课程的许多理论基于较为繁琐的数学理论和数学推导，容易使学生感到乏味，且由于该课程自身具有很强的实践性，为了加深学生对所学知识的了解，培养其对实际问题的解决能力，本书开设了相关的实验教学内容。该实验内容主要利用目前流行的MATLAB软件来仿真信号处理的基本理论，并联系实际问题强化对信号理论的综合应用环节。 MATLAB语言在矩阵运算和数值计算等方面，具有强大的功能，非常直观、简洁，且效率高，交互性好，因而它在科学与工程计算方面的优点显而易见的。基于《数字信号处理》课程的特点，本书只针对该课程中常用到的MATLAB指令进行介绍，并给出相应示例。 本书一共由11个实验组成，其中部分为基础理论实验，部分为综合实验，教师可根据情况选做相关实验。 MATLAB信号处理工具 §1-1 离散时间信号与系统 在数字信号处理中，所有信号都是离散时间信号即序列，对应的处理系统也为离散时间系统。 1-1-1、基本信号的表示 MATLAB中，可采用向量表示序列，由于MATLAB中对序列下标默认为从1开始递增，因此要表示离散信号，一般应采用两个向量分别对信号的自变量和因变量预予描述。如, n= -3～4 在MATLAB中表示为 n = [-3, -2, -1, 0, 1, 2, 3, 4]; %自变量取值 x = [ 2, 1, -1, 0, 1, 4, 3, 7]; %因变量取值 注意： 向量可用方括号[ ]表示。 当向量取值连续变化时可用冒号运算符“：”简化赋值过程，如的n值，可简化为 n=[-3：4]或n= -3:4 。 分号“；”表示不回显表达式的值。 “%”表示其后内容为注释对象。 利用MATLAB，还可对信号的波形进行描述，常采用的绘图语句有stem，plot，subplot, axis，title，xlabel，ylabel，gtext, hold on, hold off, grid 等。其中stem 绘制离散图形；plot 绘制连续图形；subplot 用于绘制子图，应在stem 或plot 语句前调用；axis 指定x 和y轴的取值范围，用在stem或plot语句后；title 标注图形名称，xlabel, ylabel 分别标注x轴和y轴名称；gtext可将标注内容放置在鼠标点击处；hold on和 hold off 用于控制对象绘制方式，是在原图上还是在新图上绘制；grid用于绘制网格。 以上函数具体使用方法可通过在命令窗口中输入help加函数名进行查阅，如“help stem”语句。绘图语句的具体应用见下面示例，对应的波形图见图1-1。 % 画子图1 subplot( 1, 2, 1); stem (n, x1); axis([-5, 5, -2, 2]); %向量的前两个值为x轴的取值范围，后两个为y轴的取值范围 title (‘例1序列’); xlabel (‘n’); ylabel (‘x1(n)’); % 画子图2 subplot( 1, 2, 2); stem (n, x2); axis([0, 20, 0, 10]); title (‘例2序列’); xlabel (‘n’); ylabel (‘x2(n)’); 图1-1 序列波形图 下面列出信号处理中常用信号的表示方法： 单位采样序列 在区间的值，MATLAB中可用下列表达式描述 n = [ n1 : n2 ]; x = [ (n - n0) = =0 ]; 为使用方便，也可将上式描述成MATLAB函数，即 function [x , n ] = impseq (n0, n1, n2) n = [ n1 : n2 ]; x = [ (n - n0) = =0 ]; 2. 单位阶跃序列 在区间的值，MATLAB中可用下列表达式描述: n = [ n1 : n2 ]; x = [ (n - n0) =0 ]; 3. 实数指数序列 例如