matlab - 如何在 Ode45 中使用额外的时变离散输入？

Question

例如如下函数所示：Ts具有初始值的变量在哪里，虽然Qu是输入，但它随时间变化，具有在矩阵中定义的特定值。

M 文件：

function Ts=temp(t,Ts,Qu)
Ts=Qu/(1500*4190)-12*10^6/(1500*4190)-11.1*3600*(Ts-20)/(1500*4190);

我以这种方式使用：

Qu=[0 0 0 0 0 0 0 0 21 41 60 75]
 for j=1:12
 Qu=Qu(t);
 [t,Ts]=ode45(@(t,Ts)temp(t,Ts,Qu),(1:1:12),45)
 end

但它给出了一个错误。请通过此示例说明如何使用随时间变化的输入。