Synchronization of a Supply Chain Model with Four Chaotic Attractors

Zheng, Jing; Zhang, Qiongxin; Xu, Qi; Xu, Fei; Shi, Victor

doi:https://doi.org/10.1155/2022/6390456

Discrete Dynamics in Nature and Society

On this page

Abstract Introduction Model Formulation Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Discrete Dynamics in Economic and Business Systems

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 6390456 | https://doi.org/10.1155/2022/6390456

Synchronization of a Supply Chain Model with Four Chaotic Attractors

Jing Zheng,^1,2Qiongxin Zhang ,³Qi Xu,¹Fei Xu,⁴and Victor Shi⁴

Academic Editor: Lei Xie

Received27 Dec 2021

Accepted15 Feb 2022

Published08 Mar 2022

Abstract

In this article, we construct a three-stage supply chain model using a system of differential equations to reveal the interplay among producers, distributors, and end customers. On the one hand, information about the products and production is exchanged between each stage in the supply chain system. Such information affects the behaviors of each stage. On the other hand, the transport of products between the stages of the supply chain affects the behavior of the system. Our findings are summarized as follows: First, we find that the supply chain model is a system with four chaotic attractors. Second, we explore the synchronization of such a chaotic system. Third, according to the characteristics of the chaotic system, we design a variety of simple control laws to realize the synchronization of two chaotic systems with the same structure. These control laws for the chaotic system are proved to realize local asymptotic synchronization or global exponential synchronization. Numerical simulations are conducted to confirm that the designed controls work well.

1. Introduction

Since the discovery of the first chaotic system, the Lorenz chaotic system [1], chaos theory has received extensive attention. Chaotic systems display complicated dynamical phenomena, sensitive dependence on initial conditions, and nonperiodic and pseudo-random behaviors. Scholars have discovered various classical chaotic systems, such as the Rössler chaotic system [2], Chua’s circuit [3], Chen’s chaotic system [4], and Lu’s chaotic system. The Lorenz system, Chen’s system, and Lu’s system belong to the family of Lorenz chaotic systems. Chaotic systems are widely used in science and industry. For example, the Lorenz chaotic attractor is generated by a weather forecast model [1], and the Rössler chaotic system can be applied to modeling chemical reactions [2]. Chua’s circuit confirmed the existence of chaos in electrical circuits [3]. Investigations have confirmed that there exists chaos in a coagulation model under the constraint of fixed density [5]. Chaos has been discovered in physics, chemistry, computer science, medicine, and engineering [6–9].

Different chaotic systems display different chaotic attractors with unique topologies. The Lorenz chaotic system exhibits an attractor with two wings, the Rössler system has a single-scroll chaotic attractor, and Chua’s circuit exhibits a double-scroll chaotic attractor. More complex chaotic systems with multiwing attractors and multiscroll attractors are designed in the literature [9–15]. The highly unpredictable behaviors of chaotic systems pose a challenge to their control and synchronization. In view of this, chaotic systems have been considered uncontrollable for a long time. Since the pioneering works of Pecora and Carroll [16] and Ott et al. [17], the control and synchronization of chaos have attracted the attention of researchers from all over the world. A variety of control and synchronization schemes are reported in the literature [13, 18, 19].

In recent years, researchers have studied various supply chain systems using mathematical modeling. Hou et al. [20] established a model to study production, inventory, and distribution processes in a supply chain system. Amorim et al. [21] applied this supply chain model to perishable products. Yuan and Hwong [22] discussed the influence of customer behavior and purchase decisions on the stability of a supply chain. Kumar and Tiwari [23] obtained the optimal factory location by studying the risk of safe inventory and operating inventory in a supply chain. Studies show that a variety of supply chains display unpredictable and chaotic behaviors [22, 24]. In this work, we study a supply chain system containing three stages, between which information is exchanged [25]. The behavior of the supply chain is influenced by the transportation of products and information exchange between these stages [26]. Based on systems of ordinary differential equations, we construct a model to study the behavior of the supply chain. We find that the interactions between producers, retailers (distributors), and consumers can lead to chaos. We found that this supply chain system has the same structure as the chaotic system proposed by Nwachioma–Pérez Cruz [27]. We also consider the chaos synchronization of the model. We design several control strategies to realize the chaos synchronization of the system. Mathematical analysis is carried out to prove the effectiveness of these designed control laws. Furthermore, we conduct numerical simulations to illustrate their effectiveness.

