Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay

Size: px
Start display at page:

Download "Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay"

Transcription

1 Throughput Optima Scheduing for Wireess Downinks with Reconfiguration Deay Vineeth Baa Sukumaran Department of Avionics Indian Institute of Space Science and Technoogy. Abstract We consider wireess downinks where the base station dynamicay switches between different users in order to transmit data intended for the respective users. When the base station switches from serving one user to another there is a reconfiguration deay. For such wireess downinks with reconfiguration deay we consider the probem of throughput optima scheduing. We propose the 1-ookahead scheduing poicy and anayticay show that it is throughput optima. We obtain the 1-ookahead poicy by using an approximate soution to a Markov decision process formuation of the scheduing probem. The approximate soution is aso used to expain the biased maxweight form of 1-ookahead as we as an existing poicy. Index Terms Wireess downink Random connectivity Reconfiguration deay Stabiity region Throughput optimaity I. INTRODUCTION We consider a wireess downink mode where the base station switches between different users in order to transmit data intended for the respective users. When switching between users the ogica inks between the user and the base station need to be configured. For exampe the user s state may need to be changed from an ide to an active state. This incurs a reconfiguration deay which is the deay between the time at which the base station scheduer decides to serve the user and the time at which the actua data transmission starts. We note that such reconfiguration deay aso arises in other cases: such as sateite systems with mechanicay steered antennae eectronic beamforg optica routers 1 and radio transceivers 3. The other important feature of such systems is that the service of a user is aso affected by the connectivity of that user to the service station. For exampe the connectivity of the user is through a wireess channe subject to fading therefore the connectivity is random over time. Motivated by such scenarios we consider the stabiity region and throughput optima scheduing for a wireess downink mode with random connectivity over time and reconfiguration deay see Fig. 1. For wireess networks with reconfiguration deay maxweight poicies 7 are not throughput optima since such poicies switch between queues very frequenty 2. Prior work in 1 2 and 4 had proposed heuristic throughput optima poicies for such systems. Ceik et a. 2 proposed the variabe frame max-weight VFMW poicy which reduces reconfiguration deay overhead by restricting switching to happen ony at the ends of scheduing frames for wireess networks with random connectivity and reconfiguration deay. Hsieh et a. 4 proposed the queue biased max weight QBMW poicy that is throughput optima for a queueing system with reconfiguration deay. The average deay performance of VFMW poicies was Fig. 1. An exampe wireess downink with two users and a base station BS. The data for the users have random packet arrivas A i t into their packet buffers at the BS. The BS can transmit at most one packet to one user from its queue in a sot. The connections between the BS and the users are randomy on or off across time sots. When the BS switches from serving one user to another there is a reconfiguration deay of sots. improved by the QBMW poicy. However why the specific biased form of the QBMW poicy was needed was not addressed in their paper. In this paper we propose the 1-ookahead poicy which is another throughput optima poicy for systems with reconfiguration deay. In contrast to VFMW and QBMW we use a forma Markov decision theoretic formuation to anayticay motivate the need for the bias term. The same approach can be used to motivate the bias term used for QBMW poicies. Outine Contributions and Notation: In Section II we discuss the queueing mode that we use for anayzing the wireess downink shown in Fig. 1 as we as its stabiity region. Our main contribution in this paper is the anayticay wemotivated throughput optima 1-ookahead poicy which we propose in Section III. We then discuss a Markov decision process formuation which is used to motivate the definition of 1-ookahead poicies in Section IV. Our first secondary contribution is that we expain the bias term appearing in the QBMW poicy which has a simiar form as 1-ookahead. Another secondary contribution is in the proof of throughput optimaity of 1-ookahead and its reation to the proof of throughput optimaity of QBMW; we provide simpifications and corrections to the proof in 4. In this paper a vectors are coumn vectors and vector transposes are denoted by. T. II. SYSTM MODL AND PROBLM FORMULATION We consider a system of N parae queues served by a singe server to mode a wireess downink for N users 1 see Fig. 1 1 We note that the mode can aso appy to a wireess upink but with additiona assumptions on the avaiabiity of queue ength information at a centraized scheduer.

