Next Article in Journal
Evolutionary Relation of Mathematical Physics Equations Evolutionary Relation as Foundation of Field Theory Interpretation of the Einstein Equation
Next Article in Special Issue
Nonlinear Control System Design of an Underactuated Robot Based on Operator Theory and Isomorphism Scheme
Previous Article in Journal
Applications of Coupled Fixed Points for Multivalued Maps in the Equilibrium in Duopoly Markets and in Aquatic Ecosystems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Weak Turnpike Property for Perturbed Dynamical Systems with a Lyapunov Function

by
Alexander J. Zaslavski
Department of Mathematics, Technion–Israel Institute of Technology, Haifa 32000, Israel
Axioms 2021, 10(2), 45; https://doi.org/10.3390/axioms10020045
Submission received: 12 February 2021 / Revised: 26 March 2021 / Accepted: 27 March 2021 / Published: 29 March 2021
(This article belongs to the Special Issue Advances in Analysis and Control of Systems with Uncertainties)

Abstract

:
In this work, we obtain a weak version of the turnpike property of trajectories of perturbed discrete disperse dynamical systems, which have a prototype in mathematical economics.

1. Introduction

In [1,2], A. M. Rubinov introduced a discrete dispersive dynamical system, which was investigated in [1,2,3,4,5,6,7]. This dynamical system is determined by a set-valued mapping and has a prototype in the economic growth theory [1,8,9]. Our dynamical system is described by a compact metric space of states and a transition operator, which is set-valued. Usually in the dynamical systems theory a transition operator is single-valued. In [1,2,3,4,5,6,7] and in the present paper we study dynamical systems with a set-valued transition operator. Such dynamical systems correspond to certain models of economic dynamics [1,8,9].
Let ( X , ρ ) be a compact metric space and let a : X 2 X { } be a set-valued mapping of which the graph
graph ( a ) = { ( x , y ) X × X : y a ( x ) }
is a closed subset of X × X . For each nonempty subset E X set
a ( E ) = { a ( x ) : x E }   and   a 0 ( E ) = E .
By induction we define a n ( E ) for any positive integer n and any nonempty subset E X as follows:
a n ( E ) = a ( a n 1 ( E ) ) .
In the present work we investigate convergence and structure of trajectories of the perturbed dynamical system determined by the set-valued mapping a. Following [1,2] this system is called a discrete dispersive dynamical system.
A sequence { x t } t = 0 X is called a trajectory of a (or just a trajectory if the mapping a is understood) if x t + 1 a ( x t ) for all nonnegative integers t.
Let T 2 > T 1 be integers. A sequence { x t } t = T 1 T 2 X is called a finite trajectory of a (or just a trajectory if the mapping a is understood) if x t + 1 a ( x t ) for all integers t { T 1 , , T 2 1 } .
Define
Ω ( a ) = { z X : for very positive number ϵ there exists a trajectory { x t } t = 0
for   which   lim inf t ρ ( z , x t ) ϵ } .
By the compactness of X, Ω ( a ) is a nonempty closed subset of the metric space ( X , ρ ) . In the dynamical systems theory the set Ω ( a ) is called a global attractor of a. In [1,2] Ω ( a ) is called a turnpike set of a. This terminology is motivated by economic growth theory [1,8,9].
For each x X and each nonempty closed subset E X put
ρ ( x , E ) = inf { ρ ( x , y ) : y E } .
Evidently, for every trajectory { x t } t = 0 ,
lim t ρ ( x t , Ω ( a ) ) = 0 .
Let ϕ : X R 1 be a continuous function such that
ϕ ( z ) 0   for all   z X ,
ϕ ( y ) ϕ ( x ) for   all   x X and   all   y a ( x ) .
Evidently, the function ϕ is a Lyapunov function for the dynamical system generated by the mapping a. In economic growth theory usually X is a subset of the finite-dimensional Euclidean space and ϕ is a linear functional on this space [1,8,9]. Our goal in [7] was to study approximate solutions of the problem
ϕ ( x T ) max ,
{ x t } t = 0 T is   a   trajectory   satisfying   x 0 = x ,
where x X and a natural number T are given.
The following theorem was obtained in [7].
Theorem 1.
The following properties are equivalent:
(1) If a sequence { x t } t = X satisfies x t + 1 a ( x t ) and ϕ ( x t + 1 ) = ϕ ( x t ) for all integers t, then
{ x t } t = Ω ( a ) .
(2) For every positive number ϵ there exists a positive integer T ( ϵ ) such that for every trajectory { x t } t = 0 X satisfying ϕ ( x t ) = ϕ ( x t + 1 ) for all integers t 0 the relation ρ ( x t , Ω ( a ) ) ϵ is valid for all integers t T ( ϵ ) .
Our results are obtained under the assumption that property (1) holds. This property indeed holds for models of economic dynamics, which are prototypes of our dynamical system [1,8,9].
For each bounded function ψ : X R 1 set
ψ = sup { | ψ ( z ) | : z X } .
We denote by Card ( A ) the cardinality of a set A and suppose that the sum over empty set is zero.
For every point ( x 1 , x 1 ) , ( y 1 , y 2 ) X × X put
ρ 1 ( ( x 1 , x 2 ) , ( y 1 , y 2 ) ) = ρ ( x 1 , y 1 ) + ρ ( x 2 , y 2 ) .
For every point ( x 1 , x 2 ) X × X and every nonempty closed subset E X × X set
ρ 1 ( ( x 1 , x 2 ) , E ) = inf { ρ 1 ( ( x 1 , x 2 ) , ( y 1 , y 2 ) ) : ( y 1 , y 2 ) E } .
In [7] we established the turnpike properties for approximate solutions of the problem
ϕ ( x T ) max ,
{ x t } t = 0 T is   a   trajectory   satisfying   x 0 = x ,
where x X and a positive integer T are given. In [10] we established a weak version of the turnpike property that holds for all finite trajectories of our dynamical system, which are of a sufficient length and which are not necessarily approximate solutions of the problem above. This turnpike result usually holds for models of economic dynamics [1,8,9].
More precisely, in [10] we prove the following result.
Theorem 2.
Assume that property (1) of Theorem 1 holds and that ϵ > 0 . Then there exists a natural number L such that for each integer T > L and each finite trajectory { x t } t = 0 T the following inequality holds:
C a r d ( { t { 0 , , T } : ρ ( x t , Ω ( a ) ) > ϵ } ) L .
In this paper we show that a weak version of the turnpike property established in Theorem 2 is stable under small perturbations.
Turnpike properties are well known in mathematical economics. See, for example, references [2,8,9,11] and the references mentioned there. Recently it was shown that the turnpike phenomenon holds for many important classes of problems arising in various areas of research [12,13,14,15,16,17,18,19,20]. For related infinite horizon problems see [9,21,22,23,24,25,26,27,28].

