Next Article in Journal
Quadrature Based Neural Network Learning of Stochastic Hamiltonian Systems
Previous Article in Journal
Pairs of Positive Solutions for a Carrier p(x)-Laplacian Type Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fold-Fold Singularity in a Piecewise Smooth Mathematical Model Describing the Dynamics of a Stockless Market

by
Oscar Emilio Molina-Díaz
1,2,
Gerard Olivar-Tost
3,* and
Deissy Milena Sotelo-Castelblanco
2
1
Faculty of Basic Sciences and Technologies, Universidad del Quindío, Armenia 630004, Colombia
2
Department of Mathematics and Statistics, Universidad Nacional de Colombia, Manizales 170004, Colombia
3
Department of Natural Sciences and Technology, Universidad de Aysén, Coyhaique 5950000, Chile
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(16), 2442; https://doi.org/10.3390/math12162442
Submission received: 24 June 2024 / Revised: 19 July 2024 / Accepted: 1 August 2024 / Published: 6 August 2024
(This article belongs to the Topic Advances in Nonlinear Dynamics: Methods and Applications)

Abstract

:
Fold-fold singularities are critical points or singularities in piecewise smooth dynamical systems (PWS) where both the stability and the structure of the system change. These singularities are of great importance in the study of specific dynamics, such as those in markets, as they indicate significant transformations in their evolution, including sudden variability in prices or changes in the behavior of offers and demand. Despite the substantial increase in the use of mathematical and computational tools applied to market dynamics, the current literature does not thoroughly address the study of the existence of fold-fold singularities in piecewise smooth systems within this context. Therefore, due to the importance of markets as economic activities, this paper proves the existence of such a singularity in a mathematical model that describes the dynamics of a stockless market, which is represented by a system of ordinary differential equations defined with piecewise smooth functions.

1. Introduction

The classical theory of dynamical systems, described by smooth differential equations, has developed in both a geometrical and qualitative sense. It has been used to analyze the dynamical behavior of different applied models. Conversely, many practical situations are defined by systems of equations that exhibit discontinuities in certain regions or surfaces. In these scenarios, the theoretical basis for the most important results remains to be fully explored.
Lately, the study of systems of piecewise smooth differential equations, which Filippov has called, in [1], discontinuous righthand sides, has increased considerably due to their applicability in different phenomena presented in science, engineering, economics, medicine, and other fields; see [2,3,4,5,6,7]. These systems are characterized by having surfaces (in the case of R 3 ) or hypersurfaces (in higher-dimensional spaces) that divide the state space, where discontinuities of the functions that define the system of differential equations occur, thus determining different vector fields in the different regions. Depending on how the system is defined, these surfaces, or so-called switching manifolds, may have crossing or sliding zones [8].
Sliding zones can be characterized by stability (attracting) or instability (repulsing or escaping). A vector field can be defined on a switching surface as a convex combination of the two vector fields adjacent to such a switching surface, which is called a Filippov vector field. The solutions are known as sliding segments, whose evolution is subject to the switching surface; see [8]. It should be clarified that this definition, called the Filippov convention, is the one implemented in this work. At this point, we highlight the efforts of some authors, such as [9], who, in their research, have defined this vector field in a more general way, including non-linear terms that generalize the definition of sliding dynamics, leading to what they have called hidden dynamics.
The aim of this paper is to prove the existence of a fold-fold singularity in a mathematical model defined by a system of piecewise smooth differential equations, which describes the dynamics of a stockless market. A fold-fold singularity is a point that satisfies Definitions 4 and 7 in the following sections. We highlight the importance of this study, as bibliographic searches in various databases and academic repositories reveal that there are few papers that combine the theory of piecewise smooth systems, fold-fold singularities, and market models. The lack of literature on this subject creates a broad field of study and an opportunity for future research to explore more deeply the implications of the existence of these singularities in market dynamics and their possible consequences for financial decision-making.
In order to ensure that this document is self-contained and to facilitate its reading, the rest of Section 1 provides a brief description of the theory framing piecewise smooth dynamical systems and some singularities. Here, we reproduce some basic definitions and concepts that can be accessed in [1,8,10,11,12,13,14,15]. In Section 2, we present the basic theory of piecewise smooth systems and tangencies. In Section 3, we show a mathematical model defined through a system of piecewise smooth differential equations that describes the dynamics of a stockless market. The main part of this paper, which is presented in Section 4, aims to prove the existence of a fold-fold singularity in the aforementioned model. To complement the study, Section 5 illustrates some numerical simulations that verify the analytical results obtained, which are of great relevance since the literature is scarce regarding the identification of such singularities or tangencies in everyday dynamics. Finally, Section 6 and Section 7 present the discussion and conclusions, respectively.

2. Review of Piecewise Smooth Systems

