Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 1485759, 8 pages
http://dx.doi.org/10.1155/2016/1485759

Research Article

A New Method to Solve Numeric Solution of Nonlinear Dynamic System

Min Hu and Fengjun Li

School of Mathematics and Statistics, Ningxia University, Yinchuan 750021, China

Received 25 July 2016; Accepted 17 October 2016

Academic Editor: Lihua Wang

Copyright © 2016 Min Hu and Fengjun Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

It is well known that the cubic spline function has advantages of simple forms, good convergence, approximation, and second-order smoothness. A particular class of cubic spline function is constructed and an effective method to solve the numerical solution of nonlinear dynamic system is proposed based on the cubic spline function. Compared with existing methods, this method not only has high approximation precision, but also avoids the Runge phenomenon. The error analysis of several methods is given via two numeric examples, which turned out that the proposed method is a much more feasible tool applied to the engineering practice.

1. Introduction

With the advancement of science and technology, the nonlinear problems have appeared in many fields. The conventional linear approach cannot meet the requirements of solving nonlinear problems. Therefore, the nonlinear dynamic has been born, which aims at discovering complexity sciences and provides an innovative way to recognize the real and complicated systems. Bifurcation and chaos are the two typical complex dynamic behaviors of nonlinear dynamic systems. In a sense, nonlinear dynamic system is the key topic of nonlinear problems, which can study the chaotic or disordered issue and excavate the complex law of them [1, 2]. For a nonlinear dynamic system, a suitable nonlinear mathematical model should be established to reflect the inherent law of the data and then obtain the characteristics of the system. However, the nonlinear dynamic is diverse, and it changes depending on the previous state in a more complex way. Consequently, there has encountered insurmountable obstruction for applying in the practical engineering. With the existence of complicated chaotic state, it is very difficult to obtain analytic solution in general case. Naturally, people give up solving the exact solution and concentrate on studying the method with characteristics of high approximation precision and easy operation in order to describe the unknown system state.

After years of accumulation and development, there are too many methods to solve the numerical solution of nonlinear dynamic systems; the main methods are as follows: perturbation method [3], averaging method [4], Runge-Kutta method [5], Euler method [6], gradient method [7] and so on. Regretfully, these methods have certain advantages in solving certain system but obtain unappealing outcomes when solving problems of general nonlinear dynamic systems, like the lower precision, the complicity and large calculation quantity, Runge phenomenon, and so forth. So now comes the question, can we find the effective method with high approximation precision as well as avoiding the Runge phenomenon to study the nonlinear dynamic systems? It is well known that the cubic spline function has advantages of good convergence, approximation, stability, and second-order smoothness. And not only that, the cubic spline function does not exist Runge phenomenon due to the restrictions of interpolation conditions.

Given this, the purpose of this paper is to construct a new method to solve the problem mentioned above which replaces the integrand used in the existing way by the constructed cubic spline function. The paper is organized as follows. In Section 2, we recall the basic concepts from approximation theory, such as modulus of continuity. In addition, some basic knowledge of dynamic system and definition of cubic spline function are introduced. In the following, the corresponding boundary conditions of general cubic spline function are given and a particular class of cubic spline function is constructed. In Section 4, we prove our conclusions. Two numerical results are provided in Section 5 by using the previously obtained theoretical results, and the result indicates that the numerical approach based on the cubic spline function has faster convergence and higher approximation precision than existing methods. In addition, the cubic spline function can approximate the analytic solution of nonlinear dynamic system very well. In Section 6, we briefly summarize our conclusions and foresee problems for the further study.

2. Preliminaries

There is a fixed rule in the dynamic system to describe how the point changes over time in the space, and a continuous dynamic system is often represented as a set of differential equations [8]:

where , a time variable function, is used to describe the state of system ; at least, it is a continuously differentiable function defined on Euclid space or a subspace

The linear combination of linearly independent solutions can be used to express the general solution of the linear dynamic system. But it is unable to work out general solution of the nonlinear dynamic system. And for all we know, the initial value problem of nonlinear dynamic system is as follows [9]:

There is at least one solution from on the interval of . However, it is impossible to give a general calculation rule like the linear dynamic system. Consequently, the numerical method is commonly used for studying the solution of nonlinear dynamic system.

This so-called numerical method is actually a dispersed method. We can obtain the approximate solutions of unknown function values on a series of discrete points . The discrete points of independent variable can be fixed before and also can select different step length along with the different nodes. , , are generally equidistant nodes, namely, , where step length , , are commonly referred to as numerical solution of initial value problem.

Considering the nonlinear dynamic system with initial value as form (2), we obtain

Definition 1. For every , , any function , the definition of th order modulus of smoothness of is as follows [10]:

where and . When , we have . That is, the first-order modulus of smoothness of is the same as modulus of continuity of .
The smooth interpolating curves of spline function are unlikely to reveal the large oscillations feature of high-dimension polynomials. It has been widely used in cartography, pyramidal, and numerical solution methods. For example, spline function may be used for solutions of initial value problems in ordinary differential equations [11, 12].
Let be a partition of the interval .

