MATLAB Answers

Making ODE45 stop at periodic solutions?

12 views (last 30 days)
Arturo Fernandez Lopez
Arturo Fernandez Lopez on 27 Dec 2019
Answered: Abhisek Pradhan on 3 Jan 2020
I am working with a Duffin oscillator:
% Recover Quadcopter variables
z = State(1,1);
w = State(2,1);
%%%%%%%%%%%%%%%%%%%%5
F(1,1) = w;
F(2,1) = epsilon*cos(omega_rho*t)-k0*z-k1*z^3-f0*w;
I wanted to find the exsistance and stability of periodic solutions. For that I wanted to evaluate the system with ODE45 and create an event that stops the integration when a periodic/quasi-periodic solution is found, but the basic event locator seems like it cnnot do it, since it cannot evaluate previous states of the solution. I have also read about Delay Differential Equation solvers and event locators, but someone said that it did not work well with quasi-periodic solutions, which I need. Is there a way to do this? If it is done using DDEs, how exactly?

  0 Comments

Sign in to comment.

Answers (1)

Abhisek Pradhan
Abhisek Pradhan on 3 Jan 2020
Event Functions can be used to determine an appropriate time to stop the solution but it does have few limitations.
Delay Differential solvers should work good for this scenario if the appropriate solver is selected.
Refer the following link to know more about Delay differential solvers, types of solvers and where they can be used.
Refer this for examples on solving Delay differential equations.

  0 Comments

Sign in to comment.


Translated by