如何在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

但是它给出了一个错误。请通过此示例说明如何使用输入更改时间。