arxiv: v1 [math.pr] 1 Aug 2014

Transcription

1 A Pollng Model wth Smart Customers M.A.A. Boon A.C.C. van Wk O.J. Boxma I.J.B.F. Adan arxv: v [math.pr] Aug 204 September, 200 Abstract In ths paper we consder a sngle-server, cyclc pollng system wth swtch-over tmes. A dstngushng feature of the model s that the rates of the Posson arrval processes at the varous queues depend on the server locaton. For ths model we study the ont queue length dstrbuton at pollng epochs and at server s departure epochs. We also study the margnal queue length dstrbuton at arrval epochs, as well as at arbtrary epochs whch s not the same n general, snce we cannot use the PASTA property. A generalsed verson of the dstrbutonal form of Lttle s law s appled to the ont queue length dstrbuton at customer s departure epochs n order to fnd the watng tme dstrbuton for each customer type. We also provde an alternatve, more effcent way to determne the mean queue lengths and mean watng tmes, usng Mean Value Analyss. Furthermore, we show that under certan condtons a Pseudo-Conservaton Law for the total amount of work n the system holds. Fnally, typcal features of the model under consderaton are demonstrated n several numercal examples. Keywords: Pollng, smart customers, varyng arrval rates, queue lengths, watng tmes, pseudo-conservaton law Introducton The classcal pollng system s a queueng system consstng of multple queues, vsted by a sngle server. Typcally, queues are served n cyclc order, and swtchng from one queue to the next queue requres a swtch-over tme, but these assumptons are not essental to the analyss. The decson at what moment the server should start swtchng to the next queue s mportant to the analyss, though. Pollng systems satsfyng a so-called branchng property The research was done n the framework of the BSIK/BRICKS proect, and of the European etwork of Excellence Euro-F. Eurandom and Department of Mathematcs and Computer Scence, Endhoven Unversty of Technology, P.O. Box 53, 5600MB Endhoven, The etherlands Eurandom, Department of Industral Engneerng & Innovaton Scences and Department of Mathematcs and Computer Scence, Endhoven Unversty of Technology, P.O. Box 53, 5600MB Endhoven, The etherlands

2 generally allow for an exact analyss, whereas pollng systems that do not satsfy ths property rarely can be analysed n an exact way. See Resng [24], or Fuhrmann [4], for more detals on ths branchng property. There s a huge lterature on pollng systems, manly because of ther practcal relevance. Applcatons are found, among others, n producton envronments, transportaton, and data communcaton. The surveys of Takag [27], Levy and Sd [2], and Vshnevsk and Semenova [29] provde a good overvew of applcatons of pollng systems. These surveys, and [30], Chapters 2.2 and 3, are also excellent references to fnd more nformaton about varous analyss technques, such as the Buffer Occupancy method, the Descendant Set approach, and Mean Value Analyss MVA for pollng systems. The vast maorty of papers on pollng models assumes that the arrval rate stays constant throughout a cycle, although t may vary per queue. The pollng model consdered n the present paper, allows the arrval rate n each queue to vary dependng on the server locaton. Ths model was frst consdered by Boxma [5], who refers to ths model as a pollng model wth smart customers, because one way to look at ths system s to regard t as a queueng system where customers choose whch queue to on, based on the current server poston. ote that Boxma s defnton of smart customers s dfferent from the defnton used by Mandelbaum and Yechal [22], who study an M/G/ queue where smart customers may decde upon arrval to on the queue, not to enter the system at all, or to wat for a whle and postpone the decson. A relevant applcaton can be found n [6], where a pollng model s used to model a dynamc order pckng system DPS. In a DPS, a worker pcks orders arrvng n real tme durng the pckng operatons and the pckng nformaton can dynamcally change n a pckng cycle. One of the challengng questons that onlne retalers now face, s how to organse the logstc fulfllment processes durng and after order recept. In tradtonal stores, purchased products can be taken home mmedately. However, n the case of onlne retalers, the customer must wat for the shpment to arrve. In order to reduce throughput tmes, an effcent enhancement to an ordnary DPS s to have products stored at multple locatons. The system can be modelled as a pollng system wth queues correspondng to each of the locatons, and customers correspondng to orders. The locaton of the worker determnes n whch of the queues an order s beng placed. In ths system arrval rates of the orders depend on the locaton of the server.e. the worker, whch makes t a typcal smart customers example. A graphcal llustraton s gven n Fgure. We focus on one specfc order type, whch s placed n two locatons, say Q and Q. Whle the pcker s on ts way to Q, say at locaton, all of these orders are routed to Q and the arrval rate at Q s zero. If the pcker s between Q and Q, say at locaton 2, the stuaton s reversed and Q receves all of these orders. Besdes practcal relevance, the smart customers model also provdes a powerful framework to analyse more complcated pollng models. For example, a pollng model where the servce dscplne swtches each cycle between gated and exhaustve, can be analysed constructng an alternatve pollng model wth twce the number of queues and arrval rates beng zero durng specfc vst perods [8]. The dea of temporarly settng an arrval rate to zero s also used n [2] for the analyss of a pollng model wth multple prorty levels. Tme varyng arrval rates also play a role n the analyss of a pollng model wth renegng at pollng nstants []. Concernng state dependent arrval rates, more lterature s avalable for systems consstng of only one queue, often assumng phase-type dstrbutons for vacatons and/or servce tmes. A system consstng of a sngle queue wth server breakdowns and arrval rates dependng on 2

3 Q Pcker Q Q 2 2 Q Q Depot Fgure : A dynamc order pckng system. Orders are placed n queues Q,..., Q. the server status s studed n [26]. A dfference wth the system studed n the present paper, besdes the number of queues, s that the machne can break down at arbtrary moments durng the servce of customers. Pollng systems wth breakdowns have been studed as well, cf. [9, 7, 20, 23]. However, only akdmon and Yechal [23] consder a model where the arrval process stops temporarly durng a breakdown. Shanthkumar [25] dscusses a stochastc decomposton for the queue length n an M/G/ queue wth server vacatons under less restrctve assumptons than Fuhrmann and Cooper [5]. One of the relaxatons s that the arrval rate of customers may be dfferent durng vst perods and vacatons. Another system, wth so-called workng vacatons and server breakdowns s studed n [8]. Durng these workng vacatons, both the servce and arrval rates are dfferent. Mean watng tmes are found usng a matrx analytcal approach. For pollng systems, a model wth arrval rates that vary dependng on the locaton of the server has not been studed n detal yet. Boxma [5] studes the ont queue length dstrbuton at the begnnng of a cycle, but no watng tmes or margnal queue lengths are dscussed. In a recent paper [], a pollng system wth Lévy-drven, possbly correlated nput s consdered. Just as n the present paper, the arrval process may depend on the locaton of the server. In [] typcal performance measures for Lévy processes are determned, such as the steady-state dstrbuton of the ont amount of flud at an arbtrary epoch, and at pollng and swtchng nstants. The present paper studes a smlar settng, but assumes Posson arrvals of ndvdual customers. Ths enables us to fnd the probablty generatng functons PGFs of the ont queue length dstrbutons at pollng nstants and customer s departure epochs, and the margnal queue length dstrbutons at customer s arrval epochs and at arbtrary epochs whch are not the same, because PASTA cannot be used. The ntroducton of customer subtypes, categorsed by ther moment of arrval, makes t possble to generalse the dstrbutonal form of Lttle s law see, e.g., [9], and apply t to the ont queue length dstrbuton at departure epochs to fnd the Laplace- Steltes Transform LST of the watng tme dstrbuton. The present paper s structured as follows: Secton 2 gves a detaled model descrpton and ntroduces the notaton used n ths paper. In Secton 3 the PGFs of the ont queue length dstrbutons of all customer types at pollng nstants are derved. The margnal queue length dstrbuton s also studed n ths secton, but we show n Secton 4 that the dervaton of 3

4 the watng tme LST for each customer type requres a more complcated analyss, based on customer subtypes. In Sectons 3 and 4 we need nformaton on the lengths of the cycle tme and all vst tmes, whch are studed n Secton 5. In Secton 6 we adapt the MVA framework for pollng systems, ntroduced n [3], to our model. Ths results n a very effcent method to compute the mean watng tme of each customer type. For pollng systems wth constant arrval rates, a Pseudo-Conservaton Law PCL s studed n [6]. In Secton 7 we show that, under certan condtons, a PCL s satsfed by our model. Fnally, we gve numercal examples that llustrate some typcal features and advantages of the model under consderaton. 2 Model descrpton and notaton The pollng model n the present paper contans queues, Q,..., Q, vsted n cyclc order by one server. Swtchng from Q to Q + =,...,, where Q + s understood to be Q, etc. requres a swtch-over tme S, wth LST σ. We assume that at least one swtch-over tme s strctly greater than zero, otherwse the mean cycle length n steady-state becomes zero and the analyss changes slghtly. See, e.g., [4] for a relaton between pollng systems wth and wthout swtch-over tmes. Swtch-over tmes are assumed to be ndependent. The cycle tme C s the tme that elapses between two successve vst begnnngs to Q, and C s the tme that elapses between two successve vst endngs to Q. The mean cycle tme does not depend on the startng pont of the cycle, so E[C ] = E[C ] = E[C]. The vst tme V of Q s the tme between the vst begnnng and vst endng of Q. The ntervst tme I of Q s the tme between a vst endng to Q and the next vst begnnng at Q. We have C = V + I, and I = S + V S +, =,...,. Customers arrvng at Q,.e. type customers, have a servce requrement B, wth LST β. We also assume ndependence of servce tmes, and frst-come-frst-served FCFS servce order. The servce dscplne of each queue determnes the moment at whch the server swtches to the next queue. In the present paper we study the two most popular servce dscplnes n pollng models, exhaustve servce the server swtches to the next queue drectly after the last customer n the current queue has been served and gated servce only vstors present at the server s arrval at the queue are served. The reason why these two servce dscplnes have become the most popular n pollng lterature, les n the fact that they are from a practcal pont of vew the most relevant servce dscplnes that allow an exact analyss. In ths respect the followng property, defned by Resng [24] and also Fuhrmann [4], s very mportant. Property 2. Branchng Property If the server arrves at Q to fnd k customers there, then durng the course of the server s vst, each of these k customers wll effectvely be replaced n an..d. manner by a random populaton havng probablty generatng functon h z,..., z, whch can be any -dmensonal probablty generatng functon. In most cases, a pollng model can only be analysed exactly, f the servce dscplne at each queue satsfes Property 2., or some slghtly weaker varant of ths property, because n ths case the ont queue length process at vst begnnngs to a fxed queue consttutes a Mult- Type Branchng Process, whch s a ncely structured and well-understood process. Gated and exhaustve servce both satsfy ths property, whereas a servce dscplne lke k-lmted servce serve at most k customers durng each vst does not. 4

5 The feature that dstngushes the model under consderaton from commonly studed pollng models, s the arrval process. Ths arrval process s a standard Posson process, but the rate depends on the locaton of the server. The arrval rate at Q s denoted by λ P, where P denotes the poston of the server, whch s ether servng a queue, or swtchng from one queue to the next: P {V, S,..., V, S }. One of the consequences s that the PASTA property does not hold for an arbtrary arrval, but as we show n Secton 3, a condtonal verson of PASTA does hold. Another dffculty that arses, s that the dstrbutonal form of Lttle s law cannot be appled to the PGF of the margnal queue length dstrbuton to obtan the LST of the watng tme dstrbuton anymore. We explan ths n Secton 4, where we also derve a generalsaton of the dstrbutonal form of Lttle s law. 3 Queue length dstrbutons 3. Jont queue length dstrbuton at vst begnnngs/endngs The two man performance measures of nterest, are the steady-state queue length dstrbuton and the watng tme dstrbuton of each customer type. In ths secton we focus on queue lengths rather than watng tmes, because the latter requres a more complex approach that s dscussed n the next secton. We restrct ourselves to branchng-type servce dscplnes,.e., servce dscplnes satsfyng Property 2.. Boxma [5] follows the approach by Resng [24], defnng offsprng and mmgraton PGFs to determne the ont queue length dstrbuton at the begnnng of a cycle. We take a slghtly dfferent approach that gves the same result, but has the advantage that t gves expressons for the ont queue length PGF at all vst begnnngs and endngs as well. Denote by LB P z,..., z the PGF of the steady-state ont queue length dstrbuton at begnnngs of perod P {V, S,..., V, S }. The relaton between these PGFs, also referred to as laws of moton n the pollng lterature, s obtaned by applcaton of Property 2. to LB V z, where z s a shorthand notaton for the vector z,..., z. Ths property states that each type customer present at the vst begnnng to Q wll be replaced durng ths vst by a random populaton havng PGF h z, whch depends on the servce dscplne. The only dfference between conventonal pollng models and the model under consderaton n the present paper, s that the arrval rates depend on the server locaton. The relatons between LB V z, LB S z, and LB V + z are gven by: LB S z = LB V z,..., z, h z, z +,..., z, 3. LB V + z = LB S z σ = λ S z, 3.2 where h z s the PGF mentoned n Property 2.. It s dscussed n the context of a pollng model wth smart customers n [5]. For gated servce, h z = β = λv z. For exhaustve servce, h z = π λv z, where π s the LST of a busy perod dstrbuton n an M/G/ system wth only type customers, so t s the root n 0, ] of the equaton π ω = β ω + λ V π ω, ω 0 cf. [2], p ow that we can relate LB V + z to LB V z, we can repeat ths and fnally obtan a recurson for LB V z. Ths 5

6 recursve expresson s suffcent to compute all moments of the ont queue length dstrbuton at a vst begnnng to Q by dfferentaton, but teraton of the expresson leads to the steadystate queue length dstrbuton at pollng epochs, wrtten as an nfnte product. We refer to [24] for more detals regardng ths approach, and for rgorous proofs of the laws of moton. Stablty condtons are studed n more detal n [], where t s shown that a necessary and suffcent condton for ergodcty s that the Perron-Frobenus egenvalue of the matrx R I should be less than 0, where I s the dentty matrx, and R s an matrx contanng elements ρ := λ V E[B ]. Ths holds under the assumpton that E[V ] > 0 for all =,...,. 3.2 Margnal queue length dstrbuton Common technques n pollng systems see, e.g. [3, 3] to determne the PGF of the steadystate margnal queue length dstrbuton of each customer type, are based on dervng the queue length dstrbuton at departure epochs. A level-crossng argument mples that the margnal queue length dstrbuton at arrval epochs must be the same as the one at departure epochs, and, fnally, because of PASTA ths dstrbuton s the same as the margnal queue length dstrbuton at an arbtrary pont n tme. In our model, the margnal queue length dstrbutons at arrval and departure epochs are also the same, but the dstrbuton at arbtrary moments s dfferent because of the varyng arrval rates durng a cycle. We can crcumvent ths problem by condtonng on the locaton P of the server P {V, S,..., V, S } and use condtonal PASTA to fnd the PGF of the margnal queue length dstrbuton at an arbtrary pont n tme. Let L denote the steady-state queue length of type customers at an arbtrary moment, and let L V and L S denote the queue length of type customers at an arbtrary tme pont durng V and S respectvely, =,...,. The followng relaton holds: E[z L ] = = [ ] E[V ] E[C] E z LV + E[S [ ] ] E[C] E z LS, =,...,. 3.3 ote that, at ths moment, E[V ] and E[C] are stll unknown. In Sectons 5 and 6 we llustrate two dfferent ways to compute them. Snce S, for =,...,, and V, for, are nonservng ntervals for customers of type, we use a standard result see, e.g., [3] to fnd the PGFs of L V and L S respectvely: [ ] E z LV = [ ] E z LS = z z E[z LBV ] E[z LBS ], =,..., ;, 3.4 E[LB S ] E[LB V ] E[z LBS ] E[z LBV + ],, =,...,, 3.5 E[LB V + ] E[LB S ] where LB P, for =,...,, are the number of type customers at the begnnng of perod V P {V, S,..., V, S }. Ther PGFs can be expressed n terms of LB z usng the relatons 3.2 and 3., and replacng argument z by the vector,...,, z,,..., where z s the element at poston. Usng branchng theory from [24], an explct expresson for 6

7 LB V V z s gven n [5]. The mean values, E[LB ] and E[LB S ], can be obtaned by dfferentaton of the correspondng PGFs and substtutng z =. ] It remans to compute E [z LV, =,...,,.e., the PGF of the number of type customers at an arbtrary pont wthn V. As far as the margnal queue length of type customers s concerned, the system can be vewed as a vacaton queue wth the ntervst tme I correspondng to the server vacaton. We can use the Fuhrmann-Cooper decomposton [5], but we have to be careful here. In a pollng system where type customers arrve wth constant arrval rate λ V, the Fuhrmann-Cooper decomposton states that ] ] [z LBS E [z LBV E[z L ] = λv E[B ] zβ λ V z V β λ z z E z E[LB V ] E[LB S ]. 3.6 The two parts n ths decomposton can be recognsed as the PGFs of the number of type customers respectvely at an arbtrary moment n an M/G/ queue, and at an arbtrary pont durng the ntervst tme I. Of course, the followng relaton also holds: Combnng 3.6 wth 3.7, results n: E[B ] λ V E[B ] E[z LV ] = λv E[z L ] = E[V V ] E[C] E[zL ] + E[I I ] E[C] E[zL ]. 3.7 z β λ V z β λ V z z E z ] [z LBS ] [z LBV E E[LB V ] E[LB S ], 3.8 for =,...,. The second part of ths decomposton s, agan, the PGF of the number of customers at an arbtrary pont durng the ntervst tme I. The frst part can be recognsed as the PGF of the queue length of an M/G/ queue wth type customers at an arbtrary pont durng a busy perod. ow we return to the model wth varyng arrval rates. The key observaton s that the behavour of the number of type customers durng a vst perod of Q, s exactly the same n ths system as n a pollng system wth constant arrval rates λ V for type customers. Equaton 3.8 no longer depends on anythng that happens durng the ntervst tme, because ths s all captured n LB V, the number of type customers at the begnnng of a vst to Q. Ths mples that, for a pollng model wth smart customers, the queue length PGF of Q at a random pont durng V s also gven by 3.8. The only dfference les n the nterpretaton of 3.8. Obvously, the frst part n 3.8 s stll the PGF of the queue length dstrbuton of an M/G/ queue at an arbtrary pont durng a busy perod. However, the last term can no longer be nterpreted as the PGF of the dstrbuton of the number of type customers at an arbtrary pont durng the ntervst tme I. Substtuton of 3.4, 3.5, and 3.8 n 3.3 gves the desred expresson for the PGF of the margnal queue length n Q. Remark 3. The margnal queue length PGF 3.3 has been obtaned by condtonng on the poston of the server at an arbtrary epoch n a cycle, whch explans the probabltes E[V ] E[C] server s servng Q and E[S ] E[C] server s swtchng to Q +. It s easy now 7

8 to obtan the margnal queue length PGF at an arrval epoch, smply by condtonng on the poston of the server at an arbtrary arrval epoch. The probablty that the server s at poston P {V, S,..., V, S } at the arrval of a type customer, s λp E[P ], wth λ E[C] λ = E[C] = λ V E[V ] + λ S E[S ]. Ths results n the followng expresson for the PGF of the dstrbuton of the number of type customers at the arrval of a type customer: E[z L arrval type ] = = λ V E[V ] λ E[C] [ ] E z LV + λs E[S ] E λ E[C] [ ] z LS, 3.9 for =,...,. A standard up-and-down crossng argument can be used to argue that 3.9 s also the PGF of the dstrbuton of the number of type customers at the departure of a type customer. As stated before, t s dfferent from the PGF of the dstrbuton of the number of type customers at an arbtrary epoch, unless λ V = λ S = λ for all, =,..., as s the case n pollng models wthout smart customers. Remark 3.2 Equatons 3.4 and 3.5 rely heavly on the PASTA property and are only vald f type arrvals take place durng the non-servng nterval. If no type arrvals take place.e. λ P = 0 for the non-servng nterval P, both the numerator and the denomnator become 0. Ths stuaton has to be analysed dfferently. ow assume that λ P = 0 for a specfc customer type =,...,, durng a non-servng nterval P {V, S,..., V, S }\V. We now dstngush between vst perods and swtch-over perods. Let us frst assume that P s a swtch-over tme, say S, =,...,. The length of a swtch-over tme s ndependent from the number of customers n the system, so the dstrbuton of the number of type customers at an arbtrary pont n tme durng S s the same as at the begnnng of S : [ ] [ ] E z LS = E z LBS,, =,...,. The case where P s a vst tme, say P = V for some, requres more attenton, because the length of V depends on the number of type customers present at the vst begnnng. Snce ths number s postvely correlated wth the number of customers n the other queues, we have to correct for the fact that t s more lkely that a random pont durng an arbtrary V, falls wthn a long vst perod wth more customers present at ts begnnng than n a short vst perod. The frst step, s to determne the probablty that the number of type customers at an arbtrary pont durng V s k. Snce we consder the case where λ V = 0, ths mples that we need the probablty that the number of customers at the begnnng of V s k. Standard renewal arguments yeld = k] = P[LBV = k] E[V LB V = k] l=0 P[LBV = l] E[V LB V = l] = E[V [LB V = k] ], E[V ] P[L V 3.0 where [A] s the ndcator functon for event A. The frst lne n 3.0 s based on the fact that the probablty s proportonal to the length of vst perods V that start wth k 8

9 type customers, and to the number of such vst perods V. The denomnator s smply a normalsaton factor. ow we can wrte down the expresson for the number of type customers at an arbtrary pont durng V f λ V = 0: E for =,..., and. [ ] z LV = k=0 = E[V ] z k P[L V = k] z k E[V [LB V = k] ] k=0 = E[V ] E[V z k [LB V = k] ] k=0 = E[V ] E[V z LB V ] = [ ] E[V ] ω E z LBV e ωv, 3. ω=0 ow we only need to determne E[z LBV e ωv ]. We use the ont queue length dstrbuton of all customers present at the begnnng of V, whch s gven mplctly by 3.2. Defne Θ as the tme that the server spends at Q due to the presence of one customer there, wth LST θ. For gated servce θ = β, and for exhaustve servce θ = π. The length of V, gven that l type customers are present at the vst begnnng, s the sum of l ndependent random varables wth the same dstrbuton as Θ, denoted by Θ,,..., Θ,l. The ont dstrbuton of the number of type customers present at the begnnng of V and the length of V s gven by: [ ] E z LBV e ωv = = = l =0 l =0 l =0 l =0 l =0 l =0 [ ] [ ] E z l e ωθ,+ +Θ,l P LB V = l, LB V = l z l E [ e ωθ ] ] [ ], E [e ωθ,l P LB V = l, LB V = l [ ] z l θ ω l P LB V = l, LB V = l = LB V,...,, z,,...,, θ ω,,...,, 3.2 where z corresponds to customers n Q, and θ ω corresponds to customers n Q. Substtuton of 3.2 n 3. gves the desred result. 4 Watng tme dstrbuton In the prevous secton we gave an expresson for the PGF of the dstrbuton of the steadystate queue length of a type customer at an arbtrary epoch, L. If the arrval rates do not 9

10 depend on the server poston,.e. λ V = λ S = λ for all, =,...,, we can use the dstrbutonal form of Lttle s law see, e.g., [9] to obtan the LST of the dstrbuton of the watng tme of a type customer, W, =,...,. Because of the varyng arrval rates, there s no λ for whch the relaton E[z L ] = E [ e λ zw +B ] holds even f we choose λ = λ. In the present secton, we ntroduce subtypes of each customer type. Each subtype s dentfed by the poston of the server at ts arrval n the system. We show that one can use a generalsed verson of the dstrbutonal form of Lttle s law that leads to the LST of the watng tme dstrbuton of a type customer, when appled to the PGF of the ont queue length dstrbuton of all subtypes of a type customer. Determnng ths PGF requres a separate treatment of exhaustve and gated servce, so results n ths secton do not apply to any arbtrary branchng-type servce dscplne. 4. Jont queue length dstrbuton at vst begnnngs/endngs for all subtypes In the present secton we dstngush between subtypes of type customers, arrvng durng dfferent vst/swtch-over perods. We defne a type P customer to be a customer arrvng at Q durng P {V, S,..., V, S }. Therefore, only n ths secton, we defne z n the followng way: z = z V,..., z S,..., z V,..., zs. ote that z has 2 2 components: customer types tmes 2 subperods wthn a cycle vst tmes plus swtch-over tmes. Let ṼBP z be the PGF of the ont queue length dstrbuton of all these customer types at the moment that the server starts servng type customers that have arrved when the server was located at poston P. ṼC P z s defned equvalently for the moment that the server completes servce of type P customers. For exhaustve servce, the vst perod V can be dvded nto the followng subperods: V = V S + V V V S + + V V. Frst the type S customers that were present at the vst begnnng are served, followed by the type V+ customers, and so on. ote that durng these servces only type V customers arrve n Q, =,...,. Vst perod V ends wth V V,.e., the exhaustve servce of all type V customers that have arrved durng V so far. The ont queue length process at pollng nstants of each of the subperods s stll a Mult-Type Branchng Process, because the servce dscplne n each subperod satsfes the Branchng Property. Hence, the laws of moton can be obtaned by applyng ths property successvely. As an example, we show the relatons for the PGFs of the ont queue length 0

11 dstrbutons at begnnngs and endngs of the subperods of V : ṼB V 2 z = ṼCS z = ṼBS z V, β λ V z V, z V 2,..., z S, ṼB S 2 z = ṼCV 2 z = ṼBV 2 ṼB V z = ṼCS z = ṼBS ṼC V. z V =,, β λ V z V, z S 2,..., z S, z V =,,...,, β λ V z V, z V 2,..., z S, z = ṼBV π λ V z V,,...,, z V 2,..., z S. Durng a swtch-over tme S only type S customers arrve,, =,...,. We can relate the PGF of the ont queue length dstrbuton at the begnnng of a vst to Q 2 startng wth the servce of type 2 S 2 customers to ṼCV z: ṼB S 2 2 z = ṼCV z σ = = λ S z S. The above expressons can be used to express ṼBS 2 2 n terms of ṼBS, and ths can be repeated to obtan a recurson for ṼBS. Remark 4. For gated servce we take smlar steps, but they are slghtly dfferent because arrvng customers wll always be served n the next cycle. Ths means that a vst to Q starts wth the servce of all type V customers present at that pollng nstant: V = V V +V S + V V V S +. The relatons for the PGF of the ont queue length dstrbuton at begnnngs and endngs of the subperods of V are: ṼB S z = ṼCV z = ṼBV ṼB V 2 z = ṼCS z = ṼBS. ṼC S z = ṼBS β λ V z V, z S,..., z S, z V =, β λ V z V, z V 2,..., z S, =,,...,, β λ V z V, z V 2,..., z S. z V The remander of ths secton s vald for any branchng-type servce dscplne treatng customers n order of arrval n each queue, such as, e.g., exhaustve, gated, globally gated and =

12 mult-stage gated [28]. Havng determned the ont queue length dstrbuton at begnnngs and endngs of all subperods wthn each vst perod, we are ready to determne the ont queue length dstrbuton at departure epochs of all customer subtypes. We follow the approach n [3, 4], whch tself s based on Esenberg s approach [3], developng a relaton between ont queue lengths at servce begnnngs/completons and vst begnnngs/endngs. In [3], for conventonal pollng systems, the ont dstrbuton of queue length vector and server poston at servce completons leads to the margnal queue length dstrbuton. Developng an equvalent for our model, requres dstngushng between customer subtypes. Frstly, the queue length vector z contans all customer subtypes. Secondly, the type of servce completon s not ust defned by the locaton of the server, but also by the subtype P of the customer that has been served. Therefore, let M P z denote the PGF of the ont dstrbuton of the subtypes of customers beng served combnaton of =,..., and P {V, S,..., V, S } and queue length vector of all customer subtypes at servce completons. Equaton 3.4 n [3], appled to our model, gves: M P z = λe[c] β = λv z V z P β = λv z V [ ṼB P z ṼCP ] z, 4. for =,..., ; P {V, S,..., V, S }, and λ = = λ. Thus, M P z s the generatng functon of the probabltes that, at an arbtrary departure epoch, the departng customer s a type P customer and the number of customers left behnd by ths departng customer s l V,..., l S. We now focus on the queue length vector of subtypes of type customers only, gven that the departure takes place at Q. The probablty that an arbtrary servce completon regardless of the subtype of the customer takes place at Q, s λ /λ. It s convenent to ntroduce the notaton z =,...,, z V,..., z S,,...,. The PGF of the ont queue length dstrbuton of all subtypes of type customers at an arbtrary departure from Q s: [ ] V E z V D S D = λ M V z + M S z 4.2 λ z S where D P s the number of type P customers left behnd at a departure from Q whch should not be confused wth L P, the number of type customers at an arbtrary moment whle the server s at poston P. Remark 4.2 Substtuton of z P = z for all P {V, S,..., V, S } n 4.2 gves the margnal queue length dstrbuton of type customers at departure epochs, whch s equal to 3.9, the margnal queue length dstrbuton at arrval epochs of a type customer. = 4.2 Watng tmes ow we present a generalsaton of the dstrbutonal form of Lttle s law that can be appled to the ont queue length dstrbuton of all subtypes of a type customer at departure epochs from Q, to obtan the watng tme LST of a type customer. 2

13 Theorem 4.3 The LST of the dstrbuton of the watng tme W of a type customer, =,...,, s gven by: E [ e ωw ] V D S = β ω E ω D λ V ω λ S. 4.3 Proof We focus on the departure of a type customer that arrved durng P A {V, S,..., V, S }. We make use of the fact that the soourn tme.e., watng tme plus servce tme of ths tagged type P A customer can be determned by studyng the subtypes of all type customers that he leaves behnd on hs departure. We need to dstngush between two cases, whch can be treated smultaneously, but requre dfferent notatons. Frstly, the case where a customer arrves n the system and departs durng another perod. In the second case, the customer departs durng the same perod n whch he arrved. Obvously, n our model the second case can only occur f a customer arrves at a queue wth exhaustve servce whle t s beng vsted by the server. Case : departure n a dfferent perod. In ths case we have that P A V, or P A = V but the cycle n whch the arrval took place s not the same as the cycle n whch the departure takes place ths stuaton cannot occur wth exhaustve servce. All type customers that are left behnd, have arrved durng the resdual perod P A, all perods between P A and V f any, and durng the elapsed part of V. Denote by P I the set of vst perods and swtch-over perods that le between P A and V. Furthermore, let P A,res be the resdual perod P A. Fnally denote by V,past the age of V at the departure nstant of the tagged type customer. Case 2: departure durng the perod of arrval. If the customer arrved durng the same vst perod n whch hs departure takes place, take P A,res = 0, P I =, and V,past s the tme that elapsed snce the arrval of the tagged type V customer. In both cases, the ont queue length dstrbuton of all customer subtypes at ths departure nstant s gven by 4.2. Snce we assume FCFS servce, at such a departure nstant no type customers are present anymore that have arrved before the arrval epoch of the tagged type customer. Ths results n: E [ z V D V z S S D ] [ = E e λp A z P A P A,res ] p P λ p I z p p λ V z V V,past 4.4 Equaton 4.3 follows from the relaton W + B = P A,res + p P I p + V,past and substtuton of z P = ω for all P {V λ P, S,..., V, S } n 4.4. Remark 4.4 Theorem 4.3 only holds f λ P > 0 for all =,...,, and P {V, S,..., V, S }. If λ P = 0 for a certan and P, we can stll fnd an expresson for E [ e ωw ], but we mght have to resort to some trcks. In Secton 8, Example 2, we show how the ntroducton of an extra vrtual customer type can help to resolve ths problem.. 3

14 5 Cycle tme, ntervst tme and vst tmes In the prevous sectons we repeatedly needed the mean cycle tme E[C] and the mean vst tmes E[V ], =,...,. In ths secton we study the LSTs of the cycle tme dstrbuton and vst tme dstrbutons, whch can be used to obtan the mean and hgher moments. The LSTs of the dstrbutons of the vst tmes V, =,...,, can easly be determned for any branchng-type servce dscplne usng the functon θ, ntroduced n Remark 3.2, and the ont queue length dstrbuton at the vst begnnng of Q not takng subtypes nto account: E[e ωv ] = LB V,...,, θ ω,,...,. 5. The cycle tme C s defned as the tme that elapses between two consecutve vst begnnngs to Q. We consder branchng-type servce dscplnes only,.e., servce dscplnes for whch Property 2. holds. The cycle tme LST for pollng models wth branchng-type servce dscplnes and arrval rates ndependent of the server poston, has been establshed n [0]. We adapt ther approach to the model wth arrval rates dependng on the server locaton. Usng θ, =,...,, we defne the followng functons n a recursve way: Smlarly, defne: ψ V ω = ω, ψ V ω = ω + ψ S ω = ω, ψ S ω = ω + k=+ k=+ λ V k θ k ψ Vk ω, =,...,. λ S k θ k ψ Vk ω, =,...,. Theorem 5. The LST of the dstrbuton of the cycle tme C s: E [ e ωc ] = LB V θ ψ V ω,..., θ ψ V ω σ ψ S ω. 5.2 Proof Smlar to the proof of Theorem 3. n [0], by gvng an expresson for the cycle tme LST condtoned on the numbers of customers n all queues at the begnnng of a cycle, and then by subsequently uncondtonng one queue at a tme. The LST of the dstrbuton of the ntervst tme I can be found n a smlar way: E [ e ωi ] = LB S, θ 2 ψ V2 ω,..., θ ψ V ω σ ψ S ω. 5.3 Equatons 5.2 and 5.3 hold for general branchng-type servce dscplnes. For gated and exhaustve servce we can gve expressons that are more compact and easer to nterpret, usng the ont queue length dstrbuton of all customer subtypes at vst begnnngs, as gven n Subsecton 4.. = = 4

15 Theorem 5.2 If Q receves exhaustve servce, the LST of the dstrbuton of the cycle tme C, startng at a vst endng to Q, and the LST of the dstrbuton of the ntervst tme I, are gven by: E [ e ωc ] = ṼB S,...,, π ω ω E [ e ωi ] = ṼBS ω,...,, λ V λ V,..., π ω ω,..., ω λ S λ S,,...,, 5.4,,...,, 5.5 provded that λ P 0 for all P {V, S,..., V, S }. In the rght hand sdes of 5.4 and 5.5, the components z P wth are. If Q receves gated servce, the LST of the dstrbuton of the cycle tme C, and the LST of the dstrbuton of the ntervst tme I, are gven by: E [ e ωc ] = ṼBV ω,...,, E [ e ωi ] = ṼBV λ V,...,,, ω λ V,..., ω λ S,..., ω λ S,,...,,,, ω λ S,..., ω λ S,,...,, 5.6 agan provded that λ P 0 for all P {V, S,..., V, S }. ote that z V = n 5.6. Proof We prove the exhaustve case only, the proof for gated servce proceeds along the same lnes. Usng I = S +V + +S + + +S +, and the fact that no type V customers are present at the begnnng of the ntervst perod and hence also at the begnnng of a cycle C, we obtan: ṼB S,...,, z V,..., z S,,..., = E [ e λs z S S λ S + z S + S + ]. 5.7 Substtuton of z P = ω for all P {V λ P, S,..., V, S } proves 5.5. Equaton 5.4 follows by usng the relaton C = I + V, and notng that V s the sum of the busy perods ntated by all type customers that have arrved durng I. In terms of LSTs: [ ] E e ωc = E [ = E e e ω+λ S π ω S ω+λ S + π ω λ S π ω ω λ S S λ S + π ω S + ] ω λ S + S +, whch, by 5.7, reduces to 5.4. Dfferentaton of the LSTs of C and C for =,...,, shows that, ust lke n pollng models wth constant arrval rates, the mean cycle tme does not depend on the startng pont of the cycle,.e. E[C ] = E[C ] = E[C]. The mean cycle tme E[C] and mean vst tmes E[V ] can be obtaned by dfferentatng the correspondng LSTs. In the next secton a more effcent method s descrbed to compute them. 5

16 6 Mean Value Analyss In ths secton we extend the Mean Value Analyss MVA framework for pollng models, orgnally developed by Wnands et al. [3], to sut the concept of smart customers. For ths purpose, we frst outlne the man deas of MVA for pollng systems. Subsequently, we determne the mean vst tmes and the mean cycle tme n a numercally more effcent way than n the prevous secton, and, fnally, we present the MVA equatons for a pollng system wth smart customers. 6. Man dea MVA For ordnary pollng models, where the arrval rates at a queue do not depend on the poston of the server, n [3] an approach s descrbed for dervng the steady-state mean watng tmes at each of the queues, E[W ] for =,...,, by settng up a system of lnear equatons, where each equaton has a probablstc and ntutve explanaton. We sketch the man deas of MVA for exhaustve servce; the cases of gated or mxed servce dscplnes requre only mnor changes. The mean watng tme E[W ] of a type customer, excludng hs servce tme, can be expressed n the followng way: upon arrval of a tagged type customer, he has to wat for the remanng tme t takes to serve all type customers already present n the system, plus possbly the tme before the server arrves at Q. By PASTA, the arrvng customer fnds n expectaton E[ˆL ] watng type customers n queue, each havng an expected servce tme E[B ]. ote that we use ˆL to denote the queue lengths excludng customers n servce. The expected tme untl the server returns to Q, s denoted by E[T ] whch depends on the servce dscplne of all queues. A fracton ρ := λ E[B ] of the tme, the server s servng Q, and hence, wth probablty ρ, an arrvng customer has to wat for a mean resdual servce tme, denoted by E[R B ]; otherwse he has to wat untl the server returns. Ths gves, for =,..., : E[W ] = E[ˆL ] E[B ] + ρ E[R B ] + ρ E[T ]. Lttle s law gves E[ˆL ] = λ E[W ], for =,...,, and so t remans to derve E[T ]. For ths, frst a system of equatons s composed for the condtonal mean queue lengths, whch can be expressed n mean resdual duratons of sums of vst and swtch-over tmes. The soluton of ths system of equatons can be used to determne E[T ], and hence E[ˆL ] and E[W ] follow. 6.2 Mean vst tmes and mean cycle tme For the case of smart customers, the vst tmes to a queue depend on all arrval rates λ V and λ S. In order to extend MVA to ths case, we frst derve the mean vst tmes to each of the queues, E[V ], for =,...,. We set up a system of lnear equatons where the mean vst tme of a queue s expressed n terms of the other mean vst tmes. We agan focus on the exhaustve servce dscplne. At the moment the server fnshes servng Q, there are no type customers present n the system any more. From ths pont on, the number of type customers bulds up at rates 6

17 λ S, λ V+,..., λ S+ dependng on the poston of the server, untl the server starts workng on Q agan. Each of these customers ntates a busy perod, wth mean E[BP ] := E[B ]/ λ V E[B ]. Ths gves: E[V ] = E[BP ] λ S ES + + =+ λ V E[V ] + λ S E[S ], for =,...,. The E[V ] follow from solvng ths set of equatons. Ths method s computatonally faster than determnng and dfferentatng the LSTs of the vst tme dstrbutons 5.. Once the mean vst tmes have been obtaned, the mean cycle tme follows from E[C] = = E[V ] + E[S ]. 6.3 MVA equatons We extend the MVA approach to pollng systems wth smart customers. Frst, we brefly ntroduce some extra notaton, then we gve expressons for the mean watng tmes, and the mean condtonal and uncondtonal queue lengths. After elmnatng varables, we end up wth a system of lnear equatons. The system can numercally be solved n order to fnd the unknowns, n partcular, the mean uncondtonal queue lengths and the mean watng tmes. Although all equatons are dscussed n the present secton, for the sake of brevty of ths secton, some of them are presented n Appendx A. The fracton of tme the system s n a gven perod P {V, S,..., V, S } s denoted by q P := E[P ] E[C]. The mean resdual duraton of a perod P, at an arbtrarly chosen pont n ths perod, s denoted by E[R P ] = E[P 2 ] 2E[P ]. The mean condtonal number of type customers n the queue durng a random pont n P s denoted by E[ˆL P number of type customers n queue s denoted by E[ˆL ]. ote that ˆL and ], and the mean uncondtonal do not nclude a potental customer n servce, whereas L and L P, ntroduced n Secton 3, denote queue lengths ncludng customers beng served. We defne an nterval, e.g. V : S, as the consecutve perods from the frst mentoned perod on, untl and ncludng the last mentoned perod, here consstng of the perods V, S, V +, S +,..., V, S. The mean resdual duraton of an nterval, e.g. V : S, s denoted by E[R V :S ]. Analogously, we defne E[R V :V ], E[R S :V ] and E[R S :S ]. An mportant concept n the remander of the analyss s the concept of condtonal duratons of a perod. Ths s an extenson of the well-known resdual duraton, or the age of a perod. It deals wth the length of a perod wthn the cycle.e., a vst tme or a swtch-over tme, gven that the system s beng observed from another perod. Before we proceed, we clarfy ths mportant concept by a smple example. Consder a vacaton system,.e., a pollng system wth =. A cycle conssts of a swtch-over tme or: vacaton S, followed by a vst tme V, durng whch the queue s served exhaustvely. ow assume that the system s beng observed at a random epoch durng the swtch-over tme S. We derve an expresson for E[ V S ], whch s the condtonal mean vst tme followng the swtch-over tme, gven that the system s beng observed durng S. Snce servce s exhaustve, the vst tme conssts of the busy perods of the customers that arrved durng the elapsed part of S, denoted by ˆL P 7

18 S,past, plus the busy perods of the customers that wll arrve durng the resdual swtch-over tme, denoted by S,res. Hence, t can be seen that n ths system E[ V S ] = λs E[B ] λ V E[B ] E[S,past] + E[S,res ] = λ S E[B ] λ V E[B ] E[S 2 ] E[S ]. Instead of studyng the mean vst tme followng the swtch-over tme durng whch the system s observed, we can also study the mean vst tme precedng ths partcular swtchover tme, denoted by E[ V S ]. ow the expresson s easer, because a swtch-over tme s ndependent of the precedng vst tme, so E[ V S ] = E[V ] = λs E[B ] λ V E[B ] E[S ]. Smlarly, we denote by E[ S V ] and E[ S V ] the mean swtch-over tmes followng, respectvely precedng, the vst perod V durng whch the system s observed. Because of the ndependence between S and the precedng V, t s mmedately clear that E[ S V ] = E[S ]. Unfortunately, E[ S V ] cannot be determned ths easly, snce V s postvely correlated wth the precedng swtch-over tme. However, t plays a role n the set of MVA equatons developed later n ths secton, relatng the condtonal mean queue lengths, the condtonal mean watng tmes and the condtonal duratons of the perods n a cycle. We leave t up to the reader to verfy that Equaton 6.8 reduces to E[ S V ] = E[S 2]/E[S ] n the exhaustve vacaton model. Ths example smply serves the purpose of llustratng the concept of these condtonal duratons. In a pollng system consstng of multple queues, these expressons become more complcated and can only be found by solvng sets of equatons, as wll be shown n the remander of ths secton. ote that, because of condtonal PASTA, an arbtrary customer arrvng durng S fnds the system n the same state as an observer who observes the system at an arbtrary epoch durng S. Hence, the condtonal duratons of perods play an mportant role n determnng the mean watng tmes. For the mean condtonal duratons of a perod, we have the followng: E[ V V ] denotes the mean duraton of the prevous perod V, seen from an arbtrary pont n V.e., V seen backwards n tme from the vewpont of V, and E[ V V ] denotes the mean duraton of the next perod V.e., V seen forwards n tme from the vewpont of V. For = they both concde, and represent the mean age, resp. the mean resdual duraton of V. Snce the dstrbuton of the age of a perod s the same as the dstrbuton of the resdual perod, we have E[ V V ] = E[ V V ] = E[R V ]. Generally, however, E[ V V ] E[ V V ] for, because of the dependences between the duratons of perods. Analogously, we defne E[ V S ], E[ V S ], E[ S V ] and E[ S V ]. ote that, e.g., E[ S V ] = E[S ], but E[ S V ] E[S ]. As swtch-over tmes are ndependent, the followng quanttes drectly smplfy: E[ S S ] = E[ S S ] = { E[S ] for, E[R S ] for =. 8

19 Havng ntroduced the requred notaton, we now present the man theorem of ths secton, whch gves a set of equatons that can be solved to fnd the mean watng tmes of customers n the system. Theorem 6. The mean watng tmes, E[W ], for =,...,, and the mean queue lengths, E[ˆL ], satsfy the followng equatons: E[W ] = qv λ V E[ˆL V ]E[B ] + E[R B ] λ + + =+ + + = E[ˆL ] = λ E[W ], E[ˆL ] = + =+ q V λ V λ q S λ S λ E[ˆL V ]E[B ] + + k= E[ˆL S ]E[B ] + E[R S ] + E[S k ] + E[ V k V ] + k=+ E[S k ] + E[ V k S ], q V E[ˆL V ] + q S E[ˆL S ], 6.3 where the condtonal mean queue lengths E[ˆL V ] and E[ˆL S ], for = +,..., +, are gven by E[ˆL V ] = E[ˆL S ] = k=+ k=+ λ V k E[ V V k ] + λ S k E[ V S k ], 6.4 λ V k E[ V k S ] + k= k= λ S k E[ S k S ], 6.5 and where all E[ P P 2 ] and E[ P P 2 ], for P, P 2 {V, S,..., V, S }, satsfy the set of equatons below, and A.2 A.7 n Appendx A. Proof In order to derve the mean watng tme E[W ], we condton on the perod n whch a type customer arrves. A fracton q V λ V /λ, and q S λ S /λ respectvely, of the type customers arrves durng perod V, and durng perod S respectvely. If a tagged type customer arrves durng perod V.e., whle hs queue s beng served, he has to wat for a resdual servce tme, plus the servce tmes of all type customers present n the system upon hs arrval, whch s by condtonal PASTA, E[ˆL V ]. As a fracton q V λ V /λ of the customers arrves durng V, ths explans the frst lne of 6.. If the customer arrves n any other perod, he has to wat untl the server returns to Q agan. For ths, we condton on the perod n whch he arrves. If the arrval perod s a vst to Q, say V for, he has to wat for the resdual duraton of V and the nterval S : S, and for the servce of the type customers present n the system upon hs arrval. Ths gves the second lne of 6.. The thrd lne, the case where the customer arrves durng the swtch-over tme from Q to Q + perod S, can be nterpreted along the same lnes as the case V. Equaton 6.3 s obtaned by uncondtonng the condtonal queue lengths E[ˆL P ]. The mean number of type customers n the queue at an arbtrary pont durng V, gven by 9

20 6.4, s the mean number of customers bult up from the last vst to Q when Q became empty untl and ncludng a resdual duraton of V as the mean resdual duraton of V s equal to the mean age of that perod, takng nto account the varyng arrval rates. The mean number of type customers queueng n the system durng perod S, gven by 6.5, can be found smlarly. Equatons 6.4 and 6.5 show one of the dffcultes n adaptng the ordnary MVA approach to that of smart customers. If the arrval rates reman constant durng a cycle, these expressons would reduce to λ multpled by the mean tme passed snce the server has left Q. However, for the smart customers case, we have to keep track of the duraton of all the ntermedate perods, from the vewpont of perod V respectvely S. As ndcated n Theorem 6., at ths pont, the number of equatons s nsuffcent to fnd all the unknowns, E[ P P 2 ] and E[ P P 2 ], for P, P 2 {V, S,..., V, S }. In the remander of the proof, we develop addtonal relatons for these quanttes to complete the set of equatons. We start by consderng E[ V V ], whch s the mean duraton of the next perod V, when observed from an arbtrary pont n V. For = ths s ust the resdual duraton of V, consstng of a busy perod nduced by a customer wth a resdual servce tme left, and the busy perods of all type customers n the queue. The cases need some more attenton. The duraton of V now conssts of the busy perod nduced by the type customers n the queue, whch are n expectaton E[ˆL V ] customers. Durng the perods V, S,..., S, however, new type customers are arrvng, each contrbutng a busy perod to the duraton of V. Hence, summng over these perods and takng nto account the varyng arrval rates, we get the mean total of newly arrvng customers n ths nterval. Ths yelds, for =,..., and = +,..., + : E[ V V ] = E[BP ] E[ˆL V ] + E[R B ]/ λ V E[B ], 6.6 E[ V V ] = E[BP ] E[ˆL V ] + + k= λ V k E[ V k V ] + λ S k E[S k ]. 6.7 Analogously E[ V S ] denotes the mean duraton of the next perod V, when observed from an arbtrary pont n S. The explanaton of ts expresson s along the same lnes as that of E[ V V ], although t should be noted that = s not a specal case. See A. n Appendx A. The last step n the proof of Theorem 6., needs the followng lemma to fnd the fnal relatons between E[ P P 2 ] and E[ P P 2 ]: 20

21 Lemma 6.2 For =,...,, and = +,..., + : E[S k ] E[ k S S k ] + E[ S S k l ] + E[ S V k l ] E[S :V ] = k= k=+ l=+ E[R Sk ] E[ V S k ] E[V k ] E[ V V k ] + E[S :V ] k= l=k E[ S V k ] E[ V k V k ] l=k+ E[S l ] + E[ V V k l ] k l=+ E[S l ] + E[ S V k l ] E[ S l V k ] + E[ V l V k ]. 6.8 Proof Equaton 6.8 can be proven by studyng all mean resdual nterval lengths E[R S :V ], E[R S :S ], E[R V :V ] and E[R V :S ]. Consder E[R S :V ], the mean resdual duraton of the nterval S, V +,..., V. We condton on the perod n whch the nterval s observed. As the mean duraton of the nterval s gven by E[S :V ], t follows that E[S k ]/E[S : V ] s the probablty that the nterval s observed n perod S k. The remanng duraton of the nterval conssts of the remanng duraton of S k plus the mean duratons of the comng perods V k+, S k+,..., V, when observed from perod S k. When observng E[S :V ] from V k, a smlar way of reasonng s used. Ths gves, for =,...,, and = +,..., + : E[S k ] E[R S :V ] = E[R Sk ] + E[ S V k ] + E[S l ] + E[ S V k l ] E[S :V ] + k=+ E[V k ] E[S :V ] E[ V V k ] + l=k l=k+ E[S l ] + E[ V l V k ]. 6.9 We now use that the dstrbuton of the resdual length of an nterval s the same as the dstrbuton of the age of ths nterval. Agan, focus on E[R S :V ], condtonng on the perod n whch the nterval s observed, but now lookng forwards n tme. Consder all the perods n S : V that have already passed when observng durng S k. The nterval has lasted for the sum of these perods, plus the age of S k. The same can be done for an arbtrary pont n V k. Ths gves, for =,...,, = +,..., + : E[R S :V ] = k= + E[S k ] E[S :V ] k=+ E[V k ] E[S :V ] E[ S S k ] + k l=+ E[ S V k ] + E[ V k V k ] + E[ S S k l ] + E[ S V k l ] k l=+ E[ S l V k ] + E[ V l V k ]. 6.0 The proof of Lemma 6.2 s completed by equatng 6.9 and 6.0 and rearrangng the terms. Smlar to the proof of Lemma 6.2, we can develop two dfferent expressons for each of the terms E[R S :S ], E[R V :V ] and E[R V :S ]. For the sake of brevty of ths secton, they are 2