2 for an exampe with N = 2. The system evoves in sotted time with the sots indexed by t { }. In each sot t a random number A i t of packets arrive to the base station BS destined for the i th user for every i { N}. For every user these packets are queued in a separate infinite ength buffers at the BS. We assume that the random process A i t t 0 is an independent and identicay distributed IID process e.g. A i t can be modeed by a Bernoui process. The arriva processes to different buffers are aso assumed to be independent. We denote the arriva rate A i 0 as λ i and the coumn vector λ 1... λ N T as λ. The BS scheduer decides which user s queue is served in a sot; S i t = 1 if the scheduer decides to serve user i in sot t and 0 otherwise. We assume that N i=1 S it = 1 i.e. at most one user can be served in a sot. We assume that there is a connection random process C i t {0 1} t 0 associated with user i that modes whether the BS is connected to user i in sot t. We assume that the the processes C i t t 0 are IID. We denote the average connection rate of queue i as µ i i.e. C i 0 = µ i. We aso assume that C i t t 0 and C j t t 0 for two different queues i and j are independent. We assume that the channe connectivity C j t j is not known to the scheduer at the start of sot t 2. We assume that if the BS is connected to user i and the scheduer decides to serve user i then at most one packet is removed from queue i in sot t. We assume that there is a reconfiguration deay of sots see Fig. 1 if the server switches from serving one queue to another. We assume that 1. We denote the number of sots to finish reconfiguration at time t by Rt. If the scheduer switches from one queue to the another at sot t Rt = and then Rt decrements by one for every sot unti Rt = 0 if the scheduer stays with the queue that it has switched to. Once reconfiguration is finished Rt stays at zero unti the scheduer switches to another queue. We note that a packet wi be removed from the i th queue in the t th sot ony if C i ts i ti {Rt=0} = 1. We define I i t = C i ts i ti {Rt=0} which is caed the service opportunity for queue i in sot t. We denote the number of packets in the i th queue at the beginning of sot t as Q i t. Then Q i t + 1 = Q i t I i t + + A i t 1 where x + = maxx 0. We note that the packets which arrive in sot t are assumed to stay in the system for at east one sot. A poicy µ is the sequence of decisions S1 S2... where St is the vector S 1 t... S N t T. The time average tota queue ength under a poicy µ is 1 im sup T T T 1 N i=1 Q it and is denoted as qµ. The queueing system with an arriva rate λ is said to be stabe under the poicy µ if the time average tota queue ength qµ is finite. The stabiity region Λ of the system is the set of arriva rate vectors λ where for each λ there exists a poicy 2 Since C i t are IID this information is not usefu. µ possiby dependent on λ such that the queueing system is stabe. From 3 we have that if λ Λ then there exists β j j such that β j 0 j β j 1 such that λ j β j I j µ j where I j is a coumn vector of ength N with 1 at the j th position and zero otherwise. It is known that the VFMW poicy from 2 can be used to achieve any point in the stabiity region Λ. We note that both VFMW and QBMW are poicies that are throughput optima for the above queueing system. In the next section we propose another throughput optima poicy 1-ookahead; the form of this poicy is anayticay motivated using a Markov decision theoretic formuation. III. TH 1-LOOKAHAD POLICY For the queueing mode any scheduing poicy at the start of a sot needs to decide whether to keep serving the current queue or to switch to any other queue. Intuitivey we shoud serve a queue with arge queue ength so that the queue ength can be reduced and argest possibe throughput since the reduction in queue ength woud be the most. This is the intuition behind the decision rue that is proposed as the 1-ookahead poicy. We define our 1-ookahead poicy in Agorithm 1. We assume that the server is serving queue i at t 1 in Agorithm 1. For the 1-ookahead poicy the weights W i t or W j t can be cacuated as the expected sum throughput that can be obtained in the immediate future consisting of + 1 sots under the decision of staying with queue i or switching to the j th queue. Then if we stay with queue i at t the expected sum throughput i.e. W i t is given by µ i. If we decide to switch to queue j at t then the expected sum throughput i.e. W j t is given by µ j. We find that a heuristic modification to the weight W i t is required to show that the 1-ookahead poicy is in fact throughput optima. We redefine W i t as F Qt where F Qt = max1 i Q it α where α 0 1. µ i Agorithm 1 1-ookahead poicy 1: Cacuate the weights W j t j { N}. 2: If W i tq i t < W j tq j t then switch to the j th queue ese stay with the i th queue. A natura question that arises in the definition of 1-ookahead poicies is why we are ooking at the expected tota throughput for + 1 sots and not + m sots for some m 1. In the next section we motivate this and the exact form for the 1-ookahead poicy by considering a Markov decision process MDP formuation for the probem of imizing time average tota queue ength for the queueing system discussed in our paper. We show that the 1-ookahead poicy arises from a heuristic approximate soution to the average cost optimaity equation for a Markov decision process formuated for the queueing system in our paper. The reason why m = 1 is reated to the need for keeping the poicy specification simpe.

3 We now show that the 1-ookahead poicy has finite average queue ength for any arriva rate vector within the stabiity region of the system and is therefore throughput optima. Theorem III.1. For IID channe connection processes for any arriva rate that is within the stabiity region the 1-ookahead poicy has finite average queue ength. The proof of this theorem is presented in Appendix A. We note that the proof of this theorem borrows ideas from the proof of the throughput optimaity of QBMW in 4. Comparison with QBMW: We note that the form of 1- ookahead poicy is simiar to that of the QBMW poicy. However for the QBMW poicy the weight W i t is cacuated as t. Since the form of 1-ookahead is amost the same as that of QBMW the derivation of the weight terms for the 1-ookahead given in the next section aso motivates the form of the QBMW poicy. We aso note that the proof of throughput optimaity of QBMW reies upon showing that the drift of a Lyapunov function defined as the sum of squared queue engths in a sot is negative when considered over frames mutipe sots. We expect that during reconfiguration deay after a switch since there is no service the Lyapunov drift woud not be negative. Hsieh et a. 4 define a frame size parameter T k in order to show that the drift is indeed negative when considered over T k sots which incude the reconfiguration deay. However the proof in 4 does not consider what happens when T k <. We address this probem in our paper. F Qt k where t k is the time of the ast switch before IV. MARKOV DCISION PROCSS FORMULATION We consider a two queue system here for ease of exposition. We note that instead of finding a scheduing poicy such that the average queue ength under that poicy is just finite we reformuate our objective to find a scheduing poicy that imizes im sup T T 1 1 T Q 1 t + Q 2 t. 2 In order to find a scheduing poicy that imizes 2 defined above we use a Markov decision process MDP formuation see 8 9. The state of the MDP is defined to be St = Q 1 t Q 2 t Rt Mt where Mt {0 1} indicates whether the first queue or second queue has been in service at sot t 1 respectivey. The state space of the MDP is the Cartesian product of the state spaces of the individua components i.e. Z + Z + {0... } {0 1}. We denote a specific state vector vaue as s = q 1 q 2 r m T. The action taken at sot t is γt which is defined as either staying with the current queue γt = 0 or switching to the other queue γt = 1. The evoution of the MDP s state from sot to sot i.e. St to St + 1 is defined in terms of its components as foows: 1 if γt = 0 and Rt > 0 then Rt + 1 = Rt 1 and since 1 there is no service from the Mt th queue; we denote this as I Mt t = 0 2 if γt = 0 and Rt = 0 then we have I Mt t = 1 3 if γt = 1 then Rt = and I Mt t = 0; Mt changes to the queue that was switched to 4 Q i t + 1 = maxq i t I i t 0 + A i t. Since we are interested in imizing the time average expected tota queue ength the singe stage cost of the MDP at sot t is chosen as Q 1 t + Q 2 t. Suppose the average cost optimaity equation ACO 9 Chapter 6 exists 3 for the MDP. The ACO is of the form 9 Chapter 6 Theorem hs = { q q 2 g + hs +1 s } 3 where hs is the reative vaue function; which is a function of the state s = q 1 q 2 r m T g is the optima imum average cost or sum of average queue engths and S +1 is the state that the MDP evoves to in one step starting from s according to the MDP evoution described above. We note that q q 2 is the singe stage cost when the state of the MDP is s. The possibe actions which can be taken when the state is s is: a 0 stay with the current queue or b 1 switch to the other queue. The optima poicy is a stationary poicy that chooses an action γ {0 1} in order to imize the expression within the imization in the RHS of 3. This optima stationary poicy prescribes an action γ as a function γs of the state s. Note that if the function hs is known then the optima stationary poicy can be competey characterized. However in 3 both hs and g are not known. In most cases the above functiona equation cannot be soved anayticay for hs and g. Vaue iteration 9 Section 6.6 is an iterative procedure that can be used to obtain a soution to the ACO. We et V 1 s = q q 2. We define V n+1 s = { } q q 2 + V n S +1 s. 4 We reca that S +1 is the state that the MDP evoves to in a singe sot starting from state s. We note that V n s caed the vaue function is the imum expected cumuative sum of queue engths when the system evoution is considered over n sots starting with state s or S0 = s i.e. n 1 N i=1 Q it S0 = s. From 9 eq for arge enough n we have that the γ which attains the imum in the above vaue iteration equation 4 for every s is an approximatey-optima stationary poicy γs for the average cost probem. Furthermore again from 9 eq we have that for arge n V n s = hs + ng. 5 The motivating idea behind the definition of the 1-ookahead poicy is that a good poicy can be obtained from the ACO 3 3 We note that one of the sufficient conditions for the existence of the ACO is that there exists a poicy under which 2 is finite. Since the VFMW poicy has a finite average queue ength if the arriva rate vector is within Λ the above sufficient condition is satisfied. Other sufficient conditions which dea with irreducibiity of the queue ength Markov chain for the existence of the ACO can be shown to hod under appropriate assumptions on the distribution of the arriva random variabes. Since this informa discussion is to motivate the form of the 1-ookahead poicy these detais are not incuded here.

