function  [f,k]=sconv(f1,f2,k1,k2,p)
%計(jì)算連續(xù)信號卷積積分f(t)=f1(t)*f2(t)
% f:  卷積積分f(t)對應(yīng)的非零樣值向量
% k： f(t)的對應(yīng)時(shí)間向量
% f1: f1(t)非零樣值向量
% f2: f2(t)的非零樣值向量
% k1: f1(t)的對應(yīng)時(shí)間向量
% k2:  序列f2(t)的對應(yīng)時(shí)間向量
% p：  取樣時(shí)間間隔
f=conv(f1,f2); 									%計(jì)算序列f1與f2的卷積和f
f=f*p;
k0=k1(1)+k2(1);									%計(jì)算序列f非零樣值的起點(diǎn)位置
k3=length(f1)+length(f2)-2;						%計(jì)算卷積和f的非零樣值的寬度
k=k0:p:k3*p;                                   %確定卷積和f非零樣值的時(shí)間向量
subplot(2,2,1)
plot(k1,f1)										%在子圖1繪f1(t)時(shí)域波形圖
grid on
title('f1(t)')
xlabel('t')
ylabel('f1(t)')
subplot(2,2,2)
plot(k2,f2)										%在子圖2繪f2(t)時(shí)波形圖
grid on

title('f2(t)')
xlabel('t')
ylabel('f2(t)')
subplot(2,2,3)
plot(k,f);										%畫卷積f(t)的時(shí)域波形
grid on

h=get(gca,'position');								
h(3)=2.5*h(3);
set(gca,'position',h)				%將第三個(gè)子圖的橫坐標(biāo)范圍擴(kuò)為原來的2.5倍
title('f(t)=f1(t)*f2(t)')
xlabel('t')
ylabel('f(t)')