How can I plot a 3D phase space for a system of differential equations?

Click For Summary

Discussion Overview

The discussion revolves around plotting a three-dimensional phase space for a system of differential equations using Mathematica. Participants explore methods for visualizing the phase space, including the representation of vector fields and critical points, as well as the implications for stability analysis.

Discussion Character

  • Exploratory
  • Technical explanation
  • Conceptual clarification
  • Debate/contested
  • Mathematical reasoning

Main Points Raised

  • One participant inquires about the necessity of solving the system numerically before plotting the phase space.
  • Another participant suggests that the phase space is essentially a subset of ##\mathbb{R}^3## and proposes plotting the vector field or typical solutions, indicating that both approaches can be combined.
  • A participant expresses a desire to plot the phase space and superimpose critical points to study their stability, but reports that their initial plot did not yield useful information.
  • There is a discussion about the meaning of "typical" solutions, with one participant clarifying that it refers to plotting orbits through various initial conditions.
  • Participants mention the calculation of eigenvalues at critical points and the potential for confirming stability results through phase space plots.
  • Alternative software options, such as MatCont, are suggested for further analysis, although familiarity with Mathematica is noted as a limitation by some participants.

Areas of Agreement / Disagreement

Participants generally agree on the methods for visualizing the phase space and the importance of critical points for stability analysis. However, there is no consensus on the best approach to achieve useful plots, and some uncertainty remains regarding the definition of "typical" solutions.

Contextual Notes

Participants express limitations in their familiarity with specific software tools, which may affect their ability to implement suggested methods. There is also an acknowledgment of the ill-defined nature of "typical" solutions in the context of the discussion.

Who May Find This Useful

Readers interested in dynamical systems, phase space analysis, and numerical methods for solving differential equations may find this discussion relevant.

Aatifa
Hi, i would like to know how can i plot a three dimentionnal phase space (mathematica), for this kind of differential equations:
x'= (z^2(x-4y+z)+y^2(x+z)+z(x*y-x^2))/(z^2+x^2-2)
y'=y(y-3)+z(4y-z)+3(1-x^2)-2(y-3z)(z^2(x-4y+z)+y^2(x+z)+z(x*y-x^2))/x(z^2+x^2-2)
z'=2y-z(z^2(x-4y+z)+y^2(x+z)+z(x*y-x^2))/x(z^2+x^2-2)
do i need to solve numerically this system and then plot its solution?
 
Physics news on Phys.org
The phase space itself would be ##\mathbb{R}^3## (or: a subset of ##\mathbb{R}^3##), so for that there is not much to plot. (Just three axes.)

You could plot the vector field defined by the right-hand side (I think Mathematica can do that?) or you could plot some "typical" solutions. For the former, you do not need to solve the ODE, for the latter you do need to solve it indeed. A combination is also possible: Plot the vector field and superimpose on this a plot of some typical solutions starting from various initial conditions.
 
  • Like
Likes   Reactions: Dr. Courtney and Greg Bernhardt
Hi, thank you for your reply. Actually what i am trying to do is to plot the phase space of my dynamical system and superimpose on this plot the critical points of this system in order to study the stability of these critical points. Here is the result of the plot, unfortunatelly i didn't get any useful information from it.
I would like to know what did you mean by "typical" solutions?
 

Attachments

  • phase-space-aatifa.jpg
    phase-space-aatifa.jpg
    18.2 KB · Views: 660
Aatifa said:
Actually what i am trying to do is to plot the phase space of my dynamical system and superimpose on this plot the critical points of this system in order to study the stability of these critical points. Here is the result of the plot, unfortunatelly i didn't get any useful information from it.
Probably you already calculated the eigenvalues of the linearization at the critical points? Could you draw any conclusions regarding (local) stability from that?
Aatifa said:
I would like to know what did you mean by "typical" solutions?
I'm sorry, "typical" is indeed an ill-defined notion here. I meant to plot a couple of orbits through various initial conditions in the phase space to get an idea of what they look like. Likely this would work in Mathematica with some experimenting, but I am not familiar with that software. Alternatively, there are also more specialized (free) toolboxes that can help you with this and subsequent analysis. One of those is MatCont (software, Scholarpedia page), which is developed academically. (It requires MATLAB, though it might also run in Octave.)
 
  • Like
Likes   Reactions: Greg Bernhardt
Krylov said:
Probably you already calculated the eigenvalues of the linearization at the critical points? Could you draw any conclusions regarding (local) stability from that?

Yes, i did, actually i found some sadlle and stable points. But i was wondering if i could confirm this result by ploting the phase space.

I'm sorry, "typical" is indeed an ill-defined notion here. I meant to plot a couple of orbits through various initial conditions in the phase space to get an idea of what they look like. Likely this would work in Mathematica with some experimenting, but I am not familiar with that software. Alternatively, there are also more specialized (free) toolboxes that can help you with this and subsequent analysis. One of those is MatCont (software, Scholarpedia page), which is developed academically. (It requires MATLAB, though it might also run in Octave.)

thank you for your clarifications about my question, i will try to use one of these softwares that you have mensionned above.
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
Replies
8
Views
2K
  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 20 ·
Replies
20
Views
3K
  • · Replies 0 ·
Replies
0
Views
4K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
7
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K