Age of Information: The Gamma Awakening

Size: px
Start display at page:

Download "Age of Information: The Gamma Awakening"

Transcription

1 Age of Information: The Gamma Awakening Eie Najm and Rajai Nasser LTHI, EPFL, Lausanne, Switzerand Emai: {eie.najm, arxiv: v [cs.it] 5 Apr 06 Abstract Status update systems is an emerging fied of study that is gaining interest in the information theory community. We consider a scenario where a monitor is interested in being up to date with respect to the status of some system which is not directy accessibe to this monitor. However, we assume a source node has access to the status and can send status updates as packets to the monitor through a communication system. We aso assume that the status updates are generated randomy as a Poisson process. The source node can manage the packet transmission to minimize the age of information at the destination node, which is defined as the time eapsed since the ast successfuy transmitted update was generated at the source. We use queuing theory to mode the source-destination ink and we assume that the time to successfuy transmit a packet is a gamma distributed service time. We consider two packet management schemes: LCFS with preemption and LCFS without preemption. We compute and anayze the average age and the average peak age of information under these assumptions. Moreover, we extend these resuts to the case where the service time is deterministic. I. INTRODUCTION In status update systems, one or severa sources send information updates to one or severa monitors at a certain effective rate. Naturay, the goa of this process is to ensure that the status updates are as timey as possibe at the receiver side. For this purpose, [] uses the term age, to refer to the time eapsed since the generation at instant ut of the newest packet avaiabe at the receiver. Formay, the age of such packet is t t ut and the timeiness requirement at the monitor corresponds to a sma average age. Indeed, rea-time status updating can be modeed as a source feeding packets at rate to a queue which deivers them to the monitor with some deay. Hence, the requirement at the destination transates into finding the optima transmission scheme and/or the optima effective update rate at the source that minimizes im τ τ τ 0 tdt. However, numerous factors affect the evauation of such as the mode of the source update process, the number of sources, the mode of the queue, the number of queues avaiabe, etc. Kau et a. in [] sove one aspect of the probem where they consider a singe source generating packets as a rate Poisson process feeding them to a singe First Come First Served FCFS queue with exponentia service time. Using Kenda notation, this is an FCFS M/M/ system. Moreover, the authors aso consider the cases of deterministic source and exponentia service time, i.e., FCFS D/M/ system, as we as random source and deterministic service time, i.e., FCFS M/D/ system. Yates and Kau in [] generaize the probem soved in [] by considering the presence of mutipe sources sending updates through the same FCFS queue to the same monitor. Aong the same generaization direction as in [], one may ask: what woud happen if we increase the number of queues avaiabe, i.e., if the source is abe to serve mutipe updates at the same time? This question is tacked in [3], where a singe Poisson process is sending updates over an infinite number of queues with exponentia service time. However, in these aforementioned works, the authors mosty consider FCFS queues. One woud focus on Last Come First Served LCFS type of queues since they are intuitivey more suitabe for the probem in hand: we are interested in deivering the newest update to the monitor, which means we gain more by sending the youngest packet in the queue first. This idea is deveoped in [4] where the authors derive an expression for by treating the foowing two modes whie assuming exponentia interarriva and service time: i LCFS queue without preemption; if the queue is busy, any new update wi have to wait in a buffer of size. This means that the new update wi repace any oder packet aready waiting to be served. ii LCFS with preemption, where contrary to the first case, any new update wi prompt the source to drop the packet being served and start transmitting the newcomer. In [4], it is shown that an LCFS queue with preemption achieves a ower average age compared to the mode without preemption. However, both modes outperform the FCFS mode presented in []. In this paper, we aso consider these ast two schemes in order to derive cosed form expressions for. However, the main novety is the assumption of a gamma distribution for the service time in age of information probems. The motivation for such a distribution is twofod: Based on the cassica appications of gamma distributions in queuing theory, these distributions can be seen as a reasonabe approximation if we want to mode reay networks. Indeed, in such network, a transmitter and a receiver are separated byk reays with each reay taking an exponentia amount of time to compete transmission to the next hop. This means that the tota transmission time is the sum of k independent exponentia random variabes which induces a gamma distribution. As we wi see ater, a deterministic random variabe can be seen as the imit of a sequence of gamma distributed random variabes. Therefore, one can study the performance of the LCFS-based schemes under deterministic service time by taking the imit of the resut obtained for a gamma distributed service time. Athough this is an indirect method of cacuating, it is simper than the direct approach. This paper is organized as foows: in Section II we present

2 the preiminary resuts that wi be used throughout our work and define the average peak age as an aternative metric. In Section III we derive the cosed form expressions for both the average age and the average peak age when assuming an LCFS scheme with preemption. On the other hand, Section IV computes the formuas for these quantities when considering an LCFS queue without preemption. In these ast two sections the service time is assumed to be gamma distributed. However, in Section V we cacuate the two ages for a deterministic service time for each of the two schemes. Finay, Section VI presents numerica simuations that vaidate our theoretica resuts. II. A. Genera definitions PRELIMINARIES As we have seen, our two schemes of interest are LCFS with preemption and LCFS without preemption. The variation of the instantaneous age for these two scenarios is given in Figure. The saw-tooth pattern depicted in those figures is due to the foowing behavior of the age. Let t i be the time the i th packet is generated and et t i be the time the ith packet is received if it is successfuy received. Moreover, without oss of generaity, we assume the beginning of observation is at time t 0 and that the queue is empty at this instant with an initia age of 0. The age t increases ineary with time and is set to a smaer vaue when a packet is received. Hence, the instantaneous age is equa to the current time minus the generation time of the newest of the received packets. It is important to note that in both schemes of interest, some packets might be dropped. Hence we ca the packets that are not dropped, and thus deivered to the receiver, as successfuy received packets or successfu packets. In addition to that, we aso define:ii i to be the true index of thei th successfuy received packet, ii Y i t I i+ t I i to be the interdeparture time between two consecutive successfuy received packets, iii X i t Ii+ t Ii to be the interarriva time between the successfuy transmitted packet and the next generated one which may or may not be successfuy transmitted, so f X x e x, iv T i to be the system time, or the time spent by the i th successfu packet in the queue and v N τ max{n;t In τ}, the number of successfuy received packets in the interva [0,τ]. B. Computing the Average Age Using these quantities and Figures a and b, the authors in [4] show that im τ τ τ 0 tdt e EQ i, where e is the effective update rate and EQ i is the expected vaue of the area Q i at steady state. Hence, we need to determine these two quantities. Computing the effective rate: As stated in [5], e P{packet is received successfuy} 3 where P{packet is received successfuy} is the probabiity that a packet in the queue wi be deivered to the receiver. Computing EQ i : Based on Figures a and b, it was shown in [4] that Y EQ i ET i Y i +E C. Computing the Average Peak Age i. 4 Another metric of interest is the average peak age. We define the peak age as P i im t, t t I i t<t I i which is the vaue of the instantaneous age just before it is reduced by the reception of the i th successfu packet. From Figures a and b, we can deduce that the peak age can be written as P i T i +Y i. Therefore, the average peak age is given by: D. Defining the service time EP i ET i +EY i. 5 A the above resuts were obtained without any assumption on the service time. However, as we have discussed before, this paper studies two modes for the service time: a gamma distributed service time with parameters k, and a deterministic service time. Here is a brief description of the gamma distribution. Definition. A random variabe S with gamma distribution Γk, has the foowing probabiity density function: f S s sk e s k Γk. The Erang distribution Ek, is a specia case of the gamma distribution where k N. Such random variabe has a mean of ES k and a variance VarS k. These quantities wi come in handy ater on. Another important property of gamma random variabes is given by the foowing emma: Lemma. Suppose S n Γk n, n is a sequence of random variabes such that ES n µ, for some µ > 0. Then the sequence S n converges in distribution to a deterministic variabe Z as k becomes very arge, i.e, where Z µ S n d Z, as k, with probabiity. The above emma obviousy sti hods if S n Ek n, n. This emma provides an additiona motivation for studying the average age and the average peak age under the assumption of a gamma distributed service time since we can easiy extend the resuts to the deterministic service time mode by etting k.

3 t Δt Age Age 0 Q Q Q 3 Q N Δ 0 Q Q4 QNτ Q Q 3 t t t 3 t 4 t 5 t 6 t 7 t' t' 5 t' t 8 7 t' t N 8 t' N T Y T Y a Age of information for LCFS with preemption scheme t N - t t t Nτ- t t 3 t 4 t 5 t t' t' 6 t' 4 t' t 6 Nτ t' Nτ T Y T Y W 3 S 3 b Age of information for LCFS without preemption scheme t Fig.. Variation of the instantaneous age for both schemes III. AGE OF INFORMATION FOR LCFS WITH PREEMPTION In this section we wi compute the average age and the average peak age EP k for the Last Come First Served LCFS scheme with preemption and a gamma distributed service time. As we have seen before, in this scenario any packet being served is preempted if a new packet arrives and the new packet is served instead. Hence, the number of packets in the queue can be modeed as a continuous-time two-state semi-markov chain depicted in Figure. The 0-state corresponds to the state where the queue is empty and no packet is being served whie the -state corresponds to the state where the queue is fu and is serving one packet. However, given that the interarriva time between packets is exponentiay distributed with rate then one spends an exponentia amount of time X in the 0-state before jumping with probabiity to the other state. Once in the -state, two independent cocks are started: the gamma distributed service time cock of the packet being served and the rate memoryess cock of the interarriva time between the current packet and the next one to be generated. We jump back to the 0-state if the service time cock happens to tick before that of the interarriva time. Given that the interarriva times between packets are i.i.d as we as the service time of each packet, then the probabiity to jump from the -state to the 0-state does not depend on the index of the current packet. Hence, the jump from the -state to the 0-state occurs with probabiity p PS < X, where S is a generic gamma distributed service time and X is a generic rate memoryess interarriva time which is independent of S. On the other hand, if the interarriva time cock happens to tick before the service time cock then the current packet being served is preempted and the new generated packet takes its pace in the queue. Therefore, we stay in the -state and the two cocks are started anew independenty from before. This expains the p probabiity seen in Figure for staying in the -state. Given that the probabiity p wi be usefu in the computation of the average age as we as the average peak age, we start by deriving its expression here: p PS < X k. 6 + Now we are ready to derive the two age metrics. A. Average age We start by deriving the expression for the average age. We need to compute two quantities for this purpose: EQ i and the effective rate e. Computing EQ i : Using 4, we obtain: Y EQ i ET i Y i +E i Y ET i EY i +E i. 7 The second equaity comes from the fact that T i and Y i are independent since the interarriva time is exponentia and hence memoryess. In fact, thei th successfu packet eaves the queue empty and hence Y i ˆX i +Z i where ˆX i X i T i is the remaining of the interarriva time between the departure of the i th successfu packet and the arriva of the next generated one and Z i is the time for a new packet to be successfuy deivered. Z i does not overap with T i and thus is independent from it. As for ˆX i, we aso get that it is independent of T i. To prove this we notice that for a successfuy received packet i the joint distribution f Xi,T i x,t can be written as { 0 if x < t f Xi,T i x,t f X,Sx,t PS<X if x > t, 8 where X and S are the generic independent interarriva time and service time respectivey. Now, using a change of variabe we get f ˆXi,T i ˆx,t f Xi T i,t i ˆx,t f Xi,T i ˆx+t,t { 0 if x < 0 f X,Sˆx+t,t PS<X if x > 0 { 0 if x < 0 hˆxgt if x > shows that ˆXi and T i are indeed independent. Moreover, one can show that ˆXi is exponentia with rate. Given that ˆX i and Z i are both independent from T i, then Y i and T i are aso independent. From now on we wi drop the subscript index since at steady state T i and T i have same the distribution, which is aso the case for Y i and Y i. The foowing emma wi be used to evauate 7:

4 Fig.. 0, p -p Semi-Markov chain representing the queue for LCFS with preemption Lemma. Let G be gamma distributed with parameters k, and F be a rate exponentia random variabe independent of G. Then, conditioned on the event {G < F}, the distribution of G becomes gamma with parameters. k, + f G/G<F t tk e t+ kγk. 0 + Proof: In order to prove this Lemma we wi compute the probabiity density function f G G<F : Pt G < t+ǫ G < F f G G<F t im ǫ 0 ǫ a Pt G < t+ǫpg < F t < G < t+ǫ im ǫ 0 ǫpg < F PF > t f G t PG < F b f G t e t p tk e e t t k Γk + tk e t+ kγk + k where a is obtained by appying Bayes rue and in b, p is given by 6. In order to appy Lemma, we first notice that for a given packet i, the event {S i < X i } is equivaent to the event {packet i was successfuy received}. Hence the probabiity P PS i < α S i < X i is the probabiity that the service time of the i th packet is ess than α given that this packet was successfuy transmitted. However, since the service times and interarriva times are i.i.d then P does not depend on the index i. Now since T is the service time of a successfu packet then this eads us to PT < α PS i < α S i < X i PS < α S < X, where S and X are the generic service and interarriva time respectivey. By repacing G by S and F by X in Lemma, we deduce that the system time T is gamma distributed with parameters. Therefore, k, + ET k +. Now we turn our attention to the distribution ofy for which we compute its moment generating function. Before going further in our anaysis, we state the foowing emma. Lemma 3. Let G be gamma distributed with parameters k, and F be a rate exponentia random variabe independent of G. If F is a random variabe such that φ F s p where p PF < α PF < α F < G, then the moment generating function of F is given by s s+ s k + k., 3 Proof: We first start by computing the probabiity density function of X. Pt F < t+ǫ F < G f X t im ǫ 0 ǫ Pt F < t+ǫpf < G t F < t+ǫ im ǫ 0 ǫpf < G e t PG > t, p where p + k. So now we can cacuate the moment generating function of F. φ F s 0 0 p f F te st dt p e t PG > te st dt s 0 PG < te st dt d Using integration by parts and the fact that dtpg < t f G t tk e t k Γk, we get φ F s p s s+ s k Lemma 4. The moment generating function of Y is given by: φ Y s k. 4 s+ s Proof: By observing Figure we notice that Y is the smaest time needed to go from the 0-state back to the 0- state. Hence Y can be written as Y X+W where X is the generic interarriva time and W is the time spent in the -state before the first jump back to the 0-state. So W can be written as: S with probabiity p X +S with probabiity pp W X +X +S with probabiity p p. M X j +S, 5 j0

5 where X j is such that PX j < α PX < α X < S, S is such that PS < α PS < α S < X and M is a geometricp random variabe which is independent of X j and S, and which gives the number of discarded packets before the first successfu reception. Appying Lemmas and 3 on S and X respectivey and using the fact that M,S and X j are a mutuay independent, it foows that φ W s E e s M j0 X j φ S s k E φ X s M + + s + φ X s j p p j + s j0 k s k. 6 s+ s Moreover, since X and W are independent and φ X s s, we get using 6 φ Y s φ X sφ W s s+ s k. Now that we have found φ Y we can compute the first two moments of Y as EY +k and EY + k + k+ k. Combining these resuts with, we obtain, EQ i +k. 7 Computing the effective rate: Using 3 we get k e p. 8 + Now we are ready to compute the average age: We concude Proposition. The average age in the LCFS with preemption scheme assuming Γk, service time is given by: e EQ i +k. 9 Proof: Using 7 and 8. B. Average peak age Proposition. The average peak age in the LCFS with preemption scheme assuming Γk, service time is given by: EP i ET+EY k + + +k. 0 Proof: Using 5, and the vaue of EY. Fig. 3. A k,/ A,/ A 0,/ A 3 3,/... A k-,/ k IV. A ',/ A,/ A 3,/ ' 3'... A k-,/ k' A k,/ eve 0 eve Markov chain representing the queue for LCFS without preemption AGE OF INFORMATION FOR LCFS WITHOUT PREEMPTION Another interesting scheme worth to study is the LCFS without preemption. In this scenario, we assume that the queue has a buffer of size and we wait for the packet being served to finish before serving a new one. If whie serving a packet a new update arrives, it repaces any packet waiting in the buffer. In this section we wi derive a cosed form expression for the average age and the average peak age EP k for LCFS without preemption whie assuming an Erang distribution for the service time with parameter k,. An Erang distribution is nothing but a specia case of the gamma distribution where k N. Moreover, an Erang distribution k, can be seen as the sum of k independent memoryess random variabes A j, each with rate. Using this observation, we mode the state of the queue as a two-eve Markov chain as shown in Figure 3. As in the previous section, we wi denote the generic rate interarriva time by X and the generic Erang distributed service time by S k j A j. Using this notation, we notice that the service time can be represented as the succession of k exponentia-time steps that need to be accompished for a successfu reception. Hence, a packet in state j {,...,k} or j {,...,k } is a packet competing his j th step out of a tota of k. Moreover, the 0-state represents an empty queue, a the states of eve 0 represent an empty buffer and those of eve represent a fu buffer. After spending an exponentia amount of time in the 0-state, we can ony jump to the -state once a new update arrives. Using the memoryess property of the exponentia distribution, we can describe the evoution of this packet in the queue as foows: at state j {,..., k}, two exponentia cocks start simutaneousy. One cock denoted A j of rate and another one denoted Λ j of rate. If cock A j ticks first then the packet jump to state j+ and the buffer stays empty. Otherwise it jumps to state j since now the buffer is fu. On the other hand, if the packet is at state j and the A j cock ticks first then the packet jump to state j + without updating the buffer. However, if the Λ j ticks first then the packet stays in state j but we update the buffer with the new arriva.

6 A. Average age Proposition 3. The average age in the LCFS without preemption scheme assuming Erang Ek, service time is: k+ +3k k q k + +k +k+ k k+k +k + + +k+ k+ + +k ++ k +k+ k Proof: As in the previous section we need to compute the effective rate given by 3 and EQ i given by 4. Computing EQ i : Foowing the same ine of thoughts as in Section III, we wi cacuate ET i Y i by expressing it as the average of two conditionay independent variabes given some set { of events. For this end we define the famiy of events Ψ i j A i j > Λi j ; k j+ Ai }, < X where j k. Hence Ψ i j is the event that during the service time of the ith successfu packet a new update arrived at the j th step of the service time i.e, state j or j and then no new update arrived for the remainder of the service time. The superscript i is used to indicate that we are deaing with the i th successfu packet. Forj 0,Ψ i 0 is the event that theith successfu packet eaves the queue empty. Note that for everyi, {Ψ i j, j k} is a partition of the probabiity space. It is sufficient to condition on the event Ψ i 0 in order to ensure conditiona independence between T i and Y i. This is due to the foowing fact: given Ψ i 0, we know that the i th successfu packet eft the queue empty and hence we have a situation identica to that of the with preemption case see Section III and T i and Y i are independent. On the other hand, given Ψ i 0, the buffer is not empty and thus a new packet wi be served directy after the departure of the i th successfu packet. In this case, the interdeparture time Y i is simpy the service time of the i th successfu packet whose vaue is independent of T i W i + S i, where W i the waiting time and S i is the service time of the i th successfu packet see Figure b. Athough conditioning on Ψ i 0 is enough to obtain independence between T i and Y i, we wi need to condition on the two independent events Ψ i j and Ψ i in order to be abe to cacuate the conditiona expectation of T i. However, it is cear that conditioning on these two events aso eads to the independence between T i and Y i. Hence we get ET i Y i E Ti Ψ i j Ψ i E Yi Ψ i j Ψ i j,0 PΨ i j PΨ i. We start by computing E T i Ψ i E W i Ψ i j Ψ i +E Si Ψ i j Ψ i. j Ψ i The waiting time of the i th successfu packet doesn t depend on Ψ i j since they are disjoint in time, but it does depend on Ψ i. In fact, given Ψ i 0, the i th successfu packet wi not wait and start service upon arriva since the i th successfu packet eft the queue empty. However, given Ψ i with 0, the i th successfu packet arrived when the i th successfu packet was at state or of its service time. In order to find the distribution of W i conditioned { on Ψ i we introduce the foowing event: Ψ i,n n g Λi,g < Ai, n+ g Λi,g > } k m Ai m, where {Λ i,g } g is the sequence of interarriva times after the i th successfu packet enters state. Notice that Ψ i,n is the event that exacty n updates arrived when the i th successfu packet was in state or and then no more updates were generated for the remainder of the service time. Hence Ψ i nψ i,n. So conditioned on Ψ i,n we have W i m A i A i m n g n g Λ i,g Λ i,g + m+ A i m 3 It can be shown that, conditioned on { n g Λi,g < Ai }, A i n g Λi,g has an exponentia distribution with rate. This means that under this condition aone, W i has the same distribution as the sum of k + independent exponentia random variabes with rate {. If we further condition on n+ g Λi,g > } k m Ai m and use Lemma, we deduce that conditioned on Ψ i,n, W i has a gamma distribution with parameters k +, + we concude that if we condition on Ψ i as Γ k +, +. Therefore,. Now since Ψ i nψ i,n,, W i is distributed E { W i Ψ i j Ψ i 0 if 0 k + + if 0. 4 Now we turn our attention to E S i Ψ i j Ψ i. One first notices that the service time S i of the i th successfu packet is independent of its arriva time given by the event Ψ i since we assumed independence between service time and interarriva time. Hence, E S i Ψ i j Ψ i E S i Ψ i j. For the case j 0, we get E S i Ψ i 0 E ma i m k + m A i m < X 5 where the ast equaity is obtained by appying Lemma with

7 G k m Ai m and F X. As for the case j 0, we get E S i Ψ i j E A i m Ai j > Λ i j, A i m < X j m m mj+ EA i m +EAi j A i j > Λ i +E mj+ A i m mj+ A i j m < X a j + + k j k +j 6 + where the third term in a is obtained by appying Lemma with G k mj+ Ai m and F X. Therefore, combining 5 and 6 we get, k + if 0,j 0 E T i Ψ i j Ψ i k++j + if 0,j > 0 k +. + if > 0,j 0 k ++j + if > 0,j > 0 7 Now we need to compute E Y i Ψ i j Ψ i. For this end, observe that Y i is independent of Ψ i given that they don t overap in time. Moreover, for j 0, the i th successfu packet eaves the queue empty and thus we wi need to wait an exponentia amount of timex of ratebefore the i th successfu packet arrives and is served directy. Hence, conditioned on Ψ i 0, Y i has same distribution as X +S with X and S independent. On the other hand, for j 0, the i th successfu packet eaves the queue with another packet waiting in the buffer ready to be served. Thus in this case, Y i is simpy the service time of the i th successfu packet. To sum up, E { Y i Ψ i j Ψ i +k if j 0 8 k if j > 0 To compute ET i Y i we sti need the probabiity PΨ i j. For j > 0, we use the fact that Ψ i j is the intersection of two independent events and find that PΨ i j. As for j 0, we have aready seen in Section III + k j+ k. that PΨ i 0 p + These probabiities are independent of the index i and thus we can find PΨ i by repacing j by in the previous expressions. Combining this resuts with 7, 8 we obtain after some tedious cacuations ET i Y i k kk + +k+qk k+k +k +q k+ qk k qk+ 9 with q +. The ast term to compute in order to obtain EQ i is EYi EY i Ψi 0 PΨ i 0 +EYi Ψi 0 PΨ i 0. Based on our previous observations we know that EYi Ψi 0 EX +S and EYi Ψi 0 ES. Using these facts we get +k EYi k +k +q k. 30 Combining 9 and 30, we finay get EQ i k+ +3k +q k+ k+k +k k +q k qk k qk+. 3 Computing the effective rate: To cacuate the effective rate we first observe that the event {packet is successfuy received} is equivaent to the event {packet passes by the - state}. Hence if we uniformize the Markov chain so that the time spent at each state is exponentia with rate +, we get e + π where π is the steady-state probabiity of the -state in the uniformized Markov chain. The anaysis of q q such chain [6], chapter 5 givesπ q k+ +k q. Therefore, e + k +k+ k. 3 Finay, repacing EQ i and e in e EQ i by their expressions in 3 and 3, we obtain our resut. B. Average peak age Proposition 4. The average peak age in the LCFS without preemption scheme assuming Erang Ek, service time is: EP i +k k + k+. 33 Proof: We know that EP i ET i + EY i. We cacuate these two terms as foows ET i j,0 E T i Ψ i j Ψ i PΨ i j PΨ i +k qk+ + +k, 34 where we used 7 for the ast equaity. For EY i we wi ony condition on Ψ i 0. Hence using 8, we get EY i EY i Ψ i 0 PΨ i 0 +EY i Ψ i 0 PΨ i 0 k + qk. 35 Thus, combining the above two resuts we obtain our resut.

8 V. AGE OF INFORMATION FOR DETERMINISTIC SERVICE TIME In order to compute the four ages of interest under a deterministic service time assumption, we use Lemma. For that, we fix the mean of the service times S n to ES n µ, for some µ > 0, and et k. It is beyond the scope of this paper to show that if S n d Z, as k then we aso have convergence in the average ages, i.e, Sn Z. Here Sn refers to the average age corresponding to service time S n. However, we wi use this resut to derive the different ages. A. LCFS with preemption Letting k in 9 and 0, we get B. LCFS without preemption e/µ EP i µ + e/µ Letting k in and 33, we get +ρ ρ e ρ +ρ+ρe ρ +3ρ +ρe ρ EP i + e ρ µ where ρ µ. VI. NUMERICAL RESULTS In this section we show that the theoretica resuts obtained in the previous sections match the simuations. We aso compare the performance of the two transmission schemes of interest as we as the effect of the parameter k on each of them. First it is worthy to specify that a simuations were done using gamma distributed service times with a having the same mean k, except for the deterministic case where the service time is fixed to. Figure 4 presents the average age under LCFS with preemption scheme and gamma distributed service time. Two observations can be made based on this pot: i the theoretica curves given by 9 and 36 coincide with the empirica curves and ii as the vaue of k increases, the average age increases for a vaues of. This means that, under LCFS with preemption, the average age assuming deterministic service time k is higher than the average age assuming a reguar gamma distributed service. In particuar, it is higher than the average age assuming memoryess time. This observation can be expained by the fact that the probabiity of packet being preempted is given k by p + refer to Section III which is an increasing function of k. Therefore, as k increases the receiver wi have to wait on average a onger time ti a new update is deivered since the preempting rate becomes higher. This anaysis is true for any vaue of, hence the phenomenon seen in Figure 4. In a parae setting, Figure 5 presents the average age under LCFS without preemption. In this case aso two observations Average age Average age, k7 Average age, k Average age, k3.5 Average age, k Fig. 4. Average age for gamma service time S with ES, different k and LCFS with preemption can be made: i the theoretica curves given by and 38 match the empirica resuts and ii as the vaue of k increases, the average age decreases for amost a except for vaues cose to 0 where a distributions behave simiary. This difference in performance is especiay seen at high. We give here a quick intuition that expains this behavior. When is high, the time where the queue is empty goes to 0 and thus the queue is aways transmitting. This aso means that on average the waiting time W i goes to 0. Given these two observations, one can say that the system time T i and the interdeparture time Y i wi have amost the same distribution as the service time, whie being amost independent. Thus EQ i ES + ES. As for the effective rate e, since the queue is amost aways busy, the average rate at which the receiver gets new update is nothing but the inverse of the average service time, i.e e ES. Therefore, ES + ES ES + 3k. This resut which is aso obtained by taking the imit over in is decreasing with k. Hence the behavior seen in Figure 5. Next, we compare the performance of the two transmission schemes in two modes: for gamma distributed and deterministic service time. Figure 6 shows the average age under LCFS with and without preemption when the service time is taken to be gamma distributed with k. In this case we notice that for sma the two schemes perform simiary. However, for s around, the LCFS with preemption scheme performs sighty better before being outperformed by the LCFS without preemption scheme at high s. Practicay, this means that if one is using a medium whose service time is modeed as a gamma random variabe, the best strategy among the considered ones is not to preempt whie increasing the update generation rate as much as possibe. This strategy aso appies when the service time is deterministic as seen in Figure 7. In fact, we observe that for deterministic service time and for a vaues of, the average age and the average peak age for the LCFS without preemption scheme are smaer than the average age and average peak age for the LCFS with preemption respectivey. VII. CONCLUSION We considered the gamma distribution as a mode for the service time in status update systems. We computed and anayzed the average and average peak age of information

9 Average age Average age, k Average age, k6 Average age, k7 Average age, k Fig. 5. Average age for gamma service time S with ES, different k and LCFS without preemption [] R. Yates and S. Kau, Rea-time status updating: Mutipe sources, in Information Theory Proceedings ISIT, 0 IEEE Internationa Symposium on, Juy 0, pp [3] C. Kam, S. Kompea, and A. Ephremides, Age of information under random updates, in Information Theory Proceedings ISIT, 03 IEEE Internationa Symposium on, Juy 03, pp [4] S. Kau, R. Yates, and M. Gruteser, Status updates through queues, in Information Sciences and Systems CISS, 0 46th Annua Conference on, March 0, pp. 6. [5] M. Costa, M. Codreanu, and A. Ephremides, On the age of information in status update systems with packet management, CoRR, vo. abs/ , 05. [Onine]. Avaiabe: [6] S. M. Ross, Stochastic Processes Wiey Series in Probabiity and Statistics, nd ed. Wiey, Feb Average age for LCFS without preemption Average age for LCFS with preemption Average age Fig. 6. Average age for gamma service time S with k and ES Average peak age for LCFS with preemption Average age for LCFS with preemption Average peak age for LCFS without preemption Average age for LCFS without preemption 7 Age Fig. 7. Average age and average peak age for deterministic service time under two schemes: LCFS with preemption and LCFS without preemption. This aowed us to evauate these metrics for deterministic service time. This suggests that considering gamma distributions for simiar probems can be a good idea since the Gamma distributions or at east Erang distributions are practicay reevant as they can be used to mode the tota service time for reay networks. ACKNOWLEDGMENT The authors woud ike to thank Emre Teatar for hepfu discussions. REFERENCES [] S. Kau, R. Yates, and M. Gruteser, Rea-time status: How often shoud one update? in INFOCOM, 0 Proceedings IEEE, March 0, pp.

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Stochastic Complement Analysis of Multi-Server Threshold Queues. with Hysteresis. Abstract

Stochastic Complement Analysis of Multi-Server Threshold Queues. with Hysteresis. Abstract Stochastic Compement Anaysis of Muti-Server Threshod Queues with Hysteresis John C.S. Lui The Dept. of Computer Science & Engineering The Chinese University of Hong Kong Leana Goubchik Dept. of Computer

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix

More information

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with? Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Explicit overall risk minimization transductive bound

Explicit overall risk minimization transductive bound 1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC (January 8, 2003) A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC DAMIAN CLANCY, University of Liverpoo PHILIP K. POLLETT, University of Queensand Abstract

More information

Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channels

Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channels Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channes arxiv:cs/060700v1 [cs.it] 6 Ju 006 Chun-Hao Hsu and Achieas Anastasopouos Eectrica Engineering and Computer Science Department University

More information

<C 2 2. λ 2 l. λ 1 l 1 < C 1

<C 2 2. λ 2 l. λ 1 l 1 < C 1 Teecommunication Network Contro and Management (EE E694) Prof. A. A. Lazar Notes for the ecture of 7/Feb/95 by Huayan Wang (this document was ast LaT E X-ed on May 9,995) Queueing Primer for Muticass Optima

More information

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness 1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea-Time Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process

More information

Multi-server queueing systems with multiple priority classes

Multi-server queueing systems with multiple priority classes Muti-server queueing systems with mutipe priority casses Mor Harcho-Bater Taayui Osogami Aan Scheer-Wof Adam Wierman Abstract We present the first near-exact anaysis of an M/PH/ queue with m > 2 preemptive-resume

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

FRIEZE GROUPS IN R 2

FRIEZE GROUPS IN R 2 FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the

More information

Reliability: Theory & Applications No.3, September 2006

Reliability: Theory & Applications No.3, September 2006 REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate

More information

Coupling of LWR and phase transition models at boundary

Coupling of LWR and phase transition models at boundary Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of

More information

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs

More information

Some Measures for Asymmetry of Distributions

Some Measures for Asymmetry of Distributions Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester

More information

8 Digifl'.11 Cth:uits and devices

8 Digifl'.11 Cth:uits and devices 8 Digif'. Cth:uits and devices 8. Introduction In anaog eectronics, votage is a continuous variabe. This is usefu because most physica quantities we encounter are continuous: sound eves, ight intensity,

More information

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness 1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea- Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process Management,

More information

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity

More information

The Binary Space Partitioning-Tree Process Supplementary Material

The Binary Space Partitioning-Tree Process Supplementary Material The inary Space Partitioning-Tree Process Suppementary Materia Xuhui Fan in Li Scott. Sisson Schoo of omputer Science Fudan University ibin@fudan.edu.cn Schoo of Mathematics and Statistics University of

More information

Stochastic Variational Inference with Gradient Linearization

Stochastic Variational Inference with Gradient Linearization Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,

More information

Akaike Information Criterion for ANOVA Model with a Simple Order Restriction

Akaike Information Criterion for ANOVA Model with a Simple Order Restriction Akaike Information Criterion for ANOVA Mode with a Simpe Order Restriction Yu Inatsu * Department of Mathematics, Graduate Schoo of Science, Hiroshima University ABSTRACT In this paper, we consider Akaike

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

On the evaluation of saving-consumption plans

On the evaluation of saving-consumption plans On the evauation of saving-consumption pans Steven Vanduffe Jan Dhaene Marc Goovaerts Juy 13, 2004 Abstract Knowedge of the distribution function of the stochasticay compounded vaue of a series of future

More information

Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay

Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay Throughput Optima Scheduing for Wireess Downinks with Reconfiguration Deay Vineeth Baa Sukumaran vineethbs@gmai.com Department of Avionics Indian Institute of Space Science and Technoogy. Abstract We consider

More information

Fitting Algorithms for MMPP ATM Traffic Models

Fitting Algorithms for MMPP ATM Traffic Models Fitting Agorithms for PP AT Traffic odes A. Nogueira, P. Savador, R. Vaadas University of Aveiro / Institute of Teecommunications, 38-93 Aveiro, Portuga; e-mai: (nogueira, savador, rv)@av.it.pt ABSTRACT

More information

VI.G Exact free energy of the Square Lattice Ising model

VI.G Exact free energy of the Square Lattice Ising model VI.G Exact free energy of the Square Lattice Ising mode As indicated in eq.(vi.35), the Ising partition function is reated to a sum S, over coections of paths on the attice. The aowed graphs for a square

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai

More information

The Group Structure on a Smooth Tropical Cubic

The Group Structure on a Smooth Tropical Cubic The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

Scalable Spectrum Allocation for Large Networks Based on Sparse Optimization

Scalable Spectrum Allocation for Large Networks Based on Sparse Optimization Scaabe Spectrum ocation for Large Networks ased on Sparse Optimization innan Zhuang Modem R&D Lab Samsung Semiconductor, Inc. San Diego, C Dongning Guo, Ermin Wei, and Michae L. Honig Department of Eectrica

More information

A Fundamental Storage-Communication Tradeoff in Distributed Computing with Straggling Nodes

A Fundamental Storage-Communication Tradeoff in Distributed Computing with Straggling Nodes A Fundamenta Storage-Communication Tradeoff in Distributed Computing with Stragging odes ifa Yan, Michèe Wigger LTCI, Téécom ParisTech 75013 Paris, France Emai: {qifa.yan, michee.wigger} @teecom-paristech.fr

More information

AST 418/518 Instrumentation and Statistics

AST 418/518 Instrumentation and Statistics AST 418/518 Instrumentation and Statistics Cass Website: http://ircamera.as.arizona.edu/astr_518 Cass Texts: Practica Statistics for Astronomers, J.V. Wa, and C.R. Jenkins, Second Edition. Measuring the

More information

The distribution of the number of nodes in the relative interior of the typical I-segment in homogeneous planar anisotropic STIT Tessellations

The distribution of the number of nodes in the relative interior of the typical I-segment in homogeneous planar anisotropic STIT Tessellations Comment.Math.Univ.Caroin. 51,3(21) 53 512 53 The distribution of the number of nodes in the reative interior of the typica I-segment in homogeneous panar anisotropic STIT Tesseations Christoph Thäe Abstract.

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

Improving the Reliability of a Series-Parallel System Using Modified Weibull Distribution

Improving the Reliability of a Series-Parallel System Using Modified Weibull Distribution Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution

More information

A proposed nonparametric mixture density estimation using B-spline functions

A proposed nonparametric mixture density estimation using B-spline functions A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),

More information

Coded Caching for Files with Distinct File Sizes

Coded Caching for Files with Distinct File Sizes Coded Caching for Fies with Distinct Fie Sizes Jinbei Zhang iaojun Lin Chih-Chun Wang inbing Wang Department of Eectronic Engineering Shanghai Jiao ong University China Schoo of Eectrica and Computer Engineering

More information

Centralized Coded Caching of Correlated Contents

Centralized Coded Caching of Correlated Contents Centraized Coded Caching of Correated Contents Qianqian Yang and Deniz Gündüz Information Processing and Communications Lab Department of Eectrica and Eectronic Engineering Imperia Coege London arxiv:1711.03798v1

More information

NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION

NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION Hsiao-Chang Chen Dept. of Systems Engineering University of Pennsyvania Phiadephia, PA 904-635, U.S.A. Chun-Hung Chen

More information

Sequential Decoding of Polar Codes with Arbitrary Binary Kernel

Sequential Decoding of Polar Codes with Arbitrary Binary Kernel Sequentia Decoding of Poar Codes with Arbitrary Binary Kerne Vera Miosavskaya, Peter Trifonov Saint-Petersburg State Poytechnic University Emai: veram,petert}@dcn.icc.spbstu.ru Abstract The probem of efficient

More information

Cryptanalysis of PKP: A New Approach

Cryptanalysis of PKP: A New Approach Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Approximated MLC shape matrix decomposition with interleaf collision constraint

Approximated MLC shape matrix decomposition with interleaf collision constraint Agorithmic Operations Research Vo.4 (29) 49 57 Approximated MLC shape matrix decomposition with intereaf coision constraint Antje Kiese and Thomas Kainowski Institut für Mathematik, Universität Rostock,

More information

Multiple Beam Interference

Multiple Beam Interference MutipeBeamInterference.nb James C. Wyant 1 Mutipe Beam Interference 1. Airy's Formua We wi first derive Airy's formua for the case of no absorption. ü 1.1 Basic refectance and transmittance Refected ight

More information

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE KATIE L. MAY AND MELISSA A. MITCHELL Abstract. We show how to identify the minima path network connecting three fixed points on

More information

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating

More information

Packet Fragmentation in Wi-Fi Ad Hoc Networks with Correlated Channel Failures

Packet Fragmentation in Wi-Fi Ad Hoc Networks with Correlated Channel Failures Packet Fragmentation in Wi-Fi Ad Hoc Networks with Correated Channe Faiures Andrey Lyakhov Vadimir Vishnevsky Institute for Information Transmission Probems of RAS B. Karetny 19, Moscow, 127994, Russia

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna

More information

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA

More information

Reichenbachian Common Cause Systems

Reichenbachian Common Cause Systems Reichenbachian Common Cause Systems G. Hofer-Szabó Department of Phiosophy Technica University of Budapest e-mai: gszabo@hps.ete.hu Mikós Rédei Department of History and Phiosophy of Science Eötvös University,

More information

Theory and implementation behind: Universal surface creation - smallest unitcell

Theory and implementation behind: Universal surface creation - smallest unitcell Teory and impementation beind: Universa surface creation - smaest unitce Bjare Brin Buus, Jaob Howat & Tomas Bigaard September 15, 218 1 Construction of surface sabs Te aim for tis part of te project is

More information

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997 Stochastic utomata etworks (S) - Modeing and Evauation Pauo Fernandes rigitte Pateau 2 May 29, 997 Institut ationa Poytechnique de Grenobe { IPG Ecoe ationae Superieure d'informatique et de Mathematiques

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

THINKING IN PYRAMIDS

THINKING IN PYRAMIDS ECS 178 Course Notes THINKING IN PYRAMIDS Kenneth I. Joy Institute for Data Anaysis and Visuaization Department of Computer Science University of Caifornia, Davis Overview It is frequenty usefu to think

More information

How many random edges make a dense hypergraph non-2-colorable?

How many random edges make a dense hypergraph non-2-colorable? How many random edges make a dense hypergraph non--coorabe? Benny Sudakov Jan Vondrák Abstract We study a mode of random uniform hypergraphs, where a random instance is obtained by adding random edges

More information

Recursive Constructions of Parallel FIFO and LIFO Queues with Switched Delay Lines

Recursive Constructions of Parallel FIFO and LIFO Queues with Switched Delay Lines Recursive Constructions of Parae FIFO and LIFO Queues with Switched Deay Lines Po-Kai Huang, Cheng-Shang Chang, Feow, IEEE, Jay Cheng, Member, IEEE, and Duan-Shin Lee, Senior Member, IEEE Abstract One

More information

A Comparison Study of the Test for Right Censored and Grouped Data

A Comparison Study of the Test for Right Censored and Grouped Data Communications for Statistica Appications and Methods 2015, Vo. 22, No. 4, 313 320 DOI: http://dx.doi.org/10.5351/csam.2015.22.4.313 Print ISSN 2287-7843 / Onine ISSN 2383-4757 A Comparison Study of the

More information

On the Goal Value of a Boolean Function

On the Goal Value of a Boolean Function On the Goa Vaue of a Booean Function Eric Bach Dept. of CS University of Wisconsin 1210 W. Dayton St. Madison, WI 53706 Lisa Heerstein Dept of CSE NYU Schoo of Engineering 2 Metrotech Center, 10th Foor

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

(Refer Slide Time: 2:34) L C V

(Refer Slide Time: 2:34) L C V Microwave Integrated Circuits Professor Jayanta Mukherjee Department of Eectrica Engineering Indian Intitute of Technoogy Bombay Modue 1 Lecture No 2 Refection Coefficient, SWR, Smith Chart. Heo wecome

More information

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation

More information

Asynchronous Control for Coupled Markov Decision Systems

Asynchronous Control for Coupled Markov Decision Systems INFORMATION THEORY WORKSHOP (ITW) 22 Asynchronous Contro for Couped Marov Decision Systems Michae J. Neey University of Southern Caifornia Abstract This paper considers optima contro for a coection of

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,

More information

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Efficient Visual-Inertial Navigation using a Rolling-Shutter Camera with Inaccurate Timestamps

Efficient Visual-Inertial Navigation using a Rolling-Shutter Camera with Inaccurate Timestamps Efficient Visua-Inertia Navigation using a Roing-Shutter Camera with Inaccurate Timestamps Chao X. Guo, Dimitrios G. Kottas, Ryan C. DuToit Ahmed Ahmed, Ruipeng Li and Stergios I. Roumeiotis Mutipe Autonomous

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, we as various

More information

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this

More information

STA 216 Project: Spline Approach to Discrete Survival Analysis

STA 216 Project: Spline Approach to Discrete Survival Analysis : Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing

More information

Efficiently Generating Random Bits from Finite State Markov Chains

Efficiently Generating Random Bits from Finite State Markov Chains 1 Efficienty Generating Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown

More information

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu

More information

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

Mat 1501 lecture notes, penultimate installment

Mat 1501 lecture notes, penultimate installment Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the

More information

Determining The Degree of Generalization Using An Incremental Learning Algorithm

Determining The Degree of Generalization Using An Incremental Learning Algorithm Determining The Degree of Generaization Using An Incrementa Learning Agorithm Pabo Zegers Facutad de Ingeniería, Universidad de os Andes San Caros de Apoquindo 22, Las Condes, Santiago, Chie pzegers@uandes.c

More information

arxiv: v2 [cond-mat.stat-mech] 14 Nov 2008

arxiv: v2 [cond-mat.stat-mech] 14 Nov 2008 Random Booean Networks Barbara Drosse Institute of Condensed Matter Physics, Darmstadt University of Technoogy, Hochschustraße 6, 64289 Darmstadt, Germany (Dated: June 27) arxiv:76.335v2 [cond-mat.stat-mech]

More information

Approximated MLC shape matrix decomposition with interleaf collision constraint

Approximated MLC shape matrix decomposition with interleaf collision constraint Approximated MLC shape matrix decomposition with intereaf coision constraint Thomas Kainowski Antje Kiese Abstract Shape matrix decomposition is a subprobem in radiation therapy panning. A given fuence

More information

From Margins to Probabilities in Multiclass Learning Problems

From Margins to Probabilities in Multiclass Learning Problems From Margins to Probabiities in Muticass Learning Probems Andrea Passerini and Massimiiano Ponti 2 and Paoo Frasconi 3 Abstract. We study the probem of muticass cassification within the framework of error

More information

Rate-Distortion Theory of Finite Point Processes

Rate-Distortion Theory of Finite Point Processes Rate-Distortion Theory of Finite Point Processes Günther Koiander, Dominic Schuhmacher, and Franz Hawatsch, Feow, IEEE Abstract We study the compression of data in the case where the usefu information

More information

Improving the Accuracy of Boolean Tomography by Exploiting Path Congestion Degrees

Improving the Accuracy of Boolean Tomography by Exploiting Path Congestion Degrees Improving the Accuracy of Booean Tomography by Expoiting Path Congestion Degrees Zhiyong Zhang, Gaoei Fei, Fucai Yu, Guangmin Hu Schoo of Communication and Information Engineering, University of Eectronic

More information

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1 Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rues 1 R.J. Marks II, S. Oh, P. Arabshahi Λ, T.P. Caude, J.J. Choi, B.G. Song Λ Λ Dept. of Eectrica Engineering Boeing Computer Services University

More information

AALBORG UNIVERSITY. The distribution of communication cost for a mobile service scenario. Jesper Møller and Man Lung Yiu. R June 2009

AALBORG UNIVERSITY. The distribution of communication cost for a mobile service scenario. Jesper Møller and Man Lung Yiu. R June 2009 AALBORG UNIVERSITY The distribution of communication cost for a mobie service scenario by Jesper Møer and Man Lung Yiu R-29-11 June 29 Department of Mathematica Sciences Aaborg University Fredrik Bajers

More information

CONSISTENT LABELING OF ROTATING MAPS

CONSISTENT LABELING OF ROTATING MAPS CONSISTENT LABELING OF ROTATING MAPS Andreas Gemsa, Martin Nöenburg, Ignaz Rutter Abstract. Dynamic maps that aow continuous map rotations, for exampe, on mobie devices, encounter new geometric abeing

More information

Two-Stage Least Squares as Minimum Distance

Two-Stage Least Squares as Minimum Distance Two-Stage Least Squares as Minimum Distance Frank Windmeijer Discussion Paper 17 / 683 7 June 2017 Department of Economics University of Bristo Priory Road Compex Bristo BS8 1TU United Kingdom Two-Stage

More information

Formulas for Angular-Momentum Barrier Factors Version II

Formulas for Angular-Momentum Barrier Factors Version II BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A

More information

Efficient Generation of Random Bits from Finite State Markov Chains

Efficient Generation of Random Bits from Finite State Markov Chains Efficient Generation of Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown

More information

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

More information