4 Reative vaue function hq q r m q 1 = 10 r = 0 m = 1 q 1 = 10 r = 2 m = 1 q 1 = 30 r = 0 m = Queue ength q 2 Fig. 2. The reative vaue function hs where s = q 1 q 2 r m potted as a function of q 2 for different q 1 r and m. The reative vaue function has been obtained from vaue iteration carried out for a system with buffer size for both queues truncated to 1000 channe connectivity parameters µ 1 = µ 2 = 0.5 and arriva rates λ 1 = 0.15 and λ 2 = 0.2. We observe that hs need to be approximated by a non-inear function of q 2. or the vaue iteration 4 if a good approximation can be found for hs or V n s. The approximations for hs or V n s can be substituted into the RHS of the ACO 3 and 4 and the imizing γ can be proposed as a candidate poicy. Suppose an approximation to hs is hs. We note that except for a constant hs is aso an approximation for V n s from 5. We first use a fuid mode 6 Chapter 10 of the queueing system in order to motivate that hs = i q2 i is a reasonabe approximation to the function hs. Under the optima poicy et µ i be the time average of the service opportunity which is given to queue i. The fuid mode of the queueing system modes each queue ength as a deteristic function q i t which evoves according to the differentia equation dq i t dt = µ i λ i. We assume that the initia state q i 0 of queue i is q i the component of our state s. We then obtain that q i t = maxq i µ i λ i t 0. The area under the q i t function or the cumuative queue ength for queue i is proportiona to qi 2. The tota cumuative queue ength which is then proportiona i q2 i is the vaue function for the fuid mode. From 6 Theorem we have that the vaue function for the fuid mode is the same as the reative vaue function hs where the queue ength components of s are q 1 and q 2 used in the fuid mode in an asymptotic regime where q 1 q 2. So we approximate the reative vaue function hs with i q2 i. The second motivation for the use of i q2 i as an approximation comes from numericay soving the ACO using vaue iteration. In Fig. 2 we pot hs as a function of q 2 for different q 1 r and m. where hs has been obtained using vaue iteration for 1000 iterations and for queue ength state space truncated to The other parameters used are indicated in Fig. 2. We find that a second order poynomia i.e. c 1 q2 2 + c 2 q 2 + c 3 is a good choice for fitting the observed function. When higher degree poynomias are used the coefficients of higher degree terms are seen to be cose to zero. We note that we excude inear terms as we as cross terms i.e. q 1 q 2 in our approximation for the sake of anaytica simpicity. We now use the sum of squares approximation to motivate the form of the 1-ookahead poicy. From 4 we have that a poicy that achieves the imum in { q q 2 + V n S + s } for arge enough n is approximatey optima for the average queue ength imization probem. Suppose we use the sum of squares approximation i q2 i for V ns motivated by 5. Then we are considering the poicy that chooses γ to achieve the imum as foows { q q 2 + Q Q s }. Here Q 1 1 and Q 2 1 are the queue engths that the MDP evoves to under an action γ according to the MDP evoution described above starting with Q 1 0 = q 1 and Q 2 0 = q 2. Let us consider the case where s is such that queue 1 is currenty being served a simiar discussion hods for the case where queue 2 is being served. For γ = 0 we have that Q Q s = q 1 C A q 2 + A 2 2 and for γ = 1 we have that Q Q s = q A q 2 + A 2 2 since we have 1. We see that γ = 1 wi never be chosen; this happens of course because our approximation does not capture the possibe service that can happen for queue 2 after sots. The approximation does not capture the possibe service since the approximation athough simpe does not incude other state variabes especiay r. In order to capture this potentia service after sots we proceed by writing V n+1 s as: { +m 1 V n+1 s = Q 1 τ + Q 2 τ s + V n +m 1S s } ++m. Here Q 1 0 = q 1 and Q 2 0 = q 2 where q 1 and q 2 are the queue ength components of s. Aso Q i τ is the queue ength that the MDP evoves to in the τ th sot after taking action γ in the first sot and then the optima actions in a sots. Then speciay note that the expectation +m 1 Q 1 τ + Q 2 τ s is computed and the state S ++m is what the MDP evoves to in + m sots under the use of γ in the first sot and the optima poicy after that. We note that the above aternate expression for V n+1 s can be written for any n > +m but we are interested in sufficienty arge vaues of n. We now expicity indicate that the optima actions are taken for the first + m 1 sots by using a imization over action variabes γ 1... γ +m 1. We have that V n+1 s is γγ 1...γ +m 1 { +m 1 Q 1 τ + Q 2 τ s V n +m 1S ++m s }. 6 We note that since n is arge the sum of squares approximation can aso be appied to V n +m 1s.

