王世兴 2013301020050
Level 1
Construct Poincare sections for various cases and compare them with Figure 3.9 in Reference [1]
Level 2
Investigate how a strange attractor is altered by small changes in one of these pendulum parameters.
Investigate the bifurcation diagrams found for the pendulum with other values of the drive frequency and damping parameter.
Chaos is an unpredictable phenomenon in a classical system domineered by the determinnistic classical mechanics. Poincare section and the bifurcation diagram are two ways to present and investigate chaos. The program calculating chaos and showing the plot of angular displacement versus time, the plot of Poincare sections and bifurcation diagram are presented. Also various parameters of the driven amplitude and frequency are tested to investigate the parameter sensitivity of the Poincare section and bifurcation diagram.
In this report we only discuss a model where the pendulum is domineered by a sinusoidal external driving force and a frinction force proportional to the angular speed, while the restoring force is also a sinuous function of angular displacement.[1]
As usual we introduce the first derivative of angular displacement and rewrite it as two first-order differential equations
$$
\begin{array}{ll}
\frac{d\omega}{dt}&=-\frac{g}{l}\text{sin}\theta-q\omega+F_D\text{sin}(\Omega_Dt)\\frac{d\theta}{dt}&=\omega
\end{array}
$$
And as we will see in the "Result" section, under some
By work of Liouville, it is well known that most differential equations could not be solved with elementary integal. Then the question of interest becomes whether we can judge the properties of the solution by the equations themselves. The French mathematicist Poincare came up with the qualitative theory and the Russian mathematicist Liapunov established the stability theory separately and contemporarily. [2]
The stability of the solution to an equation is defined as: for equations
and
Figuratively speaking, the stability means when the initial conditions deviate a little amount, the amount of the variance of the solution is also small.
Liapunov also gived methods to determine whether an equation is stable. The commonly-discussed Liapunov's second method uses a so-called Liapunov funciton
Simultaneous differential equations are called a plain autonumous system
$$\begin{matrix}\dot{x}=P(x,y)\\dot{y}=Q(x,y)\end{matrix}$$
such that the funcitons
We define the xOy plain (the
Once we get the phase trajectories, some interesting and exciting tricks could be used to find more on the properties of chaos.
If we only show the points at the same phase of the driving force, the result is called the Poincare section. Obviously if there is no chaos the Poincare section would be a single point in the phase space. However, for chaotic system the Poincare section would show a fractal behavior.
To determine the critical point where the chaotic behavior appears and changes, we observe the angular displacement on the same value of the driving force, which is called bifuraction diagram.
The source code is open on Github. Details of the codes will be discussed in the following section.
Figure 9.1 The
The definition of Poincare sections and methods to obtain it is introduced in the Background
section, and we choose four phases at
Figure 9.2 Poincare sections at four different phases.
Figure 9.3 Superposition of the four Poincare section plots.
And an interesting property that the Poincare section possesses is the fractal structure. We can zoom up a certain region of the Poincare section with phase zero and see (Figure 9.4) that there are countably infinite points in every region and they share a property of self-similarity (although not exactly the same). This is a typical fractal structure.
Figure 9.4 Fractal structure of Poincare section.
The bifurcation diagram is obtained according to the guidance from Reference [1], discussed in detail in the following
Figure 9.5 Bifurcation Diagram with
And we tried several different choices of the frequency of the driving force (in Figure 9.5) and the friction coefficient (in Figure 9.6)
Figure 9.5 Bifurcation diagram with
Figure 9.6 Bifurcation diagram with
We can see that chaos happen at a relatively large range of the choice of parameters. However, we can not find an analytical expression for the relationship between the number of points of
-
Subtle time step length.
In order to meet the requirement that we can observe the angular dislplacement in phase with the driving force we adjust the step length of the root-finding program to be an integal division of$\pi$ -
Reshape of the data with mode
$2\pi$ -
Scatter Diagram.
-
Giodano, N.J., Nakanishi, H. Computational Physics. Tsinghua University Press, December 2007
-
Differential Equation Group of Northeast University. Ordinary Differential Equation, Second Edition. Higher Education Press, April 2005.
-
Meerschaert, 《数学建模与分析方法(第一版)》,机械工业出版社,2009年5月1日