Dynamics of Fractional-Order Digital Manufacturing Supply Chain System and Its Control and Synchronization

: Digital manufacturing is widely used in the production of automobiles and aircrafts, and plays a profound role in the whole supply chain. Due to the long memory property of demand, production, and stocks, a fractional-order digital manufacturing supply chain system can describe their dynamics more precisely. In addition, their control and synchronization may have potential applications in the management of real-word supply chain systems to control uncertainties that occur within it. In this paper, a fractional-order digital manufacturing supply chain system is proposed and solved by the Adomian decomposition method (ADM). Dynamical characteristics of this system are studied by using a phase portrait, bifurcation diagram, and a maximum Lyapunov exponent diagram. The complexity of the system is also investigated by means of SE complexity and C 0 complexity. It is shown that the complexity results are consistent with the bifurcation diagrams, indicating that the complexity can reﬂect the dynamical properties of the system. Meanwhile, the importance of the fractional-order derivative in the modeling of the system is shown. Moreover, to further investigate the dynamics of the fractional-order supply chain system, we design the feedback controllers to control the chaotic supply chain system and synchronize two supply chain systems, respectively. Numerical simulations illustrate the effectiveness and applicability of the proposed methods.


Introduction
Fractional-order calculus is an extension and generalization of integer-order calculus, and fractional-order differential equations are obtained by acting the fractional-order differential operators on integer-order ones. Although the theory of fractional-order calculus had been proposed 300 years ago, and it developed slowly for a long time due to its lack of practical engineering application background and its large computational size, it now attracts extensive attention due to increases in computational power. It was not until Mandelbrot [1] first pointed out the fact that fractional dimension exists in nature and in many fields of science and technology in 1983 that the application of fractional-order calculus attracted wide attention and developed rapidly. Fractional order systems have long memory of history data [2,3]. The classical derivative in one point is only affected by the information in the local neighborhood of that point, while the fractional derivative in one point is affected by the combination of all of the information of the model in historical moments. Therefore, compared with integer-order calculus, fractional-order calculus produces more accurate and reliable results for modelling of real-world systems due to the memory effect. Recently, fractional-order systems have received special attention from researchers in various fields such as physics [4], biology [5], neural network [6], management, and economics [7][8][9][10][11], etc. On the other hand, chaos theory and applied research have developed rapidly and promoted its significant contribution in various scientific fields [12]. Nowadays, the control and synchronization of nonlinear systems are the focus of many research studies in a variety of fields [13][14][15]. Some effective controllers have been designed to control and synchronize the fractional-order chaotic systems such as the coupling controller [16], adaptive controllers [17], linear feedback controllers [18], sliding mode controllers [19], fractional order PID controllers [20,21], and so on.
Supply chain systems are complex dynamic systems with various uncertainties [22,23]. In recent years, some researchers [24][25][26][27][28] are studying nonlinear chaotic dynamics of supply chain systems and have obtained many investigations. Forrester [29] was the first to investigate the dynamics of supply chain systems and introduced the 'bullwhip effect'. The 'bullwhip effect' refers to the distortion of information in the process of transmitting information from downstream to upstream enterprises, and the distortion and gradual amplification of information in the process of transmitting information from the final customer to the original supplier in the supply chain, resulting in the phenomenon of cascading demand information. The 'bullwhip effect' can be explained by the chaos theory of dynamical systems. One of the most fundamental characteristics of chaos is that the trajectory of the system is very sensitive to the initial condition, i.e., even if the original state changes slightly, the final state of the system can be very different. Goksu et al. [30] constructed a supply chain model composed of manufacturers, distributors, and customers to achieve the synchronization and control of this chaotic supply chain management system. Gao et al. [31] proposed a new three-dimensional supply chain fractional-order difference game model composed of manufacturers, distributors, and retailers, and used the correlation theory of fractional-order difference to numerically analyze the complex dynamic behavior of this model, and discuss the effect of the output adjustment speed parameter on the dynamic behavior of the system. Recently, Yan et al. [32] integrated the computer-aided digital manufacturing process into the three-level supply chain which is composed of manufacturers, distributors, and retailers, and considered computer-aided digital design prior to the production by manufacturers, ultimately achieving synchronization and control of the system.
There are numerous attribute properties that cannot be described by the theory of integer-order calculus, so it is necessary to theoretically study the complexity of the supply chain system using the method of bifurcation and chaos of fractional nonlinear dynamics. In a supply chain system, the variables including demand, supply, and production have long memory properties, the fluctuation of which can lead to significant instability in the operation and delivery of the system. Thus, the traditional integer-order supply chain model has limitations to accurately show the operation of the system. Moreover, the prevalence of the bullwhip effect in the supply chain management increases the risks of production, supply, inventory management and marketing of suppliers, and even leads to chaos in them. However, there are few results on fractional-order supply chain systems.
Motivated by the above discussions, to be more realistic, in this study, a fractionalorder digital manufacturing supply chain system is investigated. Dynamics and complexity of this system with the variation of derivative orders and system parameter have been studied by means of bifurcation diagram and complexity measure algorithms. Furthermore, the nonlinear feedback controllers are designed to control and synchronize the chaos in this fractional order supply chain system, respectively. That is, the fractional order supply chain system has rich dynamic behaviors, and the evolution simulation is conducted for the influence of the change of fractional order and parameters on the demand order quantity, supply quantity, and digital manufacturing quantity of the supply chain enterprises. At the same time, the chaotic states appearing in the evolution process are synchronized and controlled to achieve the stability of the supply chain system. Therefore, this study contains both theoretical and practical guidance to eliminate chaotic dynamics that are unfavorable to the supply chain model. The rest of this paper is organized as follows. In Section 2, preliminaries and modeling are investigated. In Section 3, the dynamical behaviors of fractional-order digital Fractal Fract. 2021, 5, 128 3 of 17 manufacturing supply chain system are studied. In Section 4, controllers are designed to synchronize two identical systems. In Section 5, we consider the stabilization of the system. This paper ends with a conclusion in Section 6.