5 We now consider the case where m = 1. For m = 1 we note that { V n+1 s = Q 1 τ + Q 2 τ s γγ 1...γ + V n S s } Assug that s is such that queue 1 is being served the possibe action sequences can be divided into three: A1: γ = γ 1 = = γ = 0 or stay with queue 1 for sots A2: γ = 1 γ 1 = = γ = 0 or switch to queue 2 in the first sot and then stay with queue 2 for the rest of sots A3: Any other action sequence. In the foowing in order to differentiate between queue ength evoution under different action sequences we use Q j i t to denote the queue ength for queue i under action sequence Aj e.g. Q 3 1t for queue ength of queue 1 under action sequence A3. For the action sequence A1 we have for a τ { } Q 1 1τ = q 1 C 1 + A 1 and Q 1 2τ = q 2 + A 2. For the action sequence A2 we have for a τ { } and for a τ {0... } Q 2 1τ = q A 1 Q 2 2τ = q 2 + A 2. We note that there is no service under A2 for the first sots but at the + 1 th we do have a service i.e Q = q 2 C 2 + A 2. We note that the action sequence corresponding to A3 above is of the form γ = 0 γ 1 = 0... γ τs 1 = That is we switch to queue 2 from queue 1 in the τs th sot where τ s > 1. We note that for m = 1 after τ s there cannot be any more service in the system in the first sots since there is a reconfiguration deay of sots. Then we have that for A3 and any τ { } ττ s 1 Q 3 1τ = q 1 C 1 + A 1 and Q 3 2τ = q 2 + A 2. We note that for the same sampe path of arrivas A 1 A 2 and channe connectivity C 1 and C 2 and for any τ { } Q 1 1τ Q 3 1τ and Q 1 2τ = Q 3 2τ. Aso under the assumption of the sum of squares form for V n. we have that under A1 V n S ++1 s = Q 1 i s which is ess than the foowing corresponding vaue V n S ++1 s = Q 3 i s under A3. Therefore we concude that for m = 1 and a sum of squares approximation for V n s the optima sequence of actions is either A1 or A2 in the above ist but not any action sequence in A3. Thus for m = 1 with the sum of squares form for V n. we have that the optima choice is between A1 and A2; which incidentay is a choice between staying with the current queue or switching to the other queue. This choice forms the basis for the definition of the 1-ookahead note that m = 1 poicy. From 7 and the fact that the optima sequence is either A1 or A2 we have the foowing decision rue. We choose γ = 0 if Q 1 i s < Q 1 i τ + i i i Q 2 i τ + i i i Q 2 i s and γ = 1 otherwise. We have that Q 1 1τ + Q 1 2τ s = 9 q i + i A i C 1 and i Q1 i s is q A 1 τ C 1 τ 2 + q 2 + A 2 τ 2 s which can be simpified as 2 = C 1 τ C q 1 2 C 1 τ C q 1 C 1 τ C A 2 τ + 2q 2 s 8 2 A 1 τ A 1 τ A 1 τ + q1 2 + q2 2 + A 2 τ. 10

