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.

1. Introduction

Applied models defined by piecewise 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 as skew 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 -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., [58]). Moreover, the skew tent map can be used as a normal form for a so-called border 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., [1113], 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 known period adding structure (called also Arnold tongues or mode-locking tongues), which is characteristic for piecewise increasing discontinuous maps and also for certain circle maps (see, e.g., [1620]). 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-called fin 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 Section 2 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 -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., [2124]). 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 two turning 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 as a function of the income and a money supply parameter . In the real market the interest rate determines the level of investments , that is, a component of the income together with consumption . 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: where the parameter measures how strongly investment translates into new income, and the parameter includes all the autonomous components (of investment, consumption, and public expenditure). The function is typically monotonically increasing with , and the investment function 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 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 ; it is continuous but not differentiable at the points and . The dynamics are bounded in the interval . For convenience, we can normalize the interval to and denoting by the state variable income we obtain the map , defined as follows: where , , and the parameters satisfy so that . 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.

Clearly, this map (2) has a unique fixed point in the middle branch given by which is always unstable given that the slope of the middle branch () is always smaller than . The interval is trapping so in we do not have any divergent trajectory, and the attracting set may be periodic (an -cycle with ) or chaotic (a chaotic interval or an -band chaotic attractor with ). 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 maps , , as given in (2), where the parameter region of interest, which allows obtaining the required shape of the map, is the region defined as follows (as leads to while leads to ): where denotes a point in the parameter space. For the left and the right branches of map are both increasing, while the middle one is decreasing, so that is a bimodal 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, 2527], etc.). Our aim is to study the bifurcation structure of the region , that is, to describe the possible attractors of map and the parameter regions corresponding to their existence.

Let us denote by , , and the definition intervals (or, in other words, partitions) of the functions , , and , respectively. They are separated by the border points and . As already remarked, for , the fixed point is , given in (4), always exists, and is repelling, while the fixed points associated with the branches and , exist (repelling) for and . However, for the considered parameter values these fixed points do not belong to the interval given that and . For our map the absorbing interval is always associated with all three branches. In fact, given that and it follows that is an absorbing invariant interval, (see Figure 2), so that any orbit with an initial value is bounded, being trapped in . We denote by and the boundaries related to the contact of the interval with the fixed points and , given by and , respectively (although not occurring for parameters in ): Suppose that are the points of an -cycle of map . The symbolic representation of such a cycle is , obtained associating to each point the symbol depending on the partition , or which the point belongs to. In the following, to denote an -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 and called periodicity region. Clearly, the boundaries of a periodicity region can be related either to the stability loss of the cycle due to its eigenvalue crossing (recall that for a piecewise linear map such bifurcations are degenerate, see [8]) or to the appearance/disappearance of the cycle due to a border 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 a persistence border collision, while a border 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 the skew 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 such a possibility is excluded. That is, the parameter space of does not include the skew tent map structure. The period 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 and we describe it in detail in Section 3.3. The fin structure which also reveals itself in the parameter space of map (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 map . In particular, Figure 3 shows such a diagram and its enlargement in the -parameter plane for , , , and . Here different colors are related to attracting cycles of different periods , 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 , , so that all the slopes of are larger than in modulus; thus, attracting cycles cannot exist, and it is associated with chaotic attractors only (cyclic chaotic intervals).

The choice of the parameters and to be varied for fixed values of and , 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 and and vary and , 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 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 -parameter plane is shown in Figure 5 at and fixed. We will see how to get the analytic equations of the BCB curves bounding the periodicity regions evidenced in the figures.

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 attracting -cycles, , whose points only belong to the partitions and ; that is, their symbolic sequences does not include the symbol .

Following [16] all the cycles associated with the period adding structure are grouped into families according to complexity levels. The complexity level one includes two families, denoted by and , to which the so-called basic cycles belong: Note that the “central” cycle belongs to both families. To get the symbolic sequences of the cycles of families of complexity level two we apply to the families and the following symbolic replacements: This method is based on the map replacement technique (see [19, 20]). Namely, at first we substitute in each symbol by and each symbol by (replacement ), and then we substitute in each symbol by and each symbol by (replacement ). Then the index is set in order to write the two families of complexity level two, respectively, as follows: 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 and to we get the symbolic sequences of two more families: Note that the central cycle belongs to both families and , with , while the central cycle belongs to both families and , with . All the other cycles in these families are distinct. In short this procedure can be written as , , , and . So, we get 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 to the families of complexity level two we obtain families , , of complexity 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 map 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, or , 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 -cycle whose symbolic sequence does not include the symbol has multiplier , where and are the numbers of symbols and , respectively, in . Thus, a degenerate flip bifurcation (DFB for short), related to , cannot occur for such a cycle, so that the periodicity regions of the period adding structure cannot have DFB boundaries. A degenerate bifurcation (DB1 for short) associated with occurs if . 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 case at which the considered map becomes discontinuous, say , In Figures 4 and 5, when the parameters belong to the straight line then the continuous bimodal map reduces to the discontinuous one, .

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 is invertible on the absorbing interval if , that is, if , in which case is called gap map. It is noninvertible if being also called overlapping map, while for map is called circle map. The boundary defined by , that holds for is related to the changes of invertibility of . In Figure 6 we show the 2D bifurcation diagram of in the -parameter plane for , and the curve is there plotted. Note that the dynamics along the straight line defined by is related to the dynamics along the straight line in Figure 4(a). Below the curve one can observe the period adding structure. Above the curve given in (8) a generic trajectory of diverges, while in between and the existing attractors of the map are chaotic intervals.

The period adding structure in map 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 , where , , the attractor is chaotic. It can be also clearly seen that the boundaries of the periodicity regions, for example, of the -cycle, are related to the collision of the points of the cycle with the border point .

In Appendix A we give the analytic equations of the boundaries of the periodicity regions of the period adding structure for map , which holds also for map substituting . The periodicity regions of complexity level one and two of map , 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 related to the attracting -cycle of map is defined as follows: and the equations , define the two BCB curves giving the boundaries of . For fixed , , as in Figure 3, the boundaries of are just segments of the vertical and horizontal straight lines defined by and , respectively.

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 Figure 5 two -periodicity regions and two -periodicity regions attached on both sides to the period- region , as well as -, - and -periodicity regions attached on both sides to the period- regions and , and so on. These regions belong to the fin structure which is formed by the periodicity regions called -fins, , related to attracting cycles having only one point in the interval and all the other points are in and . The -periodicity region of the period adding structure to which a fin is attached is called trunk region, and its fins have the same complexity level as the complexity level of the trunk. In fact, in Figure 5, two -fins of complexity level one of the trunk regions 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 partition . As explained in [15], for an -cycle whose periodicity region has the common boundary with the region , the symbolic sequences of the cycles in the fins are , , on one side of the region , and on the other side. Interchanging and in these sequences we get the symbolic sequences of the cycles related to -fins whose trunks are 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 each -fin region for , , 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 , where is the number of symbols in ), and two other boundaries are related to two more BCBs of the cycle. Each -fin region has only three boundaries, namely, one DFB boundary and two BCB boundaries. The DFB boundary is defined by the condition , while the BCB boundaries are obtained using the skew tent map as border collision normal form.