Preliminaries
In order to solve fractional-order calculus equations, various definitions of fractionalorder calculus have been proposed, among which the most common ones are the Grunwald-Letnikov (G-L) fractional-order calculus definition, the Riemann-Liouville (R-L) fractionalorder calculus definition, and the Caputo fractional-order calculus definition [33]. The algorithm based on the Caputo definition has clear physical meaning and is beneficial to solve the actual physical system, having more practical engineering applications. In this paper, we will only use the Caputo fractional derivative. Definition 1. The Caputo fractional derivative of order q is given by where n − 1 < q < n, a and t are numbers representing the limits of the operator C a D q t , and the symbol Γ(·) is the gamma function. Lemma 1 ([34]). Consider the following fractional-order system: where The equilibrium points of the above system are solutions to the equation f (x(t)) = 0. An equilibrium is asymptotically stable if all eigenvalues λ i of the Jacobian matrix J = ∂ f ∂x = ∂( f 1 , f 2 ,..., f N ) ∂(x 1 ,x 2 ,...,x N ) evaluated at the equilibrium satisfy |arg(λ i )| > πq 2 .
As for the fractional-order continuous systems, there are several different solution algorithms such as the frequency domain method (FDM) [35], the Adams-Bashforth-Moulton algorithm (ABM) [36], and the Adomian decomposition method (ADM) [37]. The ADM has higher accuracy and smaller computational error compared to the predictioncorrection algorithm and the Runge-Kutta algorithm [38], and it is used in this paper.
For a given fractional-order chaotic system with form of is the constant in the system, and D q t 0 is the Caputo fractional derivative operator. Then it can be divided into three parts as the form where m ∈ N, m − 1 < q ≤ m. Lx(t) and Nx(t) are the linear and nonlinear terms of the fractional differential equations respectively. Then, let J q t 0 is the inverse operator of D q t 0 , thus we have m − 1, and it involves the initial condition. By applying the recursive relation The analytical solution of the fractional-order system is given by where i = 1, 2, . . . , ∞, and the nonlinear terms of the fractional differential equations Nx(t) are evaluated by Nonlinear time series complexity measure is an important technique to analyze the dynamics of a chaotic system, and is currently a hot topic in the field of nonlinear research. Complexity measure of chaotic systems is to use some algorithms to measure how close the chaotic sequence is to the random sequence, and the complexity value is larger when the sequence is closer to the random sequence. There are several methods to measure the complexity of chaotic systems including statistical complexity measure (SCM), fuzzy entropy, sample entropy, spectral entropy (SE) [39], and C 0 algorithm [40]. Among these methods, SE and C 0 algorithms are proper methods to estimate the complexity of a time series accurately. So, SE and C 0 algorithms are used in this paper.

