实验三卷积积分与卷积和运算.docVIP

a=1000; t1=-5:1/a:5; f1=stepfun(t1,0); f2=stepfun(t1,-1/a)-stepfun(t1,1/a); subplot(231); plot(t1,f1); axis([-5,5,0,1.2]); xlabel(时间(t));ylabel(幅值f1(t));title(f1); subplot(232); plot(t1,f2); ylabel(幅值f2(t)); title(f2); y=conv(f1,f2); r=2*length(t1)-1; t=-10:1/a:10; subplot(233); plot(t,y); axis([-5,5,0,1.2]); title(f1和f2的卷积); ylabel(f(t)); f11=conv(f1,f1);f22=conv(f2,f2); subplot(234); plot(t,f11); title(f1和f1的卷积); ylabel(f11(t)); axis([-5,5,0,5000]); subplot(235); plot(t,f22); title(f2和f2的卷积); ylabel(f22(t)); s=0.01; k1=0:s:2; k2=k1; f1=3*k1; f2=3*k2; f=conv(f1,f2); f=f*s; k0=k1(1)+k2(1); k3=length(f1)+length(f2)-2; k=k0:s:k3*s; subplot(3,1,1); plot(k1,f1); title(f1(t)); subplot(3,1,2); plot(k2,f2); title(f2(t)); subplot(3,1,3); plot(k,f); title(f(t)) clear all; T=0.1; t=0:T:10; f=sin(t); h=0.5*(exp(-t)+exp(-3*t)); Lf=length(f); Lh=length(h); for k=1:Lf+Lh-1 y(k)=0; for i=max(1,k-(Lh-1)):min(k,Lf) y(k)=y(k)+f(i)*h(k-i+1) ; end yzsappr(k)=T*y(k) ; end subplot(3,1,1); plot(t,f); title(f(t)); subplot(3,1,2); plot(t,h); title(h(t)); subplot(3,1,3); plot(t,yzsappr(1:length(t))); title(卷积近似计算结果); xlabel(时间); a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel(n);ylabel(幅值); k1=-5:15; f1=[zeros(1,5),ones(1,16)]; subplot(3,1,1); stem(k1,f1); title(f1(k)); k2=k1; f2=[zeros(1,5),ones(1,3),zeros(1,13)]; subplot(3,1,2); stem(k2,f2); title(f2(k)); k3=k1(1)+k2(1):k1(end)+k2(end); f3=conv(f1,f2); subplot(3,1,3); stem(k3,f3); title(f3(k)); a=1000; t1=-5:1/a:5; f1=stepfun(t1,-1/a)-stepfun(t1,1/a); f2=stepfun(t1,0); subplot(331); plot(t1,f1); axis([-5,5,0,1.2]); xlabel(时间(t));ylabel(幅值f1(t));title(单位冲激f1); subplot(332); plot(t1,f2); axis([-5,5,0,1.2]); xlabel(时间(t)); ylabel(幅值f2(t)); title(单位阶跃f2); t=0:0.0001:8; t2=0;t3=4; u1=stepfun(t,t2); u2=stepfun(t,t3); subplot(333); f3=u1-u2;plot(t,u); xlabel(时间(t));ylabel(幅值f3(t));title(u(t)-u(t-4)); y=conv(f1,f2); r=2*length(t1)-1; t=-10: