已知系統模型、初始狀态、輸入，三種方法求輸出

用matlab進行計算微分方程解結果一樣。

>> y=dsolve('Dy+y=1','y(0)=0','t')
 
y =
 
1 - exp(-t)
           

值得考慮的是在輸入u為序列，而不是t的顯式函數時，如何求輸出序列？目前思考的是用卷積方法應該可以。

目前在matlab裡面用微分方程嘗試了對任意時間序列ut求輸出序列，結果不對。matlab自帶的lsim函數很好的解決了這個問題。想用在python裡需要自己寫代碼，閱讀lsim函數改寫也并不容易。在嘗試解決，期待早點解決。