Open access peer-reviewed chapter

Control of Supply Chains

Written By

Kannan Nilakantan

Submitted: 27 August 2021 Reviewed: 19 September 2021 Published: 20 October 2021

DOI: 10.5772/intechopen.100523

From the Edited Volume

Control Systems in Engineering and Optimization Techniques

Edited by P. Balasubramaniam, Sathiyaraj Thambiayya, Kuru Ratnavelu and JinRong Wang

Chapter metrics overview

354 Chapter Downloads

View Full Metrics


This chapter aims to apply control principles in the discrete-time control of Supply Chains. The primary objective of the control is to keep the inventory levels (state variables) steady at their predetermined values and reduce any deviations to zero in the shortest possible time. The disturbances are induced by demand deviations from the planned/anticipated levels. The replenishment flows are the control variables. Thus, the control action is very similar to a “Linear Regulator with zero set-point”. A novel development in this chapter is the use of direct Operator methods to solve the system Difference Equations, thereby obviating the need for Z-Transforms, block diagrams and transfer functions of classical control theory. This chapter provides a novel application of control theory as well as an easier method of solution.


  • Supply Chain Dynamic Modeling
  • Supply Chain Control
  • Feedback Control
  • Linear Regulator Problem
  • Direct Operator Methods

1. Introduction

In the field of business, one of the most important constituents of a business system, and one which can give a business a cutting-edge over the others, is the supply chain (SC) of the business, and its effective operation and management.

A fundamental strategy associated with supply chains is that of a ‘Responsive Chain’. A responsive chain focuses on its ability to respond quickly to demand changes and meet them within the shortest possible time. A responsive chain, by its very definition, has necessarily to be able to respond swiftly to unanticipated changes in demand, and this is determined largely by the dynamics of the system. Under normal circumstances, good responsiveness can be achieved through the maintenance of adequate pipeline inventories throughout the system, and one of the factors impacting this decision in a major way, is again the dynamics of the system. Hence, given the impact of the dynamics of the SC system on its design and operation as mentioned above, it is imperative to include the tenets of dynamic analysis in its design and operation. Thus system-dynamic methods are of significant value and utility in both SC design and SC operation. The design aspects pertain to the setting of inventory levels and the operational aspects to choosing the appropriate system controls for good operational performance.

Thus, we are led naturally to the use of system-dynamic methods and concepts of control theory for effective control of a supply chain. In the following sections we explore the dynamic modeling and analysis concepts in effectively controlling a supply chain that would yield good performance characteristics.


2. Supply chain dynamics and notation

Supply chain dynamics essentially deals with the dynamic behavior and temporal variation of the inventories and flows in the system over time when subjected to demand disturbances.

We now look at a simple three stage serial supply chain, a schematic diagram of which is given below in Figure 1. Each stage represents an inventory storage facility in the chain. Stage 1 represents the raw material input storage to a manufacturing plant, and stage 2 the finished goods inventory storage at the manufacturing plant. Stage 3 represents the finished goods warehouse (W/H) at the downstream end of the chain which meets the demand for the finished product, and from which material is shipped out to customers. The inventory levels at each stage are the state variables of the system, and the material flows between the stages the control variables.

Figure 1.

A three-stage serial supply chain.

In constructing a dynamic model of a supply chain system, we adopt the standard schematic diagram (Figure 1), system notation, and model variables and equations below, as commonly found in control theory literature [1, 2, 3, 4, 5, 6, 7, 8, 9].

We use a discrete time representation of the SC system, as this is more in keeping with the prevailing practices in the SC industry, wherein the inventory levels (state variables) are recorded at the end of each period or day. In most cases the state variable records are updated at the end of each day. Hence, we take our period to be a single day. However, we need to also emphasize that this need not always be so and would be dependent entirely on the convenience of the SC practitioners. And hence the state variables are recorded at epochs corresponding to the end of each period. The flow variables however are aggregated and taken to occur within and up to the end of each period.

The inventory and flow variables are subscripted to indicate the stages in the chain. We now give the detailed representation below:

yik: inventory at stage i of the chain at time epoch k

with: i = 1, representing (the raw material) the upstream end of the production facility.

i = 2 representing (the finished goods) the downstream end of the production facility.

i = 3 representing (the finished goods) the warehouse.

qik: material flow in period (k-1,k] into stage i of the chain

