tau=1;  % 時間延遲常數
[nP1,dP1]=paderm(tau,1,3); GP1=tf(nP1,dP1)

[nP2,dP2]=paderm(tau,0,3); GP2=tf(nP2,dP2)

[nP3,dP3]=pade(tau,3); GP3=tf(nP3,dP3)

step(GP1,GP2,GP3)  % 繪制出三種近似的階躍響應曲線
h=line([0,1,1.00001,7],[0,0,1,1]); % 畫延遲模塊理論值
set(h,'LineWidth',2); % 設置其為粗線