• Sonuç bulunamadı

A sharp bound for the ergodic distribution of an inventory control model under the assumption that demands and inter-arrival times are dependent

N/A
N/A
Protected

Academic year: 2021

Share "A sharp bound for the ergodic distribution of an inventory control model under the assumption that demands and inter-arrival times are dependent"

Copied!
10
0
0

Yükleniyor.... (view fulltext now)

Tam metin

(1)

R E S E A R C H

Open Access

A sharp bound for the ergodic distribution of

an inventory control model under the

assumption that demands and inter-arrival

times are dependent

Tahir Hanalio ˘glu (Khaniyev)

1,2

and Cihan Aksop

3*

*Correspondence:

caksop@gazi.edu.tr

3Science and Society Department,

The Scientific and Technological Research Council of Turkey, Ankara, 06420, Turkey

Full list of author information is available at the end of the article

Abstract

In this study, a stochastic process which represents a single-item inventory control model with (s, S)-type policy is constructed when the demands of a costumer are dependent on the inter-arrival times between consecutive arrivals. Under the assumption that the demands can be expressed as a monotone convex function of the inter-arrival times, it is proved that this process is ergodic, and closed form of the ergodic distribution is given. Moreover, a sharp lower bound for this distribution is obtained.

Keywords: dependence; ergodic distribution; inventory model of type (s, S)

1 Introduction

Consider a single-item inventory control model as follows. Customers arrive at the depot at random times{Tn}, and the amount of their demands can be modeled by a sequence of

random variables{ηn}. If there exists enough supply in the stock, then the demanded items

of the customer are satisfied from the stock, else an immediate replenishment order takes place so that to raise the inventory level to an order-up-to level S > . In other words, if

X(t) denotes stock level just before an arrival of a customer at time t and η is the amount of his/her demand, then

X(t) = ⎧ ⎨ ⎩ X(t) – η, X(t) – η > s, S, X(t) – η≤ s.

Here s≥  is a pre-defined control level. We will assume that no product is returned or defective and the supplier of this depot is reliable so that the replenishment is not delayed. This model is known as (S, s)-type policy inventory control model and it has been exten-sively studied under various assumptions in the literature (see Scarf [], Rabta and Aissani [], Chen and Yang [], Khaniyev and Atalay [], Khaniyev and Aksop []).

The classical inventory model assumes that the inter-arrival times of the customers and the amounts of the demands are mutually independent random variables. However, real life problems are generally too complex so the assumptions of classical inventory control ©2014 Hanalio ˘glu (Khaniyev) and Aksop; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

(2)

theory are not valid. This situation shows itself mainly on the structure of the amounts of the demands and the inter-arrival times between consecutive costumers.

In most cases, these demands cannot be modeled by independent random variables. For example, demands can depend on the day of the week so that at the weekends there can be a high demand while at weekdays there can be low demand, or there can be seasonal demand for that item (see Sethi and Cheng []). Also the amount of the demand can depend on the inventory level: if the supplier can offer a wide selection of his/her items, then he/she can increase the probability of making a sale (see Urban []).

The model investigated in this study assumes that the demands are dependent on the inter-arrival times. This assumption makes sense especially in the situations when the de-mands can be met only by one supplier and reaching that supplier is not easy for the cus-tomers. In this case, if a customer did not make a demand for a long time, then he or she will probably have a need for more items. This is the situation occurring at rural districts or regions which are difficult to access.

The main purpose of this paper is to investigate the ergodicity of an inventory con-trol model with (s, S)-policy under the assumption that the amount of demands is depen-dent on inter-arrival times. Next section gives a mathematical construction of the studied stochastic process X(t). Section  gives some notations used in this study, and Section  gives the main results. In the last section, some discussion is given.

2 Mathematical construction of the process X(t)

Let{(ξi, ηi)} be an independent and identically distributed random pair, where ξiand ηi

are dependent random variables with joint distribution G(x, y) and marginal distributions

and F, respectively, that is,

G(t, x) = P{ξi≤ t, ηi≤ x}, i = , , . . .

and

(t) = P{ξi≤ t}, F(x) = P{ηi≤ x}.

Moreover, let ηibe absolutely continuous random variables.