Definition 2. A cubic interpolation function defined on is a real function with the characteristics [13]:

(a); namely, is twice continuously differentiable on .(b) coincides with a polynomial of degree three on every subinterval ,

3. Construct the Cubic Spline Function

Cubic spline interpolation function has advantages of good stability, convergence, and high approximation accuracy, which has second-order smoothness at the interpolation nodes and avoids the Runge phenomenon on account of the limitation of the interpolation condition. So we hope to replace the integrand in (3) by the constructed cubic spline function . Meanwhile satisfies the interpolation conditions, , , , and are isometric interpolation nodes. The corresponding boundary conditions are as follows [14]:

(1)The first boundary condition is as follows: , (2)The second boundary condition is as follows: , , or more special condition

The spline function is expressed by second derivative value of cubic spline function , . The polynomial is less than three order on the interval and it is a linear polynomial in particular; more specifically,

The integration constant is acquired by integrating the above equation twice and utilizing interpolation condition. And then the cubic spline function is constructed as follows:

where are unknown. Thankfully, are obtained by adopting the first and second boundary condition. As the expression of is different on the interval and , we also ensure its smoothness at nodes , namely, the smooth condition on the interpolation nodes: The tridiagonal equations about can be obtained from the above conditions. The equations have a unique solution due to the fact that coefficient matrix of the equations is diagonally dominant matrix. We can obtain based on pursing method [15] and then put into (6); the cubic spline interpolation function is acquired at last.

4. Theoretical Results

Theorem 3. If the integrand in (3) is replaced by the constructed cubic spline function (6), then the relationship of states at adjacent moment can be approximately represented as

Proof. Now, the integrand in (3) is approximated by , are interpolation nodes, and ; then

According to interpolation condition , then

In the approximate equality equation (9), is replaced by , is represented by , and then

Theorem 4. If is a kind of cubic spline function as form (6), then the nonlinear dynamic system can be approximated by with the error:

where , , , and step length

Proof. According to (6) and Peano theorem [16], we can know that if , , then is given by

Due to the fact that integral form is very complicated and does not have specific expression, can be estimated. When , , and , hence it may be

where and
And finally,

Next, the second-order smoothness property of spline function [17] is given below.

Theorem 5. The cubic spline function defined on interval meets interpolation condition , . Then the function satisfies the above condition: ,
We have that

Meanwhile, if and only if .

Proof. We prove Theorem 5 by proving

Obviously,

We use integral subsection integration such that

where are constants on every interval and at interval endpoints ,
Consequently,

Namely, .
Let in Theorem 5; we can obtain

The next section will cover two numerical examples, and the result reflects the advantages of cubic spline function in solving the nonlinear dynamic system.

5. Numerical Results

Example 1. Given a nonlinear dynamic system with an initial value (20), then we solve the numeric solution by using a variety of methods and compare them.

The numeric curve of cubic spline interpolation function is shown in Figure 1, and the approximate solution and exact solution of nonlinear dynamic system based on cubic spline interpolation function are given in Figure 2.

Figure 1: The curve of cubic spline interpolation function.

Figure 2: The approximate solution and precise curve.

In the next moment, we solve numeric solution by adopting fourth-order Runge-Kutta method, trapezoidal algorithm, and Euler method and then comparing with cubic spline interpolation function. The mean square error and the curve of numeric method are given at last; specifically see Table 1 and Figure 3.

Table 1: Mean square errors of several numeric methods for first-order nonlinear dynamic system.

Numeric method
Mean square error

Fourth-order Runge-Kutta method

Trapezoid method

Euler method

Cubic spline function

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3: The different approximation level of several numeric methods. Note: (b), (c), (d), (e), and (f) are obtained from (a), which has been magnified in different proportions.

The curves of different approximation level by several numeric methods are shown in Figure 3.

Furthermore, in order to observe approximation accuracy of nonlinear dynamic system by different numeric methods intuitively, Figure 3(a) has been magnified in different proportions, and then we can obtain Figures 3(b), 3(c), 3(d), 3(e), and 3(f). We can find out Euler method has low approximate accuracy and is great in error. The precision of the fourth-order Runge-Kutta method is also high and the calculation speed is better than trapezoidal algorithm. Although we can intuitively observe that fourth-order Runge-Kutta method has good precision compared to other methods, which is inferior to the proposed method in this paper from Figure 3.

Example 2. A second-order nonlinear dynamic system

is considered, and we solve the numeric solution of the system. Next, several different numeric methods will be used to solve this system; mean square errors and the curves of numeric approximation will be given at last, specifically shown in Table 2 and Figure 4.

Table 2: Mean square errors of several numeric methods for the second-order nonlinear dynamic system.

Numeric method
Mean square error

Fourth-order Runge-Kutta method
0.1662

Euler method
0.5479

