Matlab信号与系统仿真.doc

电子信息工程系实验报告成 绩：课程名称： 计算机仿真技术 指导教师（签名）：实验项目名称： 实验四 信号与系统仿真 实验时间： 2010-12-1 班级： xxxx 姓名： xxx 学号： xxxxxxx 一、实验目的学会运用MATLAB进行连续信号的时移、反折和尺度变换；学会运用MATLAB进行连续信号的相加、相乘运算；学会运用MATLAB数值计算方法求连续信号的卷积。二、实验环境硬件：PC机，基本配置CPU PII以上，内存256M以上；软件：Matlab 版本7.1三、实验原理 信号的时移、反折和尺度变换是针对自变量时间而言的，其数学表达式与波形变换之间存在一定的变换规律。 用MATLAB分析连续时间信号，可以通过时间间隔取足够小的离散时间信号的数值计算方法来实现。四、实验内容及结果分析1 .试用MATLAB命令绘制信号的波形图。解：新建一个项目，MATLAB源程序为：clear; clc;t = 0:0.01:3;w1=10*pi;w2=9*pi;ft = exp(-t).*sin(w1*t)+exp(-t/2).*sin(w2*t);plot(t,ft);grid on axis(0,3,-2,2); 运行后波形图如图1所示： 图12 .已知信号，画出、的波形。解：（1）先在MATLAB的工作目录下创建uCT的M文件，其MATLAB源文件为：function f = uCT(t) f = (t=0);（2）然后建立函数文件，即在MATLAB的工作目录下创建funct1.m文件，MATLAB源程序为： function f = funct1(t)f=uCT(t)-uCT(t-1)+(t-1).*(uCT(t+1)-uCT(t);（3）调用上述函数来绘制所求的信号波形。程序运行完，产生如图2所示的波形。MATLAB源程序为：clear;clc; plot(t,ft2);grid ont = -2:0.01:3; title(f(t+2);ft1 = funct1(t); axis(-2 3 -2 1.5);ft2 = funct1(t+2); subplot(223)ft3 = funct1(-t); plot(t,ft3);grid onft4 = funct1(-2*t+1); title(f(-t);subplot(221) axis(-2 3 -2 1.5);plot(t,ft1);grid on subplot(224)title(f(t); plot(t,ft4);grid onaxis(-2 3 -2 1.5); title(f(-2t+1);subplot(222) axis(-2 3 -2 1.5); 图 2 、的波形图 3 求信号与的卷积结果，并画出和的波形。解：新建一个项目，MATLAB源程序为：clear;clc; axis(-1, 3, -0.2,1.2);dt = 0.01; t = -1:dt:4; title(f1(t); xlabel(t);f1 = uCT(t)- uCT(t-2); subplot(222), plot(t,f2),grid on;f2 = uCT(t)+uCT(t-1)- uCT(t-2)-uCT(t-3); axis(-1, 4, -0.2,2.2);f = conv(f1,f2)*dt; title(f2(t); xlabel(t);n =length(f); subplot(212), plot(tt,f),grid on;tt = (0:n-1)*dt-2; axis(-2, 7, -0.2,3.5);subplot(221), plot(t,f1),grid on; title(f(t)=f1(t)*f2(t); xlabel(t);程序运行后，产生如图3 所示的波形： 图3 和的波形图4 求信号与自身的卷积结果，并画出和的波形。解：新建一个项目，MATLAB源程序为：clear;clc;dt = 0.01; t = -1:dt:2;f1 = uCT(t+0.5)- uCT(t-0.5);f = conv(f1,f1)*dt;n =length(f);tt = (0:n-1)*dt-2;subplot(221), plot(t,f1),grid on;axis(-1, 1, -0.2,1.2);title(f1(t); xlabel(t);subplot(222), plot(tt,f),grid on;axis(-2, 2, -0.2,1.2);title(f(t)=f1(t)*f1(t); xlabel(t);程序运行后，产生如图4 所示的波形： 图 4 和的波形图五、思考题MATLAB运算符中 .*和 * 的区别？可结合例子说明。 答：.* 是数组运算，例如k .*A,是k乘以A中的每个元素。而 * 是矩阵