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

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:1/a:10;

subplot(334);

plot(t,y);

axis([-5,5,0,1.2]);

title('f1和f2的卷积');

ylabel('y');

y1=conv(f2,f1);

r=2*length(t1)-1;

t=-10:1/a:10;

subplot(335);

plot(t,y1);

axis([-5,5,0,1.2]);

title('f1和f2的卷积');

ylabel('y1');