信号与信号实验.doc

信号与信号实验 MATLAB部分 实验一：基本信号在MATLAB中的表示和运算 实验目的; 1、学会用MATLAB表示常用连续信号的方法； 2、学会用MATLAB进行信号基本运算的方法； 3、学会用MATLAB实现连续时间信号的卷积的方法。 实验内容： 绘出下列信号的时域波形 (1)f(t)=(2-e-2t)u(t) (2)f(t)=cos(πt)[u(t)-u(t-1)] (3)f(t)=u(-3t+2) (4)f(t)= -(1/2)tu(t+2) 解：t1=0:0.01:5; y1=(2-exp(-2*t1)).*(t10); subplot(221);plot(t1,y1);grid; title(f(t)=(2-e-2t)u(t)); t2=0:0.01:5; y2=cos(pi*t2).*((t20)-(t21)); subplot(222);plot(t2,y2);grid; title(f(t)=cos(πt)[u(t)-u(t-1)]); t3=-2:0.01:5; y3=(-3*t3+20); subplot(223);plot(t3,y3);grid; title(f(t)=u(-3t+2)); t4=-3:0.01:5; y4=(-1/2)*t4.*(t4-2); subplot(224);plot(t4,y4);grid; title(f(t)=-(1/2)tu(t+2));  用MATLAB绘出下列信号的卷积积分f1(t)*f2(t)的时域波形 (1) f1(t)=tu(t), f2(t)=u(t) (2) f1(t)=u(t)-u(t-4), f2(t)=sin(πt)u(t) (3) f1(t)= e-2t u(t), f2(t)= e-t u(t) (4) f1(t)= e-t u(t), f2(t)=u(t) 解： （1）fs=1000; t=-1:1/fs:4; x1=stepfun(t,0); x2=x1.*t; y=conv(x1,x2)/fs; n=length(y1); tt=(0:n-1)/fs-2; subplot(311),plot(t,x1),grid; title(f1(t)=tu(t)); subplot(312),plot(t,x2),grid; title( f2(t)=u(t)); subplot(313),plot(tt,y),grid on; title(f1(t) * f2(t)); （2）fs=1000; t=-1:1/fs:4; x1=(t0)-(t4); x2=sin(pi*t).*(t0); x=conv(x1,x2)/fs; n=length(x); tt=(0:n-1)/fs-2; subplot(311);plot(t,x1);grid; title(f1(t)=u(t)-u(t-4))); subplot(312);plot(t,x2);grid; title(f2(t)=sin(πt)u(t)); subplot(313);plot(tt,x);grid; title(f1(t) * f2); （3）t=0:1/fs:4; x1=exp(-2*t).*(t0); x2=exp(-t).*(t0); x=conv(x1,x2)/fs; n=length(x); tt=(0:n-1)/fs-0; subplot(311);plot(t,x1);grid; title(f1(t)= e-2t u(t)); subplot(312);plot(t,x2);grid; title(f2(t)= e-t u(t)); subplot(313);plot(tt,x);grid; title(f1(t) * f2(t)); （4）t=0:1/fs:2; x1=exp(-2*t).*(t0); x2=(t0); x=conv(x1,x2)/fs; n=length(x); tt=(0:n-1)/fs-0; subplot(311);plot(t,x1);grid; title( f1(t)= e-t u(t))); subplot(312);plot(t,x2);grid; title(f2(t)