2

I'm studying the following dynamical system, \begin{align} \dot{x} &= y \,\, , \\ \dot{y} &= \frac{\left(-4 x^3+33 x^2-78 x+54\right) y^2+(x-3) (2 x-6)^2}{(3-2 x) (2-x) x (2 x-6)} \,\, , \end{align}

and when I plot the phase space I get the following graph: Phase space.

I don't know what is this kind of fixed point, it looks like a half saddle - half center, but there is no fixed point inside the closed orbit, but it is on this orbit. The linear stability isn't conclusive in this case (both eigenvalues are 0).

Ps: The blue and red lines are the curves of the first integral in this system. I used this first integral as Lyapunov function of the system and I checked that exactly in this fixed point, (x,y) = (3,0), there is an inflection point and the test is also inconclusive. Any clues about what is this kind of fixed point and how can I infer its stability?

Herr Schrödinger
  • 637
  • 1
  • 8
  • 20

2 Answers2

4

We can rearrange $$\dot{y}=\frac{2(x-3)^2}{(3-2x)(2-x)x} + \frac{-4 x^3 + 33 x^2 - 78 x + 54}{(3-2x)(2-x)x} \cdot \frac{y^2}{2(x-3)}$$ as $(x,y) \rightarrow (3,0)$, $\dot{y} \rightarrow 0 + 1\cdot \frac{y^2}{2(x-3)}$.

Now near $(3,0)$ you have a system $$ \dot{x}=y \\ \dot{y}=\frac{y^2}{2(x-3)} $$

$(3,0)$ doesn't seem like a fixed point. In one direction we have $\infty$, in the other we have $0$.

You can actually solve this set of equations by assuming $x(t)=A+Bt+Ct^2$.

Hope this helps you understand the equation.

  • 1
    It might be arguably useful to look at $z\equiv (x-3)/y^2$ at that point, since $\dot z=0$, so z is an undetermined constant c. Thus, $\dot y =1/2z=1/2c$, an increasing function across that point. Not a fixed point. The flow steams through. – Cosmas Zachos May 03 '20 at 20:52
  • 1
    That's a great observation! If we say $\dot{y}=1/2 c$ and $c$ is none zero, then $y^2=(x-3)/c$ would be a parabola. And we do seem to see these parabolas on the stream plot. Although we might have missed the loops on the left during this analysis. It would be great if we could do some useful change of variables to the original ODE. –  May 04 '20 at 02:19
  • Sure, it is an invitation to the OP to pick new variables, maybe (z,y), maybe not, to make the structure of his system around that critical point more transparent. A global plot in terms of those might be informative. – Cosmas Zachos May 04 '20 at 13:00
2

with:

$$\vec{R}=\left[ \begin {array}{c} y\\ {\frac { \left( -4\,{x }^{3}+33\,{x}^{2}-78\,x+54 \right) {y}^{2}+ \left( x-3 \right) \left( 2\,x-6 \right) ^{2}}{ \left( 3-2\,x \right) \left( 2-x \right) x \left( 2\,x-6 \right) }}\end {array} \right] $$

you get the matrix A to obtain the eigenvalues of your system.

$$A= \left[ \begin {array}{cc} {\frac {\partial }{\partial x}}R_{{x}} \left( x,y \right) &{\frac {\partial }{\partial y}}R_{{x}} \left( x,y \right) \\ {\frac {\partial }{\partial x}}R_{{y}} \left( x,y \right) &{\frac {\partial }{\partial y}}R_{{y}} \left( x,y \right) \end {array} \right] $$

you can now obtain the eigenvalues of the matrix A at $y=0$

$$\lambda_1\bigg|_{y=0}=\frac{1}{2}\,{\frac {\sqrt {7056\,{x}^{3}+11664\,x-2016\,{x}^{4}-3888-16\,{x}^ {6}+288\,{x}^{5}-12960\,{x}^{2}}}{ \left( -3+2\,x \right) \left( -2+x \right) x \left( x-3 \right) }} \tag 1$$

and

$$\lambda_2\bigg|_{y=0}=-\lambda_1\bigg|_{y=0}\tag 2$$

for $x_i=\frac{3}{2}\,,2\,,0\,,3$ the denominator is zero.

if you take the limit of equation (1) and (2) where the the denominator is zero, you get just for $x_i=3$ results: equation(1) and (2) are equal zero. thus

$y=0,x=3$ are semi stable points ,which you can see on your phase diagram, all other points are not stable points

Eli
  • 438