MATLAB实验报告卷积.doc

《MATLAB实验报告卷积.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验报告学院：机电班级：姓名：学号：实验名称：连续时间信号卷积运算的MATLAB实现1.实验目的：掌握卷积的概念及计算方法2.熟悉通过调用conv()函数求解连续时间信号卷积的数值分析法实验环境：MATLAB6.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）过调

2、用conv（）函数编程实现卷积计算y(t)的数值解并画图实验程序及结果：第一题：M文件(1)functionf=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.52.5-0.

3、51.5])gridonsubplot(312);plot(k2,f2)gridonaxis([-0.52.5-0.51.5])xlabel('t')ylabel('f2(t)')subplot(313)plot(k,y)axis([-0.54-0.51.5])gridonxlabel('t')ylabel('f1(t)*f2(t)')实验结果：第二题：M文件functionf=uCT(t)f=(t>=0);主程序：symstao>>t=sym('t','positive');>>f=exp(-t);>>h=t^2*exp(-2*t);>>fh_tao=subs(f,t,

4、tao)*subs(h,t,t-tao);>>yt=int(fh_tao,0,t)yt=-(-2*exp(t)+t^2+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)')gridonaxis([-0.53.5

5、-0.51.5])subplot(222)plot(k2,h);xlabel('t')ylabel('h(t)')gridonaxis([-0.511-0.050.2])subplot(223)plot(k,y);gridonxlabel('t')ylabel('f(t)*h(t)数值')axis([-0.515-0.010.1])subplot(224)yt=-(-2.*exp(k)+k.^2+2.*k+2)./exp(k).^2;plot(k,yt);gridonxlabel('t')ylabel('f(t)*h(t)理论')axis([-0.515-0.010.1]

6、)实验结果：