r3'k: finished goods demand at warehouse in period (k-1,k]

The dynamic equations of the system are written using deviation variables, as under:



yi0k is the desired or planned inventory level at time k

qi0k is the desired or planned flow in period (k-1,k]

r30k is the anticipated or forecasted demand at the warehouse in period (k-1,k], based on which the desired inventory levels in the chain have been set.

xik is the inventory deviation at time k.

qik is the material flow deviation in period (k-1,k].

r3k is the deviation in demand observed at the warehouse in (k-1, k].

Under consistent units (in equivalent units of finished goods) for all inventory and flows in the system, the dynamic equations of the system are written as:


where for i = 3, q4k=r3k, the demand outflow from the warehouse. The system behavior is controlled by the replenishment policies followed in the system.

Now, since the demand is a stochastic variable, it can be split into two components, viz., a mean demand component, and a stochastic component which represents the random variations over and above the mean demand. The mean demand is what is predicted using forecasting techniques and is the planned offtake and the demand that the system is designed to meet in the normal course. Since prediction can never be exact, the residual variation is the stochastic component.

In such stochastic systems, the demand has an additional stochastic term, which is a white noise term, given by εkWN0σ2 for all k. The sequence of stochastic disturbance components over each period is an independent and identically distributed sequence of N0σ2 random variables.

The standard initial conditions for the system are: xik=0 for all k0 and r3k=0 for all k1, i.e. the system is at zero deviation at time k = 0, and the first deviation in demand is felt at the end of the first period, at k = 1.

The demand disturbance could be of the following Types:

  1. a sudden shock demand increase, represented by a Dirac delta input function.

  2. a sudden and sustained increase in demand, represented by a step input.

  3. a demand with an increasing trend, represented by a ramp input function.

  4. a demand with second order (quadratic), third order (cubic), or higher-order trends represented by higher-order polynomial input functions.

  5. a demand with seasonality, represented by a sinusoidal input function.

  6. A random component which is represented by a White Noise process.

The first five components pertain to and define the mean demand disturbance, while the last represents the residual random variations over and above the mean demand disturbance.

The demand can have any combination of, or even all of the above components in the mean demand, and, of course, the additional stochastic component represented by a White Noise process. Such demands are often seen in supply chain warehouses.

Thus, a demand disturbance at the downstream end (at the warehouse) provides the perturbation to the system. The system is controlled through regulation/control of the replenishment flows which are the control variables in the system.

Now from the above, we can see that the system response will have two components, the Mean Response, and a stochastic component of the response. It is the mean response that characterizes the system behavior, while the stochastic component is characterized by the inventory variance.


3. Dynamic supply chain performance metrics and control action triggers

3.1 Dynamic performance metrics

The performance metrics of a supply chain that would be of interest to us from a dynamic performance point of view and control system design are explained below.

Since the demand input to the system suddenly increases, we can expect the system response to lag the demand and to fluctuate, as the system scrambles to catch up with the increased demand. Accordingly, the key performance indicators that would be of importance and interest to us from the point of view of both design and operation are the following [1, 2, 3, 4, 5, 6, 7, 8, 9]:

  1. Permanent depletion of the inventory level (the offset), if any.

  2. The trough value or the lowest dip in the inventory levels (the undershoot).

  3. The amplitude of fluctuations of the inventory.

  4. The center-line about which fluctuations occur.

  5. The fraction of time the inventory level stays depleted (in the negative region).

  6. The limiting inventory variance.

The first indicator shows whether the system is able to catch up with the demand and whether it is ultimately restored to its original level, or not. The second indicator, the trough value, represents the lowest point or value that the inventory level is likely to touch, and impacts the base stock levels that would have to be carried to maintain uninterrupted material flows in the system, (adequate pipeline inventories). The third defines the magnitude of fluctuations of the inventory levels, and high values would naturally be undesirable; and we would like to damp them down to zero as quickly as possible. The fourth indicates the center-line about which fluctuations occur and is indicative of the average inventory levels, which we would like to keep at positive levels for comfortable operation. While the fifth can be taken to be a surrogate measure of the stock-out risk, and the higher the fraction of time in the negative region (with a depleted inventory level), the higher could be taken to be the implied stock-out risk. The last, the inventory variance is the variation that we could expect even after the system has been restored to its original operating levels and is commonly taken as a measure of the robustness of the system to random demand variations.

3.2 Control action triggers