Modeling
In Ref. [32], Yan et al. designed a four-dimensional supply chain model, a demand-driven supply chain model based on digital manufacturing, focusing on the impact of digital manufacturing on manufacturers and then on the whole supply chain. The model is given by where x is the demand order quantity by retailers for the products, y is the supply quantity by distributors, z is the quantity of finished computer aided digital design, w is the quantity produced based on the digital design, m is the delivery efficiency of distributors, n is the rate of satisfying the retailer demand, r is the rate of information distortion for the demand by retailers for the products, d is the rate that the digital design is put into production, c is the conversion rate from digital design to product, k is the safety stock coefficient for manufacturers, and (θ m , θ n , θ r , θ d , θ c , θ k ) denote the parameters of the perturbations on the system. In addition, x < 0 indicates that the supply of distributors is less than retailer demand, and w < 0 indicates that there is no new production due to backlogs or returns.
In a supply chain system, the variables including demand, supply, and production have a long memory property; since the integer-order difference calculus has no longmemory effect, its theory is not suitable for studying the supply chain system with longmemory effects. Thus, the traditional integer-order supply chain model has limitations to accurately show the operation of the system. In this paper, the Caputo fractional derivative operator is applied in system (9), then the supply chain system with computer aided digital manufacturing process in the form of fractional-order differential equations are obtained as where q ∈ (0, 1) denotes the order of derivatives, in particular, model (10) degenerates to the integer-order supply chain differential equations when q = 1. Here, we take the same order of the fractional derivative in all equations, and choose the parameters m = 0.9, Since this paper focuses on the impact of digital manufacturing process on the whole supply chain system, the parameter d is taken as the control variable. Thus, the fractional-order supply chain system is proposed as follows Here, an approximate solution to the fractional-order value of the system (11) is Γ(6q+1) , j = 1, 2, 3, 4, and the detailed derivation of c 1 j , , c, · · · , c 6 j is given in Appendix A.

Dynamics Analysis of the Fractional-Order Supply Chain System
With the variation of derivative orders and system parameters, system (11) has complex dynamic behaviors such as periodic motions, period-doubling motions, and chaotic motions. If the system behaves chaotically, the supply chain system will be in a state of loss of control, which will lead to inventory, ordering, and supply chaos, and affect the decision making of supply chain enterprises at all levels, thus causing greater damage to the operation of the whole supply chain system. If the system appears in a periodic state, the supply chain system will be stable and the supply chain enterprises can make decisions based on the inventory status better. Therefore, studying the chaotic dynamic behaviors of the supply chain system can be an effective means to maintain the stability of the supply chain system and guide the scientific decision making of the supply chain enterprises.
In this section, the numerical solution of system (11) is obtained by means of ADM, and the chaotic dynamic behaviors of the system with the variation of the fractional order q and the parameter d are studied through bifurcation diagrams, maximum Lyapunov exponent diagrams, phase portraits, time series diagrams and complexity diagrams.