The rest of the manuscript is organized as follows: Section 2 considers the modeling of a supply chain. In Section 3, we introduce fundamental theorems and definitions related to chaos synchronization. We present the designed control laws in this section. Mathematical analyses are provided to prove the control schemes. We use numerical simulations to illustrate the application of the designed control laws. The main conclusions are presented in Section 4.

2. Model Formulation

In this paper, we establish a supply chain model using a system of differential equations to study a supply chain system with three layers. This supply chain includes producers, distributors, and consumers, and the relationships between them are shown in Figure 1. In this system, the producers and distributors obtain information about the demand for the product from consumers. Distributors place orders based on consumer demand and current inventory. The order is passed on to the producers. Then, the producers make decisions about production based on the orders placed by distributors and information from consumers. The producers send the products to distributors, and then, the distributors sell the products to the consumers.

Here, we use x to denote the number of products produced by the manufacturers, y the quantity demanded by retailers, and z the quantity demanded by consumers. We assume that the rate of change in the quantity of products produced by the manufacturers is positively correlated with the retailer’s demand and the demand of the consumers. The defective rate of the product is a1. We use a differential equation to model the interactions between the quantity of products produced by the manufacturers, the retailer’s demand, and the demand of the consumers. Then, we have

The rate of change of the demand of distributors has a negative correlation with the quantity of products produced by the manufacturers. That is to say, when producers produce enough products, distributors will not place large orders because they know that there is a sufficient quantity of products. When manufacturers reduce production, distributors start to place more orders due to the worry of insufficient products to sell to consumers. The rate of change of the demand of distributors has a positive correlation with the current demand for distributors. We then obtain the following differential equation:

In this model, the rate of change of the demand of consumers increases with the quantity of products produced by manufacturers and the quantity demanded by retailers. The rate of change of the demand of consumers will decrease with the increase of the demand of consumers. We then get

Summarizing the above discussion, we construct the following model:where a₁ = 1, a₂ = 1, a₃ = 2.3, a₄ = 2, a₅ = 1, a₆ = 6, a₇ = 1, and a₈ = −0.25. We notice that the model has the same structure as the chaotic system proposed by Nwachioma and Pérez-Cruz [27]. Nwachioma and Pérez-Cruz showed that system (1) can be synchronized. Here, we use numerical simulation to verify the chaotic behavior (see Figure 2). As shown in Figure 2, the dynamic behavior of the system shows a chaotic trajectory, which is sensitive to initial conditions. To demonstrate its sensitivity to initial conditions, we simulate the system with different initial conditions. For initial conditions (1.2, 1.0, and 3.0), system (4) displays the blue-colored trajectory. If we choose the initial conditions (0.1, 0.2, or 0.01), the system displays the red-colored trajectory. As can be seen from Figure 1, the two trajectories are completely different. System (4) displays complicated dynamical behavior within a wide range of parameter values. Here, we present numerical simulation results to confirm the complex behavior of the model. Figure 3 shows that system (4) has chaotic behavior for a₁ (0, 5.5) and a₁ (7.8, 8.5). Figure 4 shows that the system has chaotic behaviours for a₂ (0, 0.42) and a₁ (0.6, 1.75).

3. Synchronization of the Chaotic System with Four Attractors

The main objective of this section is to establish chaos control laws to achieve the synchronization of system (4). In view of the sensitivity of chaotic systems to initial conditions, two chaotic systems with the same structure will exhibit completely different behaviors due to different initial values. With appropriate controls, such two systems can be synchronized. Here, we name the two chaotic systems as the driving system and the receiving system, respectively. The driving system takes the form ofwhere subscript d stands for the driving system. The corresponding receiving system is then given bywhere subscript r represents the receiving system. Here, u₁, u₂, and u₃ are controls to be constructed. We thus obtain the error systemwhere e_x = x_d − x_r, e_y = y_d − y_r, and e_z = z_d − z_r. In order to make the behaviors of driving system (5) and receiving system (6) exactly the same (synchronized), we designed appropriate controls (u₁, u₂, u₃, u₄). Obviously, if the zero solution of system (7) is stabilized, then the driving system and the receiving system are synchronized.

Definition 1. x_d(0), y_d(0), z_d(0) R³ and x_r(0), y_r(0), z_r(0) R³, if the zero solution of error system (7) is locally asymptotically stable, driving system (5) and receiving system (6) are said to be locally asymptotically synchronized. Now, we construct a set of controls to realize the local asymptotic synchronization between systems (5) and (6).
Next, we will design a set of control laws and present analytic results.

