How to use extra time varying discrete inputs in Ode45?

Question

hassnain iqbal 2016-1-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/264141-how-to-use-extra-time-varying-discrete-inputs-in-ode45

评论： hassnain iqbal 2016-1-18

For example as shown in the function below: Where Ts is the variable with initial value, while Qu is input but it changes with time having specific values defined in the matrix.

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

I am using in this way: 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

But its giving error... please explain how to use inputs with changing time with this example.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2016-1-18

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/264141-how-to-use-extra-time-varying-discrete-inputs-in-ode45#answer_206511

在 MATLAB Online 中打开

function main
Qu=[0 0 0 0 0 0 0 0 21 41 60 75];
[t,Ts]=ode45(@(t,Ts)temp(t,Ts,Qu),(1:1:12),45) 
function dTdt=temp(t,Ts,Qu)
Qu_actual = interp1(1:1:12,Qu,t);
dTdt = Qu_actual/(1500*4190)-12*10^6/(1500*4190)-11.1*3600*(Ts-20)/(1500*4190);
end