# Bifurcation Structure in a Bimodal Piecewise Linear Business Cycle Model

Authors: Viktor Avrutin; Iryna Sushko; Fabio Tramontana
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

We study the bifurcation structure of the parameter space of a 1D continuous piecewise linear bimodal map which describes dynamics of a
business cycle model introduced by Day-Shafer. In particular, we obtain the analytical expression of the boundaries of several periodicity regions associated with attracting cycles of the map (principal cycles and related fin structure). By crossing these boundaries the map displays robust chaos.


---

## Body

## 1. Introduction

Applied models defined bypiecewise smooth functions appear quite often when one studies a real process characterized by some “nonsmooth” phenomena such as sharp switching between several states, impacts, friction, sliding, and the like. In economic modeling piecewise smooth systems arise, for example, taking into account that the most used economic variables have nonnegativity constraints, or when a process is studied in which the economic agents change their behavior when a relevant indicator reaches certain thresholds, and so forth. The main reason why economists still prefer to build their theoretical models avoiding piecewise smooth functions is related to the lack of knowledge and experience in the investigation of such a kind of models. In fact, a general theory for piecewise smooth dynamical systems, differently from the one for smooth systems, is not yet well established. The studies and results on these systems are growing and rapidly developing nowadays (see, e.g., the books [1, 2] and references therein). During the last decade important results have been obtained in this field, and one of the aims of the present paper is to show that such results can be successfully applied to investigate a relevant economic model, proposed by Day and Shafer in [3].We recall that among nonsmooth dynamical systems, those described by piecewise linear maps are the simplest to study due to the linearity of their components, but nevertheless they are quite rich in the outcome of the possible dynamics. In particular, a one-dimensional (1D for short) continuous piecewise linear map with one border point, known asskew tent map, depending on the parameters values can have attracting cycles of any period as well as cyclic chaotic intervals of any period, also called n-bands chaotic attractors, which have the relevant property of being robust (as introduced in [4]) with respect to parameter perturbations. The bifurcation structure of the skew tent map has been completely described (see, e.g., [5–8]). Moreover, the skew tent map can be used as a normal form for a so-calledborder collision bifurcation (BCB for short) which is characteristic in piecewise smooth maps [9, 10]. Recall that a BCB occurs when an invariant set, such as, for example, a fixed point or cycle, collides with a border separating regions of different definition of the map. This bifurcation may lead, for example, from an attracting cycle directly to chaos. The dynamic behaviors of the skew tent map are used to classify the possible dynamics which may occur after a generic BCB in a 1D continuous piecewise smooth map (see, e.g., [11–13], where the skew tent map is applied to classify BCBs in economic models).The map considered in the present paper, which represents an economic model, is described by a one-dimensional (1D for short)bimodal piecewise linear map with increasing outermost branches. Clearly, a map with two border points possesses more complicated dynamics, and all the possible outcomes are not yet fully investigated; some results can be found in [14, 15]. In particular, the bifurcation structure of its parameter space includes both regions which belong to the knownperiod adding structure (called alsoArnold tongues ormode-locking tongues), which is characteristic for piecewise increasing discontinuous maps and also for certain circle maps (see, e.g., [16–20]). The period adding structure is formed by periodicity regions related to cycles organized according to the Farey summation rule applied to the rotation numbers of the related cycles. Besides this, in bimodal piecewise linear maps the so-calledfin structure is also observed (see [15]). We will describe these regions in the parameter space of the considered map, recalling how these two structures are organized and giving formulas of the boundaries of related regions.The plan of the work is as follows. In Section2 we recall the Day-Shafer model. Its dynamics are bounded in an absorbing interval, and depending on the parameters values the system may have attracting cycles of any period as well as n-band chaotic attractors. In Section 3 we will consider the parameter space of interest, showing how the periodicity regions representing the regions in the parameter space associated with stable attracting cycles may be organized. Moreover, the boundaries of such regions (BCB curves) are obtained analytically (and reported in Appendices). Crossing these boundaries the system may either enter a different periodicity region (via the fin structure mentioned above) or it may become chaotic (and in a regime of robust chaos). Some conclusions are given in Section 4.


## 2. The Day-Shafer Model