Let us construct a sequence of integer-valued random variables{Nn} as follows:

N= , N= min{n ≥  : S – Yn< s}, Nm= min  n≥ Nm–+  : S – (Yn– YNm–) < s  , m= , , . . . , where Yn= n

i=ηi, n = , , . . . . For the sake of simplicity, we will use the following

nota-tion:

ξkn= ξNk–+n, ηkn= ηNk–+n, k= , . . . , n = , , . . . .

Let us construct a stochastic process X(t) as follows:

X(t) = S –

ν(t)



i=

(3)

Here

ν(t) = max{n ≥  : Tn≤ t}, and Tn=

n

i=ξi, n = , , . . . , T= .

Whenever the value of this process falls below to a pre-defined fixed control level s > , we kill this process and a new replica X(t) of the process X(t) is constructed with the initial value S. Let us denote this time by τ; that is,

τ= inf 

t>  : X(t) < s 

.

So the new process X(t) can be expressed as follows:

X(t) = S –

ν(t)



i=

ηi, t> ,

where, ηiand ν(t) are defined similar to ηiand ν(t), respectively.

In a similar way, let us construct a sequence of stochastic processes{Xn(t), n = , , . . .}.

By using these sequences of stochastic processes, we can define the desired stochastic process X(t) as follows: X(t) = ∞  n= Xn(t – τn–)I[τn–,τn)(t), t> ,

where IA(·) is the indicator function of set A, i.e.,

IA(t) = ⎧ ⎨ ⎩ , t∈ A, , otherwise and τ= . 3 Notation

In this section, the notations used in this article are given.

p(t, x) dx = PX(t)∈ dx  , Uθ(t) = ∞  n= R∗n (t), Rn(t) = P{τn≤ t}, n = , , . . . , pn(t, x) = p(·, x) ∗ Rn (t)tp(t – u, x)Rn(du), Yn:m= m  i=n ηi, m≥ n; Yn:m= , m< n, Fn(x) = F∗n(x), U(x)≡ ∞  n= Fn(x), F(x) = , x≥ , f(x) dx = F(dx), n(t) = ∗n(t), F(x) =  – F(x), γ = S – s.

(4)

4 Main results

In Theorem . below, the ergodicity of the constructed stochastic process is obtained un-der the assumption that the demands can be expressed as a monotone increasing function of the inter-arrival times by utilizing a theorem from Gihman and Skorohod []. Then with an additional assumption of convexity, we obtained an upper bound for the first period’s distribution function (that is, for the distribution function of X(t), ≤ t < τ). An explicit expression for the ergodic characteristics of the process X(t) is given in Corollary .. In Theorem . a lower bound for the ergodic distribution is obtained.

Proposition .(see Feller []) For all t≥ , the following equation holds true: ∞  n= n Fn–(t) – Fn(t) = U(t).

Definition .(Lehmann []) A pair of random variables (X, Y ) is positively quadrant dependent if

P{X ≤ x, Y ≤ y} ≥ P{X ≤ x}P{Y ≤ y}, x, y ∈ R. ()

Proposition .(Lehmann []) If X and Y are positively quadrant dependent, then the

following inequality holds:

P{X ≥ x, Y ≤ y} ≤ P{X ≥ x}P{Y ≤ y}.

Theorem . Let E[ξ] <∞ and E[η] > . If ξand ηare positively quadrant dependent,

then E[τ] <∞.

Proof Note that

E[τ] = E N  i= ξi  = ∞  n= nE[ξ|N= n]P{N= n} = ∞  n= n ∞  P{ξ≥ x, N= n} dx. ()

On the other hand, for n = ,  we have

P{ξ≥ x, N= n} ≤ P{ξ≥ x} ()

and for n≥  we have

P{ξ≥ x, N= n} = P{ξ≥ x, Yn–≤ γ < Yn}

= γ

(5)

= γP{Yn–≤ γ – v < Yn–}P{ξ≥ x, η∈ dv} = γγ–vF(γ – v – w)Fn–(dw)P{ξ≥ x, η∈ dv} = γFn–(γ – v) – Fn–(γ – v) P{ξ≥ x, η∈ dv}. () By applying Proposition ., we get from ()

