Abstract

Inventory management is complex nonlinear systems that are affected by various external factors, including course human action and policy. We study the inventory management model under special circumstances and analyse the equilibrium point of the system. The dynamics of the system is analysed by means of the eigenvalue trajectory, bifurcations, chaotic attractor, and largest Lyapunov exponent diagram. At the same time, according to the definition of fractional calculus, the fractional approximate entropy is used to analyse the system, and the results are consistent with those of the largest Lyapunov exponent diagram, which shows the effectiveness of this method.

1. Introduction

Over the past few decades, chaos and its applications have attracted remarkable consideration in a variety of fields, including chemistry, ecology, and economy [1]. An economics management system is a typical discrete chaotic system, which has widely concerned scientists [28]. In 1993, Stiles and Levy proposed the concept of chaotic strategic management in [3]. In 1994, Feichtinger et al. [4] studied many chaotic systems in the field of management operations research, such as queuing system, inventory system, and planning and scheduling system. These systems are mainly characterised by the chaos of queue, inventory, and planning and scheduling under different management decision rules. Murphy [5] used the chaos theory as a model to study the management of public relations, such as problems and crises. After summarising the research status of chaos management, Joseph [6] pointed out that chaos management depends on change rules, which are a set of rules based on orderly or disorderly changes, adaptability, and new orderly emergence process.

The fractional calculus theory has a history of more than 300 years. Its history is almost the same as that of traditional integral calculus theory. Compared with an integer-order differential equation, a fractional-order differential equation has a memory effect [914]; therefore, it can describe these natural phenomena more accurately [9]. Many fractional-order chaotic systems, such as the fractional-order Lorenz system [15, 16] and the fractional-order Chen system [17, 18], have been studied [19]. However, most of the studies on chaos based on the theory of fractional calculus have been mainly limited to continuous-time chaotic systems, and the study of fractional-order discrete-time chaotic systems is insufficient. Thus, researchers turned their attention to apply the theory of discrete fractional calculus to discrete chaotic maps [20, 21]. It has been proven that there are abundant dynamical behaviours in the fractional chaotic maps [2227]. All of the above studies are based on chaotic systems, and other fields of fractional order have not been studied.

Many scientists have conducted considerable research on the complexity of time series such as economic management [8, 2325]. Grassberger and Procaccia [28] used the GP algorithm to analyse the complexity of Hénon mapping, which was an early report on the complexity analysis of chaotic systems. Balasubramanian et al. [29] used the Lempel–Ziv algorithm and the ApEn algorithm to analyse the complexity of logistic mapping and used a complexity measure to distinguish among the different states of the system. In [30], the PE algorithm was used to analyse the complexity of chaotic pseudorandom sequences. In [31], the researchers analysed the randomness of discrete chaotic systems by using the SymEn algorithm and showed that discrete chaotic systems can be used as random sources. All of the above systems are based on the theory of integer-order differential information. Can we further extend them to the theory of fractional calculus?

In this paper, we propose a new special nonlinear inventory management system [32, 33]. The discretisation of the proposed system was performed using a Matlab simulation. The rest of this paper is organised as follows. In Section 2, the equilibrium point and the 0-1 test of the 2D inventory management nonlinear systems are presented. In Section 3, the bifurcation diagram, phase diagrams, and Lyapunov exponent spectrum are used to analyse the dynamics of the system. In Section 4, the fractional-order approximate entropy of 2D nonlinear inventory management systems is analysed. Finally, the obtained results are summarised in Section 5.

2. Nonlinear Inventory Management Model Description

2.1. 2D Inventory Management Model

Inventory management is an important part of enterprise organisation and management activities. A class of nonlinear inventory management systems is as follows [2]:where s, p, q, and r are the system parameters, s represents the initial sales base, p represents the inventory fund transfer rate, q represents the product resource rate, r represents the inventory efficiency, xi represents the resources for sales in the time period i, yi represents the number of customers in the time period i, and zi represents the inventory capital of the company in the time period i. After normalising the parameters of the inventory management model, we obtained the following: , , and . Here, p = 0.43, q = 0.38, s = 0.11, and r = 0.72. The attractors of system (1) are shown in Figure 1.