In this paper, we analyze, through the lens of the theory of piecewise smooth systems, the singularities presented in a mathematical model describing the dynamics of a stockless market. In particular, we focus on the definitions of fold-fold singularities and classify tangencies.
Definition 1.
In general, a piecewise smooth system (PWS) is defined as a finite set of ordinary differential equations
x ˙ = f i ( x ) ; x S i R n
where the vector fields f i are smooth and defined on the different disjoint open sets S i , i = 1 , 2 , 3 , , n [15].
The sets S i are defined in terms of smooth scalar fields H ( i ) ( x ) , with non-zero gradient H x ( i ) ( x ) on the i-th surface, for i = 1 , 2 , , n 1 , as follows:
S i = { x R n : H ( i ) ( x ) < 0 } S i + 1 = { x R n : H ( i ) ( x ) > 0 } ,
which are separated by the Σ i surfaces, i = 1 , 2 , , n 1 , called switching manifolds or switching boundaries. These are ( n 1 ) -dimensional and are defined by the 0 level set of the smooth scalar fields H ( i ) ( x ) . In particular,
Σ i = { x R n : H ( i ) ( x ) = 0 } .
Definition 2.
An equilibrium of the system x ˙ = f i ( x ) is a point x where
f i ( x ) = 0 .
In order to analyze the geometry of the solutions on or near the switching surfaces Σ i , it is necessary to introduce the following notion of derivatives. This allows us, among other aspects, to detect tangencies between solutions and boundaries. The directional derivative of H ( i ) with respect to the vector field f i is defined in terms of the Lie derivative, denoted by L f i , and through the inner product
L f i H ( i ) ( x ) = H x ( i ) ( x ) , f i ( x ) .
The m-th order Lie derivative, as noted by [13], is denoted by L f i m H ( i ) ( x ) . Specifically,
L f i 2 H ( i ) ( x ) = L f i ( L f i H ( i ) ( x ) ) = H x ( i ) ( x ) , f i ( x ) , f i ( x ) ,
where ( · ) represents the gradient.
By applying the Lie derivative, we can identify the regions or points on the switching surfaces Σ i where crossing or sliding occurs. To do this we, define
σ ( i ) ( x ) = ( L f i H ( i ) ) ( L f i + 1 H ( i ) ) ; i = 1 , 2 , , n 1 .
Now, we have the following.
  • We have crossing at the points on Σ i where σ ( i ) ( x ) > 0 . This set is referred to as the set of crossing points and is represented by Σ i c .
  • We have sliding at the points on Σ i where σ ( i ) ( x ) < 0 . This set is called the set of sliding points and is denoted by Σ i s .
As mentioned earlier, for sliding to occur, the following conditions must be met: σ ( i ) ( x ) < 0 , meaning that the expressions L f i H ( i ) and L f i + 1 H ( i ) must have opposite signs. Therefore, the following cases are considered.
  • L f i + 1 H ( i ) < 0 < L f i H ( i ) . This implies that the vector fields f i and f i + 1 are oriented towards the surface at a point x Σ i . These points constitute the set of points of stable sliding and are denoted as Σ i s s .
  • L f i + 1 H ( i ) > 0 > L f i H ( i ) . This implies that the vector fields f i and f i + 1 point away from the switching surface at a point x Σ i ; see [13]. These points form the set of unstable sliding or escape points, which is denoted by Σ i u s .
It is worth noting that Σ i s = Σ i s s Σ i u s .
For x Σ i s , it is defined a vector field g i ( x ) , called a Filippov vector field, using the convex Filippov method [15]. This vector field reflects dynamics through
g i ( x ) = λ i ( x ) f i ( x ) + ( 1 λ i ( x ) ) f i + 1 ( x ) ,
where
λ i ( x ) = H x ( i ) ( x ) , f i + 1 ( x ) H x ( i ) ( x ) , f i + 1 ( x ) f i ( x ) .
Thus, in Σ i s , we have the Filippov system
x ˙ = g i ( x ) .
As can be found in [15], points x in Σ i s , where
L f i + 1 f i H ( i ) ( x ) = H x ( i ) ( x ) , f i + 1 ( x ) f i ( x ) = 0 ; f i , i + 1 ( x ) 0 ,
are called singular sliding points. At such points, either both vectors f i ( x ) and f i + 1 ( x ) are tangent to Σ i , one of them vanishes while the other is tangent to Σ i , or they both vanish.
It is important to note that some authors have defined this vector field in a more general manner, including non-linear terms that generalize the definition of sliding dynamics [9,16].
Moreover, any smooth solution x ( t ) completely contained in one of the regions S i , defined for any t I , with I R , will be referred to as a segment. Segments confined to the switching surfaces and determined by the Filippov vector fields will be called sliding segments; the others will be referred to as non-sliding segments. Furthermore, an orbit is a continuous concatenation of sliding and non-sliding segments. See [13].
Finally, we have the following definition, which can be viewed in [11].
Definition 3.
A pseudo-equilibrium is an equilibrium of the Filippov vector field, which lies on the switching surface where the Filippov vector field has been defined.

2.1. Tangencies

The concept of tangencies is fundamental to the understanding of the dynamics in a switching manifold. They form the boundaries between the crossing, sliding, and escape regions on the Σ i surfaces. Generally, these boundaries can occur in two ways: either where the Σ i surfaces are non-smooth, which will be called a boundary intersection, or where the Σ i surfaces are smooth but one of the adjacent vector fields is tangent, satisfying the tangency condition
L f i H ( i ) = 0 .
The three simplest types of tangencies that can be observed in the smooth regions of the Σ i are a fold or quadratic tangency, a cusp or cubic tangency, and a fold-fold or double tangency. The latter is detailed below. According to [11], the fold-fold tangency or double tangency occurs when there is quadratic contact on both sides of the switching surface. More specifically, it is defined as follows.
Definition 4.
A point x Σ i is called a tangency of the fold-fold type or a double tangency if it simultaneously satisfies the following conditions:
  • L f i H ( i ) = L f i + 1 H ( i ) = 0 ,
  • L f i 2 H ( i ) 0 ; L f i + 1 2 H ( i ) 0 ;
  • the gradient vectors of H ( i ) , L f i H ( i ) , and L f i + 1 H ( i ) are linearly independent.
This tangency is of the following type:
  • Visible-visible, if L f i 2 H ( i ) > 0 > L f i + 1 2 H ( i ) ;
  • Invisible-invisible, if L f i 2 H ( i ) < 0 < L f i + 1 2 H ( i ) ;
  • Visible-invisible, if ( L f i 2 H ( i ) ) ( L f i + 1 2 H ( i ) ) > 0 .