Theorem 1. The following control laws stabilize the zero solution of error system (7) locally asymptotically, i.e., the set of control laws realizes the local asymptotic synchronization between systems (5) and (6):where k₁₂ > a₄.

Proof. By linearizing error system (7) with the control law proposed in Theorem 1 at equilibrium (0, 0, 0), we then haveIt is easy to see that the eigenvalues of J₁ are λ₁ = −a₁, λ₂ = a₄ − k₁₂, and λ₃ = −a₆. Note that all of these eigenvalues are negative. It thus follows that the zero solution of system (7) is locally asymptotically stable.
Therefore, control law (8) synchronizes systems (5) and (6). Next, numerical simulation is carried out to illustrate the application of the control law in Theorem 1. Here, we use numerical simulation to explore the behavior of the error system. As shown in Figure 5, the zero solution of the error system is stable under the designed control. Thus, the behavior of the driving system and the receiving system has achieved complete synchronization. Next, we design a control law to achieve global exponential synchronization.

(a)

(b)

(c)

Definition 2. x_d(0), y_d(0), z_d(0) R³ and x_r(0), y_r(0), z_r(0) R³, if the zero solution of error system (7) is globally exponentially stable, then driving system (5) and receiving system (6) are said to be globally exponentially synchronized.

Theorem 2. The control lawwhere k₂₂ > a₄, applied to receiving system (6), stabilizes the zero solution of (7) globally exponentially. That is to say, systems (5) and (6) are globally exponentially synchronized under control law (10).

Proof. For error system (7) with control (10), we construct the positive definite, radially unbounded Lyapunov function.We then calculate the derivative of V along the trajectory of system (7) with control (10) to getThus, the zero solution of system (7) is globally exponentially stabilized, i.e., driving system (5) and receiving system (6) are globally exponentially synchronized. Now, we conduct numerical simulations to show the applicability of the control law proposed in Theorem 2. In order to show the synchronization between the driving system and the receiving system, we plot the time histories of the two systems in the same figure. Here, the initial condition for the drive system is x_d(0) = 0.3, y_d(0) = −6, and z_d(0) = −0.1, and the initial condition of the receiving system is x_r(0) = 2, y_r(0) = 3, and z_r(0) = 5. We choose k₂₂ = 4, which satisfies the conditions of Theorem 2. The time history of driving system (5) is plotted in blue, and the time history of corresponding receiving system (6) is plotted in red. As shown in Figure 6, as t ⟶ ∞, the driving system and the receiving system are synchronized.

(a)

(b)

(c)

Figure 6

Time histories of driving system (5) (blue curve) and receiving system (6) (red curve) under control law (10) when a₁ = 1, a₂ = 1, a₃ = 2.3, a₄ = 2, a₅ = 1, a₆ = 6, a₇ = 1, and a₈ = −0.25 with initial conditions x_d(0) = 0.3, y_d(0) = −6, z_d(0) = −0.1, x_r(0) = 2, y_r(0) = 3, and z_r(0) = 5 under control law (10) with k₂₂ = 4. (a) Time histories of x_d and x_r for systems (5) and (6) under control law (10). (b) Time histories of y_d and y_r for systems (5) and (6) under control law (10). (c) Time histories of z_d and z_r for systems (3.1) and (6) under control law (10).

Theorem 3. Apply the control lawwhere k₃₃ > a₄, to receiving system (6). Then, the zero solution of (7) is globally exponentially stabilized. That is to say, we realize the global exponential synchronization of the drive system (5) and the receiving system (6) with control law (13).

Proof. Construct the positive definite, radially unbounded Lyapunov functionfor error system (7) with control (13). Calculating the derivative of V along the trajectory of system (7) with control (13) yieldsThus, the zero solution of system (7) is globally exponentially stabilized. Thus, we use the control law to link driving system (5) and receiving system (6) to achieve global exponential synchronization. We use numerical simulation to verify the correctness of the designed control law in Theorem 3. Here, we compare the time history of driving system (5) with that of the driven system by plotting them in the same figure. We choose k₃₃ = 5, which guarantees that the condition of Theorem 3 is satisfied. The initial condition for the driving system is x_d(0) = 3, y_d(0) = 2, and z_d(0) = 3, and the initial condition of the receiving system is x_r(0) = 0.3, y_r(0) = 0.5, and z_r(0) = 0.5. The time history of the driving system is plotted in blue, and the time history of the receiving system is plotted in red. As shown in Figure 7, though the initial conditions of the drive system and those of the driven systems are different, under control law (13), the behavior of the driving system is synchronized with that of the driven system.