The Day-Shafer model we are interested in dates back to 1987. Richard Day has been a pioneer in the application of nonlinear models in economics and finance (see, e.g., [21–24]). In particular, in [3] Day and Shafer argued that the trapping set of their nonlinear business cycle model in one of the most interesting cases is well approximated by a piecewise smooth map with twoturning points with dynamics bounded in an absorbing interval. Moreover, they explicitly consider a piecewise linear bimodal map as a further approximation. In the particular case in which the dynamics reduces to those on a unimodal piecewise linear map (i.e., a skew tent map), they succeeded in writing the BCB curves associated with the principal (or maximal) cycles. While in the generic case of a piecewise linear bimodal map, only a few numerical results were given. As we will show in the next sections, the results of the generic case with two turning points can be much improved and many BCB curves can be detected analytically.Let us briefly recall that in [3] Day and Shafer build a generic map defining a business cycle model of an economy with monetary and real sectors. In the monetary market, demand and supply for money are implicitly given by the interest rate r as a function of the income Y and a money supply parameter M. In the real market the interest rate determines the level of investments I
=
H
(
r
(
Y
)
,
Y
)
=
I
(
Y
), that is, a component of the income together with consumption C
=
C
(
Y
). By using the assumption that the level of income at a certain time period depends upon the one-time lagged amounts of consumption and investment, a discrete time equation of the following form is obtained:
(1)
Y
t
+
1
=
F
(
Y
t
)
=
λ
I
(
Y
t
)
+
C
(
Y
t
)
+
A
,
where the parameter λ
≥
0 measures how strongly investment translates into new income, and the parameter A
>
0 includes all the autonomous components (of investment, consumption, and public expenditure). The function C
(
Y
) is typically monotonically increasing with C
(
0
)
=
0, and the investment function I
(
Y
) usually is increasing for low level of income while it is decreasing when income is high (this is due to the fact that with high income the money market becomes too crowded, inducing the increase of the interest rate that causes a contraction in the investments). Putting together these assumptions a function F
(
Y
) is obtained that is nonlinear and has a bimodal shape. By using explicit functions, Day and Shafer found that for some parameters’ configurations the income dynamics are bounded in a trapping set like the one shown in Figure 1. The piecewise linear map shown in Figure 1 represents the function F
(
Y
); it is continuous but not differentiable at the points Y
=
Y
* and Y
=
Y
*
*. The dynamics are bounded in the interval F
(
Y
*
*
)
,
F
(
Y
*
)
=
Y
m
,
Y
M. For convenience, we can normalize the interval Y
m
,
Y
M to [
0,1
] and denoting by x the state variable income Y we obtain the map f
:
I
→
I, I
=
[
0,1
] defined as follows:
(2)
f
:
x
⟼
f
(
x
)
=
f
L
(
x
)
=
a
L
x
+
1
-
a
L
d
L
if
0
≤
x
<
d
L
,
f
M
(
x
)
=
-
x
d
R
-
d
L
+
1
+
d
L
d
R
-
d
L
if
d
L
≤
x
<
d
R
,
f
R
(
x
)
=
a
R
x
-
a
R
d
R
if
d
R
≤
x
≤
1
,
where f
L
(
d
L
)
=
f
M
(
d
L
)
=
1, f
R
(
d
R
)
=
f
M
(
d
R
)
=
0, and the parameters satisfy
(3)
a
L
>
0
,
a
R
>
0
,
0
<
d
L
<
d
R
<
1
,
f
0
>
0
,
f
1
<
1
so that f
(
I
)
=
I. More correctly, this is a version of the map in which we can investigate the dynamic behavior changing independently the slopes of the external branches and border points. In [3] the authors give the relations between the parameters here considered and those related to the economic model. Moreover, the parameters used above (slopes and border points) may all or in part depend on some economic parameters, so that as an economic parameter is varied, it is also possible that it influences several of the above parameters, leading to particular paths in the parameter space.Figure 1
The functionF
(
Y
) in the trapping interval Y
m
,
Y
M.Clearly, this map (2) has a unique fixed point in the middle branch given by
(4)
x
M
*
=
d
R
1
+
d
R
-
d
L
which is always unstable given that the slope of the middle branch (-
1
/
(
d
R
-
d
L
)) is always smaller than -
1. The interval I is trapping so in I we do not have any divergent trajectory, and the attracting set may be periodic (an n-cycle with n
>
1) or chaotic (a chaotic interval or an n-band chaotic attractor with n
>
1). Some dynamic behaviors and bifurcation structures of map (2) will be described in the next section.


## 3. Bifurcation Structure of the Parameter Space

The Day-Shafer model presented in the previous section is described by a family of 1D continuous piecewise linear mapsf
:
I
→
I, I
=
[
0,1
], as given in (2), where the parameter region of interest, which allows obtaining the required shape of the map, is the region P defined as follows (as f
(
0
)
>
0 leads to a
L
d
R
<
1 while f
(
1
)
<
1 leads to a
R
(
1
-
d
R
)
<
1):
(5)
P
=
p
:
a
L
>
0
,
a
R
>
0
,
0
<
d
L
<
d
R
<
1
,
P
=
+
a
L
d
R
<
1
,
a
R
(
1
-
d
R
)
<
1
,
where p
=
(
a
L
,
a
R
,
d
L
,
d
R
) denotes a point in the parameter space. For p
∈
P the left and the right branches of map f are both increasing, while the middle one is decreasing, so that f is abimodal map. An example is shown in Figure 2. As already mentioned in Section 1, the dynamics of 1D bimodal maps have been considered by many researchers (see, e.g., [14, 25–27], etc.). Our aim is to study the bifurcation structure of the region P, that is, to describe the possible attractors of map f and the parameter regions corresponding to their existence.Figure 2
The mapf given in (2) for a
L
=
0.4, a
R
=
0.8, d
L
=
0.3, and d
R
=
0.8.Let us denote byI
L
=
[
0
,
d
L
), I
M
=
[
d
L
,
d
R
), and I
R
=
[
d
R
,
1
] the definition intervals (or, in other words,partitions) of the functions f
L, f
M, and f
R, respectively. They are separated by theborder points
d
L and d
R. As already remarked, for p
∈
P, the fixed point is x
M
*
∈
I
M, given in (4), always exists, and isrepelling, while the fixed points
(6)
x
L
*
=
a
L
d
L
-
1
a
L
-
1
,
x
R
*
=
a
R
d
R
a
R
-
1
,
associated with the branches f
L and f
R, exist (repelling) for a
L
>
1 and a
R
>
1. However, for the considered parameter values these fixed points do not belong to the interval I given that x
L
*
<
0 and x
R
*
>
1. For our map f the absorbing interval is always associated with all three branches. In fact, given that f
(
d
L
)
=
1
>
d
R and f
(
d
R
)
=
0
<
d
L it follows that I is anabsorbing invariant interval, f
(
I
)
=
I (see Figure 2), so that any orbit with an initial value x
0
∈
I isbounded, being trapped in I. We denote by φ
L and φ
R the boundaries related to the contact of the interval I with the fixed points x
L
* and x
R
*, given by x
L
*
=
0 and x
R
*
=
1, respectively (although not occurring for parameters in P):
(7)
φ
L
:
a
L
d
L
=
1
,
(8)
φ
R
:
a
R
1
-
d
R
=
1
.
Suppose that x
i
i
=
1
n are the points of an n-cycle of map f. The symbolic representation of such a cycle is σ
=
s
1
s
2
⋯
s
n, obtained associating to each point x
i the symbol s
i
∈
L
,
M
,
R depending on the partition I
L
,
I
M, or I
R which the point x
i belongs to. In the following, to denote an n-cycle we use its symbolic representation. The region in the parameter space related to the existence and stability of a cycle with symbolic sequence σ is denoted by P
σ and calledperiodicity region. Clearly, the boundaries of a periodicity region can be related either to the stability loss of the cycle due to its eigenvalue crossing ±
1 (recall that for a piecewise linear map such bifurcations aredegenerate, see [8]) or to the appearance/disappearance of the cycle due to aborder collision bifurcation (see [9]). We recall that if some point of a cycle collides with a border point and neither the period nor the stability of the cycle changes after the collision, we say that this cycle undergoes apersistence border collision, while aborder collision bifurcation (BCB) occurs when a qualitative change in the dynamics is observed after the collision.As noticed in [15], the overall bifurcation structure of the parameter space for a generic 1D bimodal piecewise linear map is characterized by several substructures among which we recall the skew tent map structure, the period adding structure, and a particular one called fin structure (due to the shape of the periodicity regions, as it will be clear below). The simplest one is theskew tent map structure associated with absorbing intervals involving only two adjacent partitions, so that on these absorbing intervals the map is reduced to a skew tent map. However, as already mentioned, by definition for our map f such a possibility is excluded. That is, the parameter space of f does not include the skew tent map structure. Theperiod adding structure is associated with periodicity regions related to attracting cycles whose points belong to the outermost partitions only. Such periodicity regions are ordered according to the Farey summation rule applied to the rotation numbers of the related cycles. This bifurcation structure is observed for our map f and we describe it in detail in Section 3.3. Thefin structure which also reveals itself in the parameter space of map f (See Section 3.3 for a deeper analysis) is closely related to the period adding structure being formed by the periodicity regions contiguous to the regions of the period adding structure and related to attracting cycles with just one point belonging to the middle partition and all the other points belonging to the outermost partitions.


### 3.1. Two-Dimensional Bifurcation Diagrams

We first present a few 2D bifurcation diagrams in various parameter planes to illustrate the overall bifurcation structure of the parameter space of mapf. In particular, Figure 3 shows such a diagram and its enlargement in the (
a
L
,
a
R
)-parameter plane for d
L
=
0.4, d
R
=
0.8, 0
<
a
L
<
1
/
d
L
=
2.5, and 0
<
a
R
<
1
/
(
1
-
d
R
)
=
5. Here different colors are related to attracting cycles of different periods n
≤
30, where the correspondence of a color and the period is indicated in the color bar; white region corresponds either to chaotic attractors or to cycles of higher periodicity. The gray region corresponds to a
L
>
1, a
R
>
1, so that all the slopes of f are larger than 1 in modulus; thus, attracting cycles cannot exist, and it is associated with chaotic attractors only (cyclic chaotic intervals).2D bifurcation diagram of the mapf in the (
a
L
,
a
R
)-parameter plane in (a), at d
L
=
0.4, d
R
=
0.8, and an enlargement of the indicated window in (b). The color bar indicates the correspondence of a color and the period of the related cycle.
(a)
(b)The choice of the parametersa
L and a
R to be varied for fixed values of d
L and d
R, as in Figure 3, is not optimal to illustrate the bifurcation structures typical for piecewise linear bimodal maps. To this purpose it is better to fix values for a
L and a
R and vary d
L and d
R, as it is shown in Figures 4 and 5. In these figures, the characteristic shapes of the periodicity regions belonging to the period adding structure are visible. Such regions have one side on the straight line d
L
=
d
R and also are observable regions which are contiguous (or attached) to the regions of the period adding structure and are those constituting the “fin structure.” A third example in the (
d
L
,
d
R
)-parameter plane is shown in Figure 5 at a
L
=
0.7 and a
R
=
0.8 fixed. We will see how to get the analytic equations of the BCB curves bounding the periodicity regions evidenced in the figures.2D bifurcation diagram off in the (
d
L
,
d
R
)-parameter plane for a
L
=
0.5, a
R
=
1.01 in (a), and a
L
=
0.6, a
R
=
1.1 in (b). The small red region in (b) (see left bottom), defined by d
L
<
d
R
<
1
-
1
/
a
R, is related to divergent orbits and not involved in the economic model.
(a)
(b)Figure 5
2D bifurcation diagram off in the (
d
L
,
d
R
)-parameter plane for a
L
=
0.7, a
R
=
0.8.


### 3.2. Period Adding Structure

Let us briefly recall the main elements of the period adding structure. As already mentioned, the periodicity regions constituting this structure are related to attractingn-cycles, n
≥
2, whose points only belong to the partitions I
L and I
R; that is, their symbolic sequences does not include the symbol M.Following [16] all the cycles associated with the period adding structure are grouped into families according tocomplexity levels. Thecomplexity level one includes two families, denoted by Σ
1,1 and Σ
2,1, to which the so-calledbasic cycles belong:
(9)
Σ
1,1
=
L
R
n
1
n
1
=
1
∞
,
Σ
2,1
=
R
L
n
1
n
1
=
1
∞
.
Note that the “central” cycle L
R
≡
R
L belongs to both families. To get the symbolic sequences of the cycles of families ofcomplexity level two we apply to the families Σ
1,1 and Σ
2,1 the followingsymbolic replacements:
(10)
κ
m
L
∶
=
L
⟶
L
R
m
R
⟶
R
L
R
m
,
κ
m
R
∶
=
L
⟶
L
R
L
m
R
⟶
R
L
m
.
This method is based on themap replacement technique (see [19, 20]). Namely, at first we substitute in Σ
1,1 each symbol L by L
R
m and each symbol R by R
L
R
m (replacement κ
m
L), and then we substitute in Σ
1,1 each symbol L by L
R
L
m and each symbol R by R
L
m (replacement κ
m
R). Then the index m is set m
=
n
2 in order to write the two families of complexity level two, respectively, as follows:
(11)
Σ
1,2
=
L
R
n
2
R
L
R
n
2
n
1
n
1
,
n
2
=
1
∞
,
Σ
2,2
=
L
R
L
n
2
R
L
n
2
n
1
n
1
,
n
2
=
1
∞
.
We notice that the replacement technique is here used to detect the symbolic representation of the existing cycles. However, the same technique is used also to get the equations of the BCB curves of cycles of complexity level higher than one, starting from the equations of those of complexity level one.Similarly, applying the replacementsκ
m
L and κ
m
R to Σ
2,1 we get the symbolic sequences of two more families:
(12)
Σ
3,2
=
R
L
R
n
2
L
R
n
2
n
1
n
1
,
n
2
=
1
∞
,
Σ
4,2
=
R
L
n
2
L
R
L
n
2
n
1
n
1
,
n
2
=
1
∞
.
Note that the central cycle L
R
R
L
R
≡
R
L
R
L
R belongs to both families Σ
1,2 and Σ
3,2, with n
1
=
n
2
=
1, while the central cycle L
R
L
R
L
≡
R
L
L
R
L belongs to both families Σ
2,2 and Σ
4,2, with n
1
=
n
2
=
1. All the other cycles in these families are distinct. In short this procedure can be written as Σ
1,2
=
κ
n
2
L
(
Σ
1,1
), Σ
2,2
=
κ
n
2
R
(
Σ
1,1
), Σ
3,2
=
κ
n
2
L
(
Σ
2,1
), and Σ
4,2
=
κ
n
2
R
(
Σ
2,1
). So, we get 4 families of complexity level two (one more way to construct the families of the complexity level two consists in consecutive application of the concatenation rule to the “neighbour” symbolic sequences of the first complexity level. Symbolic sequences obtained in such a way are shift invariant to those obtained by symbolic replacements (10) (see [19, 20])). Further, applying the replacements (10) with m
=
n
3 to the families of complexity level two we obtain 2
3 families Σ
j
,
3, j
=
1
,
…
,
2
3, ofcomplexity level three, and so on. In this way all the symbolic sequences of cycles associated with the period adding structure are obtained.Now let us turn to the boundaries of the periodicity regions (i.e., the BCB curves) related to the cycles of mapf associated with the period adding structure. They can be confined by the boundaries related to their existence. A cycle disappears when are crossed these boundaries, which are related to the BCBs occurring when a point of the cycle close to a border point, x
=
d
L or x
=
d
R, collides with it in a saddle-node border collision bifurcation (merging with a companion unstable cycle). For the boundaries of a periodicity region associated with the stability of the cycle, first note that an n-cycle whose symbolic sequence σ does not include the symbol M has multiplier λ
σ
=
a
L
k
a
R
n
-
k
>
0, where k and n
-
k are the numbers of symbols L and R, respectively, in σ. Thus, a degenerate flip bifurcation (DFB for short), related to λ
σ
=
-
1, cannot occur for such a cycle, so that the periodicity regions of the period adding structure cannot have DFB boundaries. A degenerate +
1 bifurcation (DB1 for short) associated with λ
σ
=
1 occurs if a
L
k
a
R
n
-
k
=
1. Obviously, this condition defines a boundary of the periodicity region only if the related cycle exists.To illustrate how the periodicity regions of the period adding structure are ordered let us consider first the limit cased
L
=
d
R
≡
d at which the considered map becomes discontinuous, say f
~,
(13)
x
⟼
f
~
(
x
)
=
f
~
L
(
x
)
=
a
L
x
+
1
-
a
L
d
if
0
≤
x
<
d
,
f
~
R
(
x
)
=
a
R
x
-
a
R
d
if
d
<
x
≤
1
.
In Figures 4 and 5, when the parameters belong to the straight line d
L
=
d
R then the continuous bimodal map f reduces to the discontinuous one, f
~.The dynamics of 1D discontinuous piecewise monotone maps have been studied by many researchers (see, e.g., [17, 18, 28, 29]). In particular, the piecewise linear case has been recently reconsidered (after [16]) in [19, 20]. In the cited references, it is described as the period adding structure which is characteristic for piecewise increasing maps, when invertible on the absorbing interval. It is easy to check that map f
~ is invertible on the absorbing interval I if f
L
(
f
(
d
)
)
>
f
R
(
f
(
d
)
), that is, if f
L
(
0
)
>
f
R
(
1
), in which case f
~ is calledgap map. It is noninvertible if f
L
(
0
)
<
f
R
(
1
) being also calledoverlapping map, while for f
L
(
0
)
=
f
R
(
1
) map f
~ is calledcircle map. The boundary defined by f
L
(
0
)
=
f
R
(
1
), that holds for
(14)
κ
=
(
a
L
,
a
R
,
d
)
:
a
L
>
0
,
a
R
>
0
,
κ
+
=
0
<
d
<
1
,
1
-
d
a
L
=
a
R
(
1
-
d
)
,
is related to the changes of invertibility of f
~. In Figure 6 we show the 2D bifurcation diagram of f
~ in the (
d
,
a
R
)-parameter plane for a
L
=
0.5, and the curve κ is there plotted. Note that the dynamics along the straight line defined by a
R
=
1.1 is related to the dynamics along the straight line d
L
=
d
R in Figure 4(a). Below the curve κ one can observe the period adding structure. Above the curve φ
R given in (8) a generic trajectory of f
~ diverges, while in between φ
R and κ the existing attractors of the map are chaotic intervals.Figure 6
2D bifurcation diagram of mapf
~ in the (
d
,
a
R
)-parameter plane for a
L
=
0.5.The period adding structure in mapf
~ is illustrated by the 1D bifurcation diagram in Figure 7 corresponding to the cross section along the horizontal line with an arrow indicated in Figure 6. Note that for φ
R
<
d
<
κ, where φ
R
≈
0.167, κ
≈
0.286, the attractor is chaotic. It can be also clearly seen that the boundaries of the periodicity regions, for example, of the 2-cycle, are related to the collision of the points of the cycle with the border point x
=
d.Figure 7
1D bifurcation diagram of mapf
~ for a
L
=
0.5, a
R
=
1.2, and d
∈
(
φ
R
,
1
), where φ
R
≈
0.167. The related parameter path is indicated in Figure 6 by the horizontal line with an arrow.In AppendixA we give the analytic equations of the boundaries of the periodicity regions of the period adding structure for map f, which holds also for map f
~ substituting d
L
=
d
R
=
d. The periodicity regions of complexity level one and two of map f, corresponding to the bifurcation diagram presented in Figure 5, are shown in Figure 8 by light gray and dark gray regions, respectively. In that figure, the boundaries of the gray regions are plotted using the formulas given in Appendix A. For example, after simplifications we get that periodicity region P
R
L related to the attracting 2-cycle R
L of map f is defined as follows:
(15)
P
R
L
=
p
∈
P
:
1
-
d
L
a
R
<
d
R
<
1
-
a
L
d
L
,
and the equations d
R
=
1
-
(
d
L
/
a
R
), d
R
=
1
-
a
L
d
L define the two BCB curves giving the boundaries of P
R
L. For fixed d
L
=
0.4, d
R
=
0.8, as in Figure 3, the boundaries of P
R
L are just segments of the vertical and horizontal straight lines defined by a
L
=
(
1
-
d
R
)
/
d
L
=
0.5 and a
R
=
d
L
/
(
1
-
d
R
)
=
2, respectively.Figure 8
Periodicity regions of the cyclesL
R
n and R
L
n, n
=
1
,
…
,
11, of complexity level one are shown in light gray. A few periodicity regions of complexity level two are shown in dark gray. Here a
L
=
0.7, a
R
=
0.8 as in Figure 5.


### 3.3. Fin Structure

As already mentioned, the fin structure in the parameter space of a 1D bimodal piecewise linear map consists of periodicity regions which are attached to the regions of the period adding structure described in the previous section. For example, one can clearly see in Figure5 two 2
·
2-periodicity regions and two 2
·
3-periodicity regions attached on both sides to the period-2 region P
R
L, as well as 3
·
2-, 3
·
3- and 3
·
4-periodicity regions attached on both sides to the period-3 regions P
R
L
2 and P
L
R
2, and so on. These regions belong to the fin structure which is formed by the periodicity regions called n
·
k-fins, k
≥
1, related to attracting cycles having only one point in the interval I
M and all the other points are in I
L and I
R. The n-periodicity region of the period adding structure to which a fin is attached is calledtrunk region, and its fins have the same complexity level as the complexity level of the trunk. In fact, in Figure 5, two 2
·
1-fins of complexity level one of the trunk regions P
R
L can also be seen.A fin cycle can appear if the parameter point crosses one of two BCB boundaries of a trunk, due to which one periodic point enters the partitionI
M. As explained in [15], for an n
·
k-cycle whose periodicity region has the common boundary with the region P
L
R
n
-
1, the symbolic sequences of the cycles in the fins are (
L
R
n
-
1
)
k
-
1
M
R
n
-
1, k
=
1,2
,
…, on one side of the region P
L
R
n
-
1, and (
L
R
n
-
1
)
k
-
1
L
R
n
-
2
M on the other side. Interchanging L and R in these sequences we get the symbolic sequences of the cycles related to n
·
k-fins whose trunks are P
R
L
n
-
1 regions. The number of existing fins depends on the parameters and, in fact, some trunk regions may have fins on one side only, or have no fins at all, as can be seen, for example, in Figure 4(b).As for the boundaries of a fin, it can be shown that eachn
·
k-fin region for k
≥
2, n
≥
2, related to a cycle with symbolic sequence σ, has at most four boundaries, among which one is the common BCB boundary with the related trunk region, one boundary is related to DFB of the cycle (whose eigenvalue is λ
σ
=
a
L
m
a
R
n
k
-
m
-
1
/
(
d
L
-
d
R
)
<
0, where m is the number of symbols L in σ), and two other boundaries are related to two more BCBs of the cycle. Each n
·
1-fin region has only three boundaries, namely, one DFB boundary and two BCB boundaries. The DFB boundary is defined by the condition λ
σ
=
-
1, while the BCB boundaries are obtained using the skew tent map as border collision normal form.In AppendixB, some basic formulas related to the boundaries of a fin region are given, while below, as an example, we describe the 2
·
1- and 2
·
2-fins of the region P
R
L shown in Figure 9. There are two 2
·
1-fins, P
M
R and P
L
M, contiguous to the trunk regions P
R
L. As we mentioned above, the 2
·
1-fins are exceptional being confined by three boundaries instead of four. In particular, the fin P
M
R is confined by the BCB and DFB boundaries
(16)
B
C
L
R
=
p
∈
P
:
d
R
=
1
-
d
L
a
R
,
D
F
M
R
=
p
∈
P
:
d
R
=
d
L
+
a
R
,
respectively, and the boundary defined by d
L
=
0. The fin P
L
M is confined by the BCB and DFB boundaries
(17)
B
C
R
L
=
p
∈
P
:
d
R
=
1
-
a
L
d
L
,
D
F
L
M
=
p
∈
P
:
d
R
=
d
L
+
a
L
,
respectively, and the boundary d
R
=
1.Figure 9
The2
·
1-fins P
M
R and P
L
M, and 2
·
2-fins P
L
R
M
R and P
L
R
L
M, contiguous to the trunk regions P
L
R. Here a
L
=
0.7, a
R
=
0.8 as in Figure 5.Next, let us consider the2
·
2-fins, P
L
R
M
R and P
L
R
L
M, which are also contiguous to the trunk regions P
R
L. Applying the formulas given in Appendix B we get that the boundaries of the fin P
L
R
M
R are given by the BCB curves satisfying the equations given below:
(18)
B
C
L
R
=
p
∈
P
:
d
R
=
1
-
d
L
a
R
,
D
F
L
R
M
R
=
p
∈
P
:
d
R
=
d
L
+
a
R
2
a
L
,
B
C
L
R
M
R
1
=
p
∈
P
:
d
R
=
d
L
+
a
R
,
B
C
L
R
M
R
2
=
p
∈
P
:
d
L
=
d
R
2
-
(
a
R
+
1
)
d
R
+
a
R
a
R
a
L
+
d
R
.
Note that B
C
L
R
M
R
1
=
D
F
M
R; that is, the DFB of the cycle L
M occurs simultaneously with the BCB of the cycle L
R
M
R (see Figure 9). The fin P
L
R
L
M is confined by the following boundaries:
(19)
B
C
R
L
=
p
∈
P
:
d
R
=
1
-
a
L
d
L
,
D
F
L
R
L
M
=
p
∈
P
:
d
R
=
d
L
+
a
R
a
L
2
,
B
C
L
R
L
M
1
=
p
∈
P
:
d
R
=
d
L
+
a
L
,
B
C
L
R
L
M
2
=
p
∈
P
:
d
R
=
d
L
2
+
a
L
(
d
L
-
a
R
)
-
1
d
L
-
1
-
a
L
a
R
,
and, as already remarked in the previous case, here also we have B
C
L
R
L
M
1
=
D
F
L
M (see Figure 9).We note, from the point of view of our specific economic application, the presence of fin cycles is important because it makes it possible to reach intermediate levels of income (those belonging to the middle branch) that are excluded in period adding cycles. Moreover, the fin structure makes the bifurcation diagram even more complicated and as a consequence it becomes quite hard to foresee what would happen when a parameter of the map varies.


## 3.1. Two-Dimensional Bifurcation Diagrams

We first present a few 2D bifurcation diagrams in various parameter planes to illustrate the overall bifurcation structure of the parameter space of mapf. In particular, Figure 3 shows such a diagram and its enlargement in the (
a
L
,
a
R
)-parameter plane for d
L
=
0.4, d
R
=
0.8, 0
<
a
L
<
1
/
d
L
=
2.5, and 0
<
a
R
<
1
/
(
1
-
d
R
)
=
5. Here different colors are related to attracting cycles of different periods n
≤
30, where the correspondence of a color and the period is indicated in the color bar; white region corresponds either to chaotic attractors or to cycles of higher periodicity. The gray region corresponds to a
L
>
1, a
R
>
1, so that all the slopes of f are larger than 1 in modulus; thus, attracting cycles cannot exist, and it is associated with chaotic attractors only (cyclic chaotic intervals).2D bifurcation diagram of the mapf in the (
a
L
,
a
R
)-parameter plane in (a), at d
L
=
0.4, d
R
=
0.8, and an enlargement of the indicated window in (b). The color bar indicates the correspondence of a color and the period of the related cycle.
(a)
(b)The choice of the parametersa
L and a
R to be varied for fixed values of d
L and d
R, as in Figure 3, is not optimal to illustrate the bifurcation structures typical for piecewise linear bimodal maps. To this purpose it is better to fix values for a
L and a
R and vary d
L and d
R, as it is shown in Figures 4 and 5. In these figures, the characteristic shapes of the periodicity regions belonging to the period adding structure are visible. Such regions have one side on the straight line d
L
=
d
R and also are observable regions which are contiguous (or attached) to the regions of the period adding structure and are those constituting the “fin structure.” A third example in the (
d
L
,
d
R
)-parameter plane is shown in Figure 5 at a
L
=
0.7 and a
R
=
0.8 fixed. We will see how to get the analytic equations of the BCB curves bounding the periodicity regions evidenced in the figures.2D bifurcation diagram off in the (
d
L
,
d
R
)-parameter plane for a
L
=
0.5, a
R
=
1.01 in (a), and a
L
=
0.6, a
R
=
1.1 in (b). The small red region in (b) (see left bottom), defined by d
L
<
d
R
<
1
-
1
/
a
R, is related to divergent orbits and not involved in the economic model.
(a)
(b)Figure 5
2D bifurcation diagram off in the (
d
L
,
d
R
)-parameter plane for a
L
=
0.7, a
R
=
0.8.


## 3.2. Period Adding Structure

Let us briefly recall the main elements of the period adding structure. As already mentioned, the periodicity regions constituting this structure are related to attractingn-cycles, n
≥
2, whose points only belong to the partitions I
L and I
R; that is, their symbolic sequences does not include the symbol M.Following [16] all the cycles associated with the period adding structure are grouped into families according tocomplexity levels. Thecomplexity level one includes two families, denoted by Σ
1,1 and Σ
2,1, to which the so-calledbasic cycles belong:
(9)
Σ
1,1
=
L
R
n
1
n
1
=
1
∞
,
Σ
2,1
=
R
L
n
1
n
1
=
1
∞
.
Note that the “central” cycle L
R
≡
R
L belongs to both families. To get the symbolic sequences of the cycles of families ofcomplexity level two we apply to the families Σ
1,1 and Σ
2,1 the followingsymbolic replacements:
(10)
κ
m
L
∶
=
L
⟶
L
R
m
R
⟶
R
L
R
m
,
κ
m
R
∶
=
L
⟶
L
R
L
m
R
⟶
R
L
m
.
This method is based on themap replacement technique (see [19, 20]). Namely, at first we substitute in Σ
1,1 each symbol L by L
R
m and each symbol R by R
L
R
m (replacement κ
m
L), and then we substitute in Σ
1,1 each symbol L by L
R
L
m and each symbol R by R
L
m (replacement κ
m
R). Then the index m is set m
=
n
2 in order to write the two families of complexity level two, respectively, as follows:
(11)
Σ
1,2
=
L
R
n
2
R
L
R
n
2
n
1
n
1
,
n
2
=
1
∞
,
Σ
2,2
=
L
R
L
n
2
R
L
n
2
n
1
n
1
,
n
2
=
1
∞
.
We notice that the replacement technique is here used to detect the symbolic representation of the existing cycles. However, the same technique is used also to get the equations of the BCB curves of cycles of complexity level higher than one, starting from the equations of those of complexity level one.Similarly, applying the replacementsκ
m
L and κ
m
R to Σ
2,1 we get the symbolic sequences of two more families:
(12)
Σ
3,2
=
R
L
R
n
2
L
R
n
2
n
1
n
1
,
n
2
=
1
∞
,
Σ
4,2
=
R
L
n
2
L
R
L
n
2
n
1
n
1
,
n
2
=
1
∞
.
Note that the central cycle L
R
R
L
R
≡
R
L
R
L
R belongs to both families Σ
1,2 and Σ
3,2, with n
1
=
n
2
=
1, while the central cycle L
R
L
R
L
≡
R
L
L
R
L belongs to both families Σ
2,2 and Σ
4,2, with n
1
=
n
2
=
1. All the other cycles in these families are distinct. In short this procedure can be written as Σ
1,2
=
κ
n
2
L
(
Σ
1,1
), Σ
2,2
=
κ
n
2
R
(
Σ
1,1
), Σ
3,2
=
κ
n
2
L
(
Σ
2,1
), and Σ
4,2
=
κ
n
2
R
(
Σ
2,1
). So, we get 4 families of complexity level two (one more way to construct the families of the complexity level two consists in consecutive application of the concatenation rule to the “neighbour” symbolic sequences of the first complexity level. Symbolic sequences obtained in such a way are shift invariant to those obtained by symbolic replacements (10) (see [19, 20])). Further, applying the replacements (10) with m
=
n
3 to the families of complexity level two we obtain 2
3 families Σ
j
,
3, j
=
1
,
…
,
2
3, ofcomplexity level three, and so on. In this way all the symbolic sequences of cycles associated with the period adding structure are obtained.Now let us turn to the boundaries of the periodicity regions (i.e., the BCB curves) related to the cycles of mapf associated with the period adding structure. They can be confined by the boundaries related to their existence. A cycle disappears when are crossed these boundaries, which are related to the BCBs occurring when a point of the cycle close to a border point, x
=
d
L or x
=
d
R, collides with it in a saddle-node border collision bifurcation (merging with a companion unstable cycle). For the boundaries of a periodicity region associated with the stability of the cycle, first note that an n-cycle whose symbolic sequence σ does not include the symbol M has multiplier λ
σ
=
a
L
k
a
R
n
-
k
>
0, where k and n
-
k are the numbers of symbols L and R, respectively, in σ. Thus, a degenerate flip bifurcation (DFB for short), related to λ
σ
=
-
1, cannot occur for such a cycle, so that the periodicity regions of the period adding structure cannot have DFB boundaries. A degenerate +
1 bifurcation (DB1 for short) associated with λ
σ
=
1 occurs if a
L
k
a
R
n
-
k
=
1. Obviously, this condition defines a boundary of the periodicity region only if the related cycle exists.To illustrate how the periodicity regions of the period adding structure are ordered let us consider first the limit cased
L
=
d
R
≡
d at which the considered map becomes discontinuous, say f
~,
(13)
x
⟼
f
~
(
x
)
=
f
~
L
(
x
)
=
a
L
x
+
1
-
a
L
d
if
0
≤
x
<
d
,
f
~
R
(
x
)
=
a
R
x
-
a
R
d
if
d
<
x
≤
1
.
In Figures 4 and 5, when the parameters belong to the straight line d
L
=
d
R then the continuous bimodal map f reduces to the discontinuous one, f
~.The dynamics of 1D discontinuous piecewise monotone maps have been studied by many researchers (see, e.g., [17, 18, 28, 29]). In particular, the piecewise linear case has been recently reconsidered (after [16]) in [19, 20]. In the cited references, it is described as the period adding structure which is characteristic for piecewise increasing maps, when invertible on the absorbing interval. It is easy to check that map f
~ is invertible on the absorbing interval I if f
L
(
f
(
d
)
)
>
f
R
(
f
(
d
)
), that is, if f
L
(
0
)
>
f
R
(
1
), in which case f
~ is calledgap map. It is noninvertible if f
L
(
0
)
<
f
R
(
1
) being also calledoverlapping map, while for f
L
(
0
)
=
f
R
(
1
) map f
~ is calledcircle map. The boundary defined by f
L
(
0
)
=
f
R
(
1
), that holds for
(14)
κ
=
(
a
L
,
a
R
,
d
)
:
a
L
>
0
,
a
R
>
0
,
κ
+
=
0
<
d
<
1
,
1
-
d
a
L
=
a
R
(
1
-
d
)
,
is related to the changes of invertibility of f
~. In Figure 6 we show the 2D bifurcation diagram of f
~ in the (
d
,
a
R
)-parameter plane for a
L
=
0.5, and the curve κ is there plotted. Note that the dynamics along the straight line defined by a
R
=
1.1 is related to the dynamics along the straight line d
L
=
d
R in Figure 4(a). Below the curve κ one can observe the period adding structure. Above the curve φ
R given in (8) a generic trajectory of f
~ diverges, while in between φ
R and κ the existing attractors of the map are chaotic intervals.Figure 6
2D bifurcation diagram of mapf
~ in the (
d
,
a
R
)-parameter plane for a
L
=
0.5.The period adding structure in mapf
~ is illustrated by the 1D bifurcation diagram in Figure 7 corresponding to the cross section along the horizontal line with an arrow indicated in Figure 6. Note that for φ
R
<
d
<
κ, where φ
R
≈
0.167, κ
≈
0.286, the attractor is chaotic. It can be also clearly seen that the boundaries of the periodicity regions, for example, of the 2-cycle, are related to the collision of the points of the cycle with the border point x
=
d.Figure 7
1D bifurcation diagram of mapf
~ for a
L
=
0.5, a
R
=
1.2, and d
∈
(
φ
R
,
1
), where φ
R
≈
0.167. The related parameter path is indicated in Figure 6 by the horizontal line with an arrow.In AppendixA we give the analytic equations of the boundaries of the periodicity regions of the period adding structure for map f, which holds also for map f
~ substituting d
L
=
d
R
=
d. The periodicity regions of complexity level one and two of map f, corresponding to the bifurcation diagram presented in Figure 5, are shown in Figure 8 by light gray and dark gray regions, respectively. In that figure, the boundaries of the gray regions are plotted using the formulas given in Appendix A. For example, after simplifications we get that periodicity region P
R
L related to the attracting 2-cycle R
L of map f is defined as follows:
(15)
P
R
L
=
p
∈
P
:
1
-
d
L
a
R
<
d
R
<
1
-
a
L
d
L
,
and the equations d
R
=
1
-
(
d
L
/
a
R
), d
R
=
1
-
a
L
d
L define the two BCB curves giving the boundaries of P
R
L. For fixed d
L
=
0.4, d
R
=
0.8, as in Figure 3, the boundaries of P
R
L are just segments of the vertical and horizontal straight lines defined by a
L
=
(
1
-
d
R
)
/
d
L
=
0.5 and a
R
=
d
L
/
(
1
-
d
R
)
=
2, respectively.Figure 8
Periodicity regions of the cyclesL
R
n and R
L
n, n
=
1
,
…
,
11, of complexity level one are shown in light gray. A few periodicity regions of complexity level two are shown in dark gray. Here a
L
=
0.7, a
R
=
0.8 as in Figure 5.


## 3.3. Fin Structure