Dynamic Behaviors Analysis with Parameter d
If a manufacturing company produces digitally, it will be followed by digital design into production. Moreover, the larger the parameter d, the higher the level of digital production. In this subsection, we study the impact of the rate of digital design into production on the manufacturer and the supply chain. In system (11), let the fractional order q = 0.75, the initial conditions (x 0 , y 0 , z 0 , w 0 ) = (0.2, 0.2, 0.3, 0.32), and the parameter d is chosen as the critical variable. The bifurcation diagram and the maximum Lyapunov exponent diagram of system (11) with the parameter d varying from 0.5 to 5 are shown in Figures 1 and 2. The results indicate that the system shows the inverse period-doubling bifurcation, and as the parameter d decreases, the system goes from periodic state, after the period-doubling bifurcation, to chaotic state. When d ∈ [3. 33,5], the period one is appeared, and the period-doubling bifurcation occurs for d = 3.33. When d ∈ [2.37, 3.33), the period two is appeared, and the period-doubling bifurcation occurs for d = 2.37. When d ∈ [2.23, 2.37), the period four is appeared, and the period-doubling bifurcation occurs for d = 2.23. When d < 2.18, the system enters to chaotic state, which can be illustrated from the maximum Lyapunov exponent.  The numerical simulation results can reflect the actual situation of the supply chain system. When the rate that the digital design is put into production is high, it means that the products are digitally produced quickly, hence the replenishment demand of retailers can be met immediately, and consequently the demand of consumers can be met quickly, making the inventory system and the production of manufacturers stable. On the contrary, when the rate is low, the order demand of retailers cannot be met in time, so consumers may seek other alternatives, which will lead to overstock and chaos in production. Therefore, the strategy to increase the rate of digital design into production is feasible and it can make the whole supply chain system stable.
In order to observe the dynamic behaviors of system (11) directly, the phase portraits of the system with several different values of the parameter d are shown in Figure 3. When 0.7 d = , the system is chaotic and the chaotic attractor is presented in Figure 3a. The nu-  The numerical simulation results can reflect the actual situation of the supply chain system. When the rate that the digital design is put into production is high, it means that the products are digitally produced quickly, hence the replenishment demand of retailers can be met immediately, and consequently the demand of consumers can be met quickly, making the inventory system and the production of manufacturers stable. On the contrary, when the rate is low, the order demand of retailers cannot be met in time, so consumers may seek other alternatives, which will lead to overstock and chaos in production. Therefore, the strategy to increase the rate of digital design into production is feasible and it can make the whole supply chain system stable.
In order to observe the dynamic behaviors of system (11) directly, the phase portraits of the system with several different values of the parameter d are shown in Figure 3. When 0.7 d = , the system is chaotic and the chaotic attractor is presented in Figure 3a. The numerical analysis shows that the interaction between the quantity of demand from retailers, The numerical simulation results can reflect the actual situation of the supply chain system. When the rate that the digital design is put into production is high, it means that the products are digitally produced quickly, hence the replenishment demand of retailers can be met immediately, and consequently the demand of consumers can be met quickly, making the inventory system and the production of manufacturers stable. On the contrary, when the rate is low, the order demand of retailers cannot be met in time, so consumers may seek other alternatives, which will lead to overstock and chaos in production. Therefore, the strategy to increase the rate of digital design into production is feasible and it can make the whole supply chain system stable.
In order to observe the dynamic behaviors of system (11) directly, the phase portraits of the system with several different values of the parameter d are shown in Figure 3. When d = 0.7, the system is chaotic and the chaotic attractor is presented in Figure 3a. The numerical analysis shows that the interaction between the quantity of demand from retailers, the quantity of products available from distributors, the quantity of completed digital designs, and the quantity of production based on digital designs. When d = 2.9, the system is in the periodic four that is shown in Figure 3b. When d = 3.5, the period two is appeared, and it is shown in Figure 3c. When d = 4, the periodic one of the system is presented in Figure 3d. The phase diagrams are consistent with the bifurcation diagram and the maximum Lyapunov exponent diagram. Thus, the fractional-order supply chain system has rich dynamical properties when the parameter d of system (18)  is in the periodic four that is shown in Figure 3b. When 3.5 d = , the period two is appeared, and it is shown in Figure 3c. When 4 d = , the periodic one of the system is presented in Figure 3d. The phase diagrams are consistent with the bifurcation diagram and the maximum Lyapunov exponent diagram. Thus, the fractional-order supply chain system has rich dynamical properties when the parameter d of system (18) is varied. As shown in Figure 4, the 0 C complexity and SE complexity of system (18)  , indicating that production and inventory of manufacturers appear chaotic and difficult to predict; when 2.18 d ≥ , the value of the complexity is very small, indicating that it contributes to production planning of manufacturers and supply and inventory management of supply chain enterprises.  As shown in Figure 4, the C 0 complexity and SE complexity of system (18) have high values when 0.5 < d < 2.18, indicating that production and inventory of manufacturers appear chaotic and difficult to predict; when d ≥ 2.18, the value of the complexity is very small, indicating that it contributes to production planning of manufacturers and supply and inventory management of supply chain enterprises. is in the periodic four that is shown in Figure 3b. When 3.5 d = , the period two is appeared, and it is shown in Figure 3c. When 4 d = , the periodic one of the system is presented in Figure 3d. The phase diagrams are consistent with the bifurcation diagram and the maximum Lyapunov exponent diagram. Thus, the fractional-order supply chain system has rich dynamical properties when the parameter d of system (18) is varied. As shown in Figure 4, the 0 C complexity and SE complexity of system (18)  , indicating that production and inventory of manufacturers appear chaotic and difficult to predict; when 2.18 d ≥ , the value of the complexity is very small, indicating that it contributes to production planning of manufacturers and supply and inventory management of supply chain enterprises.

