Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

Physics for IIT JAM, UGC - NET, CSIR NET

Created by: Akhilesh Thakur

Physics : Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

The document Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev is a part of the Physics Course Physics for IIT JAM, UGC - NET, CSIR NET.
All you need of Physics at this link: Physics

Classi cation of critical points

In the previous idea we de ned the critical points and sketched how these points can be divided to stable and unstable points. We have seen that mathematical pendulum has two critical points, one is stable, the other is not. In the next example we have seen a system with two unstable critical points. The classi cation of critical points, however, is more subtle and we discuss all possibilities in this section.
Let us first recapitulate our goal. We study planar dynamical system described by equations

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

 We assume that we have found critical point of this system, i.e. point (xC, yC ) such that

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

and study the behaviour of the system near this critical point. We linearize the equations in the neighbourhood of critical point so that we obtain equations2

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

This system can be written also in the matrix form

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

where

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

Now we discuss several forms of matrix J and classify the critical points. Finally we will show how the analysis can be done for general matrix J .


Stable and unstable nodes, saddle points

Consider linear planar system of the form

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev (8.17)

 which corresponds to matrix

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev         (8.18)

 System (8.17) can be easily solved. Equations for x and y are independent; we say that these equations are decoupled which means that equation for x_ does not contain y and vice versa.

Let us solve equation

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

first. In usual mathematical notation, this equation reads

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

which is separable dierential equation. We can rewrite it as

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

This form of equation is called separated because the left hand side of the equations contains only x and the right hand side contains only time t. We can integrate the equation,

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

to obtain

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

where C is an integration constant. It is customary that if the logarithm appears in the solution, we write the constant as a logarithm as well3:

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

Exponentiating the last equation we arrive at

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

By the same procedure we solve equation for y to get

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

where L is an integration constant again. Notice that, according to the solution, we have

x(0) = K and y(0) = L:

Hence, K and L are values of x and y at time t = 0, respectively. Therefore, we can write the solution of (8.17) in the form

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev (8.19)

Clearly, the only critical point of system (8.17) is (0; 0). Having derived solution of this system, we can analyze its behaviour near the critical point. Useful function to visualise properties of the system near critical point is StreamPlot which takes the vector eld and plots tra jectories. In the following example we choose λ= λ2 = 1.

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

In this figure we can see tra jectories (8.19) for initial points (x0, y0) chosen by Mathematica. Notice that we have inserted the right hand side of (8.17) as an argument of function StreamPlot. We can see that the trajectories are straight lines emanating from the origin (critical point) and tending to in nity exponentially.

What about other choices of λ1,2? It is clear that function eλt is increasing for λ > 0 and decreasing for λ < 0. We can conclude that qualitative behaviour of the system depends on signs of λ1,2 and four possibilities are shown in gure 8.8 which was created by following commands in Mathematica. We distinguish three cases.

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

In addition to this classi cation, critical points with distinct values λ 1 = λ 2 are called singular while critical points with the same values λ 1 = λ 2 are called degenerate.
Clearly, the saddle points cannot be singular.

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

Recall that planar dynamical system (8.17) can be represented by the matrix (8.18),

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

 From elementary linear algebra we know that with matrix J we can associate a set of eigenvalues λ defined by equation

J . e = λ e

where e is called an eigenvector. It is easy to show that the eigenvalues of matrix (8.18) are λand λ2 and corresponding eigenvectors are

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

In other words, vectors e1 and e2 satisfy equations

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

We can see that tra jectories starting on lines determined by vectors ei, i = 1, 2, always remain in these lines. If the tra jectory is being repelled from the critical point along direction e, the line determined by vector e is called unstable manifold. If the tra jectory is attracted to the critical point along the vector e, the line determined by e is called stable manifold. For matrix (8.18), vectors e1 and e2 are always eigenvectors.
We can see that elies on the x axis and e2 lies on the y axis. Hence, the axes are stable or unstable manifolds of system (8.17), depending on the sign of λ1,2.

The classi cation introduced above can be reformulated in the following way. Let

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

be a matrix of general linear dynamical system

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

If matrix J has two real eigenvalues λ1 and λ2, then critical point is stable/unstable node or a saddle point, depending on the signs of these eigenvalues.

We illustrate this classi cation on the example. Consider dynamical system

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev(8.20)

with the matrix

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

This matrix is not of the form (8.18) but we can apply the second criterion. Eigenvalues and eigenvectors can be found in Mathematica using

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

which shows that eigenvalues are

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

and corresponding eigenvectors are

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

Since λ1 > 0 and λ2 < 0, vector e1 de nes the stable manifold and ede nes unstable manifold. Since both eigenvalues have different signs, the critical point is a saddle point and it is regular. Phase tra jectories together with stable and unstable manifolds can be plotted by

Dynamical systems (Part - 3) - Classical Mechanics, UGC - NET Physics Physics Notes | EduRev

The result is plotted in gure 8.9.

Dynamic Test

Content Category

Related Searches

Sample Paper

,

Free

,

Dynamical systems (Part - 3) - Classical Mechanics

,

mock tests for examination

,

Exam

,

study material

,

video lectures

,

pdf

,

Objective type Questions

,

Viva Questions

,

UGC - NET Physics Physics Notes | EduRev

,

ppt

,

Summary

,

UGC - NET Physics Physics Notes | EduRev

,

Previous Year Questions with Solutions

,

past year papers

,

Extra Questions

,

Semester Notes

,

MCQs

,

Important questions

,

shortcuts and tricks

,

Dynamical systems (Part - 3) - Classical Mechanics

,

practice quizzes

,

UGC - NET Physics Physics Notes | EduRev

,

Dynamical systems (Part - 3) - Classical Mechanics

;