6 Simiary we have that Q 2 1τ + Q 2 2τ s = 11 and q A 1 + q 2 + A 2 s Q Q s = q A 1 τ 2 + q 2 C 2 + A 2 τ 2 s = q A 1 τ + 2q 1 2 A 2 τ + 2q 2 A 1 τ + q2 2 + A 2 τ + C 2 2 C 2 0 2q 2 C 2 C C 2 C 2 0 A 2 τ. 12 We note that the comparison in 8 is equivaent to the comparison < Severa terms are common in and We cance out those common terms and keep ony those terms which have a queue ength term appearing in it. For arge vaues of queue engths these terms woud doate other constant terms. Furthermore such a choice woud resut in the decision rue for 1-ookahead to have a simpe biased max-weight form. Then we get a heuristic decision rue of choosing γ = 0 if 2q 1 C 1 τ C 1 0 < 2q 2 C 2 C 2 0 or q 1 C 1 τ C 1 0 > q 2 C 2 C 2 0 and γ = 1 otherwise. We note that this is the motivation behind the definition of the biased max-weight 1-ookahead poicy. We now consider the case where m > 1. Using a simiar sequence of steps as for m = 1 it is possibe to derive a decision rue comparing A1 with A2. This comparison between A1 and A2 eads to the rue of choosing γ = 0 if +m 1 +m 1 q 1 C 1 τ C 1 0 > q 2 C 2 τ C 2 0 τ= However for m > 1 we note that it is not possibe to concude that the optima sequence of actions is either A1 or A2 but not any in A3. We iustrate this via an exampe here. Consider 6 for = 2 and m = 3. Assug that we are serving queue 1 a possibe sequence of decisions is γ = 0 γ 1 = 1 γ 2 = 0 γ 3 = 1 γ 4 = 0 i.e. we stay with queue 1 then switch to queue 2 wait for 2 sots switch back to queue 1 and wait for 2 sots. Such a sequence of decisions might ead to service for both queue 1 and queue 2. We note that the simpe decision rue that we have proposed in this case just compares A1 and A2 and not any poicy in A3 such as the one above. If we incude other comparisons to action sequences of type A3 the simpicity of the definition of 1-ookahead poicies does not carry over. Since 1-ookahead poicies are throughput optima and have a simpe specification we restrict our attention to the case of m = 1. We note that the correction to the weight or bias using optimaity. F Qt is motivated by the proof of throughput V. CONCLUSIONS AND FUTUR WORK In this paper we considered queueing system modes with reconfiguration deay. We proposed a new 1-ookahead poicy using an approximate soution to a Markov decision process formuation. The biased max-weight form of the 1-ookahead poicy is anayticay we motivated. Furthermore since the form of 1-ookahead is simiar to that of QBMW the bias term which appears in QBMW can aso be expained using this anaytica deveopment. We aso prove that 1-ookahead is throughput optima. We note that the Markov decision process formuation can be extended to the case of correated channe connectivity by incuding channe states in the MDP s state. It is possibe to again define a 1-ookahead poicy for the correated channe case see 5; proving the throughput optimaity of such a poicy woud be part of our future work. RFRNCS 1 Güner D Ceik Sem C Borst Phiip A Whiting and ytan Modiano. Dynamic scheduing with reconfiguration deays. Queueing Systems 831-2: Güner D Ceik Long B Le and ytan Modiano. Dynamic server aocation over time-varying channes with switchover deay. I ansactions on Information Theory 589: Güner D Ceik and ytan Modiano. Scheduing in networks with timevarying channes and reconfiguration deay. I/ACM ansactions on Networking 231: Ping-Chun Hsieh I Hou and Xi Liu. Deay-optima scheduing for queueing systems with switching overhead. arxiv preprint arxiv: Amit Kumar and Vineeth B. S. Scheduing poicies for wireess downink with correated random connectivity and mutisot reconfiguration deay. Accepted to I Communications Letters S. P. Meyn. Contro Techniques for Compex Networks. Cambridge University Press Michae J Neey. Stochastic network optimization with appication to communication and queueing systems. Synthesis Lectures on Communication Networks 31: Linn I Sennott. Stochastic dynamic programg and the contro of queueing systems voume 504. John Wiey & Sons C. H. Tijms. A first course in stochastic modes. Wiey APPNDIX A PROOF OF THORM III.1 We prove the stabiity of the modified 1-ookahead poicy in this section. Suppose t k is the sot at which the k th switch happens. Let T k = t k+1 t k be the duration of the k th frame. Let us assume that in the k th frame the i th queue is being

7 served. Then we note that at t k µ i Q i t k > µ j Q j t k for every j. Furthermore we have that there exists a j such that Q i t k + T k µ i < Q j t k + T k µ j. F Qt k + T k We first prove the foowing emma which states that for arge enough queue engths the duration of the frame T k is arge. The proof is simiar to Lemma 3 in 4. Lemma A.1. The ength T k of the k th frame for every k is such that where we have used T k 1. Using this upper bound on F. in 13 we obtain that T 1+α Q i t k µ i k > A max N + Q t k µ α + A max. We aso note that Q i t k µ i Q t k µ and hence Q i t k µ i. Hence Q t k µ N T 1+α k > Q t k µ N A max N + Q t k µ α + A max. T 1+α k > Q t k µ N A max N + Q t k µ α + A max where α 0 1 is as in the definition of F Qt. Proof. From the definition of T k or t k+1 we have that there exists a j such that Q i t k + T k µ i < Q j t k + T k µ j. F Qt k + T k We note that Q i t k + T k Q i t k T k Q j t k + T k Q j t k + A max T k where we have used the fact that in T k sots the maximum amount of service for queue i is T k and the maximum number of arrivas for queue j is A max T k. We aso note that the first ower bound coud be negative. Substituting these bounds we obtain that Q i t k T k µ i < F Qt k + T k Q j t k + A max T k µ j Q i t k µ i µ j Q j t k < F Qt k + T k T k µ i + A max µ j where we have used F Qt k + T k 1. Since µ i and µ j 1 and µ i Q i t k µ j Q j t k we have that T k > Q i t k µ i F Qt k + T k + A max. 13 We note that α F Qt k + T k Q t k + A max T k µ A max T k N + Q t k µ α where we have used that there are N queues and µ 1. Then we have that F Qt k + T k T α k T α k A max N + A max N + Q t k µ α Q t k µ α. An important coroary of this emma is stated beow. Coroary A.2. Given a Z + there is a finite set Q of queue ength vectors such that T k + for Qt k Q c. Proof. From Lemma A.1 we have that T k > Q 1 t k µ N A max N + 1+α Q t k µ α. + A max We note that the ower bound on T k is an increasing function of Q t k µ. So for a given there is a finite δ such that T k + for Q t k µ > δ. There is a finite set of queue ength vectors Q such that Q t k µ δ for which the above ower bound on T k does not guarantee that T k > +. We now prove the stabiity of the 1-ookahead poicy using Lyapunov drift arguments. The Lyapunov function is chosen to be Lt = Q t 2. From Coroary A.2 we note that at t k for a queue ength vectors except for a finite set Q we have that T k +. We first consider the case of the compement set of Q. So T k + in the discussion beow. In the stabiity proof we wi consider the Lyapunov drift of the system over a frame first between t k and t k + + and then for every sot in {t k t k + T k }. We note that in 4 the drift is considered over t k and t k + T k instead of between t k and t k + +. In 4 T k is defined as T k F Qt k. It is not cear why T k woud be greater than ; this property is essentia to the drift arguments that foow in 4. In contrast we consider Tk to be a constant + but ony for those t k for which the queue ength vector ies in Q c. Then we show that whatever Lyapunov drift that we obtain for the case where Qt k Q c can be appied to the case where Qt k Q but with a different constant factor in the drift expression. The Lyapunov drift between t k and t k + + is Q t k Q t k 2. We note that Q t k + + max + 1 A t k + t. Q t k + 1 I t k + t 0 +