P{ξ≥ x, N= n} ≤ γFn–(γ – v) – Fn–(γ – v) P{η∈ dv}P{ξ≥ x} = Fn–(γ ) – Fn(γ ) P{ξ≥ x}. ()

Substituting () into () and using Proposition . yields

E[τ]≤   n= ∞  nP{ξ≥ x} dx + ∞  ∞  n= n Fn–(γ ) – Fn(γ ) P{ξ≥ x} dx = E[ξ] U(γ ) + F(γ ) – F(γ ) +  .

Since for all finite values of γ , U(γ ) is finite (Feller []), we get

E[τ]≤ E[ξ]

U(γ ) + F(γ ) – F(γ ) + 

 + U(γ )E[ξ] <∞. 

Theorem . Let ηn= h(ξn), n = , , . . . , where h∈ C(R+) is a monotone increasing

func-tion. Let h()≥  and E[η] = μ <∞.

(A) If supx∈Rh(x) > γ, then the process X(t) is ergodic. (B) If supx∈Rh(x) < γ, then, additionally, let

∞    –h(x) γ  dx<∞.

Then the process X(t) is ergodic.

Proof It is known that the following conditions are sufficient to prove that the process X(t) is ergodic (Gihman and Skorohod []):

. For a sequence of random variables{γn} such that  ≤ γ< γ<· · · , the process

Xn≡ X(γn)must form an ergodic Markov chain.

. E[γn+– γn] <∞, n = , , . . . .

Observe that Xn≡ X(τn) form an ergodic Markov chain because X(τn) = S for each n≥ .

To see that E[τn+–τn] <∞, it is enough to show only that E[τ] <∞ because τ, τ–τ, . . . have identical distribution. Note that

E[τ] = E N  i= ξi  = E N  i= E ξi|N  = E NE[ξ|N] .

(6)

On the other hand, from Proposition . and equation () we have E[τ] = ∞  ∞ h(x) U(γ – y)F(dy) dx.

Here, if a > b, then we takeabf dt= .

(A) Let us assume that there exists x∈ R+such that h(x) > γ holds. Then let us denote the infimum of such numbers with x∗(<∞); that is,

x∗= infx>  : h(x) > γ. Then we get ∞  γ h(x) U(γ – y)F(dy) dx = x∗  γ h(x) U(γ – y)F(dy) dxx∗  γU(γ – y)F(dy) dx = xU(γ ) – <∞. (B) If h(x) < γ for all x∈ R+, then

∞  γ h(x) U(γ – y)F(dy) dx≤ sup y∈[,γ ]  U(γ – y)f (y) γγ– h(x)dx ≤ sup y∈[,γ ]  U(γ – y)f (y) ∞  γ – h(x)dx = γ sup y∈[,γ ]  U(γ – y)f (y) ∞    –h(x) γ  dx<∞. Therefore E[τ] <∞.

Now, put γn= τnin . and . to see that the process X(t) is ergodic. 

Lemma . For every measurable function g, the following equation holds true:

E g X(t) = ∞  n= g∗ pn(t,·)(S). ()

Proof For every tn< t, we have

E g Xn(t – τn) χ Xn(t – τn) |τn= tn = E g Xn(t – tn) χ Xn(t – tn) = E g X(t – tn) χ X(t – tn) = ∞  g(x)PX(t – tn)∈ dx  = Sg(x)p(t – tn, S – x) dx. Here χ(x) = ⎧ ⎨ ⎩ , x≥ , , x< .

(7)

Therefore, E g Xn(t – τn) I[τn,τn+)(t) = tSg(x)p(t – tn, S – x) dxP{τn∈ dtn} and E g X(t) = ∞  n= tSg(x)p(t – tn, S – x) dxP{τn∈ dtn} = ∞  n= Sg(x)pn(t, S – x) dx.  In Theorem ., it is proved that under some assumptions the process X(t) is ergodic. Therefore, limt→∞P{X(t) ≤ x} exists for every x ∈ (s, S). Let us denote a random variable

Xwhich admits this limit as a distribution; that is,

P{X ≤ x} = lim

t→∞P



X(t)≤ x, x∈ (s, S).

Corollary . For every measurable function g, the following equation holds:

E g(X) =  E[τ] g∗ p(u, ·)(S) du.