2. The Main Results

For every pair of nonempty sets A , B X set
H ( A , B ) = max { sup { ρ ( x , B ) : x A } , sup { ρ ( y , A ) : y B } } .
We assume that the following assumption is true:
(A) for every positive number ϵ there is a positive number δ such that for every pair of points x , y X satisfying ρ ( x , y ) δ ,
H ( a ( x ) , a ( y ) ) ϵ .
We also assume that property (1) of Theorem 1 is true and obtain the following two theorems.
Theorem 3.
Let ϵ be a positive number. Then there is a positive integer L 0 such that for every natural number L > L 0 there is a positive number δ such that for each sequence { x t } t = 0 L X satisfying
ρ 1 ( ( x t , x t + 1 ) , g r a p h ( a ) ) δ , t = 0 , , T 1
the following relation holds:
C a r d ( { t { 0 , , L } : ρ ( x t , Ω ( a ) ) > ϵ } ) L 0 .
Theorem 4.
Let ϵ be a positive number. Then there is a positive integer L 1 and a positive number δ such that for every integer T > L 1 and every sequence { x t } t = 0 T X satisfying
ρ 1 ( ( x t , x t + 1 ) , g r a p h ( a ) ) δ , t = 0 , , T 1
the following relation holds:
T 1 C a r d ( { t { 0 , , T } : ρ ( x t , Ω ( a ) ) > ϵ } ) ϵ .
In Theorems 3 and 4 we deal with the structure of inexact trajectories of our dynamical system. They are important because in the real world applications computational errors and errors of measurements always take place.

3. An Auxiliary Result

