Variational Iteration Method for Solving Differential Equations with Piecewise Constant Arguments
Автор: Qi Wang, Fenglian Fu
Журнал: International Journal of Engineering and Manufacturing(IJEM) @ijem
Статья в выпуске: 2 vol.2, 2012 года.
Бесплатный доступ
In this paper, variational iteration method is applied for finding the solution of differential equations with piecewise constant arguments. A correction functional is constructed by a general Lagrange multiplier, which can be identified by variational theory. This technique provides a sequence of functions which converges to the exact solution of the problem without discretization of the variables. The flexibility and adaptation provided by the method have been verified by an example.
Variational iteration method, Piecewise constant arguments, Approximate analytical solution
Короткий адрес: https://sciup.org/15014296
IDR: 15014296
Текст научной статьи Variational Iteration Method for Solving Differential Equations with Piecewise Constant Arguments
Recently, there has been much research concerning properties of solutions of differential equations with piecewise constant arguments (EPCA) [1-4]. In these equations, the derivatives of the unknown functions depend on not just the time t before t . Since this is essentially expressing the derivatives on terms of the solution at discrete points of time before t , it is usually referred to as a hybrid system, other examples of the application of these equations to the problems of biology, cellular neural networks and mechanical systems can be found in [5-8] .
EPCA has been under intensive investigation for the last twenty years. The first studies in this field have been given in [9,10], after this, stability, contractivity and existence of periodic solutions have been treated by several authors, see [11-13] and references therein. The general theory and basic results for EPCA have been thoroughly investigated in the book of Wiener[14]. Nowadays, much research has been focused on the
* Corresponding author.
numerical solutions of EPCA, including convergence[15], stability[16] and oscillation[17]. However, we can’t find any results concerning the approximate anaytical solution of EPCA have been published. In our paper, we will apply the well-known anaytical approximate technique: the variational iteration method to the follwing EPCA:
u ’(t) = au(t) + bu([“]), t - 0, u (0) = u0, where a, b, u^ G R and [•] denotes the greatest integer function. The main goal of this paper is to extend the variational iteration method to find the solution of (1).
The variational iteration method (VIM) [18-21], which is a well-established technique with wide applications for ordinary differential equations, partial differential equations and delay differential equations, etc. More materials of the classical solution techniques that are most commonly used to solve equations and systems, such as the Adomian decomposition method (ADM) are found in [22], Differential transformation method (DTM) [23], Homotopy analysis method (HAM) [24] and the Homotopy perturbation method (HPM) [25]. The VIM is the most e ff ective and convenient one for both weakly and strongly nonlinear equations. This method has been shown to e ff ectively, easily, and accurately solve a large class of linear and nonlinear problems.
-
2. PrelimInaries
Definition 1 [14] A solution of (1) on [0, ^ ) is a function u ( t ) that satisfies the conditions:
-
(i) u ( t ) is continuous on [0, ^ ) ,
-
(ii) The derivative u '( t ) exists at each point t G [0, to ) , with the possible exception of the point t = 2 n , n = 0,1, • * • , where one-sided derivatives exist,
-
(iii) Eq. (1) is satisfied on each interval [2 n , 2( n + 1)) c [0, to ) .
-
3. VARIATIONAL ITERATION MRTHOD (VIM)
Theorem 1 [14]. If a < 0,| b | < — a , then the solution of (1) tends to zero as t ^ to for any given u0
The VIM is the general Lagrange method, in which an extremely accurate approximation at some special point can be obtained, but not an analytical solution. To illustrate the basic idea of the VIM we consider the following general differential equation
Lu + Nu = g ( x ) (2)
where L and N are linear and nonlinear operator, respectively, and g ( x ) is the inhomogeneous term. In [26–29], the author proposed the VIM where a correction functional for (2) can be written as
x
-
un + 1 ( x ) = u n ( x ) + I X(t )[ Lu n (t ) + Nu n ( t ) — g(t)]dt (3)
where X is a general Lagrange’s multiplier, which can be identified optimally via the variational theory, and u n as a restricted variation which means ^u n = 0 . It is to be noted that the Lagrange multiplier X can be a constant or a function.
The VIM should be employed by following two main steps. It is required first to determine the Lagrange multiplier X that can be identified optimally via integration by parts and by using a restricted variation. Having X determined, an iteration formula, without restricted variation, should be used for the determination of the successive approximations u(x) of the solution u(x). The zeroth approximation u can be any selective function. Consequently, the solution is given by u (x) = lim un (x) .
n ^да
Now, we apply the VIM to (1) . A correction functional is constructed as follows: ts u +i(t) = un (t) + X(s)[u (s) - aun (s) - bun ([—])]ds.
n + 1 n 00 2
Taking variational on both sides of (5), we have t'
5u,(t) = Sun (t) + 5 X(s)u (s)ds, nn0n
= 5u (t) + X(t)5u (t) - Г X'(s)5u„ (s)ds, n n0 n
= (1 + X(t))5u (t) - ft X'(s>„ (s)ds, n0 n this yields the stationary conditions:
''1 + X ( t ) = 0,
X s ) , = t = 0.
Thus
X( s) = -1, and we obtain the following iteration formula t'
u„+1 (t) = un (t) - [u (s) - aun (s) - bun ([-])]ds.
0n2
So we have u 0( t ) = u 0, t'
u (t) = uo(t) - [u (s) - au0 (s) - bu0 ([-])]ds, t'
u 2 ( t ) = ux ( t ) - [ u ( s ) - aux ( s ) - bux ([—])] ds ,
-
2 1 01 11
t'
u 3 ( t ) = u 2 ( t ) - [ u ( s ) - au 2 ( s ) - bu 2 ([—])] ds ,
-
3 2 02 22
u ( t ) = lim un ( t ) .
n ^ro
In order to overcome the difficulty from the greatest integer function [ • ] , we can consider the above iteration formula in a series of intervals : [0,2) , [2,4) , [4,6) , _, [2 n , 2( n + 1)) where n = 0,1,2, ^. Following this way, each integral in iteration formulas will be easily computed.
The above analysis yields the following theorem:
Theorem 2 . The VIM solution of (1) can be determined by (9) with the iterations (8).
In the next section, the VIM is successfully applied for solving a linear EPCA.
4. APPLICATION
In what follows, we will apply the VIM method to a physical model to illustrate the strength of the method.
Let a = 2 , b = - 3 and U o = 1 in (1), this in turn gives the successive approximations by (8).
When 1 e[0,2), u1,0( t) = 1, u1,1(1) = u1,0 (1) - C[u' 0 (5) - 2u1,0 (s) + 3u1,0 ([s])]ds , , 0 1,0 ,,
= 1- J 0 (-2+3) ds
= 1 -1, u1,2(t) = u (t) - Г[u (s) - 2и1Д (s) + 3ихл ([s])]ds
, , 0 1,1 ,,
= 1 - 1 -
Jo t [- 1 — 2(1 - s ) + 3^ , (0) ] ds
= 1 - 1 - 1 2,
U 13 ( t ) = U 12 ( t ) - f [ и ( s ) - 2 U 12 ( s ) + 3uV 2 ([ s ])] ds
, , 0 1,2 ,,
= U j 2 ( t ) - J |^- 1 - 2 s - 2(1 - s - s 2) + 3 и 2(0) ^ ds = 1 - 1 - 1 2 - ~ t 3,
U 1,4 ( 1 ) = U 1,3 ( 1 ) - f 1 [ U ' 3( s ) - 2 u 1,3 ( s ) + 3 u 1,3 ([ s ])] ds
, , 0 1,3 , ,
U\,3
( t ) - l o t - 3
t
= 1 — t — t 2-- t3
4 2 ^xx ,
+ j s + 3 ц 3(0) ds
3 ,
When 1 e [2,4) ,
U 2,o( 1 ) = 1,
s
( 5 ) - 2 u2 0 (s) + 3 u 2,o ([-D] d s
U 2,1( 1 ) = U 2,0 ( 1 ) - J0 t [ U
= 1 - Jo (-2 + 3) ds
= 1 - t ,
U 2,2 ( 1 ) = U 2,1 ( 1 ) - f t [ U' , ( s ) - 2 U 2,1 ( s ) + 3u 2,1 ([ 5 ])] ds , , 0 2,1 , ,
= 1 -t - J [-1 — 2(1 — s) + 3u21 (1)] ds
= 1 + 2 t - 1 2,
ts u2,3 (1) = u2,2 (1) -J [u2,2 (s) — 2u2,2 (s) + 3u2,2 ([^ds
= u^ 2 ( 1 ) - J 1 [ 2 - 2 s - 2(1 + 2 s - s 2) + 3 u 2 2 (1) ] ds
= 1 - 4 1 + 2 1 2 - 2 1 3,
t' s u2,4 (1) = u2,3 (1) -J [\3 (s) - 2u2,3 (s) + 3u2,3 ([^ds
= u 2,3 ( 1 ) - L
- 6 + 12 s - 6 s + — s + 3 u 23(1) ds
= 1 + 7 1 - 4 1 2 + 4 1 3 -1 1 4, 33
When 1 e[4,6), u 30( 1) = 1
u 3,1 ( t ) = u 3,0 ( t ) - Jo 1 [ u 3,0( s ) - 2 u 3,0 ( s ) + 3 u 3,0 ([ s ])] ds
= 1 - J ( - 2 + 3 ) ds = 1 - 1 ,
t' s u3,2 (1) = u3,1 (1) - In [u3 , (s) - 2u3,1 (s) + 3u3,1 ([-])]ds
, , 0 3,1 , , 2
= u 3,1 ( 1 ) - Jo [- 1 - 2(1 - s ) + 3 u 3,1 (2) ] ds = 1 + 5 1 - 1 2,
u 3,3 ( 1 ) = u 3,2 ( 1 ) - fn [ u ' 2 ( s ) - 2 u 3,2 ( s ) + 3 u 3,2 ([ s ])] ds , , 0 3,2 , ,
= u3,2(1)-Jd[5-2s-2(1 + 5s-s2) + 3u32(2)]ds = 1 -191 + 512 -113, u3,4 (1) = u3,3 (1) - J1 [u' 3 (s) - 2u3,3 (s) + 3u3,3 ([s])]ds
, , 0 3,3 , , 2
=u3,3( о - L
—21 + 48 s —12 s + — s + 3 u^ 3(2) ds
= 1 + 69 1 - 19 1 2 + 1° 1 3 - 1 1 4, 33
In a word, in the interval [2 n , 2( n + 1)) , n = 0,1,2, we have the following iteration formulas:
un +1,0( t) = 1, un +1,m (t) = un +1,m-1 (t)
-Jo [ u n + 1m - 1 ( s ) - 2 u n + 1, m - 1 ( s ) + 3 u n + 1, m - 1 ( n )] ds
In view of (9), we can obtain the approximate analytical solution. Usually, the m + 1 th approximation is used for numerical purposes.
In Table 1, we compare the absolute errors of the VIM ( m = 4,5) with the ones of the 1-Gauss-Legendre method[16] using h = 0.02 . The graphs of the exact solution and 7th approximate solutions are shown in Fig. 1.
Table. 1.Comparison of the absolute errors
t |
VIM |
1-Gauss- Legendre method |
|
m=4 |
m=5 |
||
0.7 |
3.09E-6 |
5.47E-7 |
2.33E-4 |
1.4 |
4.14E-6 |
7.72E-7 |
4.16E-4 |
2.1 |
5.16E-6 |
5.17E-6 |
3.81E-3 |
2.8 |
4.34E-5 |
6.14E-6 |
5.42E-3 |
3.5 |
5.51E-5 |
8.75E-6 |
8.15E-3 |
4.2 |
7.36E-5 |
1.08E-5 |
1.96E-2 |
4.9 |
3.84E-4 |
2.92E-5 |
2.37E-2 |
5.6 |
8.26E-4 |
6.52E-5 |
4.73E-2 |
0.8
0.6
0.4
0.2

0.8
0.6
0.4
0.2

Fig. 1. Comparison of the exact solution with the VIM solution, Upper: exact solution, Lower: 7th VIM solution
-
5. CONCLUSION
In this paper, we have demonstrated the feasibility of the VIM for solving a EPCA model. We obtain high-accuracy approximate solutions. The numerical results also show that the VIM offers a very effective and convenient approach to the approximate solution of EPCA.
Список литературы Variational Iteration Method for Solving Differential Equations with Piecewise Constant Arguments
- G.Q. Wang, “Periodic solutions of a neutral differential equation with piecewise constant arguments,” J. Math. Anal. Appl., vol. 326, pp. 736-747, 2007.
- M.U. Akhmet, “On the reduction principle for differential equations with piecewise constant argument of generalized type,” J. Math. Anal. Appl., vol. 336, pp. 646-663, 2007.
- M.U. Akhmet, “Asymptotic behavior of solutions of differential equations with piecewise constant arguments,” Appl. Math. Lett.., vol. 21, pp. 951-956, 2008.
- M.U. Akhmet, “Stability of differential equations with piecewise constant arguments of generalized type,” Nonlinear Anal., vol.68, pp. 794-803, 2008.
- D. Altintan, “Extension of the logistic equation with piecewise constant arguments and population dynamics,” M.Sc.Thesis, Middle East Technical University, 2006.
- K.L. Cooke and J. Weiner, “A survey of differential equation with piecewise continuous argument,” Lecture Notes in Mathematics, Springer, Berlin, 1991.
- L. Dai and L. Fan, “Analytical and numerical approaches to characteristics of linear and nonlinear vibratory systems under piecewise discontinuous disturbances,” Commun. Nonlinear Sci. Numer. Simul., vol. 9, pp. 417–429, 2004.
- F. Gurcan and F. Bozkurt, “Global stability in a population model with piecewise constant arguments,”J. Math. Anal. Appl., vol. 360, pp. 334–342, 2009.
- K.L. Cooke and J. Wiener, “Retarded differential equations with piecewise constant delays,” J. Math. Anal. Appl., vol. 99, pp. 265–297, 1984.
- S.M. Shah and J. Wiener, “Advanced differential equations with piecewise constant argument deviations,” Int. J. Math. Math. Sci., vol. 6, pp. 671–703, 1983.
- A.F. Ivanov, “Global dynamics of a differential equation with piecewise constant argument,” Nonlinear Anal., vol. 71, pp. e2384–e2389, 2009.
- Y. Muroya, “New contractivity condition in a population model with piecewise constant arguments,” J. Math. Anal. Appl., vol. 346, pp. 65–81, 2008.
- Y.H. Xia, Z.K. Huang, and M.A. Han, “Existence of almost periodic solutions for forced perturbed systems with piecewise constant argument,” J. Math. Anal. Appl., vol. 333, pp. 798–816, 2007.
- J. Wiener, Generalized Solutions of Functional Differential Equations, Singapore: World Scientific, 1993.
- M.H. Song, Z.W. Yang, and M.Z. Liu, “Stability of -methods for advanced differential equations with piecewise continuous arguments,”Comput.Math. Appl., vol. 49, pp. 1295-1301, 2005.
- M.Z. Liu, S.F. Ma, and Z.W. Yang, “Stability analysis of Runge-Kutta methods for unbounded retarded differential equations with piecewise continuous arguments,” Appl. Math. Comput.., vol. 191, pp. 57-66, 2007.
- Q. Wang, Q.Y. Zhu, and M.Z. Liu, “Stability and oscillations of numerical solutions of differential equations with piecewise continuous arguments of alternately advanced and retarded type,” J. Comput. Appl. Math.., vol. 235, pp. 1542-1552, 2011.
- J.H. He, “Variational iteration method-A kind of non-linear analytical technique: Some examples,” Internat. J. Non-Linear Mech., vol. 34, no. 4, pp. 699-708, 1999.
- J.H. He, “Variational iteration method-Some recent results and new interpretations,” J. Comput. Appl. Math.., vol. 207, no. 1, pp. 3-17, 2007.
- J.H. He and X.H. Wu, “Variational iteration method: New development and applications,"Comput. Math. Appl., vol. 54, no.7-8, pp. 881-894, 2007.
- J.H. He, G.C. Wu, and F. Austin, “The variational iteration method which should be followed,” Nonlinear Sci. Lett. A, vol. 1, no.1, pp. 1-30, 2010.
- M.M. Al-Sawalha, M.S.M. Noorani, and I. Hashim, “On accuracy of Adomian decompositionmethod for hyperchaotic Rössler system,” Chaos Solitons Fractals., vol. 40, no.4, pp. 1801-1807, 2009.
- M.M. Al-Sawalha and M.S.M. Noorani, “Application of the differential transformation method for the solution of the hyperchaotic Rössler system,” Comm. Non. Sci. Num. Simu., vol. 14, pp. 1509-1514, 2009.
- F.M. Allan, “Construction of analytic solution to chaotic dynamical systems using the Homotopy analysis method,” Chaos Solitons Fractals., vol. 39, no.4, pp. 1744-1752, 2009.
- J.H. He, “Homotopy perturbation method for bifurcation of nonlinear problems,” Int. J. Nonlinear Sci. Numer. Simu., vol. 6, pp. 207-208, 2005.
- J.H. He, “Some asymptotic methods for strongly nonlinear equations,” Int. J. Modern Phys.B, vol. 20, no.10, pp. 1141–1199, 2006.
- J.H. He, “Approximate analytical solution for seepage flow with fractional derivatives in porous media,” Comput. Meth. Appl. Mech. Eng., vol. 167, pp. 57–68, 1998.
- J.H. He, “Variational iteration method for autonomous ordinary differential systems,” Appl. Math. Comput.., vol. 114, pp. 115–123, 2000.
- J.H. He, “An elementary introduction to recently developed asymptotic methods and nanomechanics in textile engineering,” Int. J. Modern Phys.B, vol. 22, no.21, pp. 3487–3578, 2008.