Future scientific and technological evolution in many areas of applied mathematics and modern physics will necessarily depend on dealing with complex systems. Such systems are complex in both their composition and behavior, namely, dealing with complex dynamical systems using different types of Duffing equations, such as real Duffing equations and complex Duffing equations. In this paper, we derive an analytical solution to a complex Duffing equation. We extend the Krýlov–Bogoliúbov–Mitropólsky method for solving a coupled system of nonlinear oscillators and apply it to solve a generalized form of a complex Duffing equation.

1. Introduction

Numerous scholars have effectively used the theory of linear oscillations to analyze and model oscillatory devices. However, nonlinear behavior can be found in a wide range of real applications. Thus, scholars from several fields of science explore nonlinear systems and try to model and investigate these complicated systems in order to find solutions and explanations to some mysterious problems, whether in the manufacture of small and large machines or electronic chips. Consequently, nonlinear oscillation is one of the most popular and widely researched fields due to its diverse applications in automobiles, sensing, microscale and nanoscale, fluid and solid interaction, nonlinear oscillations in plasma physics, bioengineering, and nonlinear oscillations in optics. There are many different and various equations of motion that are used to model several nonlinear oscillations in different physical and engineering systems. The Duffing-type equation is one of the most famous and important equations that succeeded in explaining many different oscillations in different engineering, physical systems, and statistical mechanics.

The Duffing equation is a nonlinear second-order differential equation that describes an oscillator with complex, sometimes chaotic behavior. The Duffing equation was originally the result of Georg Duffing’s systematic study of nonlinear oscillations. The behavior of the solution of the Duffing equation easily changes depending on the initial value and the polynomial coefficients, and it is difficult to predict its solution. To clarify the behavior of the solution, research based mainly on numerical analysis with high-precision calculations is conducted. Interest in the equation was later revived with the advent of chaos theory. Since then, the system has come to be regarded as one of the prototype systems in chaos theory, and related equations continue to find applications today, e.g., to describe the rolling of ships. The Duffing equation readswhere is the displacement at time and the term represents a sinusoidal driving force. The cubic term describes an asymmetry in the restoring force of a spring that softens or stiffens as it is stretched. One of the most remarkable results of dynamical systems theory is the ubiquitousness of chaotic behavior in nonlinear systems. Deterministic chaos has been observed both in mathematical models and in real physical systems. Although, from the point of view of the applications, chaotic behavior can have positive effects, improving, for example, mixing processes in chemical reactions, in other situations, such behavior can have harmful consequences, as is the case in different fields of engineering: aerodynamics, electronic circuits, and magnetic confinement of plasmas.

Some recent works on complex chaos have focused on solving complex nonlinear differential equations, complex chaos control and synchronization, and so on. For example, Cveticanin developed an approximate analytic approach for solving strong nonlinear differential equations of the Duffing-type with a complex-valued function. Furthermore, excellent agreement is found between the analytic and numerical results.

In [1], authors considered the following complex Duffing equation for modelling complex signal detection:where is a complex function, , , and are the real parameters, and the dots are the time derivatives. Its dynamical behavior was analyzed. Based on the proposed (2), they constructed a complex chaotic oscillator detection system to detect complex signals in noise. They investigated the influence of noise on the detection system and the detection performance of the system for complex signals.

In their work [2], the authors considered a complex Duffing system subjected to nonstationary random excitation of the formwhere is a complex function, , , represent the natural frequency and damping coefficient, respectively, is the small perturbation parameter and nonlinearity strength, and is a random function. This equation with has connection to the complex nonlinear Schrodinger equation which appears in many important fields of physics. Authors in [2] investigated the mean square response of a complex Duffing system subjected to nonstationary random excitation using the Wiener–Hermite expansion method combining the perturbation technique.

In 2001, Mahmoud et al. [3] presented the following complex Duffing equation:

Based on the work in [3], Li et al. [4] studied the problem of chaos control for a complex Duffing oscillation system. In general, few works are devoted to the complex Duffing equation.

In this paper, we will consider the following complex Duffing equation:

To our best knowledge, no work has been devoted to seeking analytical solutions to the complex Duffing equation. This is precisely the main objective of the present paper.

2. Undamped and Unforced Complex Duffing Equation

Let us consider the i.v.p.



Assume that and obey some Duffing equations:


Equating to zero, the coefficients of and in (10) give


On the other hand, the exact solution to the i.v.p.is expressed aswhere

3. Solution to the General Complex Duffing Equation by Means of the Krýlov–Bogoliúbov–Mitropólsky Method

Let us consider the i.v.p.

Here, , , and are the real numbers. and are the real-valued functions, and . The system (16) may be written in the form

The initial conditions are

Let us consider the following -problem:

The solution is assumed to be in the ansatz form

We choose the solutions in order to avoid the presence of the so-called secularity terms. Solving the odes gives

The approximate analytical solution is obtained by letting . It reads

The expressions for

The constants , , , and are obtained from the initial conditions.

The obtained solution is valid for . Let for the sake of simplicity; we will consider only the case when . Let us change to . We are given that

In the case when and , direct calculations show that the following function will be the exact solution to :

The constants , , , and are determined from the initial conditions:

Let us solve the general case. Assume the solution in the ansatz form:


We may solve the above system using the KBM method. To this end, we consider the following -problem:

Proceeding in the same way as we did in the first part, we obtain the following first-order approximation:


The constants , , , and are determined from the initial conditions:

4. Applications

Let us check the accuracy of the obtained results in concrete examples.

Example 1. LetSee Figures 13.

Example 2. LetThe approximate analytical solution readswhereSee Figures 46.

5. Conclusions

The nonlinear complex Duffing oscillators and many related oscillators, including the unforced undamped complex Duffing oscillator (CDO), the unforced damped CDO, and the forced damped CDO, have been analyzed using the ansatz method in order to find some approximations. For the unforced undamped CDO, the exact solution of the standard Duffing oscillator (DO) with the ansatz method was used for deriving an analytical approximation in terms of the Jacobi elliptic function. Also, the unforced damped CDO has been analyzed using the ansatz method, and with the help of the approximation of the unforced damped DO, an approximation in the form of a trigonometric form was obtained. Moreover, the forced damped CDO has been examined via the Krýlov–Bogoliúbov–Mitropólsky method (KBM), and a new analytical approximation in the form of a trigonometric formula has been derived. We demonstrated the way we may use the KBM in order to solve coupled systems of nonlinear oscillators. Other works related to nonlinear oscillators may be found in [513].

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This study was financially supported by the Universidad Nacional de Colombia.