Proof Note that from Theorem . we have

E g X(t) = ∞  n= Sg(x)pn(t, S – x) dx = ∞  n= Sg(x) tp(t – u, S – x)P{τn∈ du} dx = Sg(x) tp(t – u, S – x)Uθ(du) dx = Sg(x) p(·, S – x) ∗ Uθ (t) dx.

On the other hand, it is well known from the key renewal theorem that

lim t→∞ p(·, S – x) ∗ Uθ (t) =E[τ] ∞  p(t, S – x) dt. Therefore we get E g(X) = Sg(x)E[τ] ∞  p(t, S – x) dt dx =  E[τ] Sg(x)p(t, S – x) dx dt =  E[τ] ∞  g∗ p(t, ·)(S) dt. 

(8)

The following theorem can be obtained by an application of Theorem ., Lemma . and Corollary ..

Theorem . For all x≥ , P{X ≤ x} =E[τ] ∞  xp(t, S – x) dx dt.

Remark Since X(t) is ergodic and limt→∞P{X(t) ≤ x} = P{X ≤ x}, distribution in

Theo-rem . is the ergodic distribution of the process X(t).

Lemma . In addition to the assumption in Theorem., let h(x) be a convex function.

Then the following inequality holds for t< τ:

PX(t)≤ x  ≤ ∞  n= t nh(S–xn ) (t – y)n(dy). ()

Here, if a > b, we takeab· dt = . The inequality in () will be ≥, when h is a concave

function.

Proof Note that

PX(t)≤ x  = P  Sν(t)  i= ηi≤ x  = ∞  n= P ν(t)  i= ηi≥ S – x, ν(t) = n  = ∞  n= P  n  i= ηi≥ S – x, n  i= ξi≤ t < n+  i= ξi  = ∞  n= P  n  i= h–(ξi)≥ S – x, n  i= ξi≤ t < n+  i= ξi  = ∞  n= P   n n  i= h–(ξi)≥ S– x n , n  i= ξi≤ t < n+  i= ξi  ≤ ∞  n= P  h–   n n  i= ξi  ≥S– x n , n  i= ξi≤ t < n+  i= ξi  = ∞  n= P   n n  i= ξi≥ h  S– x n  , n  i= ξi≤ t < n+  i= ξi  = ∞  n= t nh(S–xn ) (t – y)n(dy). () 

Remark Note that the series in () is convergent.

Theorem . Under the assumptions of Lemma., we have

P{X ≤ x} ≥E[τ] ∞    – ∞  n= t nh(xn) (t – y)n(dy)  dt, x≥ . ()

(9)

Proof From Theorem . we have P{X ≤ x} =E[τ] xp(t, S – x) dx dt,

which is also the ergodic distribution of X(t). Therefore, with an application of Lemma . we get P{X ≤ x} =E[τ] ∞  S S–x p(t, y) dy dt =  E[τ] ∞  PX(t)≥ S – x  dt =  E[τ] ∞   – PX(t)≤ S – x  dt ≥  E[τ] ∞    – ∞  n= t nh(xn) (t – y)n(dy)  dt. () 

Proposition . The inequality in Theorem. is sharp; that is, there exists a convex

function such that() is satisfied with equality.

Proof To prove this proposition, it is enough to show that () is sharp. Let h(x) = ax, a > . From () we have PX(t)≤ x  = ∞  n= P  n  i= ξi≥ a(S – x), n  i= ξi≤ t < n+  i= ξi  = ∞  n= t a(S–x) (t – y)n(dy). ()

Therefore, () is satisfied with equality. Moreover, note that in this case the assumptions of Theorem . are satisfied. Therefore the process X(t) is ergodic and limt→∞P{X(t) ≤ x}

is finite. The proof of this proposition follows by substituting () in ().  5 Conclusion

In this study, a stochastic model is constructed for an inventory where the customers’ de-mands are dependent on their arrival times. This assumption is important for modeling real life problems such as the models for the supply chain of items to researchers at poles or space. In Theorem ., under some assumptions, it is proved that the stochastic process

