Journal of Control Science and Engineering

Volume 2019, Article ID 6926342, 7 pages

https://doi.org/10.1155/2019/6926342

## Dynamics and Control on a Discrete Multi-Inventory System

School of Science, Tianjin University of Commerce, Tianjin 300134, China

Correspondence should be addressed to Li Xu; moc.361@lx_gnafieb

Received 12 February 2019; Revised 8 April 2019; Accepted 23 April 2019; Published 12 May 2019

Academic Editor: Michela Robba

Copyright © 2019 Yuzhen Wang and Li Xu. 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.

#### Abstract

In this paper, a discrete two item inventory model for deteriorating items with a linear stock dependent demand rate without regular supply is given first; then, regular supply is taken into account, which may be regarded as feedback control. An objective function is formed to calculate the net profit with respect to possible profits and possible loss. A necessary criterion for the steady state optimal control problem for optimizing the objective function subjected to the constraints given by the difference equations of the inventory is obtained. At last, an improved model is put forward.

#### 1. Introduction

Inventory management is one of the most important drivers for an effective enterprise. The mathematical modeling of inventory control for inventory management has caught up the mind of inventory researchers. Since the classical economic order quantity model which undoubtedly constitutes the cornerstone of inventory research [1], there are many models in the literature which can illustrate the behavior of inventory management, including not only the continuous-time framework (see [2–5], and the listed conference), but also discrete time one (for example, see [6, 7]), whose origins may be traced in the seminal work of Arrow [8].

In recent years, there are strong calls for the development of multi-item deterioration inventory system management [9]. In general, deterioration is defined as decay, damage, spoilage, evaporation, obsolescence, pilferage, loss of utility, or loss of marginal value of a commodity that results in decrease of usefulness from the original one. Because of deterioration, a certain fraction of the items is either damaged or decayed which will not satisfy the future demand of customers better. Deterioration for such items is continuous and constant or time-dependent. A number of research papers have already been published on the above type of items by some authors, the details of which can be found in the articles [10–17].

For example, some principles which occur in population dynamics can be applied to the multi -item deterioration inventory system in [17]. A two item competitive ordinary differential inventory model for deteriorating items with a linear stock dependent demand rate is given. By means of Pontryagin’s Maximum Principle, the author derived a necessary condition for the steady state optimal control problem in which objective function subjected to the constraints given by the competitive ordinary differential equations. The idea is interesting.

In most models, it is assumed that the system events are monitored continuously and that the time axis is continuous. In fact, the inventory systems, in which events occur and decisions are made during fixed length intervals, are characterized, then discrete time inventory models are relatively more reasonable. And, there is a growing research interest in discrete time queues (see [11–14]), which are mainly motivated by their applications in the computer and communication systems. In such a discrete time model, events occur at evenly-spaced points in time t =0,1,2, …N, where the integers represent any unit of time such as a day, a week, or a month. A discrete time system is a more natural manner to describe inventory dynamics.

The model of discrete system control is both a simulation model of inventory dynamics and an optimization model which can give optimal control according to the defined performance criterion [18]. Here, the key idea is to construct a function called the Hamiltonian by introducing the adjoint function to attach the difference equation to the objective function, which can be dealt with by the discrete version of Pontryagin’s Maximum Principle [19]. For the objective functional subject to the state difference equation with initial condition, the Hamiltonian can be used to find the optimal control. So far, to our knowledge, there have been very few corresponding research works focusing on discrete version of Pontryagin’s Maximum Principle on discrete inventory model for deteriorating items. It will be interesting and meaningful.

In this paper, corresponding to [17], a new discrete two item inventory model for deteriorating items with a linear stock dependent demand rate without regular supply is given first, then regular supply is taken into account. So this paper is organized as follows: in Section 2 we give some notations and assumptions. In Section 3, stability analysis is discussed for our discrete two item inventory model for deteriorating items without regular supply and with regular supply. In Section 4, an objective function is formed to calculate the net profit with respect to all possible profits and all possible loss. By means of the discrete version of Pontryagin’s Maximum Principle, a necessary condition for the steady state optimal control problem for optimizing the objective function subjected to the constraints given by the difference equations of the inventory is obtained. The next section is simulation discussion; a numerical example is used to illustrate how to find the optimal control and the corresponding state. Then an improved optimal problem with feedback control is put forward. The conclusion will be given in the last section.

#### 2. Notations and Assumptions

