信号与系统报告实验1连续时间信号的MATLAB表示与计算.doc

信号与系统 实验报告 实验一 连续时间信号的MATLAB表示与计算 连续时间信号的MATLAB表示与计算 实验目的 初步学习MATLAB语言，熟悉MATLAB软件的基本使用。 掌握用MATLAB描述连续时间信号方法，能够编写MATLAB程序，实现各种信号的时域变换和运算，并且以图形的方式再现各种信号的波形。 实验的内容与步骤 （1）、编写m文件，将实验原理中所给的单位冲激信号和单位阶跃信号的函数文件在MATLAB文件编译器中编写好，并分别以文件名chongji和heaviside存入文件夹中以便使用。 冲激函数程序源码： 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]) 阶跃函数程序源码： function y=heaviside(t) y=(t=0); （2）、编写m文件并以Q1_1为文件名存盘，产生信号f（t）=sin（2πt），时间范围在-2~2s之间，dt=0.01.若将dt改为0.2，再执行改程序，保存图形，看看所得的效果如何。 程序源码： %%sin(2*pai*t) A=1; w0=2*pi; phi=0; t1=-2:0.01:2; t2=-2:0.2:2; f1 = A*sin(w0*t1+phi); f2 = A*sin(w0*t2+phi); subplot(2,1,1); plot(t1,f1); title(f1(dt=0.01)); grid on; subplot(2,1,2); plot(t2,f2) title(f2(dt=0.2)); grid on; 图形 （3）、编写m文件，并以Q1_2为文件名存盘绘出e-2（u（t）-u（t-3））信号，给图形标出名称和坐标轴名称。 程序源码： syms t ft=sym(1*exp(-2*t)*(heaviside(t)-heaviside(t-3))); ezplot(ft); title(exp(-2*t)*((u(t)-u(t-3))); xlabel(t); ylabel(ft); grid on; 图形： （4）编写m文件，并以Q1_3为文件名存盘绘出f（t）=Sa(2 t),f(2t-2)的波形，t的范围在-2~2s。 程序源码： t=-2:0.01:2 ft=sinc(t*2); subplot(1,2,1); plot(t,ft); title(ft=Sa(2t)); grid on; f1=sinc((t*2-2)*2); subplot(1,2,2); plot(t,f1); title(f(2t-2)); grid on; 图形： （5）编写m文件，并以Q1_4为文件名存盘。实现如图1—14所示的信号f（t）。 程序源码： t=-5:0.01:5; f1=2*heaviside(t+1)-heaviside(t-1)-2*heaviside(t-2)+heaviside(t-3); plot(t,f1); grid on; title(f(t)) axis([-5 5 -3 3]); 图形： （6）编写m文件，并以Q1_5为文件名存盘绘出f(t)=e-10.1tsin(2/3t)的波形，t的范围在0~30s。 程序源码： t=0:0.01:30; ft=exp(-0.1*t).*sin(2/3*t); plot(t,ft); grid on; title(f(t)=e-10.1tsin(2/3t)); xlabel(t); ylabel(ft); 图形： 实验中遇到的问题 实验中经常出现各种错误，往往因为各种细小错误的出现使得整个程序运行错误。?第1题中要注意所定义的函数名必须与所保存的文件名一致，单位阶跃信号u(t)和单位冲激信号?(t)需要输入数据才能得到相应的信号波形。第2题中的信号是由两个信号函数相加而成，并不是两个信号波形的叠加，故所得波形中间没有一横。其中第三个图是由两个三角波脉冲信号相减而成，开始时用三角波脉冲信号减去矩形脉冲信号，运行并没有得到所要的信号波形。第4题需要增大时间t的取值范围才能得到函数的多个周期波形，以便于判断函数是否为周期性信号。????? 在这次实验中我体会到：实验就是一个发现错误并改正错误的过程。正因为有错误的出现才显示出实验的魅力。