function [f,k]=sconv(f1,f2,k1,k2,p)
%计算连续信号卷积积分 f(t)=f1(t)*f2(t)
%f:卷积积分f(t)的对应的非零样值向量
%k:f(t)对应的时间向量;
%f1:f1(t)非零样值向量;
%f2:f2(t)非零样值向量;
%k1:f(1)对应时间向量
%k2:f(2)对应时间向量
%p:取样时间间隔
f=conv(f1,f2);
f=f*p;
k0=k1(1)+k2(1);
k3=length(f1)+length(f2)-2;
k=k0:p:k3*p;
subplot 311
plot(k1,f1);
title('f1(t)');
subplot 312
plot(k2,f2);
title('f2(t)')
xlabel('t');
ylabel('f2(t)');
subplot 313
plot(k,f);
h=get(gca,'position');
h(3)=2.5*h(3);
set(gca,'position')
title('f(t)=f1(t)*f2(t)');
xlabel('t')
ylabel('f(t)')
end
评论0