实验报告
实验名称:连续时间信号的MATLAB表示与计算
实验班级:电气11-4 学号(后两位): 06
试验台号: 01
指导教师:张秋晶
试验日期: 2015年4月16日
一、实验目的
1、初步学习MATLAB 语言,熟悉MATLAB 软件的基本使用。
2、掌握用MATLAB 描述连续时间信号方法,能够编写MATLAB 程序,实现各种信号的时域变换和运算,并且以图形的方式再现各种信号的波形。
二、实验内容与步骤
1、建立以自己的汉语名字命名一个文件夹,用以存放自己所编制的M 文件和产生的结果,并设置成当前工作目录。
2、编写M 文件,将实验原理中所给的单位冲激信号和单位阶跃信号的函数文件在MATLAB 文件编辑器中编写好,并分别以文件名chongji 和ut 存入work 文件夹中以便于使用。
抄写函数文件chongji 如下: function chongji(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); x=zeros(1,n); x(1,(t0-t1)/dt+1)=1/dt; stairs(t,x); %以阶梯方式绘画 axis([t1,t2,0,1.1/dt]) 抄写函数文件ut 如下:
function f=ut(t) %阶跃信号
f=(t>0); % y = 1 for t > 0, else y = 0
3、编写M 文件,并以Q1_1为文件名存盘,产生信号()sin(2π)f t t , 时间范围在-2~2秒之间,dt = 0.01。若将dt 改为0.2,再执行该程序,保存图形,看看所得图形的效果如何? 程序如下:t1=[-2:0.01:2];
f1=sin(2*pi*t1);
figure(1);
plot(t1,f1);
title('dt=0.01');
t2=[-2:0.2:2];
f2=sin(2*pi*t2);
figure(2);
plot(t2,f2);
title('dt=0.2');
dt = 0.01时的信号波形: dt = 0.2时的信号波形
这两幅图形有什么区别,哪一幅图形看起来与实际信号波形更像?
答:第一幅图的锯齿不明显,第二幅图有很明显的锯齿。第一幅图图形看起来与实际信号
更像。
4、编写M 文件,并以Q1_2为文件名存盘绘出 ))3()((2---t u t u e t 信号,给图形标出名称和坐标轴名称。
程序如下: t=-5:0.01:5;
f=exp(-2*t).*(ut(t)-ut(t-3));
plot(t,f);
title('Q_2');xlabel('X 轴');ylabel('y 轴');
5、编写M 文件,并以Q1_3为文件名存盘绘出)π2(Sa )(t t f =、f (2t-2)的波形,t 的范围在-2 ~2s 。(结果见图一)
程序如下:
t1=-2:0.01:2 ;
f1=sinc(2*t1);
subplot(1,2,1);
plot(t1,f1);
t2=-1:0.01:2 ;
f2=sinc(2*(2*t2-2));
subplot(1,2,2);
plot(t2,f2)
;
6、编写M 文件,并以Q1_4为文件名存盘。实现如图二所示的信号f (t)。
程序如下:
a=[0 2 2 1 -1 0 ];
t=1:length(a);
stairs(t-3,a);
axis([-2 length(a) min(a)-1 max(a)+1])
7、编写M 文件,并以Q1_5为文件名存盘绘出0.12()sin()3
t f t e t -=的波形,t 的范围在0 ~30s 。 程序如下:
t=0:30;
x1=exp(-0.1*t).*sin((2/3)*t );
plot(t,x1);
图一 f (t)
t
123
-11
2
-1
图二
三.实验反思
1.实验前没有做好预习,相关函数不知道是什么意思。
2.在编写程序是遇到了很多细节问题,比如全角和半角的切换。
3.遇到问题后我会自己尝试修改,最后再问老师。
4.实验后要加强练习,理解透每个程序。
如有侵权请联系告知删除,感谢你们的配合!