8 With Γ defined as max Q t k + 1 I t k + t 0 Q t k + 1 I t k + t we have that Q t k + + Q t k + 1 I t k + t + Γ A t k + t. Then 14 can be bounded above as Q t k I t k + t + Γ + A t k + t Q t k 2 B Q A t k + t I t k + t Here B0 is a constant which is obtained by bounding I t k +t above by 1 A t k above by A max and the fact that Γ is at most + and is non-zero ony if Q +. We note that B0 is a function of. The expected Lyapunov drift which is conditioned on Qt k is then B Q + λ I {=i} µ i where I {=i} is 1 if = i and 0 otherwise since in the k th frame there is service ony for the i th queue. Furthermore in the first sots after t k we do not have service since the system is in reconfiguration. We now write the expected Lyapunov using vector notation as being B0 + 2 λ T + Ii µ i Qt k 14 where a vectors are coumn vectors and I i is a vector of zeros except with 1 at the i th position. Suppose we consider a λ in the stabiity region then we have that there exists a set of β j j β j 1 such that λ j β j I j µ j and we can write λ = 1 ɛ j β ji j µ j for some ɛ > 0. Furthermore we have that µ j I j T Qt k µ i I i T Qt k since we have switched to i at t k. Then we have that 14 is B ɛ T I i µ i Qt k + For any ɛ > 0 we choose such that 1 ɛ T = ɛ r+ 1 for some ɛ 1 > 0. Then we have that the expected Lyapunov drift between t k + + and t k is B 0 2ɛ 1 µ i I i T Qtk. Now since µ j I j T Qt k µ i I i T Qt k we have that B0 2 ɛ 1 µ Q t k. N. We note that Q t k + t Q t k + A max t or µ Q t k + t µ A max µ Q t k. Then the expected Lyapunov drift between t k and t k + + is B1 ɛ µ Q t k + t + N where B1 is a constant dependent on. We now consider the sot to sot drift of the Lyapunov functions for sots t {t k t k + T k }. We note that we are sti restricting to those t k at which Qt k Q c. The Lyapunov drift is Q t Q t 2 As in the previous case we have that Q t + 1 = Q t I t + A t + Γ and it can be shown that the Lyapunov drift is B Q t A t I t. The expected Lyapunov drift conditioned on Qt is B Q t λ I {=i} µ i. Here we note that there is ony service to queue i in frame k and we are considering sots after the reconfiguration deay. In vector notation we have T B λ I i µ i Qt. 15 Since λ is in the stabiity region we again have that 15 is B ɛ j β j I j µ j I i µ i T Qt. We note that at t µ j I j T Qt F Qt µ i I i T Qt. Substituting we get the conditiona drift as B ɛ 1 µ i I i T Qt F Qt = B ɛ + 1 ɛ µ i I i T Qt F Qt = B 2 2ɛµ i I i T Qt + 21 ɛ µ i I i T Qt. F Qt We note that µ i I i T Qt µ I T Qt. Aso µ i I i T 1 Qt N1+ µ I T Qt. There we have the conditiona drift as 2ɛ B 2 µ I T Qt + N 21 ɛ µ I T Qt. F Qt

9 We note that there woud exist a finite set Q α in the compement of which the above drift woud be negative since the doant term in the above bound is the second term. By defining a new constant B 3 we can therefore write that the expected drift is B 3 2ɛ 2 N µ I T Qt. We note that the above drifts were defined for t k such that Qt k Q. By redefining the constants in the drift expressions we have that for a Qt k the above drift conditions hod. Then the rest of the proof proceeds as in 4 from equation 127 onwards and it can be shown that the average queue ength is finite.

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

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

Delay Analysis of Maximum Weight Scheduling in Wireless Ad Hoc Networks

Delay Analysis of Maximum Weight Scheduling in Wireless Ad Hoc Networks 1 Deay Anaysis o Maximum Weight Scheduing in Wireess Ad Hoc Networks ong Bao e, Krishna Jagannathan, and Eytan Modiano Abstract This paper studies deay properties o the weknown maximum weight scheduing

More information

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

A Survey on Delay-Aware Resource Control. for Wireless Systems Large Deviation Theory, Stochastic Lyapunov Drift and Distributed Stochastic Learning

A Survey on Delay-Aware Resource Control. for Wireless Systems Large Deviation Theory, Stochastic Lyapunov Drift and Distributed Stochastic Learning A Survey on Deay-Aware Resource Contro 1 for Wireess Systems Large Deviation Theory, Stochastic Lyapunov Drift and Distributed Stochastic Learning arxiv:1110.4535v1 [cs.pf] 20 Oct 2011 Ying Cui Vincent

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

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

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

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

Age-based Scheduling: Improving Data Freshness for Wireless Real-Time Traffic

Age-based Scheduling: Improving Data Freshness for Wireless Real-Time Traffic Age-based Scheduing: Improving Data Freshness for Wireess Rea-Time Traffic Ning Lu Department of CS Thompson Rivers University amoops, BC, Canada nu@tru.ca Bo Ji Department of CIS Tempe University Phiadephia,

More information

Power Control and Transmission Scheduling for Network Utility Maximization in Wireless Networks

Power Control and Transmission Scheduling for Network Utility Maximization in Wireless Networks ower Contro and Transmission Scheduing for Network Utiity Maximization in Wireess Networks Min Cao, Vivek Raghunathan, Stephen Hany, Vinod Sharma and. R. Kumar Abstract We consider a joint power contro

More information

IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 15, NO. 2, FEBRUARY

IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 15, NO. 2, FEBRUARY IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 5, NO. 2, FEBRUARY 206 857 Optima Energy and Data Routing in Networks With Energy Cooperation Berk Gurakan, Student Member, IEEE, OmurOze,Member, IEEE,

More information

Source and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems

Source and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems Source and Reay Matrices Optimization for Mutiuser Muti-Hop MIMO Reay Systems Yue Rong Department of Eectrica and Computer Engineering, Curtin University, Bentey, WA 6102, Austraia Abstract In this paper,

More information

Arbitrary Throughput Versus Complexity Tradeoffs in Wireless Networks using Graph Partitioning