The two most common triggers for initiation of replenishment control action in the system are [1, 2, 3, 4, 5, 6, 7, 8, 9]:

  1. The inventory levels at the W/H, as found in logistic systems and warehouses,

  2. The demand variations at the W/H, as in electronic-data-interchange systems.

Thus, the control flows into the W/H are set to a function of the latest available inventory deviations and latest available Demand deviations. Thus, we have,


We discuss these points in more detail for a single stage system first.

3.3 Single stage control notation

To differentiate between the standard abbreviations in conventional control theory, we adopt the following notation used to indicate the type of control used.

P, PI, PID would represent Proportional, Proportional-integral, and Proportional-integral-derivative control as usual. Additionally, MA denotes a ‘Moving Average’ type of control (explained in detail subsequently).

An ‘I’ within parenthesis would represent ‘inventory-triggered’ control, while ‘D’ within parenthesis would denote a ‘Demand-triggered’ control. Also, ‘ID’ within parenthesis would denote a control which would have both inventory-triggered and Demand-triggered components. We give examples below.

Thus P(D) (pronounced as “P of D”) denotes Demand-triggered Proportional control, PI(I) denotes Inventory-triggered PI control. Similarly, PID(ID) (pronounced as “PID of ID”) denotes PID control with inventory-triggered and Demand-triggered components.

We also could have cases of multiple inventory triggers and multiple-demand triggers. These will be denoted as follows:

XInDm control would indicate a control of type X (P, PI, PID, MA, Composite) with n-inventory trigger terms and m-demand trigger terms. Thus, the control PDI5D3 (pronounced as “PD of I5D3”) would be Proportional-derivative control with 5 inventory triggers and 3 demand triggers.

We now first look at a single stage system below which we take to be the warehouse end of a SC.


4. Single stage control

4.1 The two response components

Since the demand has both a deterministic component as well as a stochastic component, the response of the system can also be broken down into two components, viz., a deterministic part which is the mean response, and a random or stochastic part characterized by the inventory variance. It is the mean response that portrays the behavior of the system and yields the performance indicators of the system. Whereas, the stochastic part, represents the random fluctuations that have to be accounted for even after the system is brought under full control.

From the above discussion we can see the close parallel with conventional feedback control theory. Here the information about the warehouse inventory (the state variable) is fed back to the system for initiating replenishment flow control action. Additionally, we can also have controls wherein the disturbance is also directly fed back to the system for initiation of control action.

4.2 Single stage system controls and transportation lags

In close parallel with classical feedback control theory, the control flows are functions of the state variables and the demand or input perturbation to the system. The types of functions used also closely parallel classical control theory. And hence we have P, PI, and PID controls.

Additionally, we also have Moving Average (MA) controls, wherein the control flow is set to a weighted Moving Average of the latest available inventory deviations of up to ‘r’ periods back. The parameter ‘r’ is termed the ‘Order’ of the Moving Average.

We first look at Proportional Controls.

Now, for Proportional Controls of the P(I) type, the control flows into the warehouse in stage 3 of the chain would be given by:


where l3 is the transportation lag in the flows into the warehouse from the upstream unit of the system, i.e., the finished goods inventory at the manufacturing plant in our case. And K3 is the constant of proportionality between the control flow and the latest available inventory deviation based on which the order is initiated (and hence ‘Proportional’ to the error in conventional control theory).

The development herein corresponds to the ‘Regulator Problem’ with ‘set point’ of zero. And hence control of a Responsive Chain parallels the regulator problem in conventional control theory.

In conventional modeling of SCs, the lag is taken as the number of periods strictly between the period of order initiation and the period of arrival of the consignment, not including the period of arrival of the consignment. Thus, the lag is taken to be zero if the replenishment consignment arrives in the period immediately succeeding the period of order initiation. Instantaneous replenishment is not envisaged and is very rare in SC contexts. Arrival of the consignment within the same period of order initiation is also not envisaged and is very rare in such contexts. The earliest arrival of an ordered consignment is taken to be the immediately succeeding period, for which we take the lag as zero.

This convention is based on what is normally followed in the industry and practice, as well as the literature on dynamic modeling of SCs.

Hence in our further development of dynamic models of SCs we take the lag as zero if an order initiated in period k i.e., in the interval (k – 1, k] arrives in period (k + 1), i.e., in the interval (k, k + 1].