A tangency is called visible or invisible depending on whether the orbits curve locally away from or towards the surface. This is also true for sliding segments that are tangent to the boundary of the sliding region [11]. Furthermore, it follows that if the fold-fold tangency is invisible-invisible, then it is called a T-singularity of Teixeira or T-singularity; see [17].

2.2. Bifurcations

The next definition that can be found in [13] defines the topological equivalence between piecewise smooth systems, a fundamental concept for the understanding of bifurcations.
Definition 5.
Two piecewise smooth systems are said to be topologically equivalent if there is a homeomorphism that maps the maximal segments of one system into the maximal segments of the other, preserving the direction in time and whether a segment is in the sliding zone, not sliding, or not in contact with the switching manifolds Σ i .
Based on the above, by [11], a bifurcation is defined.
Definition 6.
A bifurcation occurs when arbitrary small perturbations generate a topologically non-equivalent system. The bifurcation is said to be discontinuity-induced if it affects the phase portrait in more than one region or in the manifolds Σ i .
The following theorem, which constitutes the principal focus of study in [13], defines a sliding bifurcation and is presented below.
Theorem 1.
A generic uniparametric sliding bifurcation in R n , for any n 3 , may occur in a fold, cusp, or fold-fold tangency.
It is important to note that the objective of this paper is to prove the existence of a fold-fold singularity, but it will not study bifurcations in the model presented below. Although, according to the aforementioned theorem, which is referred to as Theorem 1.4 in [13], a uniparametric sliding bifurcation occurs at a fold-fold tangency in R n , for any n 3 .

3. Application to a System Describing Market Dynamics

The following model is analyzed, which describes the dynamics of a stockless market through a system of piecewise smooth differential equations. The state variables considered are the offer, the developing production capacity, and the demand at time t, which are denoted, respectively, by x ( t ) , y ( t ) , and z ( t ) . This model was proposed and analyzed by the authors in [18]. Here, we present a summary of the model, which is defined by the system
x ˙ = q y r x y ˙ = B ( x , z ) q y z ˙ = k · 1 B ( x , z ) + ξ I m · z .
In the aforementioned system, the function B ( x , z ) represents the investment function, which translates into the capacity to produce. This function is piecewise defined as follows:
B ( x , z ) = 0 , R ( x , z ) 0 I m , 0 < R ( x , z ) δ I h , R ( x , z ) > δ
where δ represents a positive threshold at ( 0 , 1 ) that defines the decision to invest. Hereafter, medium and high investment are denoted by I m , I h , respectively. Furthermore, the return on investment, denoted by R ( x , z ) , is defined in [19] as the difference between the profit obtained and the fixed costs of production, divided by the fixed costs of production. Therefore, the return on investment is expressed as follows:
R ( x , z ) = ( P ( x , z ) C v ) f p C f C f ,
where C v represents the costs and f p the production factor, which is interpreted as the ratio between the actual quantity of items produced during a period of time and the quantity of items that would have been produced if the same period of time had been fully worked. C f represents the operational costs. The price P ( x , z ) is defined by
P ( x , z ) = P max · ( 1 M ( x , z ) ) + 2 · P min · M ( x , z ) 1 + M ( x , z ) ,
with P max , P min as the maximum and minimum price, respectively, and M ( x , z ) representing the reserve margin. If we consider x z > 0 , then the expression for M ( x , z ) is
M ( x , z ) = x z x + m 2 ,
where m 2 is a positive constant representing the minimum offer at the beginning of the market. Therefore, it is possible to interpret x as the actual offer and the term x + m 2 as the effective offer.
The remaining parameters of the system (4) are the developing production rate, represented by q; the consumption rate, denoted by r; the minimum investment required to start the market, designated by ξ ; and the demand growth rate, represented by k.
Depending on how the system (4) is defined, it can be presented as follows.
  • Let S 1 be the set of points in R 3 for which null inversion is defined and satisfy
    z a 1 x + b 1 ,
    with
    a 1 = 2 C P min C + P max 2 · P min ; b 1 = m 2 C P max C + P max 2 · P min ,
    and
    C = C f + f p C v f p .
  • A medium investment is defined in the set of points in R 3 , denoted S 2 , with
    a 1 x + b 1 < z a 2 x + b 2 ,
    with a 1 , b 1 as defined in (8) and
    a 2 = 2 C δ P min C δ + P max 2 · P min ; b 2 = m 2 C δ P max C δ + P max 2 · P min .
    Furthermore,
    C δ = ( δ + 1 ) C f + f p C v f p .
  • The points of R 3 that exhibit a high degree of inversion are denoted as S 3 and are defined by
    z > a 2 x + b 2 .
In summary, the state space is divided into three areas denoted by S 1 , S 2 , and S 3 , which are conditioned by null, medium, and high inversion, respectively. These areas are defined by
S 1 = x R 3 : z a 1 x + b 1 , S 2 = x R 3 : a 1 x + b 1 < z a 2 x + b 2 , S 3 = x R 3 : z > a 2 x + b 2 .
These regions are separated by two surfaces or switching planes:
Σ 1 = x R 3 : z = a 1 x + b 1 , Σ 2 = x R 3 : z = a 2 x + b 2 .
In this case, we have the smooth scalar fields H ( i ) : R 3 R , which are defined by
H ( i ) ( x , y , z ) = a i x + b i z , i = 1 , 2 .
As x = ( x , y , z ) , the gradient of H ( i ) ( x ) is
H x ( i ) ( x ) = a i , 0 , 1 .
The above is evident in Figure 1, where the non-negative state variables in the first octant have been plotted. The blue plane represents the switching surface Σ 1 , the yellow plane represents Σ 2 , and the green plane is z = x .
Thus, for x = ( x , y , z ) R 3 , the system (4) is a piecewise smooth system, defined in this case as follows:
x ˙ = f 1 ( x ) ; x S 1 f 2 ( x ) ; x S 2 f 3 ( x ) ; x S 3 ,
where the vector fields f 1 , f 2 , and f 3 differ by the value taken by the function B ( x , z ) { 0 , I m , I h } in the system given by Equation (4).
The behavior of non-smooth system orbits in the regions approaching the Σ 1 and Σ 2 switching planes is studied using Filippov’s method, as detailed in [15].
For x Σ 1 , Σ 2 , the following definitions are provided:
σ ( 1 ) ( x ) = ( L f 1 H ( 1 ) ( x ) ) ( L f 2 H ( 1 ) ( x ) ) σ ( 2 ) ( x ) = ( L f 2 H ( 2 ) ( x ) ) ( L f 3 H ( 2 ) ( x ) )
Recall that, as described in the previous section, there is crossing at the points x Σ i where σ ( i ) ( x ) > 0 . Otherwise, there is sliding. Figure 2 illustrates this for each switching surface.
On the set of sliding points of each of the switching planes Σ 1 s and Σ 2 s , the following Filippov vector fields are, respectively, defined:
g 1 ( x ) = x ˙ = q y r x I m ( 1 λ 1 ( x ) ) q y a 1 ( q y r x ) ; g 2 ( x ) = x ˙ = q y r x λ 2 ( x ) ( I m I h ) + I h q y a 2 ( q y r x )
with
λ 1 ( x ) = a 1 ( q y r x ) k ( a 1 x + b 1 ) + ξ I m λ 2 ( x ) = a 2 I m ( q y r x ) k ( a 2 x + b 2 ) ( I h I m ) + ξ I h I m + 1 .
The possible pseudo-equibria for g 1 ( x ) and g 2 ( x ) are, respectively,
x ˜ 1 = ( x , y , z ) = I m ξ r , I m ξ q , a 1 I m ξ r + b 1 x ˜ 2 = ( x , y , z ) = I m ξ r , I m ξ q , a 2 I m ξ r + b 2 .
It can be proven that x ˜ 1 Σ 1 s but that x ˜ 2 Σ 2 s . Therefore, we will henceforth only analyze the dynamics on the switching plane Σ 1 . The stability of the pseudo-equilibrium x ˜ 1 was analyzed in [18].
In addition to pseudo-equilibria, in sliding zones, we can also find other points that play an important role in the analysis of smooth systems by sections. One of these points is defined below.
Definition 7.
The singular sliding points satisfy, for x Σ 1 s ,
L f 2 f 1 H ( 1 ) ( x ) = H x ( 1 ) ( x ) , f 2 ( x ) f 1 ( x ) = 0 ; f 1 , 2 ( x ) 0 .
In our case, the singular sliding point, for x Σ 1 s , is
x ^ 1 = ( x , y , z ) = b 1 a 1 , b 1 r a 1 q , 0 .
Geometrically, this point is at the intersection of the stable sliding zone designated as Σ 1 s s , and the unstable sliding zone, denoted by Σ 1 u s , and lies on the x y -plane. The following section will demonstrate that this point is a fold-fold tangency.
It should be noted that other authors have described the dynamics of a market through piecewise smooth systems, but these systems, due to the nature of the functions that define them, do not present sliding regions, nor do they study the characteristics of the tangencies [3,20]. However, in the work presented here, as has just been mentioned and analyzed, there are sliding regions. We now proceed to study and classify the tangencies or singularities.

4. Test of the Double Tangency Conditions

In this section, we will demonstrate that the singular sliding point defined in (11) satisfies the conditions to be classified as a fold-fold tangency. Recall that three conditions must be satisfied: L f 1 H ( 1 ) ( x ^ 1 ) = L f 2 H ( 1 ) ( x ^ 1 ) = 0 , L f 1 2 H ( 1 ) ( x ^ 1 ) 0 , and the gradient vectors of H ( 1 ) ( x ^ 1 ) , L f 1 H ( 1 ) ( x ^ 1 ) and L f 2 H ( 1 ) ( x ^ 1 ) are linearly independent.
Proposition 1.
In the dynamical system (4), the singular sliding point x ^ 1 given in (11) is a fold-fold tangency.
Proof. 
We shall prove that the three conditions mentioned above are satisfied.
  • Let us prove that L f 1 H ( 1 ) ( x ^ 1 ) = L f 2 H ( 1 ) ( x ^ 1 ) = 0 . It is henceforth considered x = ( x , y , z ) .
    L f 1 H ( 1 ) ( x ) = H x ( 1 ) ( x ) , f 1 ( x ) = ( a 1 , 0 , 1 ) · q y r x , q y , k z 1 ξ I m = a 1 ( q y r x ) k z 1 ξ I m .
    Replacing x ^ 1 in the above expression, it can be verified that L f 1 H ( 1 ) ( x ^ 1 ) = 0 .
    Similarly,
    L f 2 H ( 1 ) ( x ) = H x ( 1 ) ( x ) , f 2 ( x ) = ( a 1 , 0 , 1 ) · q y r x , I m q y , k z 1 I m + ξ I m = a 1 ( q y r x ) k z 1 ξ I m .
    Replacing x ^ 1 in the above expression, it can be verified that L f 2 H ( 1 ) ( x ^ 1 ) = 0 .
    For the time being, the singular sliding point x ^ 1 = ( x , y , z ) = b 1 a 1 , b 1 r a 1 q , 0 satisfies the first condition, and, as mentioned above, this point connects the stable sliding zone with the unstable one, as can be seen in Figure 3.
  • Let us now prove that L f 1 2 H ( 1 ) ( x ^ 1 ) = H x ( 1 ) ( x ^ 1 ) , f 1 ( x ^ 1 ) , f 1 ( x ^ ) 0
    L f 1 2 H ( 1 ) ( x ) = a 1 r , a 1 q , k 1 ξ I m · q y r x , q y , k z 1 ξ I m = a 1 r ( q y r x ) a 1 q 2 y k 2 z 1 ξ I m 2
    Replacing x ^ 1 in the last expression,
    a 1 r q b 1 r a 1 q + a 1 r b 1 a 1 a 1 q 2 b 1 r a 1 q = b 1 r q 0
    This is fulfilled, since b 1 , r , q 0 . Moreover, as b 1 < 0 , then L f 1 2 H ( 1 ) ( x ) < 0 .
    Similarly for L f 2 2 H ( 1 ) ( x ^ 1 ) = H x ( 1 ) ( x ^ 1 ) , f 2 ( x ^ 1 ) , f 2 ( x ^ 1 ) 0 , we see
    L f 2 2 H ( 1 ) ( x ) = a 1 r , a 1 q , k 1 I m + ξ I m · q y r x , I m q y , k z 1 I m + ξ I m = a 1 r ( q y r x ) + a 1 q ( I m q y ) k 2 z 1 I m + ξ I m 2
    In the same way, replacing x ^ 1 in the latter expression, we obtain
    a 1 r q b 1 r a 1 q + a 1 r 2 b 1 a 1 + a 1 q I m a 1 q 2 b 1 r a 1 q = a 1 I m q + b 1 r q 0 ,
    provided that b 1 a 1 I m r . Moreover, L f 2 2 H ( 1 ) ( x ) > 0 for b 1 a 1 > I m r .
    Under the above assumptions, the tangency is fold-fold invisible-invisible; therefore, it is a Teixeira singularity or Teixeira T-singularity [17].
  • Let us now check that the gradient vectors of the H ( 1 ) , L f 1 H ( 1 ) , and L f 2 H ( 1 ) are linearly independent. Recall that
    H x ( 1 ) ( x ^ 1 ) = ( a 1 , 0 , 1 ) ( L f 1 H ( 1 ) ( x ^ 1 ) ) = a 1 r , a 1 q , k 1 ξ I m ( L f 2 H ( 1 ) ( x ^ 1 ) ) = a 1 r , a 1 q , k 1 I m + ξ I m
    To check their linear independence, we calculate the determinants of these vectors.
    a 1 a 1 r a 1 r 0 a 1 q a 1 q 1 k 1 ξ I m k 1 I m + ξ I m = a 1 2 q k 0 .
    Therefore, the vectors are linearly independent.
Thus, the singular sliding point x ^ 1 = ( x , y , z ) = b 1 a 1 , b 1 r a 1 q , 0 is a fold-fold tangency; in particular, it is a T-singularity. □
As mentioned above, Theorem 1.4 is stated in [13], which indicates that, in particular, at a fold-fold tangency, a uniparametric sliding bifurcation occurs in R n , n 3 .

5. Numerical Simulations

This section presents some numerical scenarios that illustrate the behavior of trajectories near the fold-fold singularity. In order to achieve this, both positive and negative initial conditions are considered, with the latter being taken above and below the singularity. Although these scenarios are not applicable in practice, they are included here to complement the work mathematically. It should first be noted that the green plane in the graphs illustrates the z = 0 plane, the regions shown in cyan illustrate the sliding region in the switching plane Σ 1 , and the point on the sliding region represents the pseudo-equilibrium. In addition, curves of different colors represent trajectories that start from the points that indicate the initial conditions; to start the numerical computations, specific parameter values were assigned, as detailed in Table 1. The value of ξ will be defined in each numerical scenario.
Figure 4 illustrates that, by assuming negative initial conditions as starting points above the pseudo-equilibrium that is in the escape zone ( x ^ 1 Σ 1 u s ), the trajectories for the offer, developing production capacity, and demand transit the spiral-shaped crossover zone and move around the unstable sliding zone of the switching surface Σ 1 and upwards until they reach the fold-fold singularity. Moreover, in Figure 5, negative initial conditions are also considered, but, in this case, they are situated below the pseudo-equilibrium. Here, the trajectories for the offer and developing production capacity tend to zero, whereas that for the demand decreases without bounds.
The objective of Figure 4 and Figure 5 is to illustrate the dynamics of the state variables in the neighborhood of the double tangency, taking negative initial values for the variables. Moreover, depending on the parameter values, it is possible for the pseudo-equilibrium x ^ 1 to move from the stable sliding zone to the unstable sliding zone through the fold-fold tangency for some values of ξ . In the model, this represents the minimum investment required to start the market. In particular, if ξ < I m + b 1 r a 1 then x ^ 1 Σ 1 s s . Otherwise, x ^ 1 Σ 1 u s . We clarify that, according to the parameter values assigned in Table 1 for the numerical simulations, I m + b 1 r a 1 2470.7 .
Now, it is worth considering the scenario depicted in Figure 6. This scenario also assumes negative initial conditions and the pseudo-equilibrium in the stable sliding zone. As mentioned before, this can be in the stable or unstable sliding zone depending on the values of ξ . Under these conditions, the trajectories involve the escape zone through the crossing zone to a point where the offer and developing production capacity tend to zero and the demand decreases without bounds. It is similar to the situation shown in Figure 5.
In Figure 7, we take positive initial conditions and the pseudo-equilibrium in the stable sliding zone Σ 1 s s . The orbits transit through the crossing zone above the plane z = 0 until they collide with the stable sliding zone, after which the trajectories slide, seeking the pseudo-equilibrium, which is locally asymptotically stable, as mentioned before.
Finally, in the scenario shown in Figure 8, positive initial conditions and the pseudo-equilibrium in the unstable sliding zone are taken. The orbits cross the switching surface Σ 1 above the plane z = 0 until they reach the sliding zone and then slide to the fold-fold singularity.
It is noteworthy that, in this paper, based on the work and research carried out by other authors, we propose to make an additional contribution to the practical applications of the theory framing piecewise smooth systems and related topics of singularities. For example, in [12], this theory is analyzed and applied to a canonical switched feedback control model, and the classification of fold-fold singularities and their dynamics is demonstrated. Moreover, in [13], it is shown that, in general, the behavior of a uniparametric sliding bifurcation is as follows: a sliding orbit starting from the stable sliding zone can leave such a set of points either when the orbit reaches the border points of the stable sliding zone or when it passes through the fold-fold and passes into the unstable sliding region; here, the trajectory leaves the set of points by crossing the switching surface in spiral motions or when it reaches the unstable sliding boundary. The latter motion is also described in [14] and is called a non-smooth diabolo. In contrast to the above, in our work, and according to the dynamics of the model analyzed, as the z = 0 plane is invariant, the sliding orbits reach the double tangency and stop, preventing the trajectories from passing into the unstable sliding sector. We highlight that, although the pseudoqeulibrium is in this escape zone, the trajectories stop when sliding to the fold-fold tangency. The fact that the z = 0 plane is invariant makes the application studied in this work a non-generic case.

