gggfconv.m

the code of the book come form the book of marlix laboratory with the book of singal

function  [f,t]=gggfconv(f1,f2,t1,t2)
%计算连续信号卷积积分实用函数
d=input('请输入取样时间间隔d: ');
f=conv(f1,f2); 						%计算序列f1与f2的卷积和f
f=f*d;                                 %计算卷积积分信号f(t)离散样值
ts=t1(1)+t2(1)							%计算序列f非零样值的起点位置
l=length(t1)+length(t2)-2;			   %计算卷积积分f的非零样值的宽度
t=ts:d:(ts+l*d)                        %确定卷积积分f非零样值的时间向量
subplot(2,2,1)
plot(t1,f1)							    %在子图1绘制信号f1(t)的时域波形
axis([min(t1),max(t1),min(f1)-abs(min(f1)*0.2),max(f1)+max(f1)*0.2])
title('f1(t) ')
xlabel('t')
subplot(2,2,2)
plot(t2,f2)								%在子图2绘制信号f2(t)的时域波形
axis([min(t2),max(t2),min(f2)-min(f2)*0.2,max(f2)+max(f2)*0.2])
title('f2(t) ')
xlabel('t')
subplot(2,2,3)
plot(t,f);						       %在子图3绘制卷积积分f(t)的时域波形
axis([min(t),max(t),min(f)-min(f)*0.2,max(f)+max(f)*0.2])
p=get(gca, 'position');								
p(3)=2.4*p(3);
set(gca, 'position',p)			%将第三个子图的横坐标范围扩为原来的2.4倍
title('f(t)=f1(t)*f2(t)')
xlabel('t')