For a transportation lag of ‘l’ periods, an order placed in the interval (k – 1, k] would arrive in the interval (k + l, k + l + 1].

Hence for the warehouse under zero lag the control flows for P(I) control would be set as:


Now it is to be noted that the order is initiated in period k i.e., in the interval (k – 1, k] based on the latest fully observed inventory deviation which in this case would be that at the start of period k, which is the inventory recorded at the end of period (k-1) i.e., x3k1. We take it that since the inflows in period k would be the aggregate of flows in the interval (k – 1, k], and the order is initiated within the period k (and not at time point k), the latest available fully observed inventory level would be x3k1 and not x3k because x3k would not be available to the order initiator within period k. The value of x3k is recorded at the closure of period k, after the cessation of all activities including the action of ordering, of period k. Thus, replenishment orders need to be initiated within and before the closure of a period and not at the end of a period.

Thus, warehouse records would be updated at the closure of the day’s operations and would show the closing inventory and the replenishment orders placed during the day. The replenishment orders would have been placed based on the previous day’s closing stock and would arrive during the course of the next day if the lag is zero.

This is the convention that we will follow in the further development of the models.

We next look at the next type of control, which is the PI(I) control.

For the Proportional-Integral (PI(I) type) control case with zero lag, the control flows at the warehouse would be set as under:


where the second term is the integral term in our discrete-time system, and Kc is the proportionality constant (gain term) factor of the integral of the error.

Next we have the PID(I) control wherein the control flows into the warehouse under zero lag would be set to:


where the last term represents the derivative term in our discrete-time system, and Kd is the proportionality constant factor (gain term) of the derivative component of the control.

Next we have the Moving Average (MA) type of control (MA(I) type), for which the control flows into the warehouse under zero lag would be set to:


where the r is the order of the moving average, the Ks are the control parameters (the weights) of the MA terms. Thus, the control flow is set to a weighted moving average of the latest available fully observed inventory levels up to r period back.

The above controls discussed above are the conventional inventory-triggered schemes. In all these types of controls, we could additionally have demand-triggered terms also, like the PI(ID), PID(ID), MA(ID) etc.

We discuss some of them below for single stage systems.


5. Solution of single stage systems

Firstly, we note herein that in solving for the response of a supply chain system, we will not use the Z-transform nor block diagrams and transfer functions as in conventional control theory. Rather we will work directly on the system difference equation in the time domain itself. And instead, we will use direct Operator methods to obtain the system solution and response (rather than the transformed equations and inverse transforms). This is one of the advantages of this modeling paradigm.

Another advantage of the type of discrete time modeling taken up here is that transportation lags do not result in differential-delay equations as in conventional continuous-time control theory. Rather, transportation lags would increase the order of the resulting system difference equation, which is expected to be easier to solve than differential-delay equations.

We now take up the simplest form of control which is the P(I) control and illustrate the formulation of the system equation and its solution method.

5.1 P(I) control under zero lag

The control is an inventory-triggered Proportional control. And we take the Replenishment Lag = 0 in the simplest case.

The flow balance equation for the warehouse is given by:


The control flows into the warehouse are given by:


It is to be noted that the value of K3 < 0, since the control flows and inventory deviation have to be of opposite sign, i.e., when the inventory deviation is (−)ve (inventory is at a depleted level) then the flow deviation has to be (+)ve to enable more/extra flow into the warehouse to make up the inventory shortfall. Likewise, if the inventory deviation is (+)ve (extra inventory in stock) then the flow deviation has to be (−)ve (i.e., reduced flow) to reduce the inflow into the warehouse to maintain inventory levels.

Thus, we can clearly see that the controls are of the ‘feedback’ type, and they seek to keep the inventory deviation at zero level, which is just the’ Linear Regulator with zero set-point’ in standard control theory.

Thus, substituting for the control flow into the flow balance eqn. Above, yields the system-dynamic eqn. For the warehouse as:

Or equivalently,x3k+1x3kK3x3k1r3k+1E11

The above is the deterministic part of the system equation. Since demand is a stochastic variable with a stochastic component, the complete system equation is given by:


which has both parts. To solve the system equation completely, we split it into its two components and solve for each of the components separately. We hence solve the following two equations, one each for the deterministic part and the stochastic part.

x3detk+1x3detkK3x3detk1r3k+1for the deterministic part,and,E13
x3stock+1x3stockK3x3stock1εk+1for the stochastic part.E14

