信号与系统实验报告2.docVIP

实验2 连续时间信号在MATLAB中的运算 实验目的 学会运用MATLAB进行连续信号的时移、反折和尺度变换；学会运用MATLAB进行连续信号的相加、相乘运算；学会运用MATLAB数值计算方法求连续信号的卷积。 实验原理 2.1信号的时移、反折和尺度变换 信号的时移、反折和尺度变换是针对自变量时间而言的，其数学表达式与波形变换之间存在一定的变换规律。 信号的时移就是将信号数学表达式中的用替换，其中为正实数。因此，波形的时移变换是将原来的波形在时间轴上向左或者向右移动。为波形向左移动；为波形向右移动。信号的反折就是将表达式中的自变量用替换，即变换后的波形是原波形的y轴镜像。信号的尺度变换就是将表达式中的自变量用替换，其中，为正实数。对应于波形的变换，则是将原来的的波形以原点为基准压缩（）至原来的，或者扩展（）至原来的。 上述可以推广到的情况。 exp(-t).*u(t) m×n n×a real( ) imag( ) abs( ) angle( ) 2.2 MATLAB数值计算法求连续时间信号的卷积 用MATLAB分析连续时间信号，可以通过时间间隔取足够小的离散时间信号的数值计算方法来实现。可调用MATLAB中的conv( )函数近似地数值求解连续信号的卷积积分。如果对连续时间信号和进行等时间间隔均匀抽样，则和分别变为离散序列和。其中为整数。当足够小时，和即为连续时间信号和。因此连续信号的卷积积分运算转化为： 采用数值计算法，只求当时卷积积分的值，其中，为整数，即 其中，实际就是离散序列和的卷积和。当足够小时，就是卷积积分的结果，从而连续时间信号 上式表明通过MATLAB实现连续信号和的卷积，可以利用各自抽样后的离散时间序列的卷积再乘上抽样间隔。抽样间隔越小，误差也就越小。 实例分析 例1 已知信号，试用MATLAB命令画出、、、的波形图。 解：根据已知信号，先建立函数文件，即在MATLAB的工作目录下创建funct1.m文件，MATLAB源程序为： function f = funct1(t) f=uCT(t+2)-uCT(t)+(-t+1).*(uCT(t)-uCT(t-1)); 然后，可以调用上述函数来绘制所求的信号波形。程序运行完，产生如图1所示的波形。MATLAB源程序为： clear;clc; t = -2:0.01:4; ft1 = funct1(t-2); ft2 = funct1(3*t); ft3 = funct1(-t); ft4 = funct1(-3*t-2); subplot(221) plot(t,ft1);grid on title(f(t-2)); axis([-2 4 -0.5 2]); subplot(222) plot(t,ft2);grid on title(f(3t)); axis([-2 4 -0.5 2]); subplot(223) plot(t,ft3);grid on title(f(-t)); axis([-2 4 -0.5 2]); subplot(224) plot(t,ft4);grid on title(f(-3t-2)); axis([-2 4 -0.5 2]); 图1 f(t-2),f(3t),f(-t),f(-3t-2)波形图 例2 用MATLAB数值计算分析法求信号与的卷积积分。 解：因为是一个持续时间无限长的信号，而计算机数值计算不可能计算真正的无限长时间信号，所以在进行的抽样离散化时，所取的时间范围让衰减到足够小就可以了，本例中取。MATLAB源程序为： clear;clc; dt = 0.01; t = -1:dt:2.5; f1 = uCT(t)- uCT(t-2); f2 = exp(-3*t).*uCT(t); f = conv(f1,f2)*dt; n =length(f); tt = (0:n-1)*dt-2; subplot(221), plot(t,f1),grid on; axis([-1, 2.5, -0.2,1.2]); title(f1(t)); xlabel(t); subplot(222), plot(t,f2),grid on; axis([-1, 2.5, -0.2,1.2]); title(f2(t)); xlabel(t); subplot(212), plot(tt,f),grid on; title(f(t)=f1(t)*f2(t)); xlabel(t);  程序运行后，产生如图2所示的波形。 图2 例2的卷积结果 由于和的时间范围都是从开始，所以卷积结果的时间范围从开始，增量