6. Discussion

In this paper, we have detected and proved the presence of a fold-fold singularity in a mathematical model described by a system of piecewise smooth ordinary differential equations, which describes the dynamics of a stockless market.
This identification of the fold-fold singularity in our system suggests complex and potentially chaotic behavior in the dynamics of a stockless market. In contrast to previous studies on the subject, the model studied offers new insights into how discontinuities and abrupt transitions can influence market stability.
The results in this paper may have significant implications for the design and regulation of stockless markets, highlighting the need for strategies that consider possible singularities and their effects. This study contributes to the field of discontinuous dynamical systems and their application in economics, since there is little literature that combines the mathematical modeling of a market, sliding in the switching surfaces, and this type of singularity at the same time.
One of the main simplifications is the definition of the vector field in the sliding zone using the Filippov convention, which is described as a convex combination of the vector fields adjacent to the switching surface, rather than considering non-linear terms as other authors do, generalizing the Filippov definition of the vector field. This simplification, while useful for the initial analysis, may not fully capture all dynamics of the system.
Future research on this issue should explore the bifurcations in this system, using both analytical and numerical analysis to better understand how the dynamics vary under different conditions. In addition, considering non-linear terms in the vector field defined in sliding zones could provide a more complete picture of the system’s behavior.
In conclusion, this study provides a new understanding of the dynamics of stockless markets, highlighting the importance of fold-fold singularities. Although we have made significant progress, much remains to be explored, especially in the area of bifurcations and their impact on market stability.

7. Conclusions

In this paper, a brief overview of the basic concepts of piecewise smooth systems, tangencies, and the types of tangencies, and a brief description of the conditions required for bifurcations to occur in such systems, is given.
The above topics are applied to a mathematical model describing the dynamics of a stockless market, defined through a system of differential equations with piecewise smooth functions, capturing external factors such as the investment decision. This is made possible by a piecewise smooth function, which divides the state space into three regions depending on the returns obtained. Unlike other models found in the literature, the one analyzed in this paper presents sliding in the regions of the switching surfaces. Therefore, it was possible to analyze tangencies and singularities.
According to the definitions and main results, it is proven that the analyzed model has a fold-fold tangency, which could also be classified as a T-singularity or Teixeira singularity. Such a point is the boundary between the stable sliding region and the escape zone and lies on the x y -plane. The trajectories around this point indicate that, when taking negative initial conditions, below the singularity and above the pseudo-equilibrium, they pass through the crossover zone, enclosing, in a spiral shape, the escaping region, ascending until reaching the singularity. Furthermore, by taking the positive initial conditions and above the singularity, the trajectories travel through the crossover zone until they collide with the stable sliding zone. Subsequently, the trajectories slide in search of the pseudo-equilibrium, when it is in the stable sliding zone, which is locally asymptotically stable. When the pseudo-equilibrium is in the escape zone, the orbits slide through the set of sliding points, descending to the singularity, reflecting that the z = 0 plane is an invariant plane.

Author Contributions

O.E.M.-D., G.O.-T. and D.M.S.-C. participated in the conceptualization, methodology, software, formal analysis, investigation, and writing (original draft preparation). All authors have read and agreed to the published version of the manuscript.

Funding

Oscar Emilio Molina-Díaz and Deissy Milena Sotelo-Castelblanco were supported by Universidad Nacional de Colombia, Manizales-Colombia, grant number 201010035107-4063, code 51210. Gerard Olivar-Tost was supported by the project “Strengthening Research, Innovation and Technology Transfer, through a doctoral program in the fields of Nature, Intelligence, Territory, Education and Health (NITES)”, code 21991, and by the “Plan de Fortalecimiento de las Universidades Estatales 2021”, Ministry of Education of Chile.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

