%繪制矩形信號拉普拉斯變換曲面圖程序
clf;
a=-0:0.1:5;
b=-20:0.1:20;					
[a,b]=meshgrid(a,b);
c=a+i*b;						%確定繪圖區域
c=(1-exp(-2*c))./c;
c=abs(c);						%計算拉普拉斯變換
mesh(a,b,c);					%繪制曲面圖
surf(a,b,c);
view(-60,20)					%調整觀察視角
axis([-0,5,-20,20,0,2]);
title('拉普拉斯變換（S域像函數）');
colormap(hsv);