In Appendix B, some basic formulas related to the boundaries of a fin region are given, while below, as an example, we describe the - and -fins of the region shown in Figure 9. There are two -fins, and , contiguous to the trunk regions . As we mentioned above, the -fins are exceptional being confined by three boundaries instead of four. In particular, the fin is confined by the BCB and DFB boundaries respectively, and the boundary defined by . The fin is confined by the BCB and DFB boundaries respectively, and the boundary .

Next, let us consider the -fins, and , which are also contiguous to the trunk regions . Applying the formulas given in Appendix B we get that the boundaries of the fin are given by the BCB curves satisfying the equations given below: Note that ; that is, the DFB of the cycle occurs simultaneously with the BCB of the cycle (see Figure 9). The fin is confined by the following boundaries: and, as already remarked in the previous case, here also we have (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.

Appendices

A.

Consider a generic family of 1D continuous piecewise linear bimodal maps defined as where , , , .

To describe the period adding structure observed in the parameter space of map consider first the basic cycles and belonging to the families and of complexity level one defined in (9). As stated in [15], the periodicity regions and of map are defined as where with These formulas are valid for the map given in (2) substituting

The periodicity regions related to cycles belonging to the families and of complexity level two given in (11) have the form respectively, where and are defined in (10), and , , , and are defined as

To get the periodicity regions of the cycles belonging to the families and   (12) one has to interchange the indices and , as well as to change the inequality signs to the opposite ones in (A.8) and (A.7), respectively.

The periodicity regions of a complexity level are obtained from the periodicity regions of the level by a recursive algorithm described in detail in [15].

B.

Let us recall the analytic representation of the boundaries of -fins of map given in (A.1), which have the common boundary with the trunk region . Consider first the fins contiguous to , corresponding to the attracting cycles , , . It is proved in [15] that any fin , if it exists (as a fin region may be an empty set), is confined by the boundaries defined as follows: where the BCB boundary which corresponds to given in (A.2) is common with the trunk region . The boundary corresponds to the degenerate flip bifurcation, and and are related to two more BCBs of .

The -fin attached to the trunk region on the other side and related to the cycle , , is confined by the boundaries defined as follows: Here the boundary corresponds to given in (A.2).

The boundaries of the fins contiguous to the trunk regions are obtained interchanging the indexes and in the above expressions. The fin regions of the higher complexity levels are obtained using the map replacement technique (see [15] for details).

Applying these results to map given in (2), we get that the fin region , associated with the cycles , , , is confined by the following boundaries: Accordingly, the -fin related to the cycle , , is confined by the boundaries defined as

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work of V. Avrutin was supported by the Marie Curie International Fellowship within the 7th European Community Framework Programme, project “Multiple-discontinuity induced bifurcations in theory and applications.” The work of I. Sushko has been performed under the auspices of COST Action IS1104 “The EU in the new complex geography of economic systems: models, tools and policy evaluation.”