(a)

(b)

(c)

Figure 7

Time histories of driving system (5) (blue curve) and receiving system (6) (red curve) under control law (10) when a1 = 1, a2 = 1, a3 = 2.3, a4 = 2, a5 = 1, a6 = 6, a7 = 1, and a8 = −0.25 with initial conditions x_d(0) = 3, y_d(0) = 2, z_d(0) = 3, x_r(0) = 0.3, y_r(0) = 0.5, and z_r(0) = 0.5 under control law (10) with k₃₃ = 5. (a) Time histories of x_d and x_r for systems (5) and (6) under control law (10). (b) Time histories of y_d and y_r for systems (5) and (6) under control law (10). (c) Time histories of z_d and z_r for systems (5) and (6) under control law (10).

4. Conclusions

This paper established a mathematical model to study a supply chain with a three-stage structure. This model has the same structure as the chaotic system proposed by Nwachioma and Pérez-Cruz [27]. We studied the chaotic behavior of the system using numerical simulations. The chaotic behaviors of the model are the results of the interactions of various stages in the system. The unpredictability of chaotic systems poses difficulty in their application. The realization of chaos synchronization provides a theoretical basis for its further application. We design a series of controls to realize the synchronization of the chaotic system. We prove the correctness of these control laws by means of mathematical analysis. To illustrate that these control laws are applicable, we perform numerical simulations. These numerical simulation results show that our control laws are effective. Chaos in supply chains often causes loss and risk. The realization of chaos synchronization provides a theoretical basis for controlling chaos in supply chains. In particular, when there are several interconnected supply chains in an economic system, synchronizing these chains can improve their efficiency, reduce uncertainty, and minimize potential risks and disruptions.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the General Scientific Research Project of Zhejiang Education Department (Grant Number: Y202146964) and the Wenzhou Social Sciences Planning Project (Grant Number: 21wsk031), China. Financial support for this research was also provided by CN through the Centre for Supply Chain Management, Wilfrid Laurier University, Canada.

References

