Peme5220m Adv. Math. Tech. Acm Coursework-non Linear Stability

PEME5220M Adv. Math. Tech. ACM Coursework-Non linear Stability

Peishin Tan 200305149

Question: Investigate stability of given system: dx =y dt dy = −x + x 3 − y dt

Solution: dx dy and must equal to zero. Hence we have y = 0 and dt dt x ( x 2 −1) = 0 thus x = 0 or x = ±1 .

To find stationary points,

Therefore we have fixed points at (0,0) , (1,0) and (-1,0) . 

0

1 

 Jacobian of the system, J =  2  3 x − 1 − 1 2 Thus Tr ( J ) = −1 and Det ( J ) =1 − 3x 2 , λ1, 2 = − 1 ± 1 − 4(1 − 3x )

2

2

⇒ At point (0,0) 0 J ( 0, 0 ) =   −1 

1   ; Tr ( J ) = −1 ; Det ( J ) =1 ; λ1, 2 = − 1 ± −1 2 

−3 2

With conditions: Tr ( J ) < 0 , Det ( J ) > 0 , Tr ( J ) 2 − 4 Det ( J ) < 0 indicated stable focus behaviour. Eigenvector at point (0,0)

⇒ At point (1,0) 0 J (1,0 ) =  2 

1  1 3  ; Tr ( J ) = −1 ; Det ( J ) = −2 ; λ1, 2 = − ± = 1 or − 2 −1 2 2 

With conditions: Tr ( J ) < 0 , Det ( J ) < 0 , Tr ( J ) 2 − 4 Det ( J ) > 0 indicated saddle point behaviour. Eigenvector at point (1,0)

⇒ At point (-1,0) 0 J ( −1,0 ) =  2 

1  1 3  ; Tr ( J ) = −1 ; Det ( J ) = −2 ; λ1, 2 = − ± = 1 or − 2 −1 2 2 

With conditions: Tr ( J ) < 0 , Det ( J ) < 0 , Tr ( J ) 2 − 4 Det ( J ) > 0 indicated saddle

PEME5220M Adv. Math. Tech. ACM Coursework-Non linear Stability

point behaviour.

Peishin Tan 200305149