1. Introduction
Many nonlinear phenomena that appear in engineering, chemistry, physics, economics, and biology can be modeled by the nonlinear dynamical systems of the form , , where such that the system admits a Hamilton–Poisson structure (e.g., is a Hamilton–Poisson system) and is an additive term. There are two functionally independent constants of motion, (the Hamiltonian function) and (the Casimir function).
In the last decade, the dynamical properties have been examined by several researchers as bifurcation route, Poincar
map, frequency spectrum, amplitude modulation, topological horseshoe, the existence of heteroclinic orbit or homoclinic orbit, equilibria, Lyapunov exponent spectrum, a dissipative system, phase portraits, bifurcation diagrams, and Hopf bifurcation. These properties characterize the chaotic behaviors of the dynamical system. Li et al. [
1] studied a three-dimensional autonomous chaotic system that is found to possess two nonhyperbolic equilibria. Pham et al. [
2] introduced a new system with an infinite number of equilibrium points. Wang et al. [
3] presented a watermark encryption algorithm for a new memristive chaotic system. Zhang et al. [
4] proposed a numerical scheme for the study of the dislocated projective synchronization (DPS) between the fractional-order and the integer-order chaotic systems. Tong [
5] investigated the chaotic attractor for a three-dimensional (3D) chaotic system that possess invariable Lyapunov exponent spectra and controllable signal amplitude. He et al. [
6] introduced a new four-dimensional chaotic system with coexisting attractors having three quadratic nonlinearities and only one unstable fixed point. Singh et al. [
7] reported a new 4D dissipative chaotic system studying the coexistence of asymmetric hidden chaotic attractors with a curve of equilibria. Sun et al. [
8] proposed a novel kind of compound–combination antisynchronization scheme among five chaotic systems. Cicek et al. [
9] implemented in practical applications a new three-dimensional continuous time chaotic system by an electronic circuit design. Lai et al. [
10] numerically investigated a new 3D autonomous chaotic system with coexisting attractors. Varan et al. [
11] implemented a synchronization circuit model of a third-degree Malasoma system with chaotic flow. Su [
12] investigated the horseshoe chaos using the topological horseshoe theory, taking into account a three-dimensional (3D) autonomous chaotic system. Zhou et al. [
13] introduced and analyzed theoretically the basic dynamical properties of a three-dimensional chaotic system. The result shows the chaotic attractor by the realization of a circuit experiment. Akgul et al. [
14] explored a three-dimensional chaotic system with cubic nonlinearities. They applied the electronic circuit implementation for real environment application. Pham et al. [
15] introduced a three-dimensional chaotic system displaying both hidden attractors with infinite equilibria and hidden attractors without equilibrium. Zhang [
16] investigated a method for generating complex grid multiwing chaotic attractors. Kacar [
17] developed a four-dimensional chaotic system and implemented an analogue circuit and microcontroller. Tuna et al. [
18] presented numerical, analog, and digital circuit modelings by using a 3D chaotic system with a single equilibrium point. Naderi et al. [
19] explored the exponential synchronization of the chaotic system without a linear term and its application in secure communication by using the exponential stability theorem and showing the ability and effectiveness of the proposed method by numerical simulation. Li et al. [
20] studied complicated dynamical behaviors of a three-dimensional chaotic system with quadratic nonlinearities.
Recently, Liu et al. [
21] developed a new multiwing chaotic system that has an excellent effect on image encryption. Hu et al. [
22] designed a circuit implementation to verify the physical feasibility of an asymmetric memristor-based chaotic system with only one equilibrium point. Sun et al. [
23] studied a color image encryption scheme base on a 5D memristive chaotic system. Wang et al. [
24] explored the problem in image encryption on the basis of a chaotic system with time delay. Guo et al. [
25] proposed a multivortex hyperchaotic system, emphasizing its application to image encryption and outstanding anticropping and antinoise performance. Yildirim et al. [
26] used the particle swarm optimization (PSO) and ant colony optimization (ACO) to optimize the initial conditions of a continuous-time chaotic system. Ding et al. [
27] proposed a cryptosystem and its application in image encryption. Lai et al. [
28] proposed a four-dimensional multiscroll chaotic system with application to image encryption. Lu et al. [
29] proposed an encryption algorithm for 3D medical models.
Recently, Karimov et al. [
30] implemented an analog circuit and proposed a novel technique for reconstructing ordinary differential equations (ODEs) describing the circuit from data. This technique is shown for a well-studied R
ssler chaotic system. Karimov et al. [
31] studied the synchronization between a circuit modeling the R
ssler chaotic system and a computer model by using adaptive generalized synchronization.
Beyond chaotic behaviors, some systems could have nonlinear singularities. Such systems are investigated using the topological degree theory and the qualitative analysis of a Poincar
map with action angle variables [
32]. Cheng et al. [
33] established the existence of homoclinic solutions for a differential inclusion system involving the
-Laplacian by using a variational principle. Fonda et al. [
34] proved the existence and multiplicity results for periodic solutions of Hamiltonian systems using the Poincar
–Birkhoff fixed point theorem.
Many nonlinear differential problems from applied engineering are analytically solved by some methods, namely, the multiple scales technique [
35], the optimal iteration parametrization method (OIPM) [
36], the optimal homotopy asymptotic method (OHAM) [
37,
38,
39], and the optimal homotopy perturbation method (OHPM) [
40,
41,
42].
The structure of this paper is as follows: In
Section 2, we present in detail some dynamical systems involving a Hamilton–Poisson part. The steps of the mOPIM technique are the subject of
Section 3.
Section 4 presents the semianalytical solutions obtained by the mOPIM method.
Section 5 provides the numerical results and emphasizes the validation of the method. The conclusions and perspectives are highlighted in
Section 6.
3. The Basic Idea of the mOPIM Technique
Let the second-order nonlinear differential equation be
subject to the initial conditions
where
is a linear operator,
a nonlinear operator,
a boundary operator,
g a known function,
u an unknown smooth function depending on the independent variable
t, and
.
Marinca et al. [
36] proposed the following iterative scheme, namely, optimal parametric iteration method (OPIM), defined by
where
,
, and
are auxiliary continuous functions;
(obtained from Taylor series expansion of the nonlinear operator
);
is the (
n + 1)-th-order approximate solution of Equations (
56) and (
57), denoted by
; and
is the initial approximation, a solution of the linear differential problem:
The real constants , , are are unknown convergence-control parameters and can be optimally computed.
Remark 2. - (1)
In the case of nonlinear oscillators, the integration of Equation (
58)
produces secular terms of the form ,
,
,
,
,
,
and so on. The presence of has the advantage of avoiding the secular terms that appear through integration with the OPIM method, and that makes the oscillation amplitude tend toward infinity (physically, the resonance phenomenon occurs).
- (2)
The OPIM method was successfully applied in the case of ODEs with boundary conditions (see Ref) [
55],
such as - (a)
Thin film flow of a fourth-grade fluid down a vertical cylinder where
. The linear operator is chosen as .
- (b)
Thermal radiation on MHD flow over a stretching porous sheet The initial guess is chosen as
and
, with .
- (c)
The oscillator with cubic and harmonic restoring force The linear operator is chosen as .
- (d)
The Thomas–Fermi equation The linear operator is chosen as
, and the nonlinear operator yields
.
- (e)
Lotka–Volterra model with three species The initial approximations are chosen as
, , and
or
, , and
, and so on.
Next, we propose a modified version of the OPIM procedure, namely, the modified optimal parametric iteration method (mOPIM), in the following form:
where the new auxiliary continuous function
is a nonzero function and
,
, and
have the same signification. The unknown real parameters
,
,
, and
are optimally computed at least.
The
-order approximate solution of Equation (
65) is well determined if the convergence-control parameters are known.
If
is the initial approximation of Equation (
59), the nonlinear operators
,
,
, and
that appear in Equation (
65) have the form
where
is a positive integer, and
and
are known functions that depend on
.
Using the linearly independent functions
, we introduce some types of approximate solutions of Equation (
56).
Definition 1. A sequence of functions of the formis called an mOPIM sequence of Equation (56). Functions of the mOPIM sequences are called mOPIM functions of Equation (56). The mOPIM sequences with the propertyare called convergent to the solution of Equation (56), where . Definition 2. The mOPIM functions satisfying the conditionsare called ε-approximate mOPIM solutions of Equation (56). Definition 3. The mOPIM functions satisfying the conditionsare called weak ε-approximate mOPIM solutions of Equation (56) on the real interval . The existence of weak -approximate mOPIM solutions is built by the theorem presented above.
Theorem 1. Equation (56) admits a sequence of weak ε-approximate mOPIM solutions. Proof. It is similar to the theorem from [
56]. □
For
, an
-order approximate solution of Equations (
56) and (
57), the validation of this procedure is highlighted by computing the residual function given by
such that
, for all
.
4. Approximate Analytic Solutions via mOPIM
This section emphasizes the applicability of the mOPIM procedure for the nonlinear differential problems given by Equations (
8) and (
9) using only one iteration. This problem could be written in the form of Equation (
56), taking the following operators (
):
Taking into consideration the linear operator given by Equation (
71), the initial approximation
, the solution of Equation (
59) is
with
,
.
Using Equation (
71), a simple computation yields the following expressions:
Returning to Equation (
65), there are a lot of possibilities to choose the following auxiliary functions:
or
and so on.
Taking into account Equation (
74), for
, a simple computation yields
where
By the integration of Equation (
65) and using the expressions given by Equations (
71)–(
75), the first-order approximate solution
could be obtained:
with the unknown real parameters
,
,
,
, and
, depending on the parameters
,
,
,
,
,
,
,
,
,
,
,
,
,
, and
, and can be optimally identified.
Analogously, for the value , the expression is a linear combination between the elementary functions 1, , , , , , , , and .
Then the first-order approximate solution
obtained from Equation (
65) is a linear combination between the elementary functions 1,
,
,
,
,
,
,
, and
.
For an arbitrary integer number
, inductively, the expression
is a linear combination between the elementary functions 1,
and
, for
. Therefore, the the first-order approximate solution
will be of the form
where the unknown convergence-control parameters
,
,
,
, and
for
could be optimally computed.
Using the same procedure, the approximate closed-form solutions of the nonlinear problems presented in
Section 2.2 could be obtained by means of the mOPIM method.
6. Conclusions
A new analytical approach, namely, the modified optimal parametric iteration method (mOPIM), for solving second-order nonlinear differential equations is developed using only one iteration.
In this way, the closed-form analytical approximate solutions are built for a class of nonlinear dynamical systems that possess a Hamilton–Poisson structure.
The obtained results are validated by graphically comparing them with the corresponding numerical solutions. The corresponding absolute errors are tabulated.
A comparison between the approximate analytical solution obtained with mOPIM, the analytical solution obtained with the iterative method, and the corresponding numerical solution highlights the advantages of the mOPIM method.
These comparisons prove the precision of the applied method in the sense that the semianalytical solutions are approaching the exact solution; e.g., the residual functions are much smaller than 1.
The achieved results have high potential, especially given the strong alignment demonstrated between the analytical and numerical outcomes, and they encourage the study of other dynamical systems with similar properties.
The possibility of a comparison between our results and some experiments based on the dynamical systems having a Hamilton–Poisson structure could be the subject of a future work.