E. N. Lorenz, “Deterministic nonperiodic flow,” Journal of the Atmospheric Sciences, vol. 20, no. 2, pp. 130–141, 1963.
View at: Publisher Site | Google Scholar
O. E. Rössler, “An equation for continuous chaos,” Physics Letters A, vol. 57, pp. 397-398, 1976.
View at: Google Scholar
L. Chua, M. Komuro, and T. Matsumoto, “The double scroll family,” IEEE Transactions on Circuits and Systems, vol. 33, no. 11, pp. 1072–1118, 1986.
View at: Publisher Site | Google Scholar
G. Chen and T. Ueta, “Yet another chaotic attractor,” International Journal of Bifurcation and Chaos, vol. 09, no. 07, pp. 1465-1466, 1999.
View at: Publisher Site | Google Scholar
M. Escobedo and F. Pezzotti, “Propagation of chaos in a coagulation model,” Mathematical Models and Methods in Applied Sciences, vol. 23, no. 06, pp. 1143–1176, 2013.
View at: Publisher Site | Google Scholar
M. Hasegawa, K. Nanbu, and K. Iwata, “Chaotic motion of two molecules in a box,” Mathematical Models and Methods in Applied Sciences, vol. 03, no. 05, pp. 693–710, 1993.
View at: Publisher Site | Google Scholar
M. E. Yal¸cin, J. A. K. Suykens, and J. Vandewalle, “On the realization of n-scroll attractors,” Proceedings - IEEE International Symposium on Circuits and Systems, vol. 5, pp. 483–486, 1999.
View at: Google Scholar
V. H. Hoang and C. Schwab, “N-term wiener chaos approximation rates for elliptic PDEs with lognormal Gaussian random inputs,” Mathematical Models and Methods in Applied Sciences, vol. 24, no. 04, pp. 797–826, 2014.
View at: Publisher Site | Google Scholar
C. Zhang and S. Yu, “Generation of grid multi-scroll chaotic attractors via switching piecewise linear controller,” Physics Letters A, vol. 374, no. 30, pp. 3029–3037, 2010.
View at: Publisher Site | Google Scholar
J. A. K. Suykens, A. Huang, and L. O. Chua, “A family of n-scroll attractors from a generalized Chua’s circuit,” AEU Int. J. Electron. Commun., , vol. 51, pp. 131–138, 1997.
View at: Google Scholar
W. K. S. Tang, G. Q. Zhong, G. Chen, and K. F. Man, “Generation of n-scroll attractors via sine function,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 48, no. 11, pp. 1369–1372, 2001.
View at: Publisher Site | Google Scholar
F. Xu, P. Yu, and X. Liao, “Global analysis on n-scroll chaotic attractors of modified chua’s circuit,” International Journal of Bifurcation and Chaos, vol. 19, no. 01, pp. 135–157, 2009.
View at: Publisher Site | Google Scholar
F. Xu, P. Yu, and X. Liao, “Synchronization and stabilization of multi-scroll integer and fractional order chaotic attractors generated using trigonometric functions,” International Journal of Bifurcation and Chaos, vol. 23, no. 08, p. 1350145, 2013.
View at: Publisher Site | Google Scholar
C. Zhang and S. Yu, “On constructing complex grid multi-wing hyperchaotic system: Theoretical design and circuit implementation,” International Journal of Circuit Theory and Applications, vol. 41, no. 3, pp. 221–237, 2013.
View at: Publisher Site | Google Scholar
F. Xu, “Integer and fractional order multiwing chaotic attractors via the chen system and the lü system with switching controls,” International Journal of Bifurcation and Chaos, vol. 24, no. 03, p. 1450029, 2014.
View at: Publisher Site | Google Scholar
L. M. Pecora and T. L. Carroll, “Synchronization in chaotic systems,” Physical Review Letters, vol. 64, no. 8, pp. 821–824, 1990.
View at: Publisher Site | Google Scholar
E. Ott, C. Grebogi, and J. A. Yorke, “Controlling chaos,” Physical Review Letters, vol. 64, no. 11, pp. 1196–1199, 1990.
View at: Publisher Site | Google Scholar
X. Liu, K.-L. Teo, H. Zhang, and G. Chen, “Switching control of linear systems for generating chaos,” Chaos, Solitons & Fractals, vol. 30, no. 3, pp. 725–733, 2006.
View at: Publisher Site | Google Scholar
X. Wu, J. Li, and G. Chen, “Chaos in the fractional order unified system and its synchronization,” Journal of the Franklin Institute, vol. 345, no. 4, pp. 392–401, 2008.
View at: Publisher Site | Google Scholar
J. Hou, A. Z. Zeng, and L. Zhao, “Achieving better coordination through revenue sharing and bargaining in a two-stage supply chain,” Computers & Industrial Engineering, vol. 57, no. 1, pp. 383–394, 2009.
View at: Publisher Site | Google Scholar
P. Amorim, H.-O. Günther, and B. Almada-Lobo, “Multi-objective integrated production and distribution planning of perishable products,” International Journal of Production Economics, vol. 138, no. 1, pp. 89–101, 2012.
View at: Publisher Site | Google Scholar
X. Yuan and H. B. Hwarng, “Managing a service system with social interactions: stability and chaos,” Computers & Industrial Engineering, vol. 63, no. 4, pp. 1178–1188, 2012.
View at: Publisher Site | Google Scholar
S. K. Kumar and M. K. Tiwari, “Supply chain system design integrated with risk pooling,” Computers & Industrial Engineering, vol. 64, no. 2, pp. 580–588, 2013.
View at: Publisher Site | Google Scholar
E. E. Mahmoud, P. Trikha, L. S. Jahanzaib, and O. A. Almaghrabi, “Dynamical analysis and chaos control of the fractional chaotic ecological model,” Chaos, Solitons & Fractals, vol. 141, p. 110348, 2020.
View at: Publisher Site | Google Scholar
Z. Liu, K. W. Li, J. Tang, B. Gong, and J. Huang, “Optimal operations of a closed-loop supply chain under a dual regulation,” International Journal of Production Economics, vol. 233, p. 107991, 2021.
View at: Publisher Site | Google Scholar
X. Xu, C. Wang, and P. Zhou, “GVRP considered oil-gas recovery in refined oil distribution: from an environmental perspective,” International Journal of Production Economics, vol. 235, p. 108078, 2021.
View at: Publisher Site | Google Scholar
C. Nwachioma and J. H. Pérez-Cruz, “Analysis of a new chaotic system, electronic realization and use in navigation of differential drive mobile robot,” Chaos, Solitons & Fractals, vol. 144, p. 110684, 2021.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Jing Zheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

297

Downloads

481

Citations