The following lemma shows that our dynamical system has the so-called shadowing property [29,30].
Lemma 1.
Let ϵ > 0 and L be a positive integer. Then there is a positive number δ such that for every sequence { x t } t = 0 L X satisfying for all integers t = 0 , , L 1 ,
ρ 1 ( ( x t , x t + 1 ) , g r a p h ( a ) ) δ
there is a finite trajectory { y t } t = 0 L X for which
y 0 = x 0 ,
ρ ( x t , y t ) ϵ , t = 0 , , L .
Proof. 
Let
δ L = ϵ / 4 .
By induction and assumption (A), we define positive numbers δ i > 0 , i = 0 , , L 1 such that for every integer i { 1 , , L }
δ i 1 < δ i / 8
and for every pair of points x , y X for which ρ ( x , y ) δ i 1 we have
H ( a ( x ) , a ( y ) ) < δ i / 8 .
Set
δ = δ 0 .
Assume that { x t } t = 0 L X satisfies
ρ 1 ( ( x t , x t + 1 ) , graph ( a ) ) δ , t = 0 , , L 1 .
Set
x 0 = y 0 .
In view of (5) and (6),
ρ 1 ( ( y 0 , x 1 ) , graph ( a ) ) δ .
By (7), there exists
( z 0 , z 1 ) graph ( a )
such that
ρ ( y 0 , z 0 ) δ , ρ ( x 1 , z 1 ) δ .
It follows from (3), (4) and (9) that
H ( a ( y 0 ) , a ( z 0 ) ) δ 1 / 8 .
Equations (8) and (10) imply that
ρ ( z 1 , a ( y 0 ) ) δ 1 / 8 .
Equation (11) implies that there is
y 1 a ( y 0 )
for which
ρ ( y 1 , z 1 ) δ 1 / 4 .
It follows from (2), (4), (9) and (13) that
ρ ( y 1 , x 1 ) δ + δ 1 / 4 < δ 1 / 2 .
Suppose that an integer k { 1 , , L } { L } and that we have already defined a trajectory { y i } i = 1 k such that
y 0 = x 0
and that for integers i = 1 , , k ,
ρ ( x i , y i ) < δ i / 2 .
(In view of (6), (12) and (14), our assumption is valid for k = 1 ). Equation (5) implies that there is a point
( ξ k , ξ k + 1 ) graph ( a )
for which
ρ 1 ( ( x k , x k + 1 ) , ( ξ k , ξ k + 1 ) ) δ .
This implies that
ρ ( x k , ξ k ) δ , ρ ( x k + 1 , ξ k + 1 ) δ .
Equations (2), (4), (15) and (17) imply that
ρ ( ξ k , y k ) ρ ( ξ k , x k ) + ρ ( x k , y k ) δ + δ k / 2 δ k / 8 + δ k / 2 .
By (3) and (18),
H ( a ( ξ k ) , a ( y k ) ) δ k + 1 / 8 .
Equations (16) and (19) imply that
ρ ( ξ k + 1 , a ( y k ) ) δ k + 1 / 8 .
In view of (20), there exists
y k + 1 a ( y k )
such that
ρ ( ξ k + 1 , y k + 1 ) δ k + 1 / 4 .
It follows from (2), (4), (17) and (22) that
ρ ( y k + 1 , x k + 1 ) ρ ( y k + 1 , ξ k + 1 ) + ρ ( ξ k + 1 , x k + 1 ) δ + δ k + 1 / 4 < δ k + 1 / 2 .
Thus the assumption made for k is also true for k + 1 . Therefore by induction we constructed the trajectory { y t } t = 0 L X such that
y 0 = x 0 ,
ρ ( x t , y t ) ϵ , t = 0 , , L .
This completes the proof of Lemma. □
Proof of Theorem 3. 
Theorem 2 implies that there is a positive integer L 0 for which the following property is valid:
(i) for every integer T > L 0 and every finite trajectory { x t } t = 0 T ,
Card ( { t { 0 , , T } : ρ ( x t , Ω ( a ) ) > ϵ / 2 } ) L 0 .
Let L > L 0 be an integer. By Lemma 1, there is a positive number δ such that the following property is valid:
(ii) for every sequence { x t } t = 0 L X , which satisfies
ρ 1 ( ( x t , x t + 1 ) , graph h ( a ) ) δ , t = 0 , , L 1
there is a finite trajectory { y t } t = 0 L X for which
y 0 = x 0 ,
ρ ( x t , y t ) ϵ / 4 , t = 0 , , L .
Assume that { x t } t = 0 L X satisfies (23). Property (ii) and (23) imply that there exists a finite trajectory { y t } t = 0 L satisfying (24). Property (i) implies that
C a r d ( { t { 0 , , L } : ρ ( y t , Ω ( a ) ) > ϵ / 2 } ) L 0 .
In view of (24), for integers t = 0 , , L ,
ρ ( x t , Ω ( a ) ) ρ ( x t , y t ) + ρ ( y t , Ω ( a ) ) ϵ / 4 + ρ ( y t , Ω ( a ) )
and if
ρ ( x t , Ω ( a ) ) > ϵ ,
then
ρ ( y t , Ω ( a ) ) > ϵ / 2 .
Together with (25) this implies that
C a r d ( { t { 0 , , L } : ρ ( x t , Ω ( a ) ) > ϵ } )
C a r d ( { t { 0 , , L } : ρ ( y t , Ω ( a ) ) > ϵ / 2 } ) L 0 .
This completes the proof of Theorem 3. □
Proof of Theorem 4. 
We may assume without loss of generality that ϵ < 1 . Theorem 3 implies that there is a positive integer L 0 for which the following property is valid:
(i) for every integer L > L 0 there is a positive number δ such that for every sequence { x t } t = 0 L satisfying
ρ 1 ( ( x t , x t + 1 ) , graph ( a ) ) δ , t = 0 , , T 1
the relation
Card ( { t { 0 , , L } : ρ ( x t , Ω ( a ) ) > ϵ } ) L 0
is true.
Fix a natural number
L > 4 L 0 ϵ 1 .
Let a positive number δ be as guaranteed by property (i). Choose a natural number
k 0 > 4 ϵ 1 .
Set
L 1 = k 0 L .
Assume that T > L 1 is an integer and that a sequence { x t } t = 0 T satisfies
ρ 1 ( ( x t , x t + 1 ) , graph ( a ) ) δ , t = 0 , , T 1 .
There is an integer k 1 1 for which
k 1 L T < ( k 1 + 1 ) L .
By (28) and (30),
k 1 > k 0 .
Property (i), (29) and (30) imply that for integers i = 0 , , k 1 1 ,
Card ( { t { i L , , ( i + 1 ) L } : ρ ( x t , Ω ( a ) ) > ϵ } ) L 0 .
Combined with (30) the equation above implies that
Card ( { t { 0 , , T } : ρ ( x t , Ω ( a ) ) > ϵ } ) k 1 L 0 + L .
Combined with (26), (27), (30) and (31) the equation above implies that
T 1 Card ( { t { 0 , , T } : ρ ( x t , Ω ( a ) ) > ϵ } )
k 1 L 0 T 1 + L T 1 L 0 L 1 + k 0 1 < ϵ .
Theorem 4 is proved. □

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Rubinov, A.M. Turnpike sets in discrete disperse dynamical systems. Sib. Math. J. 1980, 21, 136–146. [Google Scholar]
  2. Rubinov, A.M. Multivalued Mappings and Their Applications in Economic Mathematical Problems; Nauka: Leningrad, Russia, 1980. [Google Scholar]
  3. Dzalilov, Z.; Zaslavski, A.J. Global attractors for discrete disperse dynamical systems. J. Nonlinear Convex Anal. 2009, 10, 191–198. [Google Scholar]
  4. Zaslavski, A.J. Turnpike sets of continuous transformations in compact metric spaces. Sib. Math. J. 1982, 23, 136–146. [Google Scholar] [CrossRef]
  5. Zaslavski, A.J. Uniform convergence to global attractors for discrete disperse dynamical systems. Nonlinear Dyn. Syst. Theory 2007, 4, 315–325. [Google Scholar]
  6. Zaslavski, A.J. Convergence of trajectories of discrete dispersive dynamical systems. Commun. Math. Anal. 2008, 4, 10–19. [Google Scholar]
  7. Zaslavski, A.J. Structure of trajectories of discrete dispersive dynamical systems. Commun. Math. Anal. 2009, 6, 1–9. [Google Scholar]
  8. VMakarov, V.L.; Rubinov, A.M. Mathematical theory of economic dynamics and equilibria; Springer: New York, NY, USA, 1977. [Google Scholar]
  9. Zaslavski, A.J. Turnpike Properties in the Calculus of Variations and Optimal Control; Springer Optimization and Its Applications: New York, NY, USA, 2006. [Google Scholar]
  10. Zaslavski, A.J. A turnpike property of trajectories of dynamical systems with a Lyapunov function. Games 2020, 11, 63. [Google Scholar] [CrossRef]
  11. Samuelson, P.A. A catenary turnpike theorem involving consumption and the golden rule. Am. Econ. Rev. 1965, 55, 486–496. [Google Scholar]
  12. Damm, T.; Grüne, L.; Stieler, M.; Worthmann, K. An Exponential Turnpike Theorem for Dissipative Discrete Time Optimal Control Problems. SIAM J. Control. Optim. 2014, 52, 1935–1957. [Google Scholar] [CrossRef] [Green Version]
  13. Gugat, M. A turnpike result for convex hyperbolic optimal boundary control problems. Pure Appl. Funct. Anal. 2019, 4, 849–866. [Google Scholar]
  14. Gugat, M.; Trélat, E.; Zuazua, E. Optimal Neumann control for the 1D wave equation: Finite horizon, infinite horizon, boundary tracking terms and the turnpike property. Syst. Control. Lett. 2016, 90, 61–70. [Google Scholar] [CrossRef] [Green Version]
  15. Khan, M.A.; Zaslavski, A.J. On two classical turnpike results for the Robinson-Solow-Srinivisan (RSS) model. Adv. Math. Econ. 2010, 13, 47–97. [Google Scholar]
  16. Mammadov, M.A. Turnpike Theorem for an Infinite Horizon Optimal Control Problem with Time Delay. SIAM J. Control. Optim. 2014, 52, 420–438. [Google Scholar] [CrossRef]
  17. Porretta, A.; Zuazua, E. Long Time versus Steady State Optimal Control. SIAM J. Control. Optim. 2013, 51, 4242–4273. [Google Scholar] [CrossRef]
  18. Trelat, E.; Zhang, C.; Zuazua, E. Optimal shape design for 2D heat equations in large time. Pure Appl. Funct. Anal. 2018, 3, 255–269. [Google Scholar]
  19. Trélat, E.; Zuazua, E. The turnpike property in finite-dimensional nonlinear optimal control. J. Differ. Equ. 2015, 258, 81–114. [Google Scholar] [CrossRef]
  20. Zaslavski, A.J. Turnpike Theory of Continuous-Time Linear Optimal Control Problems; Springer Optimization and Its Applications; Springer: Cham, Switzerland, 2015. [Google Scholar]
  21. Aseev, S.M.; Krastanov, M.I.; Veliov, V.M. Optimality conditions for discrete-time optimal control on infinite horizon. Pure Appl. Funct. Anal. 2017, 2, 395–409. [Google Scholar]
  22. Bachir, M.; Blot, J. Infinite Dimensional Infinite-horizon Pontryagin Principles for Discrete-time Problems. Set-Valued Var. Anal. 2014, 23, 43–54. [Google Scholar] [CrossRef]
  23. Bachir, M.; Blot, J. Infinite dimensional multipliers and Pontryagin principles for discrete-time problems. Pure Appl. Funct. Anal. 2017, 2, 411–426. [Google Scholar]
  24. Blot, J.; Hayek, N. Infinite-Horizon Optimal Control in the Discrete-Time Framework; Springer Briefs in Optimization: New York, NY, USA, 2014. [Google Scholar]
  25. Carlson, D.A.; Haurie, A.; Leizarowitz, A. Infinite Horizon Optimal Control; Springer: Berlin, Germany, 1991. [Google Scholar]
  26. Gaitsgory, V.; Grüne, L.; Thatcher, N. Stabilization with discounted optimal control. Syst. Control. Lett. 2015, 82, 91–98. [Google Scholar] [CrossRef] [Green Version]
  27. Gaitsgory, V.; Mammadov, M.; Manic, L. On stability under perturbations of long-run average optimal control problems. Pure Appl. Funct. Anal. 2017, 2, 461–476. [Google Scholar]
  28. Sagara, N. Recursive variational problems in nonreflexive Banach spaces with an infinite horizon: An existence result. Discret. Contin. Dyn. Syst. S 2018, 11, 1219–1232. [Google Scholar] [CrossRef] [Green Version]
  29. Bielecki, A.; Ombach, J. Dynamical Properties of a Perceptron Learning Process: Structural Stability Under Numerics and Shadowing. J. Nonlinear Sci. 2011, 21, 579–593. [Google Scholar] [CrossRef] [Green Version]
  30. Pilyugin, S.Y. Shadowing in Dynamical Systems; Springer: Berlin, Germany, 1999. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zaslavski, A.J. A Weak Turnpike Property for Perturbed Dynamical Systems with a Lyapunov Function. Axioms 2021, 10, 45. https://doi.org/10.3390/axioms10020045

AMA Style

Zaslavski AJ. A Weak Turnpike Property for Perturbed Dynamical Systems with a Lyapunov Function. Axioms. 2021; 10(2):45. https://doi.org/10.3390/axioms10020045

Chicago/Turabian Style

Zaslavski, Alexander J. 2021. "A Weak Turnpike Property for Perturbed Dynamical Systems with a Lyapunov Function" Axioms 10, no. 2: 45. https://doi.org/10.3390/axioms10020045

APA Style

Zaslavski, A. J. (2021). A Weak Turnpike Property for Perturbed Dynamical Systems with a Lyapunov Function. Axioms, 10(2), 45. https://doi.org/10.3390/axioms10020045

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop