ODE with a parameter that chages at each time step

2011 9 3
3 个回答
回答已采纳
4 次查看（30 天）

0 个投票

Hi all,
I have a problem solving a ODE with Matlab.
The ODE is:
y(t)'of j= [(sigma^2)/(2*k)]*(1-y(t)*k)*(y(t)-[y(t)of j-1])^2;
I've written this code:
mu = 0.20;
r = 0.05;
sigma = 0.40;
T = 10;
k = 0.2;
IC = 0.0;
TSPAN = 0:0.1:T;
TSPANlength = length(TSPAN);
%Number of crahes
N=3;
PiVektor = zeros(N,TSPANlength);
for i=1:1:N
[~, PiVektor(i, :)] = ode23(@(t, pi)myodeN(pi, PiVektor(i-1),k,mu,r,sigma),TSPAN, IC);
PiVektor(i,:) = fliplr(PiVektor(i,:));
figure(1);
hold on
plot(TSPAN,PiVektor(i,:));
end
function dpidt = myodeN(pi,piNminus1,k,mu,r,sigma)
dpidt = [(sigma^2)/(2*k)]*(1-pi*k)*(pi-piNminus1)^2;
end
The problem is that in this way the ODE is always calculated using just the first value from the previous calculation.
Is it possibile to change the value of my parameter at each time step?
Thanks a lot,
Ilaria

2 个评论
显示 无 隐藏 无

Andrew Newell
Andrew Newell 2011-9-4
I don't understand your notation. What does y(t)' of j mean?
ilaria
ilaria 2011-9-4
y(t)' of j is the differential equation for the jth crash. With an example i think that it could be clearer:
differential equation 1 crash:
y'(t)= [0.4^2/0.2*2]*(1-0.2y(t))*(y(t)- VALUE_OF_Y_FOR_0CRASH_AT_TIME_t]^2
differential equation 2 crashes:
y'(t)= [0.4^2/0.2*2]*(1-0.2y(t))*(y(t)- VALUE_OF_Y_FOR_1CRASH_AT_TIME_t]^2

请先登录，再进行评论。

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

 采纳的回答

Bjorn Gustavsson
Bjorn Gustavsson 2011-9-7

1 个投票

So you will have a sequence of ODEs to solve starting with the one for 0-crashes and then you want to plug in that solution into the solution for the 1-crash ODE?
If so I would call the ode-solver with
[tOut{i}, pOut{i}] = ode23(@(t, pi)myodeN(t,pi,{tOut{i-1},pOut{i-1},k,mu,r,sigma),TSPAN, IC);
and then modify your myodeN to something like this:
function dpidt = myodeN(t,pi,TnPprev,k,mu,r,sigma)
dpidt = [(sigma^2)/(2*k)]*(1-pi*k)*(pi-interp1(TnPprev{1},TnPprev{2},t))^2;
That way you interpolate a value for piNminus1 for times t from the "piNminus1" values (that you have values for at the choosen times). You could also use the single-output format of ode-solvers and then instead of explicitly interpolating use deval to calculate the values of piNminus1 for any given time.
HTH

更多回答（2 个）

ilaria
ilaria 2011-9-7

0 个投票

No, i think that this type of differential equation is not a d.d.e. I'm not sure but I think that ode23 is not able to solve this type of problem.. what do you think?

类别

在 帮助中心 和 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