This work was supported by the University of Quindío-Colombia, the National University of Colombia, Manizales campus—specifically, the Doctorate in Mathematical Sciences. It was also supported by the University of Aysén, Chile. Deissy M. Sotelo and Gerard Olivar-Tost thank CRM Barcelona for partially financing their stay while writing this article. Likewise, Gerard Olivar-Tost thanks the SEDA Energía consortium for its partial collaboration.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviation

The following abbreviation is used in this manuscript:
PWSPiecewise smooth system

References

  1. Filippov, A.F.; Arscott, F.M. Differential Equations with Discontinuous Righthand Sides; Kluwer Academic: New York, NY, USA, 1988. [Google Scholar]
  2. Amador, J.; Redondo, J.; Olivar-Tost, G.; Erazo, C. Cooperation-Based Modeling of Sustainable Development: An Approach from Filippov’s Systems. Complexity 2021, 2021, 4249106. [Google Scholar] [CrossRef]
  3. Valencia-Calvo, J.; Olivar-Tost, G.; Morcillo-Bastidas, J.; Franco-Cardona, C.; Dyner, I. Non-Smooth Dynamics in Energy Market Models: A Complex Approximation From System Dynamics and Dynamical Systems Approach. IEEE Access 2020, 8, 128877–128896. [Google Scholar] [CrossRef]
  4. Amador, J.; Granada, D.H.; Redondo, J.; Olivar-Tost, G. Nonlinear and Nonsmooth Dynamics in Stress-Sickness Processes. Cienc. Desarro. 2017, 8, 9–19. [Google Scholar] [CrossRef]
  5. Amador, J.; Granada, D.H.; Olivar-Tost, G. Non-linear Phenomena in Socio-Economics Systems with Natural Resource Management. Blucher Mech. Eng. Proc. 2014, 1, 3313–3325. [Google Scholar]
  6. Hou, J.; Hong, J.; Ji, B.; Dong, B.; Chen, Y.; Ward, M.P.; Zhang, Z. Changing transmission dynamics of COVID-19 in China: A nationwide population-based piecewise mathematical modelling study. Medrxiv 2020, 3. [Google Scholar] [CrossRef]
  7. Asano, T.; Yokoo, M. Chaotic dynamics of a piecewise linear model of credit cycles. J. Math. Econ. 2019, 80, 9–21. [Google Scholar] [CrossRef]
  8. Bernardo, M.; Budd, C.J.; Champneys, A.R.; Kowalczyk, P. Piecewise-Smooth Dynamical Systems: Theory and Applications; Applied Mathematical Sciences; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2008; Volume 163. [Google Scholar]
  9. Novaes, D.N.; Jeffrey, M.R. Regularization of hidden dynamics in piecewise smooth flow. J. Differ. Equ. 2015, 259, 4615–4633. [Google Scholar] [CrossRef]
  10. Bernardo, M.; Hogan, S.J. Discontinuity-induced bifurcations of piecewise smooth dynamical systems. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2010, 368, 4915–4935. [Google Scholar] [CrossRef] [PubMed]
  11. Colombo, A.; Di Bernardo, M.; Hogan, S.J.; Jeffrey, M.R. Bifurcations of piecewise smooth flows: Perspectives, methodologies and open problems. Phys. D Nonlinear Phenom. 2012, 241, 1845–1860. [Google Scholar] [CrossRef]
  12. Jeffrey, M.R. Sliding bifurcations and non-determinism in systems with switching. IFAC Proc. Vol. 2011, 44, 13275–13280. [Google Scholar] [CrossRef]
  13. Jeffrey, M.R.; Hogan, S.J. The geometry of generic sliding bifurcations. Siam Rev. 2011, 53, 505–525. [Google Scholar] [CrossRef]
  14. Jeffrey, M.R.; Colombo, A. The Two-Fold Singularity of Discontinuous Vector Fields. SIAM J. Appl. Dyn. Syst. 2009, 8, 624–640. [Google Scholar] [CrossRef]
  15. Kuznetsov, Y.A.; Rinaldi, S.; Gragnani, A. One-Parameter Bifurcations in Planar Filippov Systems. Int. J. Bifurc. Chaos 2003, 13, 2157–2188. [Google Scholar] [CrossRef]
  16. Jeffrey, M.R. Hidden dynamics in models of discontinuity and switching. Phys. Nonlinear Phenom. 2014, 273–274, 34–45. [Google Scholar] [CrossRef]
  17. Teixeira, M.A. Generic Bifurcation of Sliding Vector Fields. J. Math. Anal. Appl. 1993, 176, 436–457. [Google Scholar] [CrossRef]
  18. Molina-Díaz, O.; Olivar-Tost, G.; Sotelo-Castelblanco, D. Dynamics of a Stockless Market through Piecewise Smooth Models. Math. Comput. Simul. 2024, 226, 246–269. [Google Scholar] [CrossRef]
  19. Chiara, L.P.; Benjamin, H. Market Power in Power Markets: An Analysis of Residual Demand Curves in California’s Day-ahead Energy Market. Energy J. 2015, 36, 191–218. [Google Scholar]
  20. Redondo, J.M.; Ibarra, D.W.; Becerra, M.; Olivar-Tost, G. Making Decisions in National Energy Markets with Bifurcation Analysis. J. Phys. Conf. Ser. 2019, 1414, 012008. [Google Scholar] [CrossRef]
