[jlfocus]

首先要谢谢上面的高人，不过我的问题现在有所改变：
形如AX(导数)=BX+u的状态方程，如果A、B、u都是常系数，用ode45就可以解决问题
但是如果u=f(t)，不再是常数，ode45求解时就会报错
请问这个问题怎么解决？