Hindawi Publishing Corporation The Scientific World Journal Volume 2013, Article ID 159194, 6 pages http://dx.doi.org/10.1155/2013/159194
Research Article Dual Synchronization of Fractional-Order Chaotic Systems via a Linear Controller Jian Xiao, Zhen-zhen Ma, and Ye-hong Yang College of Mathematics and Statistics, Chongqing University, Chongqing 401331, China Correspondence should be addressed to Jian Xiao;
[email protected] Received 11 July 2013; Accepted 13 August 2013 Academic Editors: N. Herisanu, T. Li, Y. Xia, and Q. Xie Copyright Β© 2013 Jian Xiao et al. 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. The problem of the dual synchronization of two different fractional-order chaotic systems is studied. By a linear controller, we realize the dual synchronization of fractional-order chaotic systems. Finally, the proposed method is applied for dual synchronization of Van der Pol-Willis systems and Van der Pol-Duffing systems. The numerical simulation shows the accuracy of the theory.
1. Introduction In recent years, the topic of chaos synchronization has attracted increasing attention in many fields. The result of synchronization of chaotic oscillators is used in nonlinear oscillators [1], circuit experiment [2], secret communication [3], and some other fields. In 1990, the first concept of synchronization was presented by Carroll and Perora [4]. And there are many methods about chaos synchronization such as Lyapunov equation [5], Perora-Carroll (PC) [4] and backstepping control [6]. All of these methods are amid of the synchronization between one master and one slave system do not consist of the synchronization of multimaster systems and multislave systems. Dual synchronization is a special circumstance in synchronization of chaotic oscillators. The first idea of multiplexing chaos using synchronization was investigated in a small map and an electronic circuit model by Tsimring and Sushchik in 1996 in [7]; then the concept of dual synchronization was raised by Liu and Davids in 2000 in [8], which concentrates on using a scalar signal to simultaneously synchronize two different pairs chaotic oscillators, that is, the synchronization between two master systems and two slave systems. Nowadays, there are many dual synchronization methods, such as in 2000 Liu and Davids introduce the dual
synchronization of 1-D discrete chaotic systems via specific classes of piecewise-linear maps with conditional linear coupling in [8]. The dual synchronization between the Lorenz and Rossler systems by the Lyapunov stabilization theory is investigated in [9]. The output feedback strategy is used to study the dual synchronization of two different 3-D continuous chaotic systems in [10]. Then the dual synchronization in modulated time-delayed systems is investigated by designing a delay feedback controller in [11]. All of these works are amid of the dual synchronization of integer-order chaotic systems and do not consist of the dual synchronization of fractional-order chaotic systems. In this paper, a new method of dual synchronization of fractional-order chaotic systems is proposed, by a linear controller; the dual synchronization of chaos is obtained. The rest of this paper is organized as follows: in Section 2, we construct a theory frame about the dual synchronization of two different fractional-order chaotic systems. By a linear controller, we obtain dual synchronization between two different fractional-order chaotic systems in Section 3. In Section 4, the proposed method is applied to dual synchronization of Van der Pol-Willis systems and Van der Pol-Duffing systems for evaluating the performance of the method and by numerical simulation; the result shows that the controller designed by the application of this method is effective. Finally, conclusions are drawn in Section 5.
2
The Scientific World Journal
2. Problem Analysis
The error signal for dual synchronization is πβπ₯ ] = πΆπ (π β π) . πβπ¦
We define the following two systems as two master systems.
π = Vπ β Vπ = [π΄ π΅] [
Master 1: ππΌ π₯ = π (π‘, π₯) , ππ‘πΌ
(1)
Master 2:
lim βπ (π‘) β π₯ (π‘)β = 0,
π π¦ = π (π‘, π¦) , ππ‘πΌ
π
π=π
π=1
Vπ = β ππ π₯π + β ππ π¦π = [π1 , π2 , . . . , ππ ] π₯ + [π1 , π2 , . . . , ππ ] π¦ π₯ = π΄π₯ + π΅π¦ = [π΄ π΅] [ ] = πΆπ π, π¦ (3) where π΄ = [π1 , π2 , . . . , ππ ]π and π΅ = [π1 , π2 , . . . , ππ ]π are two known matrices, and ππ , ππ , π = 1, 2, . . . , π, π = 1, 2, . . . , π π
cannot be zero at the same time. So π = [π΄
(8)
π
where β β
β is the Euclidian norm.
3. Dual Synchronization Strategy Lemma 1. Considering the fractional-order system π·πΌ π§ (π‘) = ππ§,
ππΌ π = π (π‘, π) + π(1) , ππ‘πΌ
σ΅¨ πΌπ σ΅¨σ΅¨ σ΅¨σ΅¨arg (eig (πΊ (π‘) + πΎπΆπ ))σ΅¨σ΅¨σ΅¨ β₯ , σ΅¨ σ΅¨ 2
Proof. We can rewrite (1) and (2) in the following form by π defining Ξ¨ = [π π] :
(4)
ππΌ π₯ [ ππ‘πΌ ] π (π‘, π₯) [ ] ], [ ]=[ [ ππΌ π¦ ] π (π‘, π¦)
πΌ
(5)
where π = [π1 , π2 , . . . , ππ ]π and π = [π1 , π2 , . . . , ππ ]π are the state vectors of the two slave systems, π(1) = [π’11 , π’21 , . . . , 2 π ] are vectors of manipulated π’π1 ]π and π(2) = [π’12 , π’22 , . . . , π’π variables, and πΌ β (0, 1] is the order of the two slave systems. Similarly, by a linear combination of two slave systems states, a signal Vπ is generated as follows: π
π=π
π=1
Vπ = β ππ ππ + β ππ ππ = [π1 , π2 , . . . , ππ ] π + [π1 , π2 , . . . , ππ ] π
ππΌ π = Ξ¨ (π‘, π) . ππ‘πΌ
(11)
[ ππ‘πΌ ]
Similarly, (4) and (5) can be rewritten as ππΌ π [ ππ‘πΌ ] π (π‘, π) + π(1) [ ] ], [ ]=[ [ ππΌ π ] π (π‘, π) + π(2)
ππΌ π = Ξ¨ (π‘, π) + π, ππ‘πΌ
[ ππ‘πΌ ]
(12) (1)
] = where π = [(π(1) )π (π(2) )π ]π , one defines π = [ π π(2)
(6)
π = π΄π + π΅π = [π΄ π΅] [ ] = πΆπ π, π π
(10)
where πΊ(π‘) is the coefficient matrix of master systems and πΎ is a control gain vector.
Slave 2: π π = π (π‘, π) + π(2) , ππ‘πΌ
(9)
Theorem 2. The dual synchronization of fractional-order chaotic systems between the master systems and the slave systems is achieved if and only if the following condition satisfies
π΅ ] is a known
Slave 1:
π§ (0) = π§0 ,
where 0 < πΌ β€ 1, π§ β π
π , and π β π
πΓπ ; then system (9) is stable if and only if | arg(π π (π))| β₯ (πΌπ/2), π = 1, 2, . . ., where arg(π π (π)) denotes the argument of the eigenvalue π π of π.
π π
matrix. π = [π₯π π¦π ] is a combination of the two master systems states. The corresponding two slave systems are as follows:
π
σ΅© σ΅© lim σ΅©σ΅©π (π‘) β π¦ (π‘)σ΅©σ΅©σ΅© = 0,
π‘ββ σ΅©
(2)
where π₯ = [π₯1 , π₯2 , . . . , π₯π ]π and π¦ = [π¦1 , π¦2 , . . . , π¦π ]π are the state vectors of the two master systems. π β πΆ[π
+ Γ π
π , π
π ] and π β πΆ[π
+ Γ π
π , π
π ] are two known functions. πΌ β (0, 1] is the order of the two master systems. By a linear combination of the two master systems states, a signal Vπ is give as π
The main goal is to synchronize the master systems and the slave systems is equivalent to π‘ββ
πΌ
(7)
where π = [ππ ππ ] is a combination of the two slave systems states.
πΎπ
πβπ₯
[ πΎ12 π ] and πΈ = [ πβπ¦ ]. Equation (12) is transformed into
ππΌ π [ ππ‘πΌ ] π (π‘, π) + πΎ1 π [ ] [ ππΌ π ] = [ π (π‘, π) + πΎ π ] , 2 [ ππ‘πΌ ]
(13)
The Scientific World Journal
3
so the error system is transformed into ππΌ π ππΌ π₯ β πΌ] πΌ πΌ π πΈ [ ππ‘πΌ ππ‘ ] π π ππΌ π [ = β . = ππ‘πΌ [ ππΌ π ππΌ π¦ ] ππ‘πΌ ππ‘πΌ β [ ππ‘πΌ ππ‘πΌ ]
Master 2: Willis system ππΌ π¦1 = π¦2 , ππ‘πΌ
(14)
The error system is obtained as
(21)
ππΌ π¦2 = ππ¦1 + ππ¦12 + ππ¦13 + ππ¦2 + π2 cos π‘. ππ‘πΌ So the corresponding slave systems are as follows:
ππΌ π ππΌ π β ππ‘πΌ ππ‘πΌ
Slave 1: ππΌ π1 = π1 β πΎπ13 β π½π2 + π1 cos π‘ + π1 π, ππ‘πΌ
= Ξ¨ (π‘, π) + πΎπ β Ξ¨ (π‘, π) = Ξ¨ (π‘, π) β Ξ¨ (π‘, π) + πΎπ
ππΌ π2 = π (π1 β ππ2 + π) + π2 π, ππ‘πΌ
= Ξ¨ (π‘, π + πΈ) β Ξ¨ (π‘, π) + πΎ (π β π) = Ξ¨ (π‘, π + πΈ) β Ξ¨ (π‘, π) + πΎπΆπ πΈ.
Slave 2: (15)
Using the first-order Taylor expansion, the function Ξ¨(β
) is rewritten as Ξ¨ (π‘, π + πΈ) β Ξ¨ (π‘, π) =
πΞ¨ (π‘, π) πΈ + h.o.t = πΊ (π‘) πΈ + h.o.t, ππ
(16)
where h.o.t denotes the higher order terms of the series. We substitute (16) into (15) and yield ππΌ πΈ = πΊ (π‘) πΈ + h.o.t + πΎπΆπ πΈ = [πΊ (π‘) + πΎπΆπ ] πΈ + h.o.t. ππ‘πΌ (17) We can transfer the (17) into ππΌ πΈ = [πΊ (π‘) + πΎπΆπ ] πΈ ππ‘πΌ
(18)
according to Lemma 1, we can know that the error system is asymptotically stable at zero if and only if the following condition is satisfied σ΅¨ πΌπ σ΅¨σ΅¨ σ΅¨σ΅¨arg (eig (πΊ (π‘) + πΎπΆπ ))σ΅¨σ΅¨σ΅¨ β₯ . σ΅¨ σ΅¨ 2
(19)
ππΌ π2 = ππ1 + ππ12 + ππ13 + ππ2 + π2 cos π‘ + π4 π, ππ‘πΌ
Master 1: Van der Pol system πΌ
(20)
(23)
where π = π1 π1 + π2 π2 + π1 π3 + π2 π4 , π1 = π1 β π₯1 , π2 = π2 β π₯2 , π3 = π1 β π¦1 , and π4 = π2 β π¦2 . The πΊ(π‘) matrix of the master systems is achieved as [ πΊ (π‘) = [ [ [ [ [
1 β 3πΎπ₯12 π 0 0
βπ½ 0 ππ 0 0 0 0 π + 2ππ¦1 + 3ππ¦12
0 0] ], 1] ] π] ]
(24)
so the corresponding error matrix are as follows: ππΌ π1 ππ‘πΌ ( ππΌ π2 ) ( ) ( ππ‘πΌ ) ( ) ( πΌ ) (π π ) 3) ( ( πΌ ) ππ‘ ππΌ π4 ( ππ‘πΌ )
=(
Example 3 (dual synchronization of Van der Pol-Willis systems). In the first example, we can use the proposed method to achieve the dual synchronization of the Van der Pol system and the Willis system.
ππΌ π₯2 = π (π₯1 β ππ₯2 + π) . ππ‘πΌ
ππΌ π1 = π2 + π3 π, ππ‘πΌ
1β3πΎπ₯12 + π1 π1 βπ½ + π2 π1
4. The Example Analysis and Numerical Simulations
π π₯1 = π₯1 β πΎπ₯13 β π½π₯2 + π1 cos π‘, ππ‘πΌ
(22)
π1 π1
π2 π1
π + π1 π2
ππ + π2 π2
π1 π2
π2 π2
π1 π3
π2 π3
π1 π3
1 + π2 π3
π1 π4
π2 π4
)
π+2ππ¦1 +3ππ¦12 +π1 π4 π+π2 π4
π1 π2 Γ( ). π3 π4
(25)
We should choose the appropriate parameters so that all the eigenvalues of the Jacobian matrix of (25) satisfy Matignon condition; that is, the eigenvalues evaluated at the equilibrium point are satisfied: σ΅¨ πΌπ σ΅¨σ΅¨ σ΅¨σ΅¨arg (eig (πΊ (π‘) + πΎπΆπ ))σ΅¨σ΅¨σ΅¨ > . σ΅¨ σ΅¨ 2
(26)
4
The Scientific World Journal
0.1
0.3
0
0.2
β0.1
0.1 e3 e4
e1 e2
0.4
β0.2
β0.1
β0.3
β0.2
β0.4 β0.5
0
β0.3
0
10
20
30
40
50
60
70
80
90
100
β0.4
t
e1 e2
0
10
20
30
40
50 t
60
70
80
90
100
e3 e4
Figure 1: Error signals between the pair of Van der Pol system.
Figure 2: Error signals between the pair of Willis system.
The eigenvalue equation of the equilibrium point is locally asymptotically stable. From what we have discussed above, we can know that π΄ and π΅ are two known matrices; the parameter πΎ can be appropriately selected for satisfying the Matignon condition. Dual synchronization of the Van der Pol system and the Willis system is simulated. The system parameters are set to be πΎ = 1/3, π½ = 1, π1 = 0.74, π = 0.1, π = 0.8, π = 0.7, π = β0.9, π = 3, π = β2, π = β0.1, π2 = 0.1, π΄ = [1, 1, 1], π΅ = [1, 1, 1], and πΌ = 1, so
Example 4 (dual synchronization of Van der Pol and Duffing systems). For Example 4, the dual synchronization of Van der Pol and Duffing systems is investigated.
1 β π₯12 + π1
β1 + π1
π1
π1
0.1 + π2
β0.08 + π2
π2
π2
π3
π3
π3
1 + π3
π4
ππΌ π₯1 = π₯1 β πΎπ₯13 β π½π₯2 + π1 cos π‘, ππ‘πΌ ππΌ π₯2 = π (π₯1 β ππ₯2 + π) . ππ‘πΌ
(28)
Master 2: Duffing system
πΊ (π‘) + πΎπΆπ
=(
Master 1: Van der Pol system
π4
β0.9 +
6π¦1 β 6π¦12
ππΌ π¦1 = π¦2 , ππ‘πΌ ππΌ π¦2 = ππ¦1 + ππ¦13 + ππ¦2 + π2 cos π‘. ππ‘πΌ
).
+ π4 β0.1 + π4
(29)
So the corresponding slave systems are (27)
If β295 < π1 < β130, π2 = β0.1, π3 = β1, and π4 = β400, which satisfy (18), the eigenvalue equation of the equilibrium point is locally asymptotically stable. We choose π1 = β210, π2 = β0.1, π3 = β1, and π4 = β400. The initial conditions of the master system 1 and the master system 2 are taken as π₯1 (0) = 0.1, π₯2 (0) = 0.2 and π¦1 (0) = 0.2, π¦2 (0) = 0.3; the initial conditions of the slave system 1 and the slave system 2 are taken as π1 (0) = 0.3, π2 (0) = 0.4 and π1 (0) = 0.5, π2 (0) = 0.6, so the initial conditions of the error system are set to be π1 (0) = 0.2, π2 (0) = 0.2, π3 (0) = 0.3, and π4 (0) = 0.3. In Figures 1 and 2, we can see that all error variables have converged to zero; that is, we achieve the dual synchronization between the Van der Pol and the Willis systems.
Slave 1: ππΌ π1 = π1 β πΎπ13 β π½π2 + π1 cos π‘ + π1 π, ππ‘πΌ ππΌ π2 = π (π1 β ππ2 + π) + π2 π, ππ‘πΌ
(30)
Slave 2: ππΌ π1 = π2 + π3 π, ππ‘πΌ ππΌ π2 = ππ1 + ππ13 + ππ2 + π2 cos π‘ + π4 π, ππ‘πΌ
(31)
where π = π1 π1 +π2 π2 +π1 π3 +π2 π4 , π1 = π1 βπ₯1 , π2 = π2 βπ₯2 , π3 = π1 β π¦1 , and π4 = π2 β π¦2 .
5
0.2
0.2
0.1
0.1
0
0
e3 e4
e1 e2
The Scientific World Journal
β0.1
β0.1
β0.2
β0.2
β0.3
β0.3
β0.4
0
50
100
150
β0.4
0
50
100
150
t
t
e1 e2
e3 e4
Figure 3: Error signals between the pair of Van der Pol system.
Figure 4: Error signals between the pair of Duffing system.
The πΊ(π‘) matrix of the master systems is achieved as
matrices, the parameter πΎ can be appropriately selected for satisfying the Matignon condition. According to what we have studied above, parameters are set to πΎ = 1/3, π½ = 1, π1 = 0.74, π = 0.1, π = 0.8, π = 0.7, π = 1, π = β1, π = β0.15, π2 = 0.3, π΄ = [1, 1, 1], π΅ = [1, 1, 1], and πΌ = 0.98, so
[ [ πΊ (π‘) = [ [ [ [ [
1β
3πΎπ₯12
π 0
βπ½ ππ 0
0
0 0 0
0 0] ] ]. 1] ] 0 π + 3ππ¦12 π ] ]
(32)
πΊ (π‘) + πΎπΆπ
So the corresponding error matrix are as follows: ππΌ π1 ππ‘πΌ
=(
( ππΌ π2 ) ( ) ( ππ‘πΌ ) ( ) ( πΌ ) (π π ) 3) ( ( πΌ ) ππ‘ ππΌ π4 ( ππ‘πΌ ) 1 β 3πΎπ₯12 + π1 π1 βπ½ + π2 π1 =(
π1 π1
π2 π1
π + π1 π2
ππ + π2 π2
π1 π2
π2 π2
π1 π3
π2 π3
π1 π3
1 + π2 π3
π1 π4
π2 π4
π1 π Γ (π2 ) . 3 π4
)
π + 3ππ¦12 + π1 π4 π + π2 π4
(33)
We should choose the appropriate parameters so that all the eigenvalues of the Jacobian matrix of (33) satisfy Matignon condition; that is, the eigenvalues evaluated at the equilibrium point are satisfied: σ΅¨σ΅¨σ΅¨arg (eig (πΊ (π‘) + πΎπΆπ ))σ΅¨σ΅¨σ΅¨ > πΌπ . (34) σ΅¨σ΅¨ σ΅¨σ΅¨ 2 The eigenvalue equation of the equilibrium point is locally asymptotically stable. Because π΄ and π΅ are two known
1 β π₯12 + π1
β1 + π1
π1
π1
0.1 + π2
β0.08 + π2
π2
π2
π3
π3
π3
1 + π3
π4
π4
1 β 3π¦12 + π4 β0.15 + π4
).
(35)
If β275 < π1 < β117, π2 = β0.1, π3 = β1, and π4 = β400, which satisfy (34), the eigenvalue equation of the equilibrium point is locally asymptotically stable. We choose π1 = β200, π2 = β0.1, π3 = β1, and π4 = β400. The initial conditions of the master system 1 and the master system 2 are taken as π₯1 (0) = 0.1, π₯2 (0) = 0.2 and π¦1 (0) = 0.2, π¦2 (0) = 0.3, the initial conditions of the slave system 1 and the slave system 2 are taken as π1 (0) = 0.3, π2 (0) = 0.4 and π1 (0) = 0.5, π2 (0) = 0.6, so the initial conditions of the error system are set to be π1 (0) = 0.2, π2 (0) = 0.2, π3 (0) = 0.3, and π4 (0) = 0.3. In Figures 3 and 4, we can see that all error variables have converged to zero; that is, we achieve the dual synchronization between the Van der Pol and the Duffing systems.
5. Conclusions In this work, we construct a theory frame about dual synchronization of two different fractional-order chaotic systems and propose a method of dual synchronization. In addition, this method is used for designing a synchronization controller to
6 achieve the dual synchronization of two different fractionalorder chaotic systems. Finally, the proposed method is applied for dual synchronization of the Van der Pol-Willis systems and the Van der Pol-Duffing systems. The numerical simulations proves the accuracy of the theory.
Acknowledgment This work is supported by the Fundamental Research Funds for the Central Universities of China under Grant no. CQDXWL-2012-007.
References [1] G. M. Mahmoud, E. E. Mahmoud, A. A. Farghaly, and S. A. Aly, βChaotic synchronization of two complex nonlinear oscillators,β Chaos, Solitons & Fractals, vol. 42, no. 5, pp. 2858β2864, 2009. [2] Z. Xu, C.-X. Liu, and T. Yao, βStudy on a new chaotic system with analysis and circuit experiment,β Acta Physica Sinica, vol. 59, no. 1, pp. 131β139, 2010. [3] D. Matignon, βStability results for fractional differential equations with application to control processing,β in Proceedings of the International IMACS IEEE-SMC Multiconference on Computational Engineering in Systems Applications, vol. 2, pp. 963β968, Lille, France, 1996. [4] T. Carroll and L. M. Perora, βSynchronization of chaotic systems,β Physical Review Letters, vol. 64, pp. 821β824, 1990. [5] J.-B. Hu, Y. Han, and L.-D. Zhao, βSynchronizing fractional chaotic systems based on Lyapunov equation,β Acta Physica Sinica, vol. 57, no. 12, pp. 7522β7526, 2008. [6] J.-B. Hu, Y. Han, and L.-D. Zhao, βA novel stablility theorem for fractional systems and its applying in synchronizing fractional chaotic system based on back-stepping approach,β Acta Physica Sinica, vol. 58, no. 4, pp. 2235β2239, 2009. [7] L. S. Tsimring and M. M. Sushchik, βMultiplexing chaotic signals using synchronization,β Physics Letters A, vol. 213, no. 3-4, pp. 155β166, 1996. [8] Y. Liu and P. Davids, βDual synchronization of chaos,β Physical Review E, vol. 61, pp. 2176β2184, 2000. [9] D. Ning, J.-A. Lu, and X. Han, βDual synchronization based on two different chaotic systems: Lorenz systems and RΒ¨ossler systems,β Journal of Computational and Applied Mathematics, vol. 206, no. 2, pp. 1046β1050, 2007. [10] H. Salarieh and M. Shahrokhi, βDual synchronization of chaotic systems via time-varying gain proportional feedback,β Chaos, Solitons & Fractals, vol. 38, no. 5, pp. 1342β1348, 2008. [11] D. Ghosh and A. R. Chowdhury, βDual-anticipating, dual and dual-lag synchronization in modulated time-delayed systems,β Physics Letters A, vol. 374, no. 34, pp. 3425β3436, 2010.
The Scientific World Journal
Advances in
Operations Research Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Advances in
Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Applied Mathematics
Algebra
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Probability and Statistics Volume 2014
The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Differential Equations Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Submit your manuscripts at http://www.hindawi.com International Journal of
Advances in
Combinatorics Hindawi Publishing Corporation http://www.hindawi.com
Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of Mathematics and Mathematical Sciences
Mathematical Problems in Engineering
Journal of
Mathematics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Discrete Mathematics
Journal of
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Discrete Dynamics in Nature and Society
Journal of
Function Spaces Hindawi Publishing Corporation http://www.hindawi.com
Abstract and Applied Analysis
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Journal of
Stochastic Analysis
Optimization
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014