Dynamical Behaviors with q
In system (11), let the parameter d = 3.5, the initial conditions (x 0 , y 0 , z 0 , w 0 ) = (0.2, 0.2, 0.3, 0.32) and the parameter q is chosen as the critical variable to show the effect of fractional order to the behavior of chaotic system results. Figure 5 shows the bifurcation diagram with the order q as the bifurcation parameter. When q ∈ [0.47, 0.62], the system is mostly in a chaotic state. When q > 0.62, the chaos disappears and the system appears in a periodic state. In order to study the effect of different fractional orders on the supply chain system when manufacturers have a high rate of digital design into production, we consider choosing the fractional orders that appear in periodic state, so the orders are chosen as q = 1, q = 0.75, and q = 0.7, respectively.

Dynamical Behaviors with q
In system (11), let the parameter 3.5 d = , the initial conditions 0.62 , the system is mostly in a chaotic state. When 0.62 q > , the chaos disappears and the system appears in a periodic state. In order to study the effect of different fractional orders on the supply chain system when manufacturers have a high rate of digital design into production, we consider choosing the fractional orders that appear in periodic state, so the orders are chosen as 1 q = , 0.75 q = , and 0.7 q = , respectively.  Figure 6 shows the phase diagram between the demand by retailers and the digital production by manufacturers of system (11) with different fractional orders. As shown, when 1 q = , w fluctuates more at 0 x > and less at 0 x < , as the order q decreases to 0.75 and then to 0.7, w fluctuates gradually less at 0 x > and more at 0 x < . In practice, 0 x < indicates that the supply of distributors cannot meet the order demand of retailers. The analysis shows that with a small order q , then the production of the manufacturers fluctuates more when the supply from the distributors to the retailers is insufficient, while the production of the manufacturers fluctuates less when the supply from the distributors is sufficient. Thus, it indicates that the fractional-order system more accurately reflects the real-world supply chain system compare to the integer-order one. Figure 7 shows the corresponding time series diagram of the digital production of the manufacturers. As shown in the figure, the smaller the order q , the shorter the fluctuation period of the production, indicating that goods are transferred faster at all levels of the supply chain and can meet customer demand faster, resulting in a more stable inventory system and supply chain system.  Figure 6 shows the phase diagram between the demand by retailers and the digital production by manufacturers of system (11) with different fractional orders. As shown, when q = 1, w fluctuates more at x > 0 and less at x < 0, as the order q decreases to 0.75 and then to 0.7, w fluctuates gradually less at x > 0 and more at x < 0. In practice, x < 0 indicates that the supply of distributors cannot meet the order demand of retailers. The analysis shows that with a small order q, then the production of the manufacturers fluctuates more when the supply from the distributors to the retailers is insufficient, while the production of the manufacturers fluctuates less when the supply from the distributors is sufficient. Thus, it indicates that the fractional-order system more accurately reflects the real-world supply chain system compare to the integer-order one. Figure 7 shows the corresponding time series diagram of the digital production of the manufacturers. As shown in the figure, the smaller the order q, the shorter the fluctuation period of the production, indicating that goods are transferred faster at all levels of the supply chain and can meet customer demand faster, resulting in a more stable inventory system and supply chain system.
As shown in Figure 8, the C 0 complexity and SE complexity of system (11) are higher when 0.5 ≤ q ≤ 0.62, and the values of complexity both oscillate at lower values when q > 0.62. It shows that the fractional-order system has an order of magnitude higher complexity compared to the integer-order system.
From Figures 4 and 8, we find that the complexity results are consistent with the bifurcation diagrams, indicating that the complexity can reflect the dynamic characteristics of the fractional-order digital manufacturing supply chain system.  As shown in Figure 8, the . It shows that the fractional-order system has an order of magnitude higher complexity compared to the integer-order system.  As shown in Figure 8, the . It shows that the fractional-order system has an order of magnitude higher complexity compared to the integer-order system.  From Figures 4 and 8, we find that the complexity results are consistent with the bifurcation diagrams, indicating that the complexity can reflect the dynamic characteristics of the fractional-order digital manufacturing supply chain system.