Arbitrary Throughput Versus Complexity Tradeoffs in Wireless Networks using Graph Partitioning University of Pennsyvania SchoaryCommons Departmenta Papers (ESE) Department of Eectrica & Systems Engineering November 2006 Arbitrary Throughput Versus Compexity Tradeoffs in Wireess Networks using Graph

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

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

Optimal Control of Assembly Systems with Multiple Stages and Multiple Demand Classes 1

Optimal Control of Assembly Systems with Multiple Stages and Multiple Demand Classes 1 Optima Contro of Assemby Systems with Mutipe Stages and Mutipe Demand Casses Saif Benjaafar Mohsen EHafsi 2 Chung-Yee Lee 3 Weihua Zhou 3 Industria & Systems Engineering, Department of Mechanica Engineering,

More information

An Adaptive Opportunistic Routing Scheme for Wireless Ad-hoc Networks

An Adaptive Opportunistic Routing Scheme for Wireless Ad-hoc Networks An Adaptive Opportunistic Routing Scheme for Wireess Ad-hoc Networks A.A. Bhorkar, M. Naghshvar, T. Javidi, and B.D. Rao Department of Eectrica Engineering, University of Caifornia San Diego, CA, 9093

More information

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR 1 Maximizing Sum Rate and Minimizing MSE on Mutiuser Downink: Optimaity, Fast Agorithms and Equivaence via Max-min SIR Chee Wei Tan 1,2, Mung Chiang 2 and R. Srikant 3 1 Caifornia Institute of Technoogy,

More information

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing

More information

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c) A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,

More information

Asymptotic Properties of a Generalized Cross Entropy Optimization Algorithm

Asymptotic Properties of a Generalized Cross Entropy Optimization Algorithm 1 Asymptotic Properties of a Generaized Cross Entropy Optimization Agorithm Zijun Wu, Michae Koonko, Institute for Appied Stochastics and Operations Research, Caustha Technica University Abstract The discrete

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

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

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

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

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

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

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

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA) 1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using

More information

A simple reliability block diagram method for safety integrity verification

A simple reliability block diagram method for safety integrity verification Reiabiity Engineering and System Safety 92 (2007) 1267 1273 www.esevier.com/ocate/ress A simpe reiabiity bock diagram method for safety integrity verification Haitao Guo, Xianhui Yang epartment of Automation,

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

Exploring the Throughput Boundaries of Randomized Schedulers in Wireless Networks

Exploring the Throughput Boundaries of Randomized Schedulers in Wireless Networks Exporing the Throughput Boundaries of Randomized Scheduers in Wireess Networks Bin Li and Atia Eryimaz Abstract Randomization is a powerfu and pervasive strategy for deveoping efficient and practica transmission

More information

Age of Information: The Gamma Awakening

Age of Information: The Gamma Awakening Age of Information: The Gamma Awakening Eie Najm and Rajai Nasser LTHI, EPFL, Lausanne, Switzerand Emai: {eie.najm, rajai.nasser}@epf.ch arxiv:604.086v [cs.it] 5 Apr 06 Abstract Status update systems is

More information

Provisions estimation for portfolio of CDO in Gaussian financial environment

Provisions estimation for portfolio of CDO in Gaussian financial environment Technica report, IDE1123, October 27, 2011 Provisions estimation for portfoio of CDO in Gaussian financia environment Master s Thesis in Financia Mathematics Oeg Maximchuk and Yury Vokov Schoo of Information

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

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

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

Seung Jun Baek 1 and Joon-Sang Park Introduction. example, autonomous vehicles, remote surgery, and automated

Seung Jun Baek 1 and Joon-Sang Park Introduction. example, autonomous vehicles, remote surgery, and automated Hindawi Mathematica Probems in Engineering Voume 2017, Artice ID 4362652, 15 pages https://doi.org/10.1155/2017/4362652 Research Artice Deay-Optima Scheduing for Two-Hop Reay Networks with Randomy Varying

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

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

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

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

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

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

Fast Blind Recognition of Channel Codes

Fast Blind Recognition of Channel Codes Fast Bind Recognition of Channe Codes Reza Moosavi and Erik G. Larsson Linköping University Post Print N.B.: When citing this work, cite the origina artice. 213 IEEE. Persona use of this materia is permitted.

More information

Target Location Estimation in Wireless Sensor Networks Using Binary Data

Target Location Estimation in Wireless Sensor Networks Using Binary Data Target Location stimation in Wireess Sensor Networks Using Binary Data Ruixin Niu and Pramod K. Varshney Department of ectrica ngineering and Computer Science Link Ha Syracuse University Syracuse, NY 344

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

Approximate Bandwidth Allocation for Fixed-Priority-Scheduled Periodic Resources (WSU-CS Technical Report Version)

Approximate Bandwidth Allocation for Fixed-Priority-Scheduled Periodic Resources (WSU-CS Technical Report Version) Approximate Bandwidth Aocation for Fixed-Priority-Schedued Periodic Resources WSU-CS Technica Report Version) Farhana Dewan Nathan Fisher Abstract Recent research in compositiona rea-time systems has focused

More information

Discrete Applied Mathematics

Discrete Applied Mathematics Discrete Appied Mathematics 159 (2011) 812 825 Contents ists avaiabe at ScienceDirect Discrete Appied Mathematics journa homepage: www.esevier.com/ocate/dam A direct barter mode for course add/drop process

More information

Heavy-traffic Delay Optimality in Pull-based Load Balancing Systems: Necessary and Sufficient Conditions

Heavy-traffic Delay Optimality in Pull-based Load Balancing Systems: Necessary and Sufficient Conditions Heavy-traffic Deay Optimaity in Pu-based Load Baancing Systems: Necessary and Sufficient Conditions XINGYU ZHOU, The Ohio State University JIAN TAN, The Ohio State University NSS SHROFF, The Ohio State

More information

BALANCING REGULAR MATRIX PENCILS

BALANCING REGULAR MATRIX PENCILS BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

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

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,

More information

High Spectral Resolution Infrared Radiance Modeling Using Optimal Spectral Sampling (OSS) Method