The first is a deterministic Linear Difference Equation, and the second, a Stochastic Linear Difference Equation (SDE).

Both are second-order Linear Difference Equations (LDEs) in the state variable x3k and can be solved by direct Operator methods for any type of input or disturbance term on the RHS of the system LDE.

An excellent treatment of difference calculus and solution methods for LDEs is given in [10]. We follow the methods given therein.

In order to solve the LDE, we first introduce the Forward Shift Operator E as under: Ex3kx3k+1, with the important property that: EEx3k=E2x3k=x3k+2. Before using the Operator, we first write the LDE in standard form as under: (the lowest time value is taken as k):


Now using the forward Shift Operator E, we can write the LDE in Operator form as:

E2EK3x3detkr3k+2for the deterministic part,and,E16
E2EK3x3stockεk+2for the stochastic part.E17

5.1.1 The mean response: solution of the deterministic LDE

We first look at the deterministic part of the solution below, which will yield the mean response. For notational convenience, we drop the superscript on x3k.


Now this is an LDE of order two (the order being the highest power of the Operator E).

We write the Characteristic Equation of the LDE as [10]:


to determine the characteristic roots of the LHS Operator.

From which we get the rootsas:α1/2=1±1+4K32.E19

The stability of the system is entirely controlled by the roots of the LHS Operator. And elementary analysis leads to the following stability conditions:

K30: instability

1/4K3<0: stability with non-oscillatory response

1K3<1/4: stability with oscillatory behavior

K3<1: instability with oscillatory behavior

Now we take up the solution of the system LDE for a unit step increase in demand, i.e.,


Substituting for the demand disturbance in the system equation yields the LDE:


which we can call the “Original Non-Homogeneous Eqn.” (O-NHE).

We first look at the solution of the homogeneous LDE (i.e., with RHS = 0). The homogeneous LDE is:


which upon factoring the LHS Operator can be written as: (λ1,λ2 being the distinct roots of the LHS Operator):


which has the solution as: x3kC1λ1k+C2λ2k for the case of distinct roots.

For a repeated root, the solution is given by: x3kC0+C1kλk for a repeated root of algebraic multiplicity two.

Now that we have the solution of the homogeneous LDE, we next look for a particular solution of the Original Non-Homogeneous LDE (O-NHE).

A standard method of solution of the Non-homogeneous eqn. is by the ‘Annihilator Method’ [10].

We look for the Operator that annihilates the RHS terms of the O-NHE, say A(E). Then operating by the Annihilator on both sides of the O-NHE yields:


which is a homogeneous LDE albeit of a higher order, but which can be solved by factorizing the LHS Operator. As an example, in our case of a unit step disturbance, r3k1,k1.

And the Annihilator is given by:


And hence the equivalent Homogeneous LDE is given by:


which has the solution:

x3kD1k+C1λ1k+C2λ2kD+C1λ1k+C2λ2kfor the case of distinct roots.E25
x3kD1k+C0+C1kλkfor the case of repeated rootsE26

Now we note that the O-NHE being of order two, will admit only two undetermined constants. The above solution, however, has three undetermined constants. The third was introduced by us due to the Annihilator. Hence to determine the extra constant D in the solution, we substitute the solution into the O-NHE to determine D. Thus, and noting that the terms involving the roots of the LHS operator of the O-NHE are precisely the homogeneous solution terms of the O-NHE, we only need substitute the extra terms introduced by the Annihilator into the O-NHE. Hence, we have:


Noting that E(D) = D itself (EDDk+1D,k), we obtain Dk+2Dk+1K3DDDK3D=1 from which we readily obtain

D=1/K3which isvedueto thevesignofK3.E27

Thus, the full solution is given by:

x3k1/K3+C1λ1k+C2λ2kfor the case of distinct roots,andE28
x3k1/K3+C0+C1kλkfor the case of repeated roots.E29

Now, for stable solutions, we will have λ1<1,λ2<1, and hence the terms involving the roots of the LHS Operator decay to zero for large k, leaving the ‘Offset’ term as 1/K3 which is the Offset value. The Offset value is (−)ve because of the (−)ve sign of K3.

The Damping rate, which is the rate at which the fluctuations decay to zero are given by the magnitudes of the roots of the LHS Operator λ1kλ2k, which can clearly be seen to decrease in magnitude with decrease in magnitude of K3. However, the Offset value increases in magnitude with decrease of K3.

