《离散时间信号的表示及运算资料》由会员分享，可在线阅读，更多相关《离散时间信号的表示及运算资料（14页珍藏版）》请在金锄头文库上搜索。
1、第2章 离散时间信号的表示及运算2.1 实验目的l 学会运用MATLAB表示的常用离散时间信号；l 学会运用MATLAB实现离散时间信号的基本运算。2.2 实验原理及实例分析2.2.1 离散时间信号在MATLAB中的表示离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用来表示，自变量必须是整数。离散时间信号的波形绘制在MATLAB中一般用stem函数。stem函数的基本用法和plot函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill”、“filled”，或者参数“.”。由于MATLAB中矩阵元素的个数有限，所以MAT
2、LAB只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。1. 单位取样序列单位取样序列，也称为单位冲激序列，定义为 （12-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:3;x=impDT(n);stem(n,x,fill),xlabel(n),grid ontitle(单位冲激序列)axis(-3 3 -0.1 1.1)图2-1 单位冲激序列程序运行结果如图12-1所示。2. 单位阶跃序列单位阶跃序列定义为 （12-2）在MATLAB中，冲激序列可以通过编写uDT.m文件来实现，即function y=uDT(n)y=n=0; %当参数为非负时输出1调用该函数时n也同样必须为整数或整数向量。【实例2-2】 利用MATLAB的uDT函数绘出单位阶跃序列的波形图。解：MATLAB源程序为n=-3:5;x=uDT(n);stem(n,x,fill),xla
4、bel(n),grid ontitle(单位阶跃序列)axis(-3 5 -0.1 1.1)图2-2 单位阶跃序列程序运行结果如图12-2所示。3. 矩形序列矩形序列定义为 （12-3）矩形序列有一个重要的参数，就是序列宽度N。与之间的关系为 因此，用MATLAB表示矩形序列可利用上面所讲的uDT函数。【实例2-3】 利用MATLAB命令绘出矩形序列的波形图。解：MATLAB源程序为n=-3:8;x=uDT(n)-uDT(n-5);stem(n,x,fill),xlabel(n),grid ontitle(矩形序列)axis(-3 8 -0.1 1.1)程序运行结果如图2-3所示。图2-3 矩
5、形序列4. 单边指数序列单边指数序列定义为 （12-4）【实例2-4】 试用MATLAB命令分别绘制单边指数序列、的波形图。解：MATLAB源程序为n=0:10;a1=1.2;a2=-1.2;a3=0.8;a4=-0.8;x1=a1.n;x2=a2.n;x3=a3.n;x4=a4.n;subplot(221)stem(n,x1,fill),grid onxlabel(n),title(x(n)=1.2n)subplot(222)stem(n,x2,fill),grid onxlabel(n),title(x(n)=(-1.2)n)subplot(223)stem(n,x3,fill),grid
6、 onxlabel(n),title(x(n)=0.8n)subplot(224)stem(n,x4,fill),grid onxlabel(n),title(x(n)=(-0.8)n)图2-4 单边指数序列单边指数序列的取值范围为。程序运行结果如图12-4所示。从图可知，当时，单边指数序列发散；当时，该序列收敛。当时，该序列均取正值；当时，序列在正负摆动。5. 正弦序列正弦序列定义为 （12-5）其中，是正弦序列的数字域频率；为初相。与连续的正弦信号不同，正弦序列的自变量n必须为整数。可以证明，只有当为有理数时，正弦序列具有周期性。【实例2-5】 试用MATLAB命令绘制正弦序列的波形图。解
7、：MATLAB源程序为n=0:39; x=sin(pi/6*n);stem(n,x,fill),xlabel(n),grid ontitle(正弦序列)axis(0,40,-1.5,1.5);程序运行结果如图2-5所示。图2-5 正弦序列6. 复指数序列复指数序列定义为 （2-6）当时，得到虚指数序列，式中是正弦序列的数字域频率。由欧拉公式知，复指数序列可进一步表示为 （2-7）与连续复指数信号一样，我们将复指数序列实部和虚部的波形分开讨论，得出如下结论：（1）当时，复指数序列的实部和虚部分别是按指数规律增长的正弦振荡序列；（2）当时，复指数序列的实部和虚部分别是按指数规律衰减的正弦振荡序列；
8、（3）当时，复指数序列即为虚指数序列，其实部和虚部分别是等幅的正弦振荡序列。【实例2-6】 用MATLAB命令画出复指数序列的实部、虚部、模及相角随时间变化的曲线，并观察其时域特性。解：MATLAB源程序为n=0:30;A=2;a=-1/10;b=pi/6;x=A*exp(a+i*b)*n);subplot(2,2,1)stem(n,real(x),fill),grid ontitle(实部),axis(0,30,-2,2),xlabel(n)subplot(2,2,2)stem(n,imag(x),fill),grid ontitle(虚部),axis(0,30,-2,2) ,xlabel(
9、n)subplot(2,2,3)stem(n,abs(x),fill),grid ontitle(模),axis(0,30,0,2) ,xlabel(n)subplot(2,2,4)stem(n,angle(x),fill),grid ontitle(相角),axis(0,30,-4,4) ,xlabel(n)图2-6 复指数序列程序运行后，产生如图2-6所示的波形。2.2.2 离散时间信号的基本运算 对离散时间序列实行基本运算可得到新的序列，这些基本运算主要包括加、减、乘、除、移位、反折等。两个序列的加减乘除是对应离散样点值的加减乘除，因此，可通过MATLAB的点乘和点除、序列移位和反折来实
10、现，与连续时间信号处理方法基本一样。【实例2-7】 用MATLAB命令画出下列离散时间信号的波形图。（1）；（2）（3）；（4）解：设，MATLAB源程序为a=0.8;N=8;n=-12:12;x=a.n.*(uDT(n)-uDT(n-N);n1=n;n2=n1-3;n3=n1+2;n4=-n1;subplot(411)stem(n1,x,fill),grid ontitle(x1(n),axis(-15 15 0 1)subplot(412)stem(n2,x,fill),grid ontitle(x2(n),axis(-15 15 0 1)subplot(413)stem(n3,x,fil
11、l),grid ontitle(x3(n),axis(-15 15 0 1)subplot(414)stem(n4,x,fill),grid ontitle(x4(n),axis(-15 15 0 1)图2-7 离散时间信号的基本运算及波形图其波形如图2-7所示。2.3 编程练习1. 试用MATLAB命令分别绘出下列各序列的波形图。（1） （2）（3） （4）（5） （6）2. 试用MATLAB分别绘出下列各序列的波形图。（1） （2）（3） （4）第3章 离散时间LTI系统的时域分析3.1 实验目的l 学会运用MATLAB求解离散时间系统的零状态响应；l 学会运用MATLAB求解离散时间系统
12、的单位取样响应；l 学会运用MATLAB求解离散时间系统的卷积和。3.2 实验原理及实例分析3.2.1 离散时间系统的响应离散时间LTI系统可用线性常系数差分方程来描述，即 （3-1）其中，（，1，N）和（，1，M）为实常数。MATLAB中函数filter可对式（13-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。函数filter的语句格式为y=filter(b,a,x)其中，x为输入的离散序列；y为输出的离散序列；y的长度与x的长度一样；b与a分别为差分方程右端与左端的系数向量。【实例3-1】 已知某LTI系统的差分方程为试用MATLAB命令绘出当激励信号为时，该系统的零状态
13、响应。解：MATLAB源程序为a=3 -4 2;b=1 2;n=0:30;x=(1/2).n;y=filter(b,a,x);stem(n,y,fill),grid onxlabel(n),title(系统响应y(n)程序运行结果如图3-1所示。图3-1 实例3-1系统的零状态响应3.2.2 离散时间系统的单位取样响应系统的单位取样响应定义为系统在激励下系统的零状态响应，用表示。MATLAB求解单位取样响应可利用函数filter，并将激励设为前面所定义的impDT函数。例如，求解实例13-1中系统的单位取样响应时，MATLAB源程序为a=3 -4 2;b=1 2;n=0:30;x=impDT(n);h=filter(b,a,x);stem(n,h,fill),grid onxlabel(n),title(系统单位取样响应h(n)程序运行结果如图3-2所示。图3-2 实例