Synchronization of the Chaotic Fractional-Order Digital Manufacturing Supply Chain System
From the analysis in Section 3, it is understood that the fractional-order digital manufacturing supply chain system is chaotic when the rate at which the digital design of the manufacturer is put into production is low. Without applying appropriate synchronization controllers to the system, the trajectories of the system with different initial values will exhibit different behaviors. In order to achieve synchronization between two supply chain systems with different initial conditions, the controllers can be designed to synchronize the two chaotic systems.
We establish the fractional-order error system by two chaotic fractional-order supply chain systems with different initial conditions, which are called the driving supply chain system and the responding supply chain system, respectively. Let d = 0.7, then the drive supply chain system is defined as follows And the response supply chain system is given as where u 11 (t), u 12 (t), u 13 (t), u 14 (t) are the controllers to be determined later.
Denote the error variables as where, e 1 , e 2 , e 3 , e 4 denote errors between the drive system and the response system. The fractional-order error system is obtained by subtracting the drive supply chain system from the response supply chain system To synchronize the drive system (12) and the response system (13), suitable controllers are chosen so that the fractional-order error system (15) is asymptotically stable at the origin. In this paper, we consider the controllers as The fractional-order error system (17) is asymptotically stable based on Lemma 1, which implies that synchronization between (12) and (13) will be realized. This completes the proof.
Furthermore, the effect of the fractional orders on the synchronization of two chaotic fractional-order supply chain systems with different initial conditions is studied. Let the parameter d = 0.7, then both the drive system and the response system are chaotic. Let the initial values of the drive system (x 10 , y 10 , z 10 , w 10 ) = (0.2, 0.2, 0.3, 0.32), and the initial values of the response system (x 20 , y 20 , z 20 , w 20 ) = (20,20,30,32). The fractional order q is chosen as q = 1, q = 0.75, and q = 0.7, respectively. The numerical solution of system (17) is obtained by the ADM, and the time series diagrams of system (17) with different fractional orders are shown in Figure 9. As shown in the figure, for two identical fractional-order chaotic supply chain systems (although there are differences in their initial states), the errors gradually converge to zero as time grows. This indicates that the two chaotic systems reach synchronization under the designed controllers, and the smaller the fractional order, the faster the synchronization. This implies that despite the factors such as the bullwhip effect, inventory can be synchronized faster when the order of the fractional-order supply chain systems is decreased, allowing the enterprises to reduce risk more effectively.  In order to see the comparison of the time histories of the drive system and the response system more visually and to further verify the synchronization, the fractional order q is chosen to be 0.75 and the time series of systems (12) and system (13) are plotted in the same figure, as shown in Figure 10. The green line in the figure represents the time series of the drive supply chain system under its initial conditions, and the yellow line represents the time series of the response supply chain system under its initial conditions, and the comparison reveals that the time series of the two systems overlap after 1.6 t = as time grows. Again, it shows that the designed synchronization controllers are effective and can synchronize the two systems. In order to see the comparison of the time histories of the drive system and the response system more visually and to further verify the synchronization, the fractional order q is chosen to be 0.75 and the time series of systems (12) and system (13) are plotted in the same figure, as shown in Figure 10. The green line in the figure represents the time series of the drive supply chain system under its initial conditions, and the yellow line represents the time series of the response supply chain system under its initial conditions, and the comparison reveals that the time series of the two systems overlap after t = 1.6 as time grows. Again, it shows that the designed synchronization controllers are effective and can synchronize the two systems.
in the same figure, as shown in Figure 10. The green line in the figure represents the time series of the drive supply chain system under its initial conditions, and the yellow line represents the time series of the response supply chain system under its initial conditions, and the comparison reveals that the time series of the two systems overlap after 1.6 t = as time grows. Again, it shows that the designed synchronization controllers are effective and can synchronize the two systems.
From Lemma 1, the equilibrium point E 0 is unstable, whereas the stability of E 1 and E 2 is determined by the value of the fractional order q.

