MATLAB实验报告卷积.doc

实 验 报 告学院： 机电 班级： 姓名： 学号： 实验名称：连续时间信号卷积运算的MATLAB实现1. 实验目的：掌握卷积的概念及计算方法2. 熟悉通过调用conv()函数求解连续时间信号卷积的数值分析法实验环境：MATLAB 6.5.1软件实验要求：1、已知信号f1（t）=t/2*(t)- (t-2), f2（t）= (t)- (t-1)，通过调用conv（）函数编程实现卷积计算y(t)= f1（t）* f2（t）,画出波形。2、已知信号f（t）=e-t *(t), h（t）= t2 *e-2t *(t),y(t)= f（t）* h（t）（1）用符号分析法编程实现计算y（t）的理论解；（2）过调用conv（）函数编程实现卷积计算y(t)的数值解并画图实验程序及结果：第一题：M文件 (1) function f=uCT(t)f=(t=0);主程序：p=0.01;k1=0:p:2;k2=0:p:1;f1=k1/2.*uCT(k1)-uCT(k1-2);f2=uCT(k2)-uCT(k2-1);y=conv(f1,f2)*p;k0=k1(1)+k2(1);k3=length(f1)+length(f2)-2;k=k0:p:k3*p+k0;subplot(311)plot(k1,f1);xlabel(t)ylabel(f1(t)axis(-0.5 2.5 -0.5 1.5)grid onsubplot(312);plot(k2,f2)grid onaxis(-0.5 2.5 -0.5 1.5)xlabel(t)ylabel(f2(t)subplot(313)plot(k,y)axis(-0.5 4 -0.5 1.5)grid onxlabel(t)ylabel(f1(t)*f2(t)实验结果：第二题：M文件function f=uCT(t)f=(t=0);主程序：syms tao t=sym(t,positive); f=exp(-t); h=t2*exp(-2*t); fh_tao=subs(f,t,tao)*subs(h,t,t-tao); yt=int(fh_tao,0,t) yt = -(-2*exp(t)+t2+2*t+2)/exp(t)2p=0.01;k1=0:p:3;k2=0:p:12;f=exp(-k1).*uCT(k1);h=(k2).2.*exp(-2.*k2).*uCT(k2);y=conv(f,h)*p;k0=k1(1)+k2(1);k3=length(f)+length(h)-2;k=k0:p:k3*p+k0;subplot(221)plot(k1,f);xlabel(t)ylabel(f(t)grid onaxis(-0.5 3.5 -0.5 1.5)subplot(222)plot(k2,h);xlabel(t)ylabel(h(t)grid onaxis(-0.5 11 -0.05 0.2)subplot(223)plot(k,y);grid onxlabel(t)ylabel(f(t)*h(t)数值)axis(-0.5 15 -0.01 0.1)subplot(224)yt =-(-2.*exp(k)+k.2+2.*k+2)./exp(k).2;plo