High Spectral Resolution Infrared Radiance Modeling Using Optimal Spectral Sampling (OSS) Method High Spectra Resoution Infrared Radiance Modeing Using Optima Spectra Samping (OSS) Method J.-L. Moncet and G. Uymin Background Optima Spectra Samping (OSS) method is a fast and accurate monochromatic

More information

V.B The Cluster Expansion

V.B The Cluster Expansion V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f( q ) = exp ( βv( q )), which is obtained by summing over

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

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia

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

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

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

A Statistical Framework for Real-time Event Detection in Power Systems

A Statistical Framework for Real-time Event Detection in Power Systems 1 A Statistica Framework for Rea-time Event Detection in Power Systems Noan Uhrich, Tim Christman, Phiip Swisher, and Xichen Jiang Abstract A quickest change detection (QCD) agorithm is appied to the probem

More information

Algorithms to solve massively under-defined systems of multivariate quadratic equations

Algorithms to solve massively under-defined systems of multivariate quadratic equations Agorithms to sove massivey under-defined systems of mutivariate quadratic equations Yasufumi Hashimoto Abstract It is we known that the probem to sove a set of randomy chosen mutivariate quadratic equations

More information

Multiple Loop Self-Triggered Model Predictive Control for Network Scheduling and Control

Multiple Loop Self-Triggered Model Predictive Control for Network Scheduling and Control Mutipe Loop Sef-Triggered Mode Predictive Contro for Networ Scheduing and Contro Eri Henrisson, Danie E. Quevedo, Edwin G.W. Peters, Henri Sandberg, Kar Henri Johansson arxiv:5.38v [cs.sy] Feb 5 Abstract

More information

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................

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

Optimal Distributed Scheduling under Time-varying Conditions: A Fast-CSMA Algorithm with Applications

Optimal Distributed Scheduling under Time-varying Conditions: A Fast-CSMA Algorithm with Applications Optima Distributed Scheduing under Time-varying Conditions: A Fast-CSMA Agorithm with Appications Bin Li and Atia Eryimaz Abstract Recenty, ow-compexity and distributed Carrier Sense Mutipe Access (CSMA)-based

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

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

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS

More information

Delay Asymptotics with Retransmissions and Fixed Rate Codes over Erasure Channels

Delay Asymptotics with Retransmissions and Fixed Rate Codes over Erasure Channels Deay Asymptotics with Retransmissions and Fixed Rate Codes over Erasure Channes Jian Tan, Yang Yang, Ness B. Shroff, Hesham E Gama Department of Eectrica and Computer Engineering The Ohio State University,

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

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

Joint Congestion Control and Routing Optimization: An Efficient Second-Order Distributed Approach

Joint Congestion Control and Routing Optimization: An Efficient Second-Order Distributed Approach IEEE/ACM TRANSACTIONS ON NETWORKING Joint Congestion Contro and Routing Optimization: An Efficient Second-Order Distributed Approach Jia Liu, Member, IEEE, Ness B. Shroff, Feow, IEEE, Cathy H. Xia, and

More information

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University Turbo Codes Coding and Communication Laboratory Dept. of Eectrica Engineering, Nationa Chung Hsing University Turbo codes 1 Chapter 12: Turbo Codes 1. Introduction 2. Turbo code encoder 3. Design of intereaver

More information

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential Lecture 6 Povh Krane Enge Wiiams Properties of -nuceon potentia 16.1 4.4 3.6 9.9 Meson Theory of Nucear potentia 4.5 3.11 9.10 I recommend Eisberg and Resnik notes as distributed Probems, Lecture 6 1 Consider

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

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

A Graphical Approach for Solving Single Machine Scheduling Problems Approximately

A Graphical Approach for Solving Single Machine Scheduling Problems Approximately A Graphica Approach for Soving Singe Machine Scheduing Probems Approximatey Evgeny R Gafarov Aexandre Dogui Aexander A Lazarev Frank Werner Institute of Contro Sciences of the Russian Academy of Sciences,

More information

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP

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

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. 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

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

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized

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

Methods for Ordinary Differential Equations. Jacob White

Methods for Ordinary Differential Equations. Jacob White Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa

More information

Tracking Control of Multiple Mobile Robots

Tracking Control of Multiple Mobile Robots Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart

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

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

Inductive Bias: How to generalize on novel data. CS Inductive Bias 1

Inductive Bias: How to generalize on novel data. CS Inductive Bias 1 Inductive Bias: How to generaize on nove data CS 478 - Inductive Bias 1 Overfitting Noise vs. Exceptions CS 478 - Inductive Bias 2 Non-Linear Tasks Linear Regression wi not generaize we to the task beow

More information

Spring Gravity Compensation Using the Noncircular Pulley and Cable For the Less-Spring Design

Spring Gravity Compensation Using the Noncircular Pulley and Cable For the Less-Spring Design The 14th IFToMM Word Congress, Taipei, Taiwan, October 25-30, 2015 DOI Number: 10.6567/IFToMM.14TH.WC.PS3.010 Spring Gravity Compensation Using the Noncircuar Puey and Cabe For the Less-Spring Design M.C.

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

V.B The Cluster Expansion

V.B The Cluster Expansion V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f(q ) = exp ( βv( q )) 1, which is obtained by summing over

More information

On Stochastic Sensor Network Scheduling for Multiple Processes

On Stochastic Sensor Network Scheduling for Multiple Processes On Stochastic Sensor Network Scheduing for Mutipe Processes Duo Han, Junfeng Wu, Yiin Mo, Lihua Xie arxiv:6.08v [cs.sy 8 Dec 06 Abstract We consider the probem of mutipe sensor scheduing for remote state

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

Fractional Power Control for Decentralized Wireless Networks

Fractional Power Control for Decentralized Wireless Networks Fractiona Power Contro for Decentraized Wireess Networks Nihar Jinda, Steven Weber, Jeffrey G. Andrews Abstract We consider a new approach to power contro in decentraized wireess networks, termed fractiona

More information