Feedback Controllers
Considering the equilibrium point E 0 , to control chaos in fractional-order system (10), the feedback controllers u 11 (t), u 12 (t), u 13 (t), u 14 (t) are considered, then the controlled system is defined as For the fractional-order system (10), the following feedback controllers are designed in order to control the system asymptotically stably at the equilibrium point E 0 (0,0,0,0).
The simulation result shows the effectiveness of the designed control functions. The simulation result from Figure 11 shows the effectiveness of the designed control functions. This suggests that the designed controllers can control chaos in fractional-order system (10), and the smaller the fractional order, the faster the system is controlled to the equilibrium point. This implies that when the order of the fractional-order supply chain system is decreased, the system can reach a stable state more quickly.
In the above numerical results from Figure 5, we find that the derivative orders can change the bifurcation types and dynamics of the system. Figures 6 and 7 show that with different orders, supply chain companies have different decisions for inventory and management. Furthermore, Figure 9 indicates that the rate of synchronization is affected by the value of order. Figure 11 also shows that the rate of the system is controlled to the equilibrium point which is affected by the value of fractional order. It again indicates that the derivative order is very important in the fractional-order supply chain chaotic system.

Conclusions
In the present research, a fractional-order digital manufacturing supply chain system model was established. Through well-known tools and methods, including the phase portrait, bifurcation diagram, and maximum Lyapunov exponent diagram, the characteristics The simulation result shows the effectiveness of the designed control functions. The simulation result from Figure 11 shows the effectiveness of the designed control functions. This suggests that the designed controllers can control chaos in fractional-order system (10), and the smaller the fractional order, the faster the system is controlled to the equilibrium point. This implies that when the order of the fractional-order supply chain system is decreased, the system can reach a stable state more quickly.
In the above numerical results from Figure 5, we find that the derivative orders can change the bifurcation types and dynamics of the system. Figures 6 and 7 show that with different orders, supply chain companies have different decisions for inventory and management. Furthermore, Figure 9 indicates that the rate of synchronization is affected by the value of order. Figure 11 also shows that the rate of the system is controlled to the equilibrium point which is affected by the value of fractional order. It again indicates that the derivative order is very important in the fractional-order supply chain chaotic system.

Conclusions
In the present research, a fractional-order digital manufacturing supply chain system model was established. Through well-known tools and methods, including the phase portrait, bifurcation diagram, and maximum Lyapunov exponent diagram, the characteristics of the system were explored. The behavior of the system and the effects of the fractional-order derivative on the results of the system were displayed. It was demonstrated that the fractional-order system more accurately reflects the real-world supply chain system compared to the integer-order one. Furthermore, controllers were designed to synchronize two systems with different initial conditions. In addition, the equilibrium point analysis was carried out, and to suppress the chaotic behavior, feedback controllers were designed to stabilize the supply chain system. Our research results can help achieve more stable inventory systems and supply chain systems. In the future, we aim to further study the complexity evolution of fractional-order chaotic supply chain systems. Also, future research can be devoted to the control and synchronization of the proposed model through some more simple controllers with less parameters, such as fractional-order fuzzy controllers, fractional-order PID controllers, and fractional-order PI λ D controllers.

Data Availability Statement:
The data used during the study appear in the submitted article.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
According to the ADM, the system (11) can be represented as The nonlinear terms in the above equation can be decomposed as follows: A 0 xy = x 0 y 0 A 1 xy = x 1 y 0 + x 0 y 1 A 2 xy = x 2 y 0 + x 1 y 1 + x 0 y 2 A 3 xy = x 3 y 0 + x 2 y 1 + x 1 y 2 + x 0 y 3 A 4 xy = x 4 y 0 + x 3 y 1 + x 2 y 2 + x 1 y 3 + x 0 y 4 A 5 xy = x 5 y 0 + x 4 y 1 + x 3 y 2 + x 2 y 3 + x 1 y 4 + x 0 y 5 (A3) where the superscript in the decomposition formula is the number of ADM decompositions. The initial condition is According to the property of fractional-order calculus, then we obtain