time dependent parameter in ode-function matlab

1 次查看（过去 30 天）
sejo
sejo 2014-7-23
评论： kiran gk 2016-11-20
I want to solve an ode with a time dependent parameter. cA sould be 10000 if t is >=10 and <=11 else it should have the value of 0. cA is then used in an differential equation to calculate cB. See the code:
function dcB = myode(t,y)
cB=y(1,:);
if t>=10 && t<=11
cA=10000
else
cA=0
end
dcB=(cA-cB)*100/1750;
[t,y]=ode45(@myode,[tdown tup],0);
Fallowing problems show up:
  • if i print cA it has not the correct values at the specified times.
  • if tup is e.g. 20 cB has values, if tup is e.g. 100 cB is zero.

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

采纳的回答

Mischa Kim
Mischa Kim 2014-7-23
Sejo, one fairly clean approach would be to split the integration interval: [tdown 10),[10 11],(11 tup] and call ode45 three times with two separate ODE functions. You'd just have to make sure that the final conditions of the respective previous integration are passed on as initial conditions to the next integration.
  1 个评论
kiran gk
kiran gk 2016-11-20
What if cA has multiple values with respect to t? In that case it will be difficult to create multiple ode function please suggest the possible answer

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by