The company had surplus inventory in stage i, and its main goal was to deal with this surplus inventory. At this time, the company shifted resources from production to sales. The new customers absorbed by the company in stage i + 1 mainly depended on the direct investment in sales in stage i. In addition, some old customers had to be added to the total customer volume. If the company’s inventory was insufficient and could not meet the requirements of the regular customers for timely pickup, some regular customers would disappear from the company’s customer list. indicates that the company had no resources available for sale. The absolute value of implies that it had to borrow funds and promote production to supplement the sales resources, denoted losing customers, and denoted that the sales were very strong and the products were in short supply.

The problems of convenient calculation and original speed, in this paper, considers a special case of inventory management system, P = 0, that is, there is no transfer of resources. The mathematical model of inventory management (1) is rewritten aswhere , , and . The definition of the parameters is the same as in formula (1). Here, q = 0.051, s = 0.58, and r = 0.72. The attractors of system (2) are shown in Figure 2.

2.2. Equilibrium Stability Analysis

The fixed point of system (2) can be obtained by solving the following equations:

Apparently, the fixed point can be obtained by setting . The Jacobian matrix J at the equilibrium set can be expressed as follows:

The systems’ characteristic equation corresponding to matrix J can be written as follows:

From equation (7), we inferred that when , at the fixed point p0, system (2) is unstable. When, at the fixed point p0, system (2) is stable. When , system (2) is at the bifurcation point.

2.3. 0–1 Test

We implemented the 0-1 test method. The 0-1 test is a novel test approach to determine whether a given deterministic nonlinear dynamic system is chaotic. It was proposed by Sun et al. [31] and has already been successfully tested for integer-order chaotic systems [31] and fractional order chaotic systems with delay time [8]. If a set of discrete data φ(n) (n = 1, 2, 3, ...) represents a one-dimensional observable dataset obtained from the iterative, then the following two real-valued functions can be defined:where , plotting the trajectories in the p-s diagram. In any case, the trajectory in the p-s graph can provide a simple visual measurement. If the bounded trajectory in the p-s plots is a regular cloud slave, then the unbounded trajectory similar to Brownian motion indicates that the system is a chaotic system. This method was used to study the y and z sequence of the 2D system (2). Its parameters were the same as those shown in Figure 2. The obtained p-s diagram is shown in Figure 3; it can be seen that the system is chaotic.

3. Analysis of Dynamic Characteristics

3.1. Dynamics with q Varying

We fixed s = 0.58 and r = 0.72 and let the derivative order q vary from 0.05 to 0.08 with a step size of 0.0002 and the initial values of state variables [y0, z0] = [1.1, 0.1]. Figure 4(a) shows that, upon the variation of q, it was obvious that the system chaos and periodic variation were repeated and crossed and were not a process from period to chaos directly. Obviously, we observed that was in a chaos state and the other regions were in a chaotic state with periodic windows. Secondly, when , the largest Lyapunov exponent (LLE) was positive; LLE was equal to 0 in the other regions, as shown in Figure 4(b).

3.2. Dynamics with r Varying

The chaotic system (2) had two system parameters besides q. In this study, the dynamic properties of the system parameters with r varied were analysed. Firstly, the bifurcation diagrams of the 2D system with r[0.5, 0.75] were investigated, as shown in Figure 5(a). Obviously, we observed that was in a chaos state and the other regions were in a chaotic state with periodic windows. Secondly, when , the LLE was positive; the LLE was equal to 0 in the other regions, as shown in Figure 5(b). Meanwhile, in eigenvalue trait of system (2) with r varying, as shown in Figure 6, r ranges from 0.5 to 0.75, and the eigenvalues are unstable, stable, and unstable processes, which are consistent with Figure 5. Eigenvalues are stable in the unit circle and unstable outside the unit circle.