Figure 1. Switching surfaces and regions where there is zero, medium, and high inversion, denoted by S 1 , S 2 , and S 3 , respectively.
Figure 1. Switching surfaces and regions where there is zero, medium, and high inversion, denoted by S 1 , S 2 , and S 3 , respectively.
Mathematics 12 02442 g001
Figure 2. Regions in the Σ 1 and Σ 2 surfaces where crossing or sliding occurs. View from the top. The set of unstable sliding points or escaping is below the z = 0 plane, and the stable sliding zone is above the z = 0 plane.
Figure 2. Regions in the Σ 1 and Σ 2 surfaces where crossing or sliding occurs. View from the top. The set of unstable sliding points or escaping is below the z = 0 plane, and the stable sliding zone is above the z = 0 plane.
Mathematics 12 02442 g002
Figure 3. Set of points on Σ 1 where the conditions are satisfied: L f 1 H ( 1 ) = L f 2 H ( 1 ) = 0 . At the green point, both conditions are satisfied simultaneously.
Figure 3. Set of points on Σ 1 where the conditions are satisfied: L f 1 H ( 1 ) = L f 2 H ( 1 ) = 0 . At the green point, both conditions are satisfied simultaneously.
Mathematics 12 02442 g003
Figure 4. Behavior of the trajectories for offer, developing production capacity, and demand, assuming negative initial conditions above the pseudo-equilibrium and in the unstable sliding zone. Here, ξ = 2480 , implying that x ^ 1 Σ 1 u s .
Figure 4. Behavior of the trajectories for offer, developing production capacity, and demand, assuming negative initial conditions above the pseudo-equilibrium and in the unstable sliding zone. Here, ξ = 2480 , implying that x ^ 1 Σ 1 u s .
Mathematics 12 02442 g004
Figure 5. Behavior of the trajectories for offer, developing production capacity, and demand, taking negative initial conditions—in this case, below the pseudo-equilibrium, which is in the unstable sliding zone. Assuming ξ = 2480 implies that x ^ 1 Σ 1 u s .
Figure 5. Behavior of the trajectories for offer, developing production capacity, and demand, taking negative initial conditions—in this case, below the pseudo-equilibrium, which is in the unstable sliding zone. Assuming ξ = 2480 implies that x ^ 1 Σ 1 u s .
Mathematics 12 02442 g005
Figure 6. Behavior of the trajectories for offer, developing production capacity, and demand, assuming negative initial conditions and the pseudo-equilibrium in the stable sliding zone. Taking ξ = 2400 , i.e., x ^ 1 Σ 1 s s .
Figure 6. Behavior of the trajectories for offer, developing production capacity, and demand, assuming negative initial conditions and the pseudo-equilibrium in the stable sliding zone. Taking ξ = 2400 , i.e., x ^ 1 Σ 1 s s .
Mathematics 12 02442 g006
Figure 7. The behavior of the trajectories for offer, developing production capacity, and demand is analyzed, taking positive initial conditions and the pseudo-equilibrium in the stable sliding zone into account. Here, ξ = 2300 , implying that x ^ 1 Σ 1 s s .
Figure 7. The behavior of the trajectories for offer, developing production capacity, and demand is analyzed, taking positive initial conditions and the pseudo-equilibrium in the stable sliding zone into account. Here, ξ = 2300 , implying that x ^ 1 Σ 1 s s .
Mathematics 12 02442 g007
Figure 8. Behavior of the trajectories for offer, developing production capacity, and demand, taking positive initial conditions and the pseudo-equilibrium in the unstable sliding zone. At ξ = 2515 , then x ^ 1 Σ 1 u s .
Figure 8. Behavior of the trajectories for offer, developing production capacity, and demand, taking positive initial conditions and the pseudo-equilibrium in the unstable sliding zone. At ξ = 2515 , then x ^ 1 Σ 1 u s .
Mathematics 12 02442 g008
Table 1. Model parameter values.
Table 1. Model parameter values.
SymbolValue/UnitName
M 2 1000 U.A.Minimum quantity of articles to start the market
ξ U.A.Minimum investment quantity to start the market
δ 0.5Investment threshold
r1/35 d 1 Consumption rate
q1/5 d 1 Production rate in development
C f 60 M.U.Fixed production costs
C v 150 M.U.Variable costs
f p 0.7Production factor
k0.03 d 1 Demand growth rate
P max 350 U.M.Maximum price
P min 180 U.M.Minimum price
I m 2500 U.A.Medium investment
I h 4200 U.A.High investment
U.A.: units of articles, M.U.: monetary units.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Molina-Díaz, O.E.; Olivar-Tost, G.; Sotelo-Castelblanco, D.M. Fold-Fold Singularity in a Piecewise Smooth Mathematical Model Describing the Dynamics of a Stockless Market. Mathematics 2024, 12, 2442. https://doi.org/10.3390/math12162442

AMA Style

Molina-Díaz OE, Olivar-Tost G, Sotelo-Castelblanco DM. Fold-Fold Singularity in a Piecewise Smooth Mathematical Model Describing the Dynamics of a Stockless Market. Mathematics. 2024; 12(16):2442. https://doi.org/10.3390/math12162442

Chicago/Turabian Style

Molina-Díaz, Oscar Emilio, Gerard Olivar-Tost, and Deissy Milena Sotelo-Castelblanco. 2024. "Fold-Fold Singularity in a Piecewise Smooth Mathematical Model Describing the Dynamics of a Stockless Market" Mathematics 12, no. 16: 2442. https://doi.org/10.3390/math12162442

APA Style

Molina-Díaz, O. E., Olivar-Tost, G., & Sotelo-Castelblanco, D. M. (2024). Fold-Fold Singularity in a Piecewise Smooth Mathematical Model Describing the Dynamics of a Stockless Market. Mathematics, 12(16), 2442. https://doi.org/10.3390/math12162442

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