We take the case of distinct roots first, say for a value of K3=1/2, in the stable region.

Hence, we have the system equation as:


The roots of the LHS Operator are readily obtained as: λ=1/2±1/2j,λ1<1,λ2<1,j=1, which can be written in the form.

λ=ρe±1/2e±/4, which hence yields the homogeneous solution as:

ACos/4+BSin/41/2k, and hence the full solution as:


which shows a Damping Rate of the Order of O1/2k=O0.707k, with an Offset of - 2 units. The response if plotted in Figure 2, and the ‘Undershoot’ is obtained as −2.5 units.

Figure 2.

P(I) control.

We next take a value of K3=1/4 for the case of repeated roots, again in the stable region. The roots of the LHS Operator are obtained as: λ1=λ2=λ=1/2, and the full solution as:

x3k1/K3+C0+C1k1/2kwhich forK3=1/4yieldsE33

which shows a Damping Rate of the Order of O1/2k=O0.5k, but with a high Offset value of - 4 units, and an Undershoot of - 4 units.

Thus, while we would like the oscillations to be damped out rapidly, this would compromise on the Offset value, which would impact the base stock requirements of the system.

Thus, for the practitioner, the implications are quite clear: there is a trade-off between stability and rapid damping on the one hand, and the base stock requirements of the system (for low stock-out risk) on the other. The higher the damping (stability) required, the higher would be the base-stock requirements to keep stock-out risk low; and alternatively, the higher the (stability) damping achieved, the higher would be the stock-out risk at fixed base-stock levels.

We next obtain the undetermined constants in the general solutions above, using the Initial Conditions (ICs) of the system, as under:

The standard ICs of the system are as: x3k0r3k0k0.

Now the system LDE: E2EK3x3krk+21,k0 is valid for k0.

And hence the ICs for our warehouse system can be obtained from the system equation itself using the standard system ICs and yields: x30=0x31=1. Substituting these ICs into the solutions above yields the full solutions as:


We additionally examine the case for K3=1 (the maximum possible magnitude for stability), the response is given by

x3k1+ACos/3+BSin/31k,k0, which using the LDE ICs yields:


which can be simplified to.


where H(.) is the unit Heaviside step function and yields a sinusoidal pattern with a center-line of – 1 and constant amplitude of 2/3 and the maximum negative deviation in inventory, the undershoot, equal to −2.

This last case is that of marginal stability characterized by constant amplitude perpetual oscillations which never die down to zero.

The responses for the three cases above are plotted in Figure 2.

5.1.2 The limiting inventory variance: solution of the SDE

We next look at the determination of the limiting inventory variance, which is a measure of the variation that we could expect even after the system (mean inventory level) has been restored to its original value.

Also, for the behavior of the mean response as well as our inferences from it to be meaningful, it is necessary that the limiting inventory variance be bounded and finite. We can then expect the inventory levels to be within the band given by: x3kdet±3limtvarx3k, where x3kdet is the mean response that has been determined above from the deterministic system LDE.

In order to determine the limiting inventory variance, we make use of the stochastic component of the system equation and determine the stochastic part of the response. For our system the Stochastic LDE (or SDE) is as under:


where the term on the RHS of the SDE is the random variation represented by a White Noise Process, with εkWN0σ2.

We can note that the LHS Operator is again the same as in the deterministic part of the system LDE that we have solved for above. Hence the roots of the LHS Operator remain unaltered in the SDE also.

Now following the method used in [11], we can note that if unity is not a root of the LHS Operator, then the SDE admits as a solution, an infinite weighted Moving Average representation in terms of the White Noise disturbance terms as under:


where the βl s are the weights.

And hence the stochastic part of the solution of the system eqn. can be written as an infinite linear combination of the white noise disturbance terms.

Now since the individual white noise terms are Uncorrelated and Normal, i.e., with εkN0σ2, and, with Covεiεj=δijσ2,δij=0,ij1,i=j, the limiting variance of x(k) is given by:


In order to solve for the weighting terms, the βs, we substitute the solution into the system SDE above, as under:

E2EK3x3kstocεk+2,valid forallk0,E42

which is


Another and more convenient way to write the SDE is to use the Backward Shift Operator L, defined by:


The SDE for the βs becomes:


which is


Now comparing coefficients of ε(k) for each k, yields the system of equations as below:

LDE termCoefft of ε(k)Coefft of ε(k-1)Coefft of ε(k-2)Coefft of ε(k-3)Coefft of ε(k-4)Coefft of ε(k-5)….Coefft of ε(k)
- L0β0β1β2β3β4…..βk1

The above set of equations yields: β0=1,β1β0=0,&,βkβk1K3βk2=0,k2, which is an LDE for the βs, as.


We can note that the LHS Operator is the same as for the system LDE, and hence has the same characteristics and form of solution.

We first illustrate the computation for the stable case K3=1/2, for which the solution form has already been obtained as: ACos/4+BSin/41/2k, and hence plugging in the ICs, we have the solution for βs as:


To obtain the limiting inventory variance, we firstly note that:


And hence,


Hence, we have: Limkvarx3k3.5σ2, showing that the inventory variance is bounded and finite for this case.

We next illustrate the computation for the marginally stable case K3=1, for which the solution form has already been obtained as: ACos/3+BSin/3, and hence plugging in the ICs, we have the solution for the βs as:


And we can see from the above that the βs oscillate in value, from −1 to 0 to +1 infinitely often, i.e., the sequence {…0, −1, −1, 0, +1,+1, 0…..} repeats infinitely often. And hence the series k=0βk2 diverges to infinity.

Hence in this case the limiting inventory variance is not bounded and not finite.

It can similarly be shown that the limiting inventory variance is not bounded for unstable solutions also.

Thus, the limiting inventory variance will be bounded only for stable cases.

5.1.3 The complete solution

We can also note from the above discussion and work up that we have also obtained the stochastic part of the solution, as: x3kstocl=0k1βlεkl, where the βs are as has been obtained above. Thus, the complete solution can be written as:


where x3kdet is the mean response derived from the solution of the deterministic part of the LDE, and the βs in the stochastic part are as obtained above by solution of the stochastic part of the system equation, the SDE.

The above representation proves useful for simulation purposes.

We next take up the P(ID) control. We discuss only the deterministic system LDE hereafter, since the stochastic part of the solution and the limiting inventory can be obtained by methods similar to that discussed above in all cases to follow.

5.2 P(ID) control under zero lag

In this type of control an additional demand-triggered component is also added to the control thereby making it more proactive. The control initiates corrective replenishment action no sooner than a demand deviation is observed. It does not wait for an inventory deviation to take place before initiating replenishment action though it does have an inventory-triggered component also.

The replenishment control flow is given by:


where the first term is the inventory-triggered component and the second the demand-triggered component. Substituting for the control flow into the system equation yields:


which can be written as


We can note that the addition of the demand-triggered component has left the LHS of the LDE unaltered. Thus, the LHS Operator of the LDE remains the same and is unaffected by addition of demand-triggered components to the control. The system eqn. can hence be written in Operator form as:

E2EK3x3kK03r3kr3k+2valid ink0E55

Since r3k1,k1, the system LDE can be written as:

E2EK3x3kK031valid ink1,E56
with theICsnowas:x31=1x32=2.E57

The ICs for the LDE have been obtained from the system Eq. (55) above using the standard system ICs; x3kr3k00k0 applied to the Eq. (55) above.

Since the LHS Operator is the same as for the earlier P(I) control, the stability analysis remains the same as earlier, as also the roots of the LHS Operator for various values of the inventory-trigger parameter discussed earlier, i.e.,K3=1/41/21.

The solutions are the same as given earlier in Eqs. (25) and (26).

And substituting the solution back into the O-NHE yields the value of the extra constant D as


which hence yields the solutions for the two cases as:

x3kK031/K3+C1λ1k+C2λ2kfor the case of distinct rootsE59
x3kK031/K3+C0+C1kλkfor the case of repeated roots(K3=1/4)E60

where the offset term is now given by K031/K3 for both cases.

The important point to note in the above solution is that the offset can be made zero by choice of the demand-trigger parameter as K03=1.

And hence we can observe the enhanced response of the P(ID) control over the earlier P(I) control, in that the Offset can now be controlled by us by choice of K03.

In fact, we can also achieve a (+)ve value of the offset by choosing K031.

We can now obtain the full solutions for the three cases above, using the LDE ICs x31=1x32=2. We take K03=1 to be able to obtain zero offset. Hence, we have:

ForK3=1/4the repeated roots case:x3k46k1/2k,k1E61
ForK3=1the marginal stability case:x3k2Cosk+1π/3,k1E63

The solution curves are plotted in Figure 3, from which we can see that the response in all cases has zero offset.

Figure 3.

P(ID) control.

We can similarly extend the modeling and analysis to PI(I), PID(I), and MA(ID) controls.

We could also have different types of input disturbances as indicated earlier in Section 3.1. Additionally, the third dimension of our analysis could be to have non-zero lags, i.e., lags of one, two periods, and so on. Cases with non-zero and higher lags will result in higher order system LDEs, and the LHS Operator would be of a higher order.

Further details of the above can be obtained in [1, 2, 3, 4, 5, 6, 7, 8, 9].


6. Controls for multi-stage supply chains

We look at the serial supply chain system as given in Figure 1.

We can see from Figure 1 that the immediately succeeding downstream stage in a supply chain will provide the “demand perturbation” for the immediately preceding stage. Thus, the demand perturbation at the warehouse at the downstream end will successively be felt up the chain. And the single-stage analysis described above can be used in turn for each stage of the chain.

For non-serial supply chains, the arguments are similar and single-stage analysis can be used as described above.

Details of some of these analyses can be found in [1, 2, 3, 4, 5, 6, 7, 8, 9].


7. Conclusion

This chapter has presented the application of control concepts to the control of supply chains. The state variables have been taken to be the inventory levels, while the control variables are the replenishment flows into the various stages of the system. The conventional P, PI, PID controls have been discussed, as also some newer forms of control which are especially applicable to supply chains and warehouses. The performance of P(I) and P(ID) controls have been derived in detail, and their performance analyzed.

A significant feature of this chapter is that the conventional block diagrams and transfer functions of conventional control theory have not been used. Rather direct Operator Methods have been used to good advantage to solve the system equations.


  1. 1. Nilakantan, Kannan; Enhancing Supply Chain performance with improved Order-control Policies, International Journal of System Science Sep2010, Vol 41-9, p 1099-1113, Taylor and Francis
  2. 2. Nilakantan, Kannan; Replenishment Policies for Warehouse Systems under Cyclic demand, International Journal of Business Performance and Supply Chain Modelling, 2013,Inderscience, Vol 5-2, pp 148-176
  3. 3. Nilakantan, Kannan; Supply Chains subject to Demand Shocks, International Journal of Logistics Systems and Management, 2015, Inderscience, Vol 21-2, pp133-159
  4. 4. Nilakantan, Kannan; Supply Chains under Demands with Trend Components, and Robust Controls, Annals of Management Science, International Centre for Business and Management Excellence, Franklin, Tennessee, USA.Vol 3-2, Dec 2014, pp 27-64
  5. 5. Nilakantan, Kannan; Warehouse Control Systems for General Forms of Demand, International Journal of Operations and Quantitative Management, 2014, INFOMS, USA, Vol 20–4, pp 273–299,Dec 2014
  6. 6. Nilakantan, Kannan; Replenishment systems for responsive supply chains under dynamic and sudden lead-time disturbances, International Journal of Systems Science-Operations and Logistics, 2019, Elsevier, Vol 6-4, pp 320-345
  7. 7. Nilakantan, Kannan; Supply Chain Resilience to Sudden and Simultaneous Lead-time and Demand Disturbances International Journal of Supply Chain and Operations Resilience, 2019, Inderscience (USA), Vol 3 - 4, pp. 292-329
  8. 8. Nilakantan Kannan, Responsiveness and Recovery Performance Analysis, and Replenishment System Design in Supply Chains under Large Replenishment Lead-times and Uncertain Demand, International Journal of Supply Chain Management 2021,(UK), Vol 10-3, pp 1-17
  9. 9. Nilakantan Kannan The Ripple Effect in a Supply Chain: A Sudden Demand Increase and with a Sinusoidal Component, International Journal of Operations and Quantitative Management March 2021, Vol 27-1, pp. 1-37
  10. 10. Kelley, W. G., & Peterson, A. C., Difference Equations, 2/e, Academic Press, Elsevier, USA
  11. 11. Gourieroux, G, & Monfort, A., Time Series and Dynamic Models, Cambridge University Press,(1990), Cambridge, UK

Written By

Kannan Nilakantan

Submitted: 27 August 2021 Reviewed: 19 September 2021 Published: 20 October 2021