X(t) is ergodic. Moreover, an explicit expression for the ergodic characteristics of the pro-cess X(t) is obtained. A bound for the ergodic distribution is given in Theorem . which is sharp.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TH formulated the considered problem and gave the construction of the process X(t). Additionally, TH participated in the process of the proofs of Theorem 4.1, Theorem 4.2, Theorem 4.3 and Theorem 4.4. CA carried out the proofs of Lemma 4.1, Corollary 4.1, Lemma 4.2 and Proposition 4.3. Additionally, CA participated in the process of the proofs of Theorem 4.1, Theorem 4.2, Theorem 4.3 and Theorem 4.4. Both authors read and approved the final manuscript.

(10)

Author details

1Department of Industrial Engineering, TOBB University of Economics and Technology, Sö ˘gütözü, Ankara, 06560, Turkey. 2Institute of Cybernetics, Azerbaijan National Academy of Sciences, Baku, Azerbaijan.3Science and Society Department,

The Scientific and Technological Research Council of Turkey, Ankara, 06420, Turkey.

Acknowledgements

This study was partially supported by TÜB˙ITAK 110T559 coded project.

Received: 15 October 2013 Accepted: 17 January 2014 Published:13 Feb 2014

References

1. Scarf, H: The optimality of (s, S) policies in the dynamic inventory problem. In: Sheshinski, E, Weiss, Y (eds.) Mathematical Methods in the Social Sciences 1959. Proceedings of the First Stanford Symposium, pp. 196-202. Stanford University Press, Stanford (1959)

2. Rabta, B, Aissani, D: Strong stability in an (R, s, S) inventory model. Int. J. Prod. Econ. 97, 159-171 (2005)

3. Chen, Z, Yang, Y: Optimality of (s, S, p) policy in a general inventory-pricing model with uniform demands. Oper. Res. Lett. 38, 256-260 (2010)

4. Khaniyev, T, Atalay, KD: On the weak convergence of the ergodic distribution for an inventory model of type (s, S). Hacet. J. Math. Stat. 39(4), 599-611 (2010)

5. Khaniyev, T, Aksop, C: Asymptotic results for an inventory model of type (s, S) with a generalized beta interference of chance. TWMS J. Appl. Eng. Math. 1(2), 223-236 (2011)

6. Sethi, SP, Cheng, F: Optimality of (s, S) polices in inventory models with Markovian demand. Oper. Res. 45(6), 931-939 (1997)

7. Urban, TL: Inventory models with inventory-level-dependent demand: a comprehensive review and unifying theory. Eur. J. Oper. Res. 162, 792-804 (2005)

8. Gihman, II, Skorohod, AV: Theory of Stochastic Processes, vol. II, 1st edn. Springer, Berlin (1975)

9. Feller, W: An Introduction to Probability Theory and Its Applications, vol. II, 1st edn. Wiley, New York (1971) 10. Lehmann, EL: Some concepts of dependence. Ann. Math. Stat. 37, 1137-1153 (1966)

10.1186/1029-242X-2014-75

Cite this article as: Hanalio ˘glu (Khaniyev) and Aksop: A sharp bound for the ergodic distribution of an inventory control model under the assumption that demands and inter-arrival times are dependent. Journal of Inequalities and

Referanslar

Benzer Belgeler

Birinci Dünya Savaşı’na giriş karann- da, zamanlama ve tertiplenmelerde, 2 Ağustos 1914 gizli Alman ittifak anlaş­ masında, aynı gün genel seferberlik ila­

The turning range of the indicator to be selected must include the vertical region of the titration curve, not the horizontal region.. Thus, the color change

In this example, total sales remained at the original figure, but prime costs expressed as a percentage of sales increased because of lowered menu prices and increased labor cost

In order to determine the number of moles of reactant B remaining at time t, for every mole of A that consumed, b/a moles of B must consumed; therefore, the number of moles of B

Extensive property is the one that is dependent on the mass of the system such as volume, kinetic energy and potential energy.. Specific properties are

In the new public management, accountability contains all the legal, political and financial dimensions, unlike traditional public administration, it takes on managerial

Shirley Jackson’s famous story “The Lottery” takes place in an American town and like many of her works it includes elements of horror and mystery.. Name symbolism,

Clearly, second-language learners who have no extensive access to native speakers are likely to make slower progress, particularly in the oral/aural aspects of language