数字信号处理实验(1-4)..docVIP

《数字信号处理》实验报告 学院：___________________ 专业：__________________ 班级：___________________ 学号：___________________ 姓名：___________________ 湖南工业大学理学院 2014年12月7日 实验一 离散时间信号的表示及运算 1.1 实验目的 学会运用MATLAB表示的常用离散时间信号； 学会运用MATLAB实现离散时间信号的基本运算。 1.2 实验原理及实例分析 1.2.1 离散时间信号在MATLAB中的表示 离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用来表示，自变量必须是整数。 离散时间信号的波形绘制在MATLAB中一般用stem函数。stem函数的基本用法和plot函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill”、“filled”，或者参数“.”。由于MATLAB中矩阵元素的个数有限，所以MATLAB只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。 例1. 单位取样序列 单位取样序列，也称为单位冲激序列，定义为 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n=0处是取确定的值1。在MATLAB中，冲激序列可以通过编写以下的impDT.m文件来实现，即 function y=impDT(n) y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n必须为整数或整数向量。 【实例2-1】 利用MATLAB的impDT函数绘出单位冲激序列的波形图。 解：MATLAB源程序为 n=-3:3; x=impDT(n); stem(n,x,fill),xlabel(n),grid on axis([-3 3 -0.1 1.1]) title(单位阶跃序列) 其波形如下图 图2-2 单位阶跃序列 1.2.2 离散时间信号的基本运算 图2-2 单位阶跃序列 对离散时间序列实行基本运算可得到新的序列，这些基本运算主要包括加、减、乘、除、移位、反折等。两个序列的加减乘除是对应离散样点值的加减乘除，因此，可通过MATLAB的点乘和点除、序列移位和反折来实现，与连续时间信号处理方法基本一样。 例2 用MATLAB命令画出下列离散时间信号的波形图。 （1）； （2） 解：设，，MATLAB源程序为 a=0.8;N=8;n=-12:12; x=a.^n.*(uDT(n)-uDT(n-N)); n1=n;n2=n1-3; subplot(121) stem(n1,x,fill),grid on title(x1(n)),axis([-15 15 0 1]) subplot(122) stem(n2,x,fill),grid on title(x2(n)),axis([-15 15 0 1] 其波形如下图 1.3 编程练习 试用MATLAB命令分别绘出下列各序列的波形图。 （1） （2） （3） （4） （5） （6） 解： 附： n=0:10; a1=1/2;a2=-1/2;a3=2;a4=-2;a5=2;a6=1/2; x1=a1.^n;x2=a2.^n;x3=a3.^n;x4=a4.^n; x5=a5.^(n-1).*uDT(n-1);x6=a6.^(n-1); subplot(231) stem(n,x1,fill),grid on xlabel(n),title(x(n)=1/2^{n}) subplot(232) stem(n,x2,fill),grid on xlabel(n),title(x(n)=(-1/2)^{n}) subplot(233) stem(n,x3,fill),grid on xlabel(n),title(x(n)=2^{n}) subplot(234) stem(n,x4,fill),grid on xlabel(n),title(x(n)=(-2)^{n}) subplot(235) stem(n,x5,fill),grid on xlabel(n),title(x(n)=2^{n-1}*u(n-1)) subplot(236)