As already mentioned, the fin structure in the parameter space of a 1D bimodal piecewise linear map consists of periodicity regions which are attached to the regions of the period adding structure described in the previous section. For example, one can clearly see in Figure5 two 2
·
2-periodicity regions and two 2
·
3-periodicity regions attached on both sides to the period-2 region P
R
L, as well as 3
·
2-, 3
·
3- and 3
·
4-periodicity regions attached on both sides to the period-3 regions P
R
L
2 and P
L
R
2, and so on. These regions belong to the fin structure which is formed by the periodicity regions called n
·
k-fins, k
≥
1, related to attracting cycles having only one point in the interval I
M and all the other points are in I
L and I
R. The n-periodicity region of the period adding structure to which a fin is attached is calledtrunk region, and its fins have the same complexity level as the complexity level of the trunk. In fact, in Figure 5, two 2
·
1-fins of complexity level one of the trunk regions P
R
L can also be seen.A fin cycle can appear if the parameter point crosses one of two BCB boundaries of a trunk, due to which one periodic point enters the partitionI
M. As explained in [15], for an n
·
k-cycle whose periodicity region has the common boundary with the region P
L
R
n
-
1, the symbolic sequences of the cycles in the fins are (
L
R
n
-
1
)
k
-
1
M
R
n
-
1, k
=
1,2
,
…, on one side of the region P
L
R
n
-
1, and (
L
R
n
-
1
)
k
-
1
L
R
n
-
2
M on the other side. Interchanging L and R in these sequences we get the symbolic sequences of the cycles related to n
·
k-fins whose trunks are P
R
L
n
-
1 regions. The number of existing fins depends on the parameters and, in fact, some trunk regions may have fins on one side only, or have no fins at all, as can be seen, for example, in Figure 4(b).As for the boundaries of a fin, it can be shown that eachn
·
k-fin region for k
≥
2, n
≥
2, related to a cycle with symbolic sequence σ, has at most four boundaries, among which one is the common BCB boundary with the related trunk region, one boundary is related to DFB of the cycle (whose eigenvalue is λ
σ
=
a
L
m
a
R
n
k
-
m
-
1
/
(
d
L
-
d
R
)
<
0, where m is the number of symbols L in σ), and two other boundaries are related to two more BCBs of the cycle. Each n
·
1-fin region has only three boundaries, namely, one DFB boundary and two BCB boundaries. The DFB boundary is defined by the condition λ
σ
=
-
1, while the BCB boundaries are obtained using the skew tent map as border collision normal form.In AppendixB, some basic formulas related to the boundaries of a fin region are given, while below, as an example, we describe the 2
·
1- and 2
·
2-fins of the region P
R
L shown in Figure 9. There are two 2
·
1-fins, P
M
R and P
L
M, contiguous to the trunk regions P
R
L. As we mentioned above, the 2
·
1-fins are exceptional being confined by three boundaries instead of four. In particular, the fin P
M
R is confined by the BCB and DFB boundaries
(16)
B
C
L
R
=
p
∈
P
:
d
R
=
1
-
d
L
a
R
,
D
F
M
R
=
p
∈
P
:
d
R
=
d
L
+
a
R
,
respectively, and the boundary defined by d
L
=
0. The fin P
L
M is confined by the BCB and DFB boundaries
(17)
B
C
R
L
=
p
∈
P
:
d
R
=
1
-
a
L
d
L
,
D
F
L
M
=
p
∈
P
:
d
R
=
d
L
+
a
L
,
respectively, and the boundary d
R
=
1.Figure 9
The2
·
1-fins P
M
R and P
L
M, and 2
·
2-fins P
L
R
M
R and P
L
R
L
M, contiguous to the trunk regions P
L
R. Here a
L
=
0.7, a
R
=
0.8 as in Figure 5.Next, let us consider the2
·
2-fins, P
L
R
M
R and P
L
R
L
M, which are also contiguous to the trunk regions P
R
L. Applying the formulas given in Appendix B we get that the boundaries of the fin P
L
R
M
R are given by the BCB curves satisfying the equations given below:
(18)
B
C
L
R
=
p
∈
P
:
d
R
=
1
-
d
L
a
R
,
D
F
L
R
M
R
=
p
∈
P
:
d
R
=
d
L
+
a
R
2
a
L
,
B
C
L
R
M
R
1
=
p
∈
P
:
d
R
=
d
L
+
a
R
,
B
C
L
R
M
R
2
=
p
∈
P
:
d
L
=
d
R
2
-
(
a
R
+
1
)
d
R
+
a
R
a
R
a
L
+
d
R
.
Note that B
C
L
R
M
R
1
=
D
F
M
R; that is, the DFB of the cycle L
M occurs simultaneously with the BCB of the cycle L
R
M
R (see Figure 9). The fin P
L
R
L
M is confined by the following boundaries:
(19)
B
C
R
L
=
p
∈
P
:
d
R
=
1
-
a
L
d
L
,
D
F
L
R
L
M
=
p
∈
P
:
d
R
=
d
L
+
a
R
a
L
2
,
B
C
L
R
L
M
1
=
p
∈
P
:
d
R
=
d
L
+
a
L
,
B
C
L
R
L
M
2
=
p
∈
P
:
d
R
=
d
L
2
+
a
L
(
d
L
-
a
R
)
-
1
d
L
-
1
-
a
L
a
R
,
and, as already remarked in the previous case, here also we have B
C
L
R
L
M
1
=
D
F
L
M (see Figure 9).We note, from the point of view of our specific economic application, the presence of fin cycles is important because it makes it possible to reach intermediate levels of income (those belonging to the middle branch) that are excluded in period adding cycles. Moreover, the fin structure makes the bifurcation diagram even more complicated and as a consequence it becomes quite hard to foresee what would happen when a parameter of the map varies.


## 4. Conclusions

In this work we have considered a pioneering model by Day and Shafer [3] which describes a business cycle by using a bimodal piecewise linear map. Our investigation shows how rich are the dynamic behaviors of the system, going from attracting cycles of any period to robust chaotic intervals, depending on the parameters values. By using recently developed techniques, we have studied a typical bifurcation diagram and obtained analytically the border collision bifurcation curves that separate different periodicity regions and degenerate flip bifurcation curves. The transition to chaos was also studied.Complex dynamics are usually the consequence of the introduction of some nonlinearity in a system. Moreover, specially in economics, it is quite natural to have a nonlinearity coming simply by some constraint. That is, an economic model is often characterized by different functional definitions depending on some threshold reached by the dynamic variables (income, prices, etc.). In this way, the models are described by piecewise smooth systems (in place of smooth ones), whose theoretical results are still under study among scholars in dynamical systems. Thus, results similar to those of our model can be obtained in the study of other continuous piecewise linear systems.



---

*Source: 401319-2014-11-19.xml*