Strange fixed point in state space

131 Views Asked by At

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?

2

There are 2 best solutions below

0
On BEST ANSWER

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

3
On

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.