For convenience, we use the similar assumptions and notation used in [17] which are as follows:

##### 2.1. Notations

: the index of products : economic order quantity which enters the inventory at time point , : inventory level at discrete time point : number of discrete time periods (inventory cycle length) : rate of supply to which is inversely proportional to and may be , are scale parameter : linear demand rate at the instantaneous level of inventory at discrete time point are scale parameter : deterioration coefficient due to self-contact of (self-deterioration coefficient) : deterioration coefficient of due to presence of unit of at discrete time point (cross-deterioration coefficient) : order cost for one order : variable cost to maintain supply to at the rate are scale parameter : price per unit of : the perishable rate of product

##### 2.2. Assumptions

A product may have a lifetime which ends until its utility reaches zero.

Shortages are not allowed.

The percentage of profit does not depend on the initial inventory.

The rate of changing of holding cost is proportional to the square of ; it is given by , are constants, .

#### 3. Model and Stability Analysis

In this section, we will show the discrete model which can express the multi-item deterioration inventory. It can also be seen as the discrete version of the model in [17], which can be explained by some principles occurring in population dynamics. Furthermore, stability conditions will be discussed for our discrete two item inventory model without regular supply and with regular supply.

*Case 1 (without regular supply). *We give the following two item discrete inventory modelwhere are nonnegative, and the above system has four fixed points , , and , , is the only nonnegative one. The linearization of (1) about has the Jacobian matrixwhich has two eigenvaluesif or , then steady state is locally asymptotically stable; this means that the inventory is empty at last. Then we discuss the case with supply.

*Case 2 (with regular supply). *If the regular supply is taken into account, the following model can be givenClearly, system (4) has four possible steady states, i.e., , exclusion points , , and nontrivial coexistence point , whereand or

At the fixed point , the Jacobian matrix has the formand the corresponding eigenvalues of (6) areIf or , then steady state is locally asymptotically stable. That is to say, if rate of supply exceeds the sum of the natural deterioration rate and linear demand rate at the instantaneous level of inventory, the inventory will not be empty.

The linearization of (4) about has the Jacobian matrix and the corresponding eigenvalues areIf or and , then steady state is locally asymptotically stable. From this, we can conclude that if the above conditions are met, the number of will go to zero, and the number of will go to .

At the fixed point , the corresponding stable conditions can be given as follows:Similarly, the number of will go to zero, and the number of will go to

For the fixed point , the Jacobian matrix has the formwhere

The characteristic equation of is given bywhere

According to the Jury’s criteria [20], is asymptotically stable if the following conditions hold:then we have the following conclusion that if holds, the fixed point is asymptotically stable. We can say that and can coexist.

#### 4. Necessary Conditions for Dynamic Optimality

From Case 1 in Section 2, it is a fact that the inventory will be empty at last if . On the other hand, it is unreasonable if is too large. The adjustment for the system is necessary which will result in Case 2. From the control theory, the supply action may be regarded as a feedback control one. In other words, one may wish to move the undesired trajectory towards the desired orbit. Here we will alter the positions of positive equilibrium, but to obtain its stability. This motivates us to adjust the model form as follows.where and , and are control variables.

Under the condition the inventory will not be empty; we will calculate the net profit of the inventory. The objective function for the inventory problem should first be formed. Note that the percentage of profit does not depend on the initial inventory and the nonnegativity constraints hold for inventory variables, and percentage of profit maybe is found to be more for a lower initial stock than that for a higher initial stock [17]. Naturally we can formulate the optimization problem to maximize the net profit (performing index) taking the form

andwhere

For simplicity, we note

We can easily obtain the next results.

Theorem 1. *There exists an optimal control with corresponding states from (14) such that .*

Now we consider the Hamiltonian are the costate variables, then, from discrete Pontryagin’s Maximum Principle [19], we can derive the necessary conditions in the next theorem.

Theorem 2. *Given an optimal control and corresponding states and from (14), there exist adjoint functions and satisfying*

andwhere

*Proof. *Using the discrete version of Pontryagin’s Maximum PrincipleAdditionally, from the adjoint final conditions, we can getFurthermore,from which we can get

#### 5. Simulation and Discussion

In this section, we will show a simple example in order to intuitively illustrate how to find the optimal control firstly. We fix some of the constants in the parameter region , and let . The initial values will be set to , and let , by means of forward-backward sweep, the optimal solutions and optimal control can be displayed in Figures 1 and 2.