已知系统模型、初始状态、输入，三种方法求输出

已知系统模型、初始状态、输入，三种方法求输出

用matlab进行计算微分方程解结果一样。

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

值得考虑的是在输入u为序列，而不是t的显式函数时，如何求输出序列？目前思考的是用卷积方法应该可以。
目前在matlab里面用微分方程尝试了对任意时间序列ut求输出序列，结果不对。matlab自带的lsim函数很好的解决了这个问题。想用在python里需要自己写代码，阅读lsim函数改写也并不容易。在尝试解决，期待早点解决。