%dm08102
%求LTI連續系統響應的數值解
a=[1 2 1];
b=[1 2];
sys=tf(b,a);                 %定義系統函數對象
p=0.5;                       %定義采樣時間間隔
t=0:p:5;                     %定義時間范圍向量
f=exp(-2*t);                 %定義輸入信號
y=lsim(sys,f,t)              %求出系統響應數值解