0

I have a problem in the form of 2 equations and known initial conditions. The equations are:

dx/dt = (-a1*sin(y) + a2 + a3*sin(y-x)) / ((dy/dt)*a4*cos(y-x))

dy/dt = (a1*sin(x) -a5 + a6*x + a7*sin(y-x)) / ((dx/dt)*a8*cos(y-x))

where a1 to a8 are variables.

I am trying to plot x vs t and y vs t on MATLAB, but I am not sure how to solve this numerically or analytically. Any help would be appreciated!!

OriolAbril
  • 7,315
  • 4
  • 29
  • 40
Bob Smith
  • 1
  • 3
  • You should check the following likns, from generic to more specific https://es.mathworks.com/help/matlab/math/choose-an-ode-solver.html about ODEs, https://es.mathworks.com/help/matlab/math/solve-differential-algebraic-equations-daes.html about your specific kind of ODE, which is called DAE and https://es.mathworks.com/help/matlab/ref/ode15i.html – OriolAbril Apr 02 '18 at 22:13
  • Do you have initial conditions on x and y only or also on their derivatives? – OriolAbril Apr 03 '18 at 19:03
  • Note that you can rewrite your equations in the form `f(x,y) = cos(y-x)*x'*y' = g(x,y)`. This shows explicitly that it is a DAE. – Lutz Lehmann Apr 03 '18 at 20:59
  • Is the question solved? – OriolAbril Apr 06 '18 at 11:51

1 Answers1

1

The ODE of your problem cannot be written as dy/dt=f(t,y) nor M(t,y)dy/dt=f(t,y). That means it is a Differential Algebraic Equation which has to be solved numerically in the form:

f(t, y, dy/dt)=0

In matlab this can be done with the command ode15i.

Therefore, the first step is to write the function in a proper way, in this case, one option is:

function f = cp_ode(t,y,yp,a)
    f1 = (-a(1)*sin(y(2))+a(2)+a(3)*sin(y(2)-y(1)))/yp(2)*a(4)*cos(y(2)-y(1)) - yp(1);
    f2 = (a(1)*sin(y(1))-a(5)+a(6)*y(1)+a(7)*sin(y(2)-y(1)))/yp(1)*a(8)*cos(y(2)-y(1)) - yp(2);
    f = [f1 ; f2] ; 
end

Then, the initial conditions and integration timespan can be set in order to call ode15i:

tspan = [0 10];
y0 = [1; 1] ; 
yp0 = fsolve(@(yp)cp_ode(0,y0,yp,a),<yp0_guess>);
a = ones(1,8) ;
[t,y] = ode15i(@(t,y,yp)cp_ode(t,y,yp,a), tspan, y0, yp0);

It has to be noted that the initial conditions t, y, yp should fulfill the equation f(t, y, yp)=0. Here an fsolve is used to satisfy this condition.

The use of annonymus functions is to define the function with the parameters as inputs and then execute the ODE solver having defined the parameters in the main code.

OriolAbril
  • 7,315
  • 4
  • 29
  • 40