For the case of s = 0.58, q = 0.051, and r = 0.6, the periodic state phase diagram is given in Figure 7(a) and its corresponding p-s plot is presented in Figure 7(c). For the case of s = 0.58, q = 0.051, and r = 0.73, the chaos state is shown in Figure 7(b) and the corresponding p-s plot is presented in Figure 7(d). Thus, the system had different states with different parameters, and chaos was verified.

4. Complexity Analysis of the System

Entropy measures the complexity of a time series by measuring the probability of generating new patterns in signals [3436]. The greater the probability of signal generation is, the greater is the value of the complexity of the sequence [29, 30].

4.1. Fractional Entropy Calculation

Information entropy can be calculated for both continuous and discrete signals and is defined as follows:

Let c = 1; then, according to the definition of a derivative,

Comparing equation (7) with equation (8), we obtain the following:

According to [7], a fractional-order derivative can be defined as follows:

Equation (10) was used to first derive h, and then, we set h⟶0 to seek the limit. On the basis of this fractional derivative, KARCI proposed the first type of fractional derivative [1, 3]:

When , equation (11) cannot be degenerated into equation (7), so there is an error. The fractional entropy is thus redefined as follows:

When , equation (12) can be degenerated into equation (7), but and ɑ is a complex number. Thus, if , ; then, . If , then is a real number. In practical applications, fractional entropy is redefined as follows:

4.2. Fractional-Order Approximate Entropy Calculation

Fractional-order entropy is combined with approximate entropy to obtain fractional-order approximate entropy (FAE). The specific steps are as follows.

In general, for the specific calculation method of the sample entropy of the time series {x(n)} composed of N data items, the specific steps are as follows:S1: according to the serial number to form a one-dimensional m vector sequence , where , . These vectors are m consecutive x values starting from point 1.S2: define the distance between Xm(i) and Xm(j), i.e., d[Xm(i) Xm(j)], as the difference between the maximum absolute values of the corresponding elements. S3: define the criterion of similarity r. Calculate the value of (, ) and its ratio to the total distance Nm, denoted as . S4: on the basis of equation (12), the average value of can be calculated as follows: S5: similarly, change m to m + 1, and repeat S1 to S3 to obtain . S6: theoretically, according to equation (13), the FAE complexity can be calculated as follows:

4.3. Complexity Analysis

The FAE analysis results of the inventory management system with varying q and r values are shown in Figure 8. Here, the same parameter settings and step size are used as those in the corresponding bifurcation diagrams and LLE. As demonstrated in these FAE plots, the results matched the corresponding bifurcation diagrams well. When the system was chaotic, higher complexity measurement results were obtained, whereas when the system was nonchaotic, the measured results were relatively small. Obviously, FAE provided different information on the dynamics of the system from that obtained from the bifurcation diagrams. Thus, we could determine the complexity variation trend of the system clearly.

As it is well known and quite obvious, we do not need chaos in the inventory management system, at least in most cases. This implies that the system is unstable and the inventory management situation has become unpredictable. We diagnosed the status of the inventory management system by means of the complexity method using the time series.

5. Conclusion

In the present study, the dynamic behaviour of a nonlinear inventory management system was investigated. The different aspects of the dynamic behaviour of the system, including the phase diagrams and their corresponding p-s plots, bifurcation diagrams, LLE with parameters and FAE complexity with parameters, were studied. We found that chaos and different states could be observed with different parameter settings. Meanwhile, the AEP complexity showed the changes in the complexity of the system. It also showed that the system had a wide region of high complexity in the corresponding parameter planes.

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 work was supported the Major Scientific and Technological Innovation Projects of Shandong Province (Grant no. 2019JZZY010111), Natural Science Foundation of Shandong Province (Grant no. ZR2017PA008), Key Research and Development Plan of Shandong Province (Grant no. 2019GGX104092), and Science and Technology Plan Projects of Universities of Shandong Province (Grant no. J18KA381).