The article searched on mathematics and numerical solutions for the nonlinear pendulum (Chaotic pendulum). The numerical solution that was used for our research suitably the pseudospectral methods. With these equations, we studied and calculated on the interval −1; 1, with boundary conditions already known. We used the software Mathematica 10.4 to calculate the results of the problems.
Trang 1УДК 517.9
Pseudospectral Methods for Nonlinear Pendulum
Equations
Le Anh Nhat∗
RUDN University
6 Miklukho-Maklaya str., Moscow, 117198
Russia Tan Trao University Tuyen Quang 22227
Vietnam
Received 29.05.2018, received in revised form 10.08.2018, accepted 20.10.2018
The article searched on mathematics and numerical solutions for the nonlinear pendulum (Chaotic pen-dulum) The numerical solution that was used for our research suitably the pseudospectral methods With these equations, we studied and calculated on the interval [−1, 1], with boundary conditions already known We used the software Mathematica 10.4 to calculate the results of the problems.
Keywords: chaotic pendulum, Chebyshev, pseudospectral methods, differentiation matrices, collocation method, nonlinear equations.
DOI: 10.17516/1997-1397-2019-12-1-79-84
The differential equation of a single pendulum is
d2
here v is the time, θ is value of the angular displacement, the number value v2= g/l, the clarify the value g is the acceleration due to gravity and l is the length of the pendulum.
The fact that the pendulum is placed in a friction environment, the drag force on the pendu-lum is calculated accordance with the Stokes law and is propertied to the instantaneous velocity
of the pendulum because of the environment has value friction the equation of motion (1) becomes
d2
dt2θ(t) = −v2sin θ(t) − αd
and here by α is the coefficient of friction of the environment and show that is α > 0.
To resist the friction, we need to add some external force, to simplified we choose the external periodic force and so the equation of motion becomes the equation:
d2
dt2θ(t) = −v2sin θ(t) − αd
here β is amplitude and w is frequency.
When α = 0 and then (3) becomes
d2
∗leanhnhat@mail.ru
c
⃝ Siberian Federal University All rights reserved
Trang 2all the above equations are nonlinear second order differential equations The nonlinear
equa-tion (3) has no analytic soluequa-tion However, if θ → 0 then sin(θ) ≈ θ, so the nonlinear differential
equations above can be transformed into linear differential equations This has got a lot of solutions, so we do not study
With these nonlinear equations, we will be study and compute by the pseudospectral method
Let p(x) is a polynomial of degree n, and we know that it is values at the points p(x0), p(x1),
, p(x n ), then the first and second derivatives p(x) at the same points are expressed in matrix
p ′ (x0)
p ′ (x n)
= D
p(x0)
p(x n)
,
p ′′ (x0)
p ′′ (x n)
= D2
p(x0)
p(x n)
where D = {d i,j } is the so-called differentiation matrix [4].
In case when the Chebyshev-Gauss-Lobatto points are chose as the collocation points y k =
= cos(kπ/n) [5],
D i,j=
2n2+ 1
c i
2c j
(−1) i+j
sin[π(i + j)/(2n)] sin[π(i − j)/(2n)] i ̸= j
cos(πj/n)
2 sin(πj/n) 0 < i = j < n
− 2n2+ 1
(6)
here c k =
{
2 if k = 0, n
1 if k = 1, 2, , n − 1.
The application of differential algebra in ordinary differential equations can also extend to
nonlinear differential equations, so we transformed the matrix D into matrices [3]:
E(1)={d i,j }, 1 6 i, j 6 n − 1,
e0={d i,0 }, e n={d i,n }, 0 < i < n (7)
for a first-order differential element, the form u ′ (x i ) = E(1)u(x i ) + be0+ ae n here a and b are
the two-point boundary-value on the range [–1; 1] of problem
For a second-order differential element, we use D2={d(2)
i,j } and define the matrices:
E(2)={d(2)
i,j }, 1 6 i, j 6 n − 1,
e2={d(2)
i,0 }, e2
n={d(2)
has the form u ′′ (x
i ) = E(2)u(x i ) + be(2)0 + ae(2)n When we applied Pseudospectral method to solve nonlinear differential equations, nonlinear equations may or may not have a unique solution
If we have identified a solution, we will proceed with an iterative procedure Therefore, it
is important to determine the iterative equation The iterative procedure is simple, we assume
u(0)= const, then find u(1), u(2), , stop it until the error ε = |u (k) − u (k −1) | < ε0
Trang 32 Results
Equation (1), we shall consider the two-point boundary-value of problem on the range [ −1, 1], θ( −1) = a, θ(1) = b We can find and transfer to form of method:
E(2)θ(t j) =−v2sin θ(t j)− be(2)
0 − ae(2)
Equation (9) are iterative equation We give the error ε6 10−12 In the Fig 1 and Tab 1
are the result are calculated based on the program by pseudospectral method, and then the solid line show the result calculated by Mathematica 10.4
Fig 1 Graphic of equation with v = 1.7, n = 64 in the case a = b = 0.9
Table 1 Competition the numerical results and error with Mathematica’s calculations of the
equation (1) with v = 1.7, n = 64 in the case a = b = 0.9
1 0.998795 0.903167 0.903167 1.53307×10 −7
10 0.881921 1.19371 1.19371 1.33523×10 −7
20 0.55557 1.80554 1.80554 7.05996×10 −8
30 0.0980171 2.17294 2.17294 4.01061×10 −8
40 –0.382683 2.00793 2.00793 2.20135×10 −8
50 –0.77301 1.43174 1.43174 1.19427×10 −8
60 –0.980785 0.950125 0.950125 1.65847×10 −9
Remasks: by this case, when 0 < v < 1.8, the results are correct and convergence When
v < 1.8, means that the length of the string is at least g/v2 When v > 1.8, equation (1) non-compliance boundary conditions θ( −1) = a, θ(1) = b.
Equation (2), we will be consider on the range [ −1, 1] and the boundary conditions θ(−1)=a, θ(1) = b We can find iterative equation and transfer to form of method:
(E(2)+ αE(1))θ = f − b(e(2)
0 + αe0)− a(e(2)
here θ and f denote the vectors with the elements {t j } and {−v2sin θ(t j)}.
In this problem, we give the error ε < 10 −8 In the Fig 2, points of value are the calculated
by result of the program by pseudospectral method, and then the solid line show the result
Trang 4Fig 2 Graphic of equation (2) with α = 10, v = 1.7, n = 96 in the case a = b = −0.3
Table 2 Competition the numerical results and error with Mathematica’s calculations of the
equation (2) with α = 10, v = 1.7, n = 96 in the case a = b = −0.3
5 0.986643 –0.301648 –0.301648 6.69452×10 −10
20 0.793353 –0.326501 –0.326501 3.08531×10 −9
35 0.412707 –0.381289 –0.381289 5.35524×10 −9
50 –0.0654031 –0.462363 –0.462363 6.22867×10 −9
65 –0.528068 –0.551825 –0.551825 6.90365×10 −9
80 –0.866025 –0.533131 –0.533131 1.19349×10 −7
95 –0.999465 –0.301757 –0.301757 7.68573×10 −7
calculated by Mathematica In the Tab 2 are the competition the numerical results and error
with Mathematica’s calculations with α = 10, v = 1.7, n = 96 in the case a = b = −0.3 Remasks: The fixed point, θ(t) = dθ/dt = 0 is linearly stable, i.e., small perturbations
from this point will decay in time The fixed point θ(t) = π, dθ/dt = 0 is linearly unstable,
which means that small perturbations from this point will grow exponentially in time [2] When
0 < v < 2.7, the results are correct and convergence; when v > 2.7 method for big error and not convergence When v fixed, α increase, then graph are left-leaning.
Equation (3), we will be consider on the range [ −1, 1] and in the case of boundary conditions θ( −1) = a, θ(1) = b Equation (3) is non-Hamiltonian and it does not have an analytical
solution This is a nonlinear equation and has three dynamic variables We found iterative equation and transferred to form of the method:
(E(2)+ αE(1))θ = f − b(e(2)
0 + αe0)− a(e(2)
here θ and f denote the vectors with the elements {t j } and {β cos wt j − v2sin θ(t j)}.
With v, α, β, w > 0, we give the error ε6 10−12 In the Fig 3, points of value are calculated
the result of the program by the PSM, and then the solid line show the result calculated by Mathematica Numerical results are given in Tab 3
Remasks: When 0 < v < 1.7 the results are correct and convergence; when β cos(wt) descended
then θ(t) ascended; when α ascended then θ(t) descended; the complex motion one would expect
when the three forces are comparable; from here we can orient the pendulum control
Trang 5Fig 3 Graphic of equation (3) with α = 0.5, v = 1.1, w = 2π, β = 1.3, n = 128 in the case
a = b = 0
Table 3 Competition the numerical results and error with Mathematica’s calculations of the
equation (3) with α = 0.5, v = 1.1, w = 2π, β = 1.3, n = 128 in the case a = b = 0
1 0.999699 0.0000121827 0.0000121743 8.4617×10 −10
25 0.817585 0.027557 0.027557 1.50768×10 −8
50 0.33689 0.0854559 0.0854558 2.98413×10 −8
75 –0.266713 0.0722918 0.0722917 3.26543×10 −8
100 –0.77301 0.0487747 0.0487747 4.7381×10 −9
125 –0.99729 0.000283445 0.000283453 7.51629×10 −9
Equation (4), we will be consider on the range [ −1, 1] and in the case of boundary conditions θ( −1) = a, θ(1) = b We found iterative equation and transferred to form of method:
E(2)θ = f − be(2)
0 − ae(2)
here θ and f denote the vectors with the elements {t j } and {β cos wt j − v2sin θ(t j)} with
v, β, w > 0, we give the error ε 6 10−12 In the Fig 4, the value are the calculated base on
the program by the PSM, and then the solid line show the result calculated by Mathematica The calculation results are shown in Tab 4
Remasks: When 0 < v < 1.9 the results is show correct and convergence; when β cos(wt)
ascended then θ(t) ascended; when w ascended then θ(t) descended.
Fig 4 Graphic of equation (4) with v = 1.0, β = 20, w = 20, n = 128 in the case a = b = 0
Trang 6Table 4 Competition the numerical results and error with Mathematica’s calculations of the
equation (4) with v = 1.0, β = 20, w = 20, n = 128 in the case a = b = 0
1 0.999699 –0.000265697 –0.000265549 1.47924×10 −7
30 0.740951 0.0595138 0.059514 1.47715×10 −7
60 0.0980171 0.0566888 0.0566889 1.26391×10 −7
90 –0.595699 –0.0085085 –0.00850845 5.6271×10 −8
127 –0.999699 –0.000265697 –0.000265696 9.03764×10 −10
We have monitored four equations of the nonlinear pendulum with the range [−1, 1] and the
know boundary conditions by the pseudospectral methods The method has issued good results And then we can survey the variation of the function through our parameters and design of force acting on the pendulum to control the operation the pendulum
The publication was prepared with the support of the “RUDN University Program 5-100”.
References
[1] A.Belendez, E.Arribas, Approximate solutions for the nonlinear pendulum equation using a
rational harmonic representation, Article in Computers & Mathematics with Applications,
9(2012), no 64, 1602–1611
[2] G.Moshe, The Chaotic Pendulum, World Scientific Publishing Co Pte Ltd, 2010
[3] G.L.Baker, J.A.Blackburn, The Pendulum — a case study in physics, New York, Oxford University Press Inc., 2005
[4] J.C.Mason, D.C.Handscomb, Chebyshev Polynomials, CRC Press LLC, 2003
[5] W.S.Don, A.Solomonoff, Accuracy and speed in computing the Chebyshev collocation
de-vivative, SIAM Juarnal of Scientific Computing, 16(1991), no 6, 1253–1268.
Псевдоспектральные методы нелинейных уравнений
маятника
Лe Ань Ньат
Российский университет дружбы народов Миклухо-Маклая, 6, Москва, 117198
Россия
В статье исследованы математические и численные решения для нелинейного маятника (хаоти-ческий маятник) Решение, которое использовалось для наших исследований, соответствовало псевдоспектральным методам Вычисления проводили на интервале [–1, 1] с уже известными граничными условиями Для расчета использовалось программное обеспечение Mathematica 10.4 Ключевые слова: хаотический маятник, Чебышев, псевдоспектральные методы, дифференцирую-щие матрицы, метод коллокации, нелинейные уравнения.