how to constraint the input variables of ode45 with the odeset function (nonnegative option)

Question

yosra welhazi 2015-4-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/196804-how-to-constraint-the-input-variables-of-ode45-with-the-odeset-function-nonnegative-option

评论： Antonio Garcia-Martinez 2016-3-28

Dear friends; I have the differential equations given as follows:

    x1'=x1*(1-x2^2)-x2
    x2'=x1
I have created an m-file which contains these differential equations and I have constrained the variables x1 and x2 to their upper and lower limits (-limit and +limit)
    function xdot=fun(t,x)
limit1=2;
limit2=1;
if (abs(x(1))>limit1)
x(1)=sign(x(1))*limit1;
end
if (abs(x(2))>limit2)
  x(2)=sign(x(2))*limit2;
end
 xdot=zeros(2,1);
 xdot(1)=x(1)*(1-x(2)^2)-x(2);
 xdot(2)=x(1);
end

then I have simulated the differential equation defined in the function fun over the interval 0<=t<=20;

    x0=[0;0.25];
    [t,x]=ode45('fun',[0:0.01:20],x0);
    plot(t,x)
My problem consists in how to limit the state variable x because I have the condition
    -2<=x(1)<=2
      -1<=x(2)<=1
So, how can I simulate the differential equation over the interval 0<=t<=20 with satisfying this condition, I will be very grateful if someone can help me, because I have tried but it not works

Thanks

Can I solve this problem with the odeset function (nonnegative option), and how ??

Thanks

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Antonio Garcia-Martinez 2016-3-28

http://www.mathworks.com/help/matlab/math/parameterizing-functions.html

请先登录，再进行评论。

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

Answer 1

Torsten 2015-4-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/196804-how-to-constraint-the-input-variables-of-ode45-with-the-odeset-function-nonnegative-option#answer_174530

You can't limit the solution of an ODE.

Once you have defined the ODE and given its initial conditions, the solution over your interval of integration [0:20] is entirely determined.

So the only thing you can do is to vary the initial conditions to satisfy your constraints.

Best wishes

Torsten.

2 个评论
显示无隐藏无

yosra welhazi 2015-4-9

Can I solve this problem with the odeset function (nonnegative option)

Torsten 2015-4-9

在 MATLAB Online 中打开

The nonnegative option is a tool that helps to keep the solution of your ODE system positive if you know that the solution of the ODE system is positive.

So if e.g. you solve the ODE

y'(x)=-1, y(0)=0

with solution y(x)=-x and you activate the nonnegative option, the solver will fail.

I guess this is true for your problem - the solution exceeds the limits you want to impose by