Adams algorithm
0.3792

Cubic spline interpolation function
0.0179

(a)

(b)

(c)

(d)

Figure 4: The numeric approximation curves by different methods. Note: (b), (c), and (d) are obtained from (a), which has been magnified to some extent.

It can be seen from Figure 4 that the approximation ratio of second-order nonlinear dynamic system based on cubic spline interpolation function is better than other numeric methods, which avoids the Runge phenomenon on account of the limitation of the interpolation condition. It can approach the analytic solution of second-order nonlinear dynamic system with a high degree and is simple in calculation.

6. Conclusions and Prospects

We have discussed the numerical method based on cubic spline interpolation function to solve the numeric solution of nonlinear dynamic system in this paper. The basic theories and knowledge are introduced primarily. In the next moment, the cubic spline function is constructed according to the boundary conditions and some theory results are also given. Finally, we also demonstrate two numerical examples to reveal the effectiveness of the method proposed in this paper. The results indicate that the proposed method based on cubic spline interpolation function is obviously advantageous compared with other methods, which has quick calculation speed and avoids the Runge phenomenon.

We finish this paper with the following prospects:

(a)We only give the upper approximation of the constructed cubic spline function in this paper. Therefore, we hope to study the lower approximation in the further study.(b)The main theories only applied to first- and second-order systems, but there are many of the more complex systems in actual engineering. Consequently, it is interesting and significant to extend the main theories in this paper to higher order nonlinear dynamic system.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 11261042 and 61662060).

References

G. Chen and J. L. Moiola, “An overview of bifurcation, chaos and nonlinear dynamics in control systems,” Journal of the Franklin Institute, vol. 331, no. 6, pp. 819–858, 1994.

A. Farshidianfar and A. Saghafi, “Global bifurcation and chaos analysis in nonlinear vibration of spur gear systems,” Nonlinear Dynamics, vol. 75, no. 4, pp. 783–806, 2014.

L. V. Stepanova and S. A. Igonin, “Perturbation method for solving the nonlinear eigenvalue problem arising from fatigue crack growth problem in a damaged medium,” Applied Mathematical Modelling, vol. 38, no. 14, pp. 3436–3455, 2014.

J. Giné, J. Llibre, K. Wu, and X. Zhang, “Averaging methods of arbitrary order periodic solutions and integrability,” Journal of Differential Equations, vol. 260, no. 5, pp. 4130–4156, 2016.

H. de la Cruz, R. J. Biscay, J. C. Jimenez, and F. Carbonell, “Local linearization—Runge–Kutta methods: a class of A-stable explicit integrators for dynamical systems,” Mathematical and Computer Modelling, vol. 57, no. 3-4, pp. 720–740, 2013.

C. M. Duque and M. P. Almeida, “The Euler-Galerkin finite element method for a non-local coupled system of reaction diffusion type,” Applied Mathematics and Computation, vol. 296, pp. 116–126, 2016.

Y. H. Dai and Y. Yuan, “A nonlinear conjugate gradient method with a strong global convergence property,” SIAM Journal on Optimization, vol. 10, no. 1, pp. 177–182, 1999.

D. A. Mâceş and M. A. Stadtherr, “Computing fuzzy trajectories for nonlinear dynamic systems,” Computers and Chemical Engineering, vol. 52, pp. 10–25, 2013.

J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields, Springer, New York, NY, USA, 2013.

O. Blasco, “Modulus of continuity with respect to semigroups of analytic functions and applications,” Journal of Mathematical Analysis and Applications, vol. 435, no. 1, pp. 616–626, 2016.

F. R. Loscalzo and T. D. Talbot, “Spline function approximations for solutions of ordinary differential equations,” SIAM Journal on Numerical Analysis, vol. 4, no. 3, pp. 433–445, 1967.

A. Bellour, D. Sbibih, and A. Zidna, “Two cubic spline methods for solving Fredholm integral equations,” Applied Mathematics and Computation, vol. 276, no. 5, pp. 1–11, 2016.

B. Sepehrian and M. K. Radpoor, “Numerical solution of non-linear Fokker-Planck equation using finite differences method and the cubic spline functions,” Applied Mathematics and Computation, vol. 262, pp. 187–190, 2015.

V. Baramidze, “Smooth bivariate shape-preserving cubic spline approximation,” Computer Aided Geometric Design, vol. 44, pp. 36–55, 2016.

D. Wang, G. G. Wang, and G. F. Naterer, “Collaboration pursuing method for multidisciplinary design optimization problems,” AIAA Journal, vol. 45, no. 5, pp. 1091–1103, 2007.

P. Hájek and M. Johanis, “On Peano's theorem in Banach spaces,” Journal of Differential Equations, vol. 249, no. 12, pp. 3342–3351, 2010.

Z. Chungang and W. Renhong, “Functional splines with differential degree of smoothness and their applications,” Computer Aided Design, vol. 40, no. 5, pp. 616–624, 2008.