arxiv: v1 [math.pr] 19 Oct 2016

Transcription

1 Performance analyss of pollng systems wth retrals and glue perods Murtuza Al Abdn, Onno Boxma, Bara Km, Jeongsm Km, Jacques Resng EURADOM and Department of Mathematcs and Computer Scence arxv: v1 [math.pr] 19 Oct 2016 Endhoven Unversty of Technology P.O. Box 513, 5600 MB Endhoven, The etherlands e-mal: Department of Mathematcs, Korea Unversty 145 Anam-ro, Seongbuk-gu, Seoul, 02841, Korea e-mal: Department of Mathematcs Educaton, Chungbuk atonal Unversty 1 Chungdae-ro, Seowon-gu, Cheongju, Chungbuk, 28644, Korea e-mal: jeongsmkm@chungbuk.ac.kr Abstract We consder gated pollng systems wth two specal features: () retrals, and() glue or reservaton perods. When a type- customer arrves, or retres, durng a glue perod of staton, t wll be served n the next vst perod of the server to that staton. Customers arrvng at staton n any other perod jon the orbt of that staton and wll retry after an exponentally dstrbuted tme. Such pollng systems can be used to study the performance of certan swtches n optcal communcaton systems. For the case of exponentally dstrbuted glue perods, we present an algorthm to obtan the moments of the number of customers n each staton. For generally dstrbuted glue perods, we consder the dstrbuton of the total workload n the system, usng t to derve a pseudo conservaton law whch n ts turn s used to obtan accurate approxmatons of the ndvdual mean watng tmes. We also consder the problem of choosng the lengths of the glue perods, under a constrant on the total glue perod per cycle, so as to mnmze a weghted sum of the mean watng tmes. Keywords: Pollng system, Retrals, Glue perods 1

2 1 Introducton Ths paper s devoted to the performance analyss of a class of sngle server queueng systems wth multple customer types. Our motvaton s twofold: () to obtan nsght nto the performance of certan swtches n optcal communcaton systems, and () to obtan nsght nto the effect of havng partcular reservaton perods, wndows of opportunty durng whch a customer can make a reservaton for servce. Our class of queueng systems combnes several features, vz., pollng, retrals, and the new feature of so-called glue perods or reservaton perods. These wll frst be dscussed separately, whle ther relaton to optcal swtchng wll also be outlned. Pollng systems are queueng models n whch a sngle server, alternatngly, vsts a fnte number of, say, queues (or statons) n some prescrbed order. Pollng systems have been extensvely studed n the lterature. For example, varous dfferent servce dscplnes (rules whch descrbe the server s behavour whle vstng a queue) have been consdered, both for models wth and wthout swtchover tmes between queues. We refer to Takag [18, 19] and Vshnevsk and Semenova [20] for lterature revews and to Boon, van der Me and Wnands [4], Levy and Sd [13] and Takag [17] for overvews of the applcablty of pollng systems. Swtches n communcaton systems form an mportant applcaton area of pollng systems. Here, packets must be routed from source to destnaton, passng through a seres of lnks and nodes. In copper-based transmsson lnks, packets from varous sources are tme-multplexed, and ths may be modelled by a pollng system. In recent years optcal networkng has become very mportant, because optcal fbers offer major advantages wth respect to copper cables: huge bandwdth, ultra-low losses and an extra dmenson, vz., a choce of wavelengths. When one wants to model the performance of an optcal swtch by a pollng system [14, 16], one s faced wth the followng dffculty. Bufferng of optcal packets s not easy, as photons can not wat. Whenever there s a need to buffer photons, they are sent nto a local fber loop, thus provdng a small delay to the photons wthout losng or dsplacng them. If, at the completon of the loop, a photon stll needs to be buffered, t s agan sent nto the fber delay loop, etc. From a queueng theoretc perspectve, 2

3 ths rases the need to add the feature of retral queue to a pollng system: nstead of havng a queueng systemwth oneserverand ordnaryqueues, t hasoneserverand retralqueues. Retralqueueshave receved much attenton n the lterature, see, e.g., the books by Faln and Templeton [9] and by Artalejo and Gomez-Corral [3], but they have hardly been studed n the settng of pollng models. Langars [10, 11, 12] has poneered the study of pollng models wth retral queues. However, for our purpose - the performance analyss of optcal swtches - hs assumptons about the servce dscplne of the server at the varous queues are not sutable. A thrd mportant feature n the present paper s that of so-called glue or reservaton perods. Just before the server arrves at a staton there s some glue perod. Customers (both new arrvals and retryng customers) arrvng at the staton durng ths glue perod stck and wll be served durng the vst of the server. Customers arrvng at the staton n any other perod jon the orbt of that staton and wll retry after an exponentally dstrbuted tme. One motvaton for studyng glue perods s the followng. A sophstcated technology that one mght try to add to the use of fber delay loops n optcal networkng s varyng the speed of lght by changng the refractve ndex of the fber loop, cf. [15]. Usng a hgher refractve ndex n a small part of the loop one can acheve slow lght, whch mples slowng down the packets. Ths feature s n our model ncorporated as glue perods, where we slow down the packets arrvng at the end of the fber loop just before the server arrves, so that they do not have to retry but get served durng the subsequent vst perod. ot restrctng ourselves to optcal networks, one can also nterpret a glue perod as a reservaton perod,.e., a perod n whch customers can make a reservaton at a staton for servce n the subsequent vst perod of that staton. In our model, the reservaton perod mmedately precedes the vst perod, and could be seen as the last part of a swtchover perod. A frst attempt to study a pollng model whch combnes retrals and glue perods s [7], whch manly focuses on the case of a sngle server and a sngle staton, but also outlnes how that analyss can be extended to the case of two statons. In [1] an -staton pollng model wth retrals and wth constant glue perods s consdered, for the case of gated servce dscplne at all statons. The gated dscplne s an mportant dscplne n pollng systems; t mplcates that the server, when vstng a staton, serves exactly those customers whch were present upon hs arrval. The steady-state jont staton sze (.e., 3

4 the number of customers n each staton) dstrbuton was derved n [1], both at an arbtrary epoch and at begnnngs of swtchover, glue and vst perods. In the current paper we present an algorthm to obtan the moments of the staton sze for the case of exponentally dstrbuted glue perods. Usng Lttle s law, that also gves mean sojourn tmes. Thereafter for each ndvdual staton we allow generally dstrbuted glue perods and we focus our attenton on other performance measures next to staton szes. In partcular, we consder the steady-state dstrbuton of the total workload n the system, whch leads us to a pseudo conservaton law,.e., an exact expresson for a weghted sum of the mean watng tmes at all statons. We use that pseudo conservaton law to derve an accurate approxmaton for the ndvdual mean watng tmes. We further consder the problem of choosng the lengths of the glue perods, gven the total glue perod n a cycle, so as to mnmze a weghted sum of the mean watng tmes. The rest of the paper s organzed as follows. In Secton 2 we present a model descrpton. Secton 3 contans a detaled analyss for generatng functons and moments of staton szes at dfferent tme epochs when the glue perods are exponentally dstrbuted. We also present a numercal example, whch n partcular provdes nsght nto the behavor of the pollng system n the case of long glue perods. In Secton 4 we derve a pseudo conservaton law for a system wth generally dstrbuted glue perods. Subsequently we use ths pseudo conservaton law for dervng an approxmaton for the mean watng tmes at all statons. In ts turn, ths approxmaton s used to mnmze weghted sums of the mean watng tmes by optmally choosng the lengths of the glue perods, gven the total glue perod per cycle. Fnally, Secton 5 lsts some topcs for further research. 2 The model We consder a sngle server cyclc pollng system wth retrals and so-called glue perods. Ths model was frst ntroduced n [7] for a sngle staton vacaton model and a two-staton model wth swtchover tmes. Further n [1] ths was extended to an -staton model wth swtchover tmes. In both papers, the model was studed for determnstc glue perods. We ndex the statons by, = 1,...,, n the order of server movement. For ease of presentaton, all references to staton ndces greater than or 4

5 less than 1 are mplctly assumed to be modulo. Customers arrve at staton accordng to a Posson process wth rate λ ; they are called type- customers, = 1,...,. The overall arrval rate s denoted by λ = λ 1 + +λ. The servce tmes at staton are ndependent and dentcally dstrbuted (..d.) random varables wth a generc random varable B, = 1,...,. Let B (s) = E[e sb ] be the Laplace- Steltjes transform (LST) of the servce tme dstrbuton at staton. The swtchover tmes from staton to staton +1 are..d. random varables wth a generc random varable S. Let S (s) = E[e ss ] be the LST of the swtchover tme from staton to staton +1, = 1,...,. The nterarrval tmes, the servce tmes, and the swtchover tmes are assumed to be mutually ndependent. After a swtch of the server to staton, there s a glue perod for collectng retryng customers (whch wll be followed by the vst perod of the server to staton ). We assume that the successve glue perods at staton are..d. random varables wth a generc random varable G. Let G (s) = E[e sg ] be the LST of the glue perod dstrbuton at staton. Each staton conssts of an orbt and a queue. When customers (both new arrvals and retryng customers) arrve at staton durng a glue perod, they stck and wat n the queue to get served durng the vst of the server to that staton. When customers arrve at staton n any other perod, they jon the orbt of staton and wll retry after a random amount of tme. The nter-retral tme of each customer n the orbt of staton s exponentally dstrbuted wth mean ν 1 and s ndependent of all other processes. A sngle server cyclcally moves from one staton to another servng the glued customers at each of the statons. The servce dscplne at all statons s gated. Durng the servce perod of staton, the server serves all glued customers n the queue of staton,.e., all type- customers watng at the end of the glue perod (but none of those n orbt, and nether any new arrvals). Let (X () 1,X() 2,...,X() ) denote the vector of numbers of customers of type 1 to type n the system (hence n the orbt) at the start of a glue perod of staton, = 1,...,, n steady state. Further, let (Y () 1 (),Y 2 (),...,Y ) denote the vector of numbers of customers of type 1 to type n the system at the start of a vst perod at staton, = 1,...,, n steady state. We dstngush between those who are queueng (glued) and those who are n the orbt of staton : We wrte Y () = Y (q) +Y (o), = 1,...,, 5

6 where q denotes n the queue and o denotes n the orbt. Fnally, let (Z () 1,Z() 2,...,Z() ) denote the vector of numbers of customers of type 1 to type n the system (hence n the orbt) at the start of a swtchover from staton to staton +1, = 1,...,, n steady state. The utlzaton of the server at staton, ρ, s defned by ρ = λ E[B ] and the total utlzaton of the server ρ s gven by ρ = =1 ρ. It can be shown that a necessary and suffcent condton for stablty of ths pollng system s ρ < 1. We hence assume that ρ < 1. The cycle length of staton, = 1,..., s defned as the tme between two successve arrvals of the server at ths staton. The mean cycle length, E[C], s ndependent of the staton nvolved (and the servce dscplne) and s gven by E[C] = =1 (E[G ]+E[S ]), (2.1) 1 ρ whch can be derved as follows: Snce the probablty of the server beng dle (n steady state) s 1 ρ, and ths equals =1 (E[G]+E[S]) E[C] by the theory of regeneratve processes, we have (2.1). 3 The pollng system wth retrals and exponental glue perods In [1] the authors calculated the generatng functons and the mean values of the number of customers at dfferent tme epochs when the glue perods are determnstc. In ths secton we assume that the glue perods are exponentally dstrbuted wth mean E[G ] = 1/γ, = 1,...,. We wll derve a set ofpartal dfferental equatons for the jont generatng functon of the staton sze (.e., the number of customers n each staton) and then obtan a system of lnear equatons for the frst and the second moments of the staton sze. We also provde an teratve algorthm for solvng the system of lnear equatons. Observe that the generatng functon for the vector of numbers of arrvals at staton 1 to staton durng the servce tme of a type- customer, B, s β (z) := B ( j=1 λ j(1 z j )) for z = (z 1,z 2,...,z ). Smlarly, the generatng functon for the vector of numbers of arrvals at staton 1 to staton durng a swtchover tme from staton to staton +1, S, s σ (z) := S ( j=1 λ j(1 z j )). 6

7 3.1 Staton sze analyss at embedded tme ponts In ths subsecton we study the steady-state jont dstrbuton and the mean of the numbers of customers n the system at the start of a glue perod, vst perod and swtchover perod. Let us defne the followng jont generatng functons of the number of customers n each staton at the start of a glue perod, vst perod and swtchover perod: R () R () g 1 (z) = E[zX() 1 z X() 2 2 z X() ], v (z,w) = E[z Y () 1 1 z Y () 2 2 z Y (o) R () s (z) = E[z Z() 1 1 z Z() 2 2 z Z() ], z Y () (q) wy ], for z = (z 1,z 2,...,z ) wth z 1, = 1,...,, and w 1. Let M o (t) represent the number of customers n the orbt of staton, = 1,, and Υ(t) the number of glued customers, at tme t. Further, let τ j be the tme at whch an arbtrary glue perod starts at staton j, j = 1,...,. ote that M o (τ j) = X (j). We defne φ (z;w;t) = E[z Mo 1 (τ+t) φ (z,w) = 1 z Mo (τ+t) 0 φ (z;w;t)dt. w Υ(τ+t) ½ {G>t}], Then all the generatng functons for the numbers of customers n steady state descrbed above, can be expressed n terms of φ (z,w), as shown below n Proposton 1. Proposton 1. The generatng functons () () R v (z,w), R s (z) and R () g (z) satsfy the followng: R () v (z,w) = γ φ (z,w), (3.1) R () s (z) = γ φ (z,β (z)), (3.2) R () g (z) = γ 1 σ 1 (z)φ 1 (z,β 1 (z)). (3.3) 7

8 Proof. Equaton (3.1) s obtaned as follows: By the law of total expectaton, R () v (z,w) = = E[z Mo 1 (τ+t) 1 z Mo (τ+t) 0 E[z Mo 1 (τ+t) 1 z Mo (τ+t) 0 = γ φ (z,w). w Υ(τ+t) G > t]γ e γt dt w Υ(τ+t) ½ {G>t}]γ dt To obtan (3.2), observe that the customers at the end of a vst perod are the customers n the orbt at the begnnng of that vst plus the customers who arrve durng the servce tmes of the glued customers at the begnnng of that vst. Hence R () s () 1 (z) = E[zY 1 z Y () 2 2 z Y (o) z Y () [β (z)] Y (q) ] = R () v (z,β (z)) = γ φ (z,β (z)). Also, to obtan (3.3), observe that the customers at the end of a swtchover from staton 1 to staton are the customers n the orbt at the begnnng of that swtchover plus the customers who arrved durng that swtchover perod. Hence R () g (z) = ( 1) R s (z)σ 1 (z), from whch and (3.2) we get (3.3). We have the followng result for the generatng functons φ (z,w), = 1,...,. Theorem 1. The generatng functons φ (z,w), = 1,...,, satsfy the followng equaton: ν (w z ) ( φ (z,w) (λ j (1 z j ))+λ (1 w)+γ )φ (z,w) z j=1,j +γ 1 φ 1 (z,β 1 (z))σ 1 (z) = 0. (3.4) 8

9 Proof. ote that φ (z;w;t+ t) = E [ z Mo 1 (τ+t+ t) = = 1 z Mo (τ+t+ t) n 1=0 n =0k=0 E [ z Mo 1 (τ+t+ t) n 1=0 w Υ(τ+t+ t) ] ½ {G>t+ t} P(M1 o (τ +t) = n 1,...,M o (τ +t) = n,υ(τ +t) = k,g > t) 1 z Mo (τ+t+ t) M o (τ +t) = n,υ(τ +t) = k,g > t ] n =0k=0 w Υ(τ+t+ t) ½ {G>t+ t} M1 o (τ +t) = n 1,..., P(M1(τ o +t) = n 1,...,M(τ o +t) = n,υ(τ +t) = k,g > t) z n1 1 zn 1 zn+1 zn wk ((1 e ν t )w +e ν t z ) n e ( j=1,j (λj(1 zj))+λ(1 w)) t e γ t = e ( j=1,j (λj(1 zj))+λ(1 w)+γ ) t φ (z 1,...,z 1,z +(1 e ν t )(w z ),z +1,...,z ;w;t). Thus, we have t φ (z;w;t) = ν (w z ) ( φ (z;w;t) z Snce φ (z;w;0) = E[z X() 1 1 z X() 2 2 z X() j=1,j ntegratng the above equaton wth respect to t from 0 to yelds (λ j (1 z j ))+λ (1 w)+γ )φ (z;w;t). ] = R () g (z) = γ 1 σ 1 (z)φ 1 (z,β 1 (z)) and φ (z;w; ) = 0, γ 1 σ 1 (z)φ 1 (z,β 1 (z)) =ν (w z ) z φ (z,w) Ths completes the proof. ( j=1,j λ j (1 z j )+λ (1 w)+γ )φ (z,w). We now calculate the mean value of the staton szes at embedded tme ponts usng the dfferental equaton (3.4). For an -tuple l = (l 1,...,l ) of nonnegatve ntegers, we defne l = l 1 + +l, l! = l 1!l 2! l!, and z l = z l1 1 zl2 2 zl. Wth ths notaton, we defne the followng scaled moment: Φ (l,m) = 1 l!m! l +m z l w mφ (z,w), z=1,w=1 9

10 where z l = z l1 1 zl, and 1 s the -dmensonal row vector wth all ts components equal to one. The frst scaled moments of φ (z,w), = 1,2,...,, can be obtaned from the followng theorem. Theorem 2. We have () Φ (0,0) = 1 γ, = 1,...,. () Φ (1j,0) and Φ (0,1), 0,j, are gven by the followng recurson: for j = 1,...,, Φ (0,1) j Φ (1j,0) j = λ j γ j E[C], (3.5) = λ j ν j ( E[C] 1 γ j ), (3.6) Φ (1j,0) = γ +1 Φ (1j,0) +1 + λ ( j (δ,j 1 ρ )E[C] 1 ) E[S ], = j 1,j 2,...,j +1, γ γ γ +1 where 0 s the -dmensonal row vector wth all ts elements equal to zero, 1 j s the -dmensonal row vector whose jth element s one and all other elements are zero, and δ j s the Kronecker delta. ote that f s nonpostve n (3.7), then t s nterpreted as +. (3.7) Proof. Takng the partal dervatve of Equaton (3.4) wth respect to z j and puttng z = 1,w = 1, we have ν δ j Φ (1,0) + (1 δ j)λ j γ γ Φ (1j,0) +γ 1 Φ (1j,0) 1 +γ 1 λ j E[B 1 ]Φ (0,1) 1 +λ je[s 1 ] = 0. (3.8) Takng the partal dervatve of Equaton (3.4) wth respect to w and puttng z = 1,w = 1 yelds ν Φ (1,0) Summng (3.8) over = 1,...,, we have ν j Φ (1j,0) j +λ j j 1 γ +λ j + λ γ γ Φ (0,1) = 0. (3.9) =1 γ E[B ]Φ (0,1) +λ j Addng (3.9) and (3.10) and multplyng the resultng equaton by E[B j ] yelds ρ j =1 ( 1 +E[S ]) γ j E[B j ]Φ (0,1) j +ρ j γ =1 E[S ] = 0. (3.10) =1 γ E[B ]Φ (0,1) = 0, 10

11 and summng ths over j = 1,..., gves where we have used (2.1). Pluggng (3.11) nto (3.10) leads to =1 Φ (1j,0) j γ E[B ]Φ (0,1) = ρe[c], (3.11) = λ j ν j ( E[C] 1 γ j ), whch s (3.6). Insertng ths equaton nto (3.9) yelds (3.5). When = j n Equaton (3.8), we have Φ (1j,0) j 1 = γ j Φ (1j,0) j + λ ( j (1 ρ j 1 )E[C] 1 ) E[S j 1 ]. (3.12) γ j 1 γ j 1 γ j On the other hand, when j,.e., = j 1,j 2,...,j +1, n Equaton (3.8), we have Φ (1j,0) 1 = γ Φ (1j,0) + λ ( j ρ 1 E[C] 1 ) E[S 1 ]. (3.13) γ 1 γ 1 γ Fnally, (3.7) follows from (3.12) and (3.13). ext, we calculate Φ (l,m) for l +m 2. Equaton (3.4) can be wrtten as ( ν (w 1) ν (z 1) ) ( φ (z,w)+ λ (w 1)+ z +γ 1 φ 1 (z,β 1 (z))σ 1 (z) = 0. j=1,j λ j (z j 1) γ )φ (z,w) From ths we get (γ +l ν )Φ (l,m) =½ {m 1} (l +1)ν Φ (l+1,m 1) + j ½ {lj 1}λ j Φ (l 1j,m) + ½ {m 1} λ Φ (l,m 1) + ½ {m=0} γ 1 l l l l k=0 Φ (l,k) 1 Γ(l l ) 1,k, (3.14) where Γ (l),m = ( 1 l l! (β (z) 1) m σ z l (z) ) z=1 and the nequalty l l s nterpreted componentwse. Therefore, from (3.14) we have the followng proposton. Proposton 2. For l +m 2, Φ (l,m) = ½ {m 1}λ Φ (l,m 1) + ½ {lj 1}λ j Φ (l 1j,m) + ½ {m=0}γ 1 γ +l ν γ +l ν γ +l ν j + ½ {m 1}(l +1)ν Φ (l+1,m 1) + ½ {m=0}γ 1 γ +l ν γ +l ν l l l l l l 1 k=0 Φ (l,k) 1 Γ(l l ) 1,k Φ (l, l l ) 1 Γ (l l ) 1, l l. (3.15) 11

12 We note that (3.15) s a system of lnear equatons for Φ (l,m). Ths system of lnear equatons can be solved by the Gaussan elmnaton method. However, we wll use an teratve method to solve the system of lnear equatons (3.15). In the followng theorem the teratve algorthm s presented and the convergence of teraton s guaranteed. Theorem 3. For,l,m,n wth = 1,...,, l +m = k, n = 0,1,..., defne Φ (l,m) (n) as follows: Φ (l,m) (0) =0, Φ (l,m) (n) = ½ {m 1}λ Φ (l,m 1) + ½ {lj 1}λ j Φ (l 1j,m) γ +l ν γ +l ν j Then we have that + ½ {m=0}γ 1 γ +l ν + ½ {m=0}γ 1 γ +l ν l l l l l l 1 k=0 Φ (l,k) 1 Γ(l l ) 1,k + ½ {m 1}(l +1)ν γ +l ν Φ (l, l l ) 1 (n 1)Γ (l l ) 1, l l, n 1. Φ (l+1,m 1) (n 1) () Φ (l,m) (n) s nondecreasng n n. () lm n Φ (l,m) (n) = Φ (l,m). Proof. By nducton on n, we have that Φ (l,m) (n) s ncreasng n n and Φ (l,m) (n) Φ (l,m) for all n. Thus () s proved. Moreover, lm n Φ (l,m) that {lm n Φ (l,m) (n) : = 1,...,, l + m = k} and {Φ (l,m) (n) exsts and lm n Φ (l,m) (n) Φ (l,m). Suppose : = 1,...,, l + m = k} are dfferent solutons of the system of equatons (3.15). Then {alm n Φ (l,m) (n) + (1 a)φ (l,m) : = 1,...,, l +m = k} s a soluton for any a R. Snce lm n Φ (l,m) that (n) Φ (l,m), there exsts a such a lm n Φ(l,m) for all = 1,..., and (l,m) wth l +m = k and (n)+(1 a)φ (l,m) 0 a lm n Φ(l,m) (n)+(1 a)φ (l,m) = 0 for some = 1,..., and (l,m) wth l +m = k. Hence there exsts a nonnegatve (vector) soluton of (3.15)wthazerocomponent,whchsacontradcton,because ½ {m 1}λ 12 γ +l ν Φ (l,m 1) + j ½ {lj 1}λ j γ +l ν Φ (l 1j,m) +

13 ½ {m=0} γ 1 l l γ +l ν l 1 l k=0 Φ (l,k) ) 1 Γ(l l 1,k s postve for all = 1,, and (l,m) wth l +m = k. Therefore, lm n Φ (l,m) (n) = Φ (l,m) for all = 1,..., and (l,m) wth l +m = k. 3.2 Staton sze analyss at arbtrary tme ponts In the prevous subsecton we have found the generatng functons of the number of customers at the begnnng of glue perods, vst perods, and swtchover perods n terms of φ (z,w). We now represent the generatng functon of the number of customers at arbtrarytme ponts n terms of φ (z,w), as shown below n Theorem 4. Ths wll allow us to obtan the moments of the staton sze dstrbuton at arbtrary tme ponts. Theorem 4. (a) The jont generatng functon, R () s (z), of the number of customers n the orbt at an arbtrary tme pont n a swtchover perod from staton s gven by R s () (z) = γ E[S ] φ 1 σ (z) (z,β (z)) j=1 λ j(1 z j ). (3.16) (b) The jont generatng functon, R () g (z,w), of the number of customers n the queue and n the orbt at an arbtrary tme pont n a glue perod of staton s gven by R () g (z,w) = γ φ (z,w). (3.17) (c) The jont generatng functon, R () v (z,w), of the number of customers n the queue and n the orbt at an arbtrary tme pont n a vst perod of staton s gven by R () v (z,w) = γ ρ E[C] φ (z,w) φ (z,β (z)) w β (z) 1 β (z) j=1 λ j(1 z j ). (3.18) Proof. (a) otce that the number of customers n the orbt at an arbtrary tme pont n a swtchover perod from staton s the sum of two ndependent terms: the number of customers at the begnnng of the swtchover perod and the number of customers who arrved durng the elapsed swtchover perod. The generatng functon of the former s R () s (z) and the generatng functon of the latter s gven by 1 σ (z) E[S ]( j=1 λj(1 zj) ). Thus R () s (z) = 1 σ (z) s (z) E[S ] ( j=1 λ j(1 z j ) ), R () 13

14 from whch and (3.2) we get (3.16). (b) By the theory of Markov regeneratve processes, Ths yelds (3.17). R () g (z,w) = γ 0 φ (z;w;t)dt. (c) otce that the number of customers n the system at an arbtrary tme pont n a vst perod conssts of two parts: the number of customers n the system at the begnnng of the servce of the customer currently n servce and the number of customers who arrved durng the elapsed tme of the current servce. The generatng functon of the former s gven by (see Remark 3 of [1] for a detaled proof) R v () () (z,w) R v (z,β (z)) E[Y (q) ](w β (z)) and the generatng functon of the latter s gven by From (3.19) and (3.20) we have R () v (z,w) = γ E[Y (q) ]E[B ] = γ E[Y (q) ] φ (z,w) φ (z,β (z)), (3.19) w β (z) 1 β (z) ( ). (3.20) E[B ] j=1 λ j(1 z j ) φ (z,w) φ (z,β (z)) w β (z) 1 β (z) j=1 λ j(1 z j ). Snce ρ = (q) E[Y ]E[B ] E[C], (3.18) follows from the above equaton. We ntroduce the followng scaled moments: Ψ (l,m) g, = 1 l!m! Ψ (l,m) v, = 1 l!m! Ψ (l) s, = 1 l! l +m z l w mr() g (z,w) z=1,w=1, l +m z l w mr() v (z,w) z=1,w=1, l z lr() s (z) z=1. These moments satsfy the followng theorem, whch can be derved by usng Equatons (3.16), (3.17) and (3.18). Theorem 5. We have 14

15 () Ψ (l,m) g, = γ Φ (l,m). () Ψ (l,m) v, = γ l l,m+k+1) ρ E[C] l l k=0 Φ(l () Ψ (l) s, = γ ) E[S ] l l Θ(l ζ (l l ) Moreover, Θ (l) s gven by η (l l ),k, where ζ (l) where (l),m = 1 l l! (β z l (z) 1) m z=1. ( ), where η (l),m = 1 l (β(z) 1) m+1 z=1 l! z l. j=1 λj(1 zj) ( ) = 1 l 1 σ (z) z=1 l! z l and Θ (l) j=1 λj(1 zj) = 1 l! Θ (l) = l l Φ (l,k) (l l ),k, l l k=0 l z l φ (z,β (z)) z=1. From now on we obtan the frst and second moments of the staton szes of each type of customers n steady state. Let M o and Υ be the steady state random varables correspondng to M o (t) and Υ(t), respectvely. That s, M o s the number of customers n the orbt of staton n steady state and Υ s the number of glued customers n steady state. Let M oq be the number of customers n the orbt of staton plus the glued customers n the queue of staton n steady state, and M be the number of customers n staton (ncludng the customer n servce at staton ) n steady state. Moreover, we defne the followng ndcator random varables: for = 1,...,, 1 f the server s servng at staton n steady state, I v, = 0 otherwse, 1 f the server s n the glue perod of staton n steady state, I g, = 0 otherwse, 1 f the server s swtchng from staton to staton +1 n steady state, I s, = 0 otherwse. Then we have that for = 1,...,, M o = M o (I v,k +I g,k +I s,k ), k=1 M oq = M o +Υ(I v, +I g, ), M = M oq +I v,. 15

16 Therefore, the mean staton szes, E[M o ],E[Moq ], and E[M ], = 1,...,, are gven by E[M o ] = k=1 ( ρ k Ψ (1,0) v,k + E[G k] E[C] Ψ(1,0) g,k + E[S ) k] E[C] Ψ(1) s,k, (3.21) E[M oq ] = E[M o ]+ρ Ψ (0,1) v, + E[G ] E[C] Ψ(0,1) g,, (3.22) E[M ] = E[M oq ]+ρ. (3.23) ow, n ordertoobtan the secondmomentsofthe statonszes, E[M omo j ],E[Moq M oq j ], ande[m M j ],,j = 1,...,, note that M o Mo j = M o Mo j (I v,k +I g,k +I s,k ), k=1 M oq M oq j = M o Mo j +Mo Υ(I v,j +I g,j )+M o j Υ(I v, +I g, )+Υ 2 (I v, +I g, )½ {=j}, M M j = M oq M oq j +M oq I v,j +M oq j I v, +I v, ½ {=j}. Therefore, the second moments of the staton szes are gven by ) E[M o Mo j ] = k=1(ρ k Ψ (1+1j,0) v,k + E[G k] E[C] Ψ(1+1j,0) g,k + E[S k] E[C] Ψ(1+1j) s,k 2 ( ) k=1 ρ k Ψ (21,0) v,k + E[G k] E[C] Ψ(21,0) g,k + E[S k] E[C] Ψ(21) s,k E[M oq M oq j ] = E[M omo j ]+ρ Ψ (1j,1) v, ( E[(M o)2 ]+2 E[M M j ] = ρ Ψ (1,1) v, +ρ j Ψ (1,1) v,j f j, f = j, + E[G] E[C] Ψ(1j,1) g, + E[Gj] E[C] Ψ(1,1) g,j f j, + E[G] E[C] Ψ(1,1) g, +ρ Ψ (0,2) v, + E[G] E[C] Ψ(0,2) g, E[M oq M oq j ]+ρ Ψ (1j,0) v, +ρ j Ψ (1,0) v,j f j, E[(M oq ) 2 ]+2ρ Ψ (1,0) v, +2ρ Ψ (0,1) v, +ρ f = j. ) f = j, (3.24) (3.25) (3.26) 3.3 A numercal example In ths subsecton we present numercal results for the frst and second moments of the number of customers n each staton. The expresson for the mean number of customers n each staton s gven by (3.23), together wth (3.21) and (3.22). By usng the formulas (3.24)-(3.26), we can obtan an expresson for the varance of the number of customers n each staton and an expresson for the covarance of the numbers of customers n two dfferent statons. ote that these moments are expressed n terms of Φ (l,m), refer to Theorem 5. Therefore, these moments can be obtaned by usng Theorems 2 and 3. In the followng numercal example we consder a sngle server pollng model wth fve statons (.e., = 5). 16

17 E[M ] 4.5 E[M ] =1 =2 =3 =4 = =1 =2 =3 =4 = E[G] E[G] (a) 0 E[G] 10. (b) 0 E[G] Fgure 1: The mean number of customers n staton, E[M ], = 1,...,5, varyng E[G] =1 =2 =3 =4 =5 SCV[M ] 0.6 SCV[M ] =1 =2 =3 =4 = E[G] E[G] (a) 0 E[G] 10. (b) 0 E[G] Fgure 2: The squared coeffcent of varaton for the number of customers n staton, SCV[M ], = 1,...,5, varyng E[G]. Example 1. We assume that the arrval rate of type- customers s λ = 0.025for all, = 1,...,5. The servce tmes of type- customers are exponentally dstrbuted wth means E[B 1 ] = 1,E[B 2 ] = 2,E[B 3 ] = 4,E[B 4 ] = 8 and E[B 5 ] = 16, respectvely. Hence the total utlzaton of the server s ρ = 5 =1 ρ = < 1. The swtchover tmes from staton to staton +1 are determnstc wth E[S ] = 1 for all, = 1,...,5. The retral rate of customers n the orbt of staton s ν = 1 for all, = 1,...,5. The glue perods at staton are exponentally dstrbuted wth parameters γ, = 1,...,5. We assume that γ s the same for all,.e., E[G ] = E[G] for all, = 1,...,5. In Fgure 1 we plot the mean number of customers n staton, E[M ], = 1,...,5, varyng the 17

18 =1,j=2 =1,j=3 =1,j=4 =3,j= =1,j=2 =1,j=3 =1,j=4 =3,j=5 Cor(M,M j ) Cor(M,M j ) E[G] E[G] (a) 0 E[G] 10. (b) 0 E[G] Fgure 3: The correlaton coeffcent of the numbers of customers n staton and staton j, Cor(M,M j ), (,j) = (1,2), (,j) = (1,3), (,j) = (1,4) and (,j) = (3,5), varyng E[G]. mean glue perod E[G]. In Fgure 2 we plot the squared coeffcent of varaton (SCV) for the number of customers n staton, SCV[M ], = 1,...,5, varyng E[G]. In Fgure 3 we plot the correlaton coeffcent of the numbers of customers n two dfferent statons, Cor[M 1,M 2 ], Cor[M 1,M 3 ], Cor[M 1,M 4 ] and Cor[M 3,M 5 ], varyng E[G]. In Fgures 1(a), 2(a) and 3(a) we vary E[G] from 0 to 10 n order to better reveal the behavor of the system for small E[G]. In Fgures 1(b), 2(b) and 3(b) we vary E[G] from 0 to 1000 n order to examne the behavor of the system for large E[G]. We can draw the followng conclusons from these plots: For small glue perod lengths, the chances for a customer to retry are very low, hence the staton sze s large. If the glue perod s very large, the customers face a long delay before gettng served. There exsts an optmal glue length at whch each staton has a mnmum mean staton sze. The fgures suggest that the followng happens when the mean glue perod grows large: () The mean numbers of customers grow lnearly n E[G]. () The squared coeffcent of varaton tends to a lmt when E[G]. 18

19 () The correlaton coeffcents between the numbers of customers n dfferent statons tend to some lmt when E[G]. In [21] the author consders classcal pollng systems wth a branchng-type servce dscplne lke exhaustve or gated servce, and wthout glue perods, for the case that swtchover tmes become large. It s readly seen that our pollng model starts to behave very smlarly as such a pollng model, when the glue perods grow large; ndeed, every type- customer wll now almost surely become glued durng the frst glue perod of staton that t experences durng ts stay n the system, and hence wll be served durng the frst vst perod of staton after ts arrval to the system - just as n an ordnary gated pollng system. However, we cannot mmedately apply the asymptotc results of [21] where swtchover tmes become large, because t consders determnstc swtchover tmes, whle the focus s on the watng tme dstrbuton. In a future paper we ntend to study the asymptotc behavour of pollng systems wth large swtchover tmes, thus also obtanng the asymptotc behavour of pollng systems wth large glue perods. We shall, among others, derve asymptotc expressons for the kth moment of the staton sze. Our prelmnary fndngs are n agreement wth the lmtng behavours of the mean staton sze, the squared coeffcent of varaton of the staton szes, and the correlaton coeffcent of the staton szes. The mean staton sze s asymptotcally lnear n the mean swtchover tmes, and the squared coeffcent of varaton and the correlaton coeffcent of the staton szes converge as the mean swtchover tmes go to nfnty (wth the ratos of the swchover tmes beng constant), as dsplayed n the fgures. 4 The pollng system wth retrals and general glue perods In [1] and n Secton 3 of the current paper we have presented the dstrbuton and mean of the number of customers at dfferent tme epochs for a gated pollng model wth retrals and glue perods, where the glue perods are determnstc and exponentally dstrbuted, respectvely. In ths secton, we assume that glue perods have general dstrbutons. We frst consder the dstrbuton of the total workload n the system and present a workload decomposton. Subsequently we use ths to obtan a pseudo conservaton law,.e., an exact expresson for a weghted sum of the mean watng tmes. In ts turn, the 19

20 pseudo conservaton law s used to obtan an approxmaton for the mean watng tmes of all customer types. We present numercal results that ndcate that the approxmaton s very accurate. Fnally we use ths approxmaton to optmze a weghted sum of the mean watng tmes, =1 c E[W ], where c, = 1,..., are postve constants and E[W ] s the mean watng tme of a type- customer untl the start of ts servce, by choosng the glue perod lengths, gven the total glue perod n a cycle. 4.1 Workload dstrbuton and decomposton DefneV astheamountofworknthesystemnsteadystate. Furthermore,let B(s) = =1 λ (1 B (s)). The LST of the amount of work at an arbtrary tme can be wrtten as E[e sv ] = 1 E[C] ( E[S ]E[e =1 (S) sv ]+E[G ]E[e (G) sv (D) sv ]+ρ E[C]E[e ] ), (4.1) where V (S), V (G) and V (D) are the amount of work n the system durng the swtchovertme from staton, glue perod of staton and vst perod of staton, respectvely. Let V (X), V (Y) and V (Z) be the work n the system at the start of glue perod of staton, vst perod of staton and swtchover perod from staton, respectvely. We know that Therefore, (E[e =1 = = = = sv (S) ( E[e =1 (E[e =1 (E[e =1 =1 (Z) sv E[e E[e (S) sv (G) sv ] = E[e ] = E[e (Z) sv (X) sv ) (G) sv ]E[S ]+E[e ]E[G ] sv (Z) sv (Z) ] 1 S ( B(s)) B(s) ] E[e (Z) sv +E[e (X) sv ] S ( B(s))+E[e B(s) ] 1 S ( B(s)) E[S ] B(s), ] 1 G ( B(s)) E[G ] B(s). ] 1 G ) ( B(s)) B(s) (X) sv ] E[e (X) sv ] E[e +1 ]+E[e sv (X) (Y ) sv ] E[e ] B(s) (Z) (Y ) sv (E[e sv ] E[e ] B(s) ) (X) sv ) ) ] G ( B(s)). (4.2) 20

21 Furthermore, usng the last formula of the proof of Theorem 2 n Boxma et al. [6], but wth our notatons, we have (D) sv ρ E[C]E[e (Y ) ] = E[e sv B(s) s (Z) sv ] E[e ]. (4.3) Substtutng (4.2) and (4.3) n (4.1), we have E[e sv s ] = E[C](s B(s)) ) (Z) (Y ) sv (E[e sv ] E[e ]. (4.4) B(s) =1 Defne the dle tme as the tme the server s not servng customers (.e., the sum of all the swtchover and glue perods). Let V (Idle) be the amount of work n the system at an arbtrary moment n the dle tme. We have, by (4.2) and (2.1), ( E[e sv (Idle) 1 ] = E [ =1 (S +G ) ] E[e =1 1 = E [ =1 (S +G ) ] = 1 (1 ρ)e[c] =1 =1 (S) sv (G) sv ]E[S ]+E[e ]E[G ] (Z) (Y ) sv (E[e sv ] E[e ] B(s) (Z) (Y ) sv (E[e sv ] E[e ] B(s) ) ) ). (4.5) We know that the LST of the amount of work at steady state, V M/G/1, n the standard M/G/1 queue where the arrval rate s =1 λ and the LST of the servce tme dstrbuton s λ =1 B (s), s j=1 λj gven by E[e sv M/G/1 ] = (1 ρ)s. (4.6) s B(s) From Equatons (4.4), (4.5) and (4.6) we have E[e sv ] = E[e sv M/G/1 ]E[e sv (Idle) ]. In Theorem 2.1 of [5], a workload decomposton property has been proved for a large class of sngleserver mult-class queueng systems wth servce nterruptons (lke swtchover perods or breakdowns). It amounts to the statement that, under certan condtons, the steady-state workload s n dstrbuton equal to the sum of two ndependent quanttes: () the steady-state workload n the correspondng queueng model wthout those nterruptons, and () the steady-state workload at an arbtrary nterrupton epoch. The gated pollng model wth glue perods and retrals of the present paper satsfes all the assumptons of 21

22 Theorem 2.1 of [5], and hence, n agreement wth what we have seen above, the workload decomposton ndeed holds. 4.2 Pseudo conservaton law By the workload decomposton, t s shown n [5] that =1 =1 ρ E[W ] = ρ λ E[B 2] 2(1 ρ) +ρ E[( =1 (S +G ) ) 2] 2E [ =1 (S +G ) ] + E [ =1 (S +G ) ] ( ρ 2 2(1 ρ) =1 ρ 2 ) + E[F ], where F s the work left n staton at the end of a vst perod of staton (and hence at the start of a swtchover from staton ). Other than E[F ], Equaton (4.7) s ndependent of the servce dscplne. ote that E[F ] = E[Z () ]E[B ]. To fnd E[Z () ] we wll derve a relaton between E[Z () ] and E[Y (q) ]. E[Y (q) ] conssts of the followng three parts: =1 (4.7) () Mean number of type- customers who were already present at the end of the prevous vst to staton and who are glued durng the glue perod just before the current vst to staton. () Mean number of type- customers who have arrved durng the tme nterval from the end of the prevous vst to staton to the start of the glue perod of staton just before the vst to staton, and who are glued durng that glue perod. () Mean number of type- customers who arrve durng the glue perod of staton just before the vst to staton. ote that () equals (1 G (ν ))E[Z () ] because the mean number of type- customers who were present at the end of the prevous vst to staton s E[Z () ], and the probablty that a customer who was present at the end of the prevous vst to staton s glued durng the glue perod just before the current vst to staton, s 1 G (ν ). () equals (1 G (ν ))λ ( (1 ρ )E[C] E[G ] ). Here λ ( (1 ρ )E[C] E[G ] ) s the mean number of type- customers who have arrved durng the tme nterval from the end of the prevous vst to staton to the start of the glue perod of staton just before the vst to staton. 22

23 Fnally, () equals λ E[G ]. Therefore, E[Y (q) ] =(1 G (ν ))E[Z () ]+(1 G (ν ))(1 ρ )λ E[C]+ G (ν )E[G ]λ. (4.8) Snce ρ = (q) E[Y ]E[B ] E[C], we have E[Y (q) ] = λ E[C]. Hence, by (4.8), we get E[Z () ] = λ ρ E[C]+ λ G (ν ) 1 G (ν ) (E[C] E[G ]). Therefore, E[F ] s gven by E[F ] = ρ 2 E[C]+ ρ G (ν ) 1 G (ν ) (E[C] E[G ]). (4.9) The frst term on the rght-hand sde equals the mean amount of work for type- customers who arrved at statondurngavstperodofstaton. Thesecondtermsnterpretedasfollows: Snceλ (E[C] E[G ]) s the mean number of type- customers who arrve durng one cycle excludng the glue perod of staton n that cycle, ( G (ν)) k ρ (E[C] E[G ]) s the mean amount ofworkfor type- customerswho arrvedurng the kth prevous cycle excludng the glue perod of staton n that cycle, and who are present at the end of the current vst perod of staton. Hence, the second term, whch s k=1 ( G (ν)) k ρ (E[C] E[G ]), s the mean amount of work for type- customers who were present n the orbt of staton at the begnnng of the vst perod of staton. From Equatons (4.7) and (4.9) together wth (2.1), we obtan the followng pseudo conservaton law: ( =1 ρ EW =ρ λ E[B 2] + E[( =1 (S +G ) ) 2] ) ( ) [ E 2(1 ρ) 2E [ =1 (S +G ) ] + ρ 2 + ρ 2 =1 (S +G ) ] 2(1 ρ) =1 ( [ ρ G (ν ) E j=1 + (S j +G j ) ] ) 1 G E[G ]. (4.10) (ν ) 1 ρ =1 =1 4.3 Approxmaton of the mean watng tmes We now use the pseudo conservaton law to fnd an approxmaton for the mean watng tmes of all customer types. Below we brefly sketch the dea behnd the approxmaton. Evertt [8] has developed a method to approxmate the mean watng tmes n an ordnary gated pollng system (wthout retrals and glue perods). The dea n ths approxmaton s that an arrvng customer frst has to wat for the resdual cycle tme, untl the server begns a new vst to ts staton. Subsequently, t has to wat for the 23

24 servce tmes of all customers of the same type, who arrved before t, n the elapsed cycle tme. Ths leads to E[W j ] = (1+ρ j )E[R cj ], where E[R cj ] s the mean of the resdual tme of a cycle startng wth a vst to staton j, whch s the same as the mean of the elapsed tme of a cycle startng wth a vst to staton j. ext, Evertt assumed that for all j, ths mean resdual cycle tme s ndependent of j,.e., E[R cj ] E[R c ], leadng to the approxmaton E[W j ] (1+ρ j )E[R c ] for the model wthout retrals and glue perods. In ths paper, we ntroduce a smlar type of approxmaton, ncludng one extra term, for the mean watng tmes n the model wth retrals and glue perods: E[W j ] (1+ρ j )E[R c ]+ G j (ν j ) 1 G j (ν j ) (E[C] E[G j]). (4.11) In the appendx we provde a detaled dervaton of (4.11). The frst term on the rght-hand sde of (4.11) s the same as the term n [8]. The second term on the rght-hand sde s added because not every customer who arrves n a partcular cycle receves servce n that cycle. The type-j customers arrvng durng any perod other than the glue perod of staton j receve servce n the followng vst perod wth probablty 1 G j (ν j ). Furthermore, the type-j customers arrvng durng any perod other than a glue perod of staton j, have to wat for a geometrc number (wth parameter 1 G j (ν j )) of cycles before recevng servce. Snce a type-j customer arrves durng a perod other than a glue perod of staton j wth probablty E[C] E[Gj] E[C], the mean number of cycles untl an arbtrary type-j customer receves ts servce, s E[C] E[Gj] E[C] G j(ν j) 1 G. The second term s obtaned by multplyng ths mean number of cycles j(ν j) and the mean cycle tme. It should be noted that, n realty, the mean resdual cycle tmes for staton and staton j (j ) are not equal. A key element of our approxmaton s to assume that they are equal. We can now use the pseudo conservaton law to determne the one unknown term E[R c ]: By substtutng (4.11) n (4.10) and usng (2.1), we get ( ρ E[R c ] ρ+ =1 λ E[B 2] + E[( =1 (S +G ) ) 2] [ =1 ρ2 2(1 ρ) 2E [ =1 (S +G ) ] +ρe =1 (S +G ) ] ) 2(1 ρ) ( [ =1 + ρ2 E ρ+ =1 (S +G ) ] ). (4.12) =1 ρ2 2(1 ρ) 24

25 Substtuton of (4.12) nto (4.11) yelds the approxmaton for the mean watng tmes of all customer types: { ( 1+ρ j E[W j ] ρ+ =1 ρ λ E[B 2] =1 ρ2 2(1 ρ) ( [ E + ρ 2 =1 (S +G ) ] )} 2(1 ρ) =1 + E[( =1 (S +G ) ) 2] 2E [ + G j (ν j ) 1 G j (ν j ) [ =1 (S +G ) ] ) =1 (S +G ) ] +ρe 2(1 ρ) ( [ E =1 (S +G ) ] 1 ρ E[G j ] ). (4.13) We now consder varous examples to compare the above approxmaton results wth the exact analyss from[1] for determnstc glue perods and from Secton 3 of the present paper for exponentally dstrbuted glue perods. Further, we wll compare the results of ths approxmaton wth smulaton results for the case that the glue perods follow a gamma dstrbuton. Determnstc glue perods In the numercal example of Table 1 we consder a two-staton pollng system. The swtchover tmes and servce tmes are exponentally dstrbuted. We keep the parameters of staton 1 fxed, λ 1 = 1,E[B 1 ] = 0.45,E[S 1 ] = 1,G 1 = 0.5,ν 1 = 1, and vary the parameters of staton 2. Exponental glue perods In the numercal example of Table 2 we consder a three-staton pollng system. The swtchover tmes are determnstc, and the servce tmes are exponentally dstrbuted. We keep the parameters of staton 1 fxed, λ 1 = 1,E[B 1 ] = 0.45,E[G 1 ] = 0.5. Further, the swtchover tmes and exponental retral rates of all three statons are fxed, S 1 = S 2 = S 3 = 1 and ν 1 = ν 2 = ν 3 = 1. Gamma dstrbuted glue perods In the above two examples we can get the exact mean watng tmes usng the method n [1] and Secton 3 of ths paper. In the numercal examples of Table 3 we compare the approxmate mean watng tmes wth smulaton results, for a pollng system where the lengths of glue perods are gamma dstrbuted. We consder a fve-staton pollng system n whch the glue perods, swtchover tmes and servce tmes are all gamma dstrbuted. We smulate such a system to fnd the mean watng tmes. We also gve a 25

26 λ 2 E[B 2] E[S 2] G 2 ν 2 Exact (E[W 1],E[W 2]) Approx (E[W 1],E[W 2]) (71.61, 71.61) (71.61, 71.61) (21.44, 20.34) (21.49, 20.24) (15.18, 13.96) (15.21, 13.83) (20.52, 18.82) (20.55, 18.71) (23.01, 11.48) (22.99, 11.67) (22.97, 20.31) (22.99, 20.20) Table 1: Comparson of exact and approxmate mean watng tmes for a pollng system wth determnstc glue perods. λ 2 E[B 2] E[G 2] λ 3 E[B 3] E[G 3] Exact (E[W 1],E[W 2],E[W 3]) Approx (E[W 1],E[W 2],E[W 3]) (121.0, 121.0, 121.0) (121.0, 121.0, 121.0) (47.59, 47.58, 46.74) (47.71, 47.71, 46.24) (33.65, 33.64, 32.54) (33.69, 33.69, 31.97) (246.8, 246.6, 242.3) (242.4, 257.1, 230.2) (33.52, 33.51, 32.42) (33.56, 33.56, 31.86) (44.88, 19.71, 43.64) (45.22, 19.50, 42.92) (48.66, 21.42, 28.75) (49.03, 21.17, 27.98) Table 2: Comparson of exact and approxmate mean watng tmes for a pollng system wth exponentally dstrbuted glue perods. 95% confdence nterval for the mean watng tmes obtaned usng smulatons. We have generated one mllon cycles, splttng ths nto ten perods of 10 5 cycles, and usng the results of these ten perods to obtan confdence ntervals. Then we compare the smulaton results wth the results obtaned usng the approxmaton formula. Here, k and θ are, respectvely, the shape and the scale parameters of the gamma dstrbuton wth probablty densty functon 1 Γ(k)θ k x k 1 e x θ. The values of the parameters are lsted n Table 3(a). Table 3(b) shows the mean watng tmes by smulaton, along wth 95% lower and upper confdence bounds. Table 3(c) shows the approxmate mean watng tmes. We can draw the followng conclusons about the mean watng tme approxmaton. The mean watng tme approxmaton s very accurate. In only two cases (the fourth case n Table 2 and the case () n Table 3) the error s n the order of 5%; n all other cases, we fnd errors 26

27 Parameters () () () (v) (λ 1, λ 2, λ 3, λ 4, λ 5) (0.1, 0.1, 0.1, 0.1, 0.1) (0.1, 0.1, 0.1, 0.1, 0.2) (0.1, 0.1, 0.1, 0.3, 0.2) (0.1, 0.2, 0.1, 0.2, 0.2) (ν 1, ν 2, ν 3, ν 4, ν 5) (2.0, 2.0, 2.0, 2.0, 2.0) (2.0, 2.0, 2.0, 2.0, 3.0) (2.0, 2.0, 2.0, 1.0, 3.0) (2.0, 5.0, 2.0, 4.0, 3.0) (k B1, k B2, k B3, k B4, k B5 ) (1.0, 1.0, 1.0, 1.0, 1.0) (1.0, 1.0, 1.0, 1.0, 0.8) (1.0, 1.0, 1.0, 0.5, 0.8) (1.0, 0.5, 1.5, 0.5, 0.8) (θ B1, θ B2, θ B3, θ B4, θ B5 ) (1.5, 1.5, 1.5, 1.5, 1.5) (1.5, 1.5, 1.5, 1.5, 1.5) (1.5, 1.5, 1.5, 1.5, 1.5) (1.5, 1.5, 1.5, 1.5, 1.5) (k S1, k S2, k S3, k S4, k S5 ) (2.0, 2.0, 2.0, 2.0, 2.0) (2.0, 2.0, 2.0, 2.0, 1.0) (2.0, 2.0, 2.0, 3.0, 1.0) (2.0, 5.0, 2.0, 3.0, 1.0) (θ S1, θ S2, θ S3, θ S4, θ S5 ) (1.0, 1.0, 1.0, 1.0, 1.0) (1.0, 1.0, 1.0, 1.0, 1.0) (1.0, 1.0, 1.0, 1.0, 1.0) (1.0, 1.0, 1.0, 1.0, 1.0) (k G1, k G2, k G3, k G4, k G5 ) (1.0, 1.0, 1.0, 1.0, 1.0) (1.0, 1.0, 1.0, 1.0, 2.0) (1.0, 1.0, 1.0, 0.5, 2.0) (1.0, 3.0, 1.0, 0.5, 2.0) (θ G1, θ G2, θ G3, θ G4, θ G5 ) (1.0, 1.0, 1.0, 1.0, 1.0) (1.0, 1.0, 1.0, 1.0, 1.0) (1.0, 1.0, 1.0, 1.0, 1.0) (1.0, 1.0, 1.0, 1.0, 1.0) (a) The parameter values for the four cases. Cases (E[W 1],E[W 2],E[W 3],E[W 4],E[W 5]) 95% lower confdence bound 95% upper confdence bound () (68.94, 68.92, 68.87, 68.91, 68.84) (68.65, 68.64, 68.60, 68.58, 68.54) (69.23, 69.21, 69.14, 69.24, 69.15) () (108.59, , , , 72.43) (107.86, , , , 71.99) (109.33, , , , 72.87) () (217.54, , , , ) (216.04, , , , ) (219.04, , , , ) (v) (276.39, , , , ) (275.30, , , , ) (277.47, , , , ) (b) Smulaton results. Cases Approx (E[W 1],E[W 2],E[W 3],E[W 4],E[W 5]) () (69.00, 69.00, 69.00, 69.00, 69.00) () (108.25, , , , 72.84) () (210.22, , , , ) (v) (274.10, , , , ) (c) Approxmaton results. Table 3: Comparson of mean watng tmes from smulatons and approxmate mean watng tmes for a pollng system wth gamma dstrbuted glue perods. whch typcally are less than 2%. The mean watng tme approxmaton at one staton s ndependent of the change n retral rates of other statons, whch s not true n realty. The mean watng tme approxmatons for two totally symmetrc statons are the same, ndependent of ther order n the system; but ths s also not qute true n realty. 27

28 4.4 Optmal choce of the glue perod dstrbutons In ths subsecton we dscuss an optmzaton problem for the choce of the dstrbutons of the glue perods, G, = 1,...,, to mnmze the weghted sum of the mean watng tmes =1 c E[W ], subject to the constrant =1 E[G ] = L, where c, = 1,...,, and L are postve constants. Because we do not have an explct formula for the mean watng tme, t s dffcult to solve exactly the constraned mnmzaton problem. Instead of fndng the exact soluton of the constraned mnmzaton problem, we wll fnd the optmal choce of the dstrbutons of G, = 1,...,, to mnmze subject to c U =1 E[G ] = L, =1 where U s the approxmaton of E[W ] gven by the rght-hand sde of (4.13). ote that under the constrant =1 E[G ] = L, the objectve functon of the mnmzaton problem becomes [ ( c (1+ρ ) c U = =1 =1 ρ+ j=1 ρ λ je[bj 2] + E[( j=1 S j) 2 ]+2L j=1 E[S j]+e[( j=1 G j) 2 ] j=1 ρ2 2(1 ρ) j 2E[ j=1 S j +L] ( [ E j=1 S j +L ] ) + E[ j=1 S j +L] 2(1 ρ) ( ρ 2 + j=1 ρ 2 j) ] + =1 c E[e νg ] 1 E[e νg ] 1 ρ E[G ] ). (4.14) By Jensen s nequalty, t can be shown that f the nondetermnstc glue perod dstrbutons wth means g, are changed to the degenerate (determnstc) ones wth the same means g, = 1,...,, then the rght-hand sde of (4.14) becomes strctly smaller. Therefore, the above optmzaton problem becomes as follows: mnmze U(g 1,...,g ) subject to g > 0, = 1,...,, (4.15) g = L, (4.16) =1 28

29 where U(g 1,...,g ) [ ( c (1+ρ ) = =1 ρ+ j=1 ρ λ je[bj 2] j=1 ρ2 2(1 ρ) j j=1 + E[S j]+l( ρ 2 + ρ 2 2(1 ρ) j) ] + j=1 j=1 E[S ) j]+l 2 ] 2E[ j=1 S j +L] 1 ) ( j=1 c ( 1+ E[S j]+l g ). (4.17) 1 e νg 1 ρ + E[( j=1 S j) 2 ]+2L =1 Snce U(g 1,...,g ) s contnuous on D {(g 1,...,g ) : g 1 > 0,...,g > 0,g g = L} and U(g 1,...,g ) as mn{g 1,...,g } 0+, U(g 1,...,g ) takes a mnmum at a pont n D. At a mnmum pont (g 1,...,g ), there exsts a Lagrange multpler κ satsfyng f (g ) = κ, = 1,...,, (4.18) where f (g ) c ( c 1 e c ν e νg j=1 E[S j]+l νg (1 e νg ) 2 g ), = 1,...,. 1 ρ For each = 1,...,, the functon f : (0,L) (,f (L)) s bjectve, contnuous and strctly ncreasng. Therefore, t has the nverse functon h : (,f (L)) (0,L), whch s also contnuous and strctly ncreasng. Therefore, Equaton (4.18) and the constrants (4.15) and (4.16) can be wrtten as h j (κ) = L, < κ < mn{f 1 (L),...,f (L)}, (4.19) j=1 g = h (κ), = 1,...,. (4.20) Snce lm κ j=1 h j(κ) = 0, lm κ (mn{f1(l),...,f (L)}) j=1 h j(κ) > L and j=1 h j(κ) s strctly ncreasng n κ, (4.19) has a unque soluton, say κ. Therefore, from (4.20), the optmal soluton (g1,...,g ) s gven by g = h (κ ), = 1,...,. We wll now consder a few numercal examples to look at the dependency of dfferent system characterstcs and the respectve optmal glue perods. In [2] a smlar system was studed wth a focus on optcal swtches, where the revenue of the system depended on dstrbutng glue perods optmally to each 29