Stable Scheduling Policies for Maximizing Throughput in Generalized Constrained Queueing Systems

Size: px
Start display at page:

Download "Stable Scheduling Policies for Maximizing Throughput in Generalized Constrained Queueing Systems"

Transcription

1 1 Sable Scheduling Policies for Maximizing Throughpu in Generalized Consrained Queueing Sysems Prasanna Chaporar, Suden Member, IEEE, Saswai Sarar, Member, IEEE Absrac We consider a class of queueing newors referred o as generalized consrained queueing newors which form he basis of several differen communicaion newors and informaion sysems. These newors consis of a collecion of queues such ha only cerain ses of queues can be concurrenly served. Whenever a queue is served, he sysem receives a cerain reward. Differen rewards are obained for serving differen queues, and furhermore, he reward obained for serving a queue depends on he se of concurrenly served queues. We demonsrae ha he dependence of he rewards on he schedules aler fundamenal relaions beween performance merics lie hroughpu and sabiliy. Specifically, maximizing he hroughpu is no longer equivalen o maximizing he sabiliy region; we herefore need o maximize one subjec o cerain consrains on he oher. Since sabiliy is criical for bounding pace delays and buffer overflow, we focus on maximizing he hroughpu subjec o sabilizing he sysem. We design provably opimal scheduling sraegies ha aain his goal by scheduling he queues for service based on he queue lenghs and he rewards provided by differen selecions. The proposed scheduling sraegies are however compuaionally complex. We subsequenly develop echniques o reduce he complexiy and ye aain he same hroughpu and sabiliy region. We demonsrae ha our framewor is general enough o accommodae random rewards and random scheduling consrains. Index Terms Consrained queueing newors, sabiliy, hroughpu, opimizaion, wireless, mulicas, randomized algorihms I. INTRODUCTION Consrained queueing newors have been exensively used o model several sysems of pracical ineres including wireless newors 35, 34, 25, 27, inpu queued swiches 23 and daabase sysems 34. A consrained queueing newor is a collecion of queues such ha only cerain ses of queues can be concurrenly served; hese schedulable ses depend on he underlying sysem. Whenever a queue is served, he sysem receives a cerain reward. In such sysems, queues need o be seleced for service such ha i) he oal reward earned by he sysem per uni ime hroughpu ) is maximized and ii) each queue is served ofen enough such ha he mean queue lengh in each queue is bounded sysem sabiliy ). The wo goals urn ou o be equivalen if he service of each queue i.e., he ransmission of each pace) feches he same reward. The performances of such newors are now reasonably well-undersood owing o several seminal conribuions 1, 3, 4, 6, 24, 25, 26, 35. We now invesigae consrained queueing newors where differen rewards are obained for ransmiing paces from differen queues, and furhermore, he reward obained for serving a queue depends on he se of concurrenly served queues. Such generalized consrained queueing newors form he basis of several communicaion and informaion sysems of pracical ineres, bu have no received adequae aenion in he research communiy. We firs provide examples of such sysems, and subsequenly demonsrae ha new resource allocaion goals and echniques are required for capuring he rade-off beween differen performance merics in hese sysems. Firs, consider one-o-many communicaions in wireless newors. Here, a sender may wish o ransmi is paces o muliple receivers in is communicaion range. Due o he broadcas propery of he wireless ransmission, a single ransmission may reach all hese receivers. Here, each sender consiues a queue, and he reward aained by a ransmission is he number of receivers who successfully receive i. Since differen mulicas groups have differen P. Chaporar is wih Indian Insiue of Technology, Mumbai, India. His address is chaporar@ee.iib.ac.in. S. Sarar is wih he Deparmen of Elecrical and Sysems Engineering a Universiy of Pennsylvania. Her address is swai@seas.upenn.edu. This wor was suppored by he Naional Science Foundaion under grans ANI , NCR and CNS

2 2 01 R3 Transmission Range of S 2 R4 01 S1 R 01 5 Transmission Range of S 1 R 2 01 R S2 R 6 Fig. 1. Figure shows an example o demonsrae he applicaion of generalized consrained queueing newors in one-o-many communicaion in wireless newors. There are wo senders S 1, S 2 and 6 receivers R 1 o R 6. The dashed circles indicae he communicaion ranges of he senders. A single ransmission from S 1 can reach all is receivers, R 1,..., R 5. Here, R 6 is S 2 s receiver. Each sender corresponds o a queue. Here, L = { l 1 = 0 0, l 2 = 1 0 l 3 = 0 1, l 4 = 1 1}. Here, r l 1) = 0, {1,2}, r 1 l 2) = 5, r 2 l 2) = 0, r 1 l 3) = 0, r 2 l 3) = 1, r 1 l 4) = 3, r 2 l 4) = 1. T1 T2 T3 T4 U1 U2 U3 Fig. 2. Figure shows a daabase sysem wih four ables T1,..., T4 ha are accessed by hree applicaions U1, U2 and U3. The arrows indicae he ables each applicaion updaes. When here are concurren requess for updaes in he same able, he reques from an applicaion wih he lowes id is honored. Noe ha if all hree applicaions ry o simulaneously updae he daabase, hen U1, U2 and U3 achieve rewards 3, 1 and 0 respecively. If only U2 and U3 ry o simulaneously updae he daabase, hen hey achieve rewards 2 and 1 respecively. number of receivers, he reward aained by serving differen queues will be differen. Furhermore, wheher a receiver can successfully decode a ransmission depends on oher ransmissions in is neighborhood. Thus, he reward associaed wih each ransmission depends on he se of queues served concurrenly. For example in Figure 1 when S 2 is ransmiing o R 6, R 1 and R 2 canno receive a ransmission from S 1 as boh he ransmissions will collide a hese receivers. Hence, S 1 receives a reward of 5 when S 1 alone is served, and i receives a reward of 3 when S 1 and S 2 are served ogeher. Thus, he reward for S 1 depends on he se of queues served. Now, consider one-o-one communicaion in wireless newors. Success of each ransmission depends upon he inerference due o concurren ransmissions in he newor and he channel sae. Le he reward for each ransmission be 1 if he ransmission is successful. Thus, differen ransmissions aain differen rewards depending on he se of queues served. Furhermore, here, he same selecion of sessions may generae differen rewards a differen imes as he inerferences randomly change due o fading - rewards may herefore be random. Nex, in many daabase sysems, a single updae operaion from an applicaion involves updaes in many ables. Here, each applicaion consiues a queue, and he reward aained by an updae operaion is he number of ables ha are successfully updaed. Since differen applicaions require o updae differen number of ables, rewards received by serving differen queues will be differen. Moreover, if many applicaions ry o updae he same able, hen only one of hem can do so, as he access o hese ables is conrolled o avoid inconsisencies due o concurren updaes. Thus, he reward for a queue depends on he se of queues served. We demonsrae his using a specific applicaion in Figure 2. Our conribuion is o provide a mahemaical framewor for modeling and opimizing ey performance aribues in generalized consrained queueing newors. Firs, we define appropriae performance merics Secion II). Nex, we demonsrae ha he fundamenal relaions beween performance merics such as hroughpu and sabiliy change due o he dependence of he rewards on he se of queues served Secion III). Specifically, maximizing he hroughpu is no longer equivalen o maximizing he sabiliy region; we herefore need o maximize one subjec o cerain consrains on he oher. Since sabiliy is criical for bounding pace delays and buffer overflow, we focus on maximizing he hroughpu subjec o sabilizing he sysem. We design provably opimal scheduling sraegies ha aain his goal by scheduling he queues for service based on he queue lenghs and he rewards provided by differen selecions Secion IV). These scheduling sraegies are however compuaionally complex. We nex develop a framewor o reduce he compuaional complexiy and ye aain he opimum performance Secion V).

3 3 Finally, we consider some possible generalizaions Secion VI) and describe he relaed wor Secion VII). II. SYSTEM MODEL We consider a queueing newor wih n queues. We assume ha ime is sloed. In each queue {1,...,n} paces arrive as per arrival process {Λ )} =1, where Λ ) is he number of arrivals in queue during slo. Arrivals for he same session in differen slos are independen and idenically disribued. The arrival processes for differen sessions are independen bu no idenically disribued. We assume ha Λ ) A max in each slo def and for any. Le λ = EΛ ) and λ = λ 1,...,λ n ) denoe he arrival rae vecor. Each pace can be served in a mos one slo, and i depars he sysem a he end of he slo in which i is served. This assumpion has been moivaed by he fac ha in wireless newors muliple ransmissions of he same pace consume addiional energy and increase he inerference for oher ransmissions. We denoe by Q ) he queue lengh of he h queue a he beginning of slo. Also, Q) = Q 1 ) Q n ). A queue can only be served if i has a pace o ransmi, and in each slo in which i is served i ransmis one pace. The indicaor l ) = 1 if he h queue is served in slo, and is 0 oherwise. The vecor l) = l 1 ) l n ) denoes he service vecor in slo. The sysem consrains may prohibi simulaneous service of cerain queues. Thus, all 2 n n-dimensional binary vecors may no consiue a valid service vecor. Le L = { l 1,..., l m } denoe he se of all valid service vecors, and l i denoe he h elemen of l i L. Clearly m 2 n. For example, Figure 1 elucidaes a consrained queueing newor wih n = 2 and m = 4. Now, if he sysem has an addiional consrain ha all he receivers { should receive every pace, hen boh S 1 and S 2 can no be served concurrenly. Thus, in his case, L = l1 = 0 0, l 2 = 1 0, } l 3 = 0 1 and m = 3. We assume he following abou L. If l L, hen every l 1 l also belongs o L, where he inequaliy is elemen-wise. In oher words, if a cerain se of queues can be served simulaneously, hen any subse of hese queues can also be served simulaneously. Noe ha his assumpion holds in wireless newors. For each l i L and queue lengh vecor Q, we define an n-dimensional vecor l i Q) as follows. The h componen of l i Q) equals l i if Q > 0, and is 0 oherwise. Clearly, for each l i L and Q, l i Q) L. The sysem receives a reward for serving each queue, and he reward obained for serving he h queue in slo, r l)), is a funcion of he service vecor l) in slo, for each. We assume ha r l)) G < for each. We iniially assume ha he reward for each queue is a deerminisic funcion of he service vecor, and laer generalize o allow he reward o randomly depend on he service vecor Secion VI). Refer o Figure 1 for some example rewards. We assume he following properies of he reward funcion. Firs, if l = 0 hen r l) = 0. Thus, if a queue is no served hen i does no receive any reward. Nex, for any l 1, l 2 L if l 1 l 2 and l 1 > 0, r l 1 ) r l 2 ). Thus, r l Q)) r l) for any Q, l L and such ha Q > 0. We jusify his assumpion in conex of one of he applicaion scenarios, wireless newors. In wireless newors, when fewer queues ransmi, he inerference is less in he sysem and herefore, usually, he queues ha ransmi receive higher reward. If his is no he case, e.g., when he probabiliy of success increases wih increase in inerference due o he use of sophisicaed decoding sraegies, hen if an empy queue is seleced, i can ransmi a signal so as o ensure ha oher queues do no receive less reward because i is empy. This may increase he overall energy consumpion, bu our focus here is o maximize he hroughpu. Join minimizaion of he energy consumpion and maximizaion of he hroughpu consiss of ineresing opics for fuure research. The assumpion can also be similarly jusified for daabase sysems. Nex, we presen some imporan definiions. Definiion 1 Scheduling Policy): A scheduling policy decides he service vecor l ) in each slo 1 such ha l ) L and l i ) = 0 if Q i) = 0. This class includes offline policies ha decide heir service vecors based on he nowledge of pace arrivals in each pas, presen and even fuure slos. Since r l u)) = 0 if l u) = 0, and l u) = 0 if Q u) = 0, r l u)) = 0 if Q u) = 0. Thus, irrespecive of he scheduling policy, no queue receives any reward in a slo in which i is empy. Transmission of a signal from an empy queue is no considered service for he empy queue.

4 4 Definiion 2 Throughpu): For an arrival rae vecor λ, he hroughpu under a scheduling policy, Ω λ), is he reward i receives per uni ime. Mahemaically, Ω λ) = lim inf Since r l u)) = 0 if l u) = 0, and l u) {0,1}, Ω 1 λ) = lim inf 1 u=1 u=1 r l u)). r l u))l u). 1) Noe ha if he reward r l) is he number of receivers of session ha receive a pace when he service vecor is l, he hroughpu under is he sum, over all receivers, of he number of paces each receiver receives per uni ime. This is consisen wih he usual definiion of hroughpu in a communicaion newor. Definiion 3 Loss): The loss under a scheduling policy a any slo is he difference beween he sum of he maximum possible rewards of he queues i serves a and he reward i obains a. The loss under a scheduling policy, L λ), is is oal loss per uni ime. Mahemaically, L 1 λ) = lim inf G r ) l u)) l u). u=1 In a communicaion newor, usually, he loss experienced by a receiver denoes he number of paces ransmied by is source ha i does no receive per uni ime, and he newor loss denoes he sum of he losses of all receivers. Again, if he reward r l) is he number of receivers of session ha receive a pace when he service vecor is l, hen he formal definiion of loss in Definiion 3 has he same connoaion as above. Definiion 4 Sysem Sabiliy): The queueing P sysem is said o be sable if he ime average of queue lenghs is finie for each queue, i.e., lim sup Qiu) u=1 < wih probabiliy w.p.) 1 for each i. A scheduling policy ha sabilizes he sysem is called a sable scheduling policy. The sabiliy region of a scheduling policy is he se of arrival rae vecors for which he sysem is sable under he policy. The sabiliy region of he sysem Θ is he union of he sabiliy regions of all scheduling policies. A scheduling policy whose sabiliy region equals Θ is said o maximize he sabiliy region. Le C denoe he convex hull of he vecors in L and C denoe he inerior of C. In heir seminal wor, Tassiulas e al Theorem 3.2, 35) showed ha C Θ C. Definiion 5 Sabilizable Arrival Rae Vecor): We denoe he arrival rae vecor λ as sabilizable if λ C. Definiion 6 Throughpu Opimaliy): A sable scheduling policy is said o be hroughpu opimal if w.p. 1 i aains he maximum hroughpu among all he sable scheduling policies. We denoe he hroughpu aained by such a policy for arrival rae vecor vecor λ Θ by Ω max λ). Definiion 7 ǫ-throughpu Opimaliy): A scheduling policy is said o be ǫ-hroughpu opimal for a ǫ > 0 if a) i is sable, and b) Ω λ) Ω max λ) ǫ w.p. 1. In he nex secion, we show ha in generalized consrained queueing newors maximizing he sabiliy region is no equivalen o maximizing he hroughpu. Since sabiliy is imperaive for guaraneeing bounded delay and for limiing pace drop due o buffer overflow, we aim o maximize he hroughpu subjec o sabilizing he sysem. Specifically, our goal is o design ǫ-hroughpu opimal policies. We now invesigae he relaion beween he hroughpu and he loss. Now, L 1 λ) = G lim inf l u) Ω λ). Noe ha if a sysem is sable under policy, lim 1 u=1 u=1 l u) = λ w.p. 1. Thus, L λ) = n λ G Ω λ) w.p. 1. Thus, if λ is in he sabiliy region of policies 1, 2, Ω 1 λ) + L 1 λ) = Ω 2 λ) + L 2 λ) w.p. 1. Thus, for any sabilizable arrival rae vecor λ, a hroughpu opimal policy mus also minimize he loss, and an ǫ-hroughpu opimal policy aains a loss which is a mos ǫ more han he loss of any sable policy. Thus, we focus on obaining ǫ-hroughpu opimal policies.

5 5 III. RELATION BETWEEN THROUGHPUT AND STABILITY Firs, we examine wha decisions policies are liely o mae if hey wan o maximize only he sabiliy region, or if hey wan o maximize only he hroughpu. A policy ha aims o maximize he sabiliy region serves as many paces as possible in a slo while giving prioriy o longer queues. If he policy aims o maximize he hroughpu, hen i may wai and ransmi only when he reward is high so ha each pace feches he maximum possible reward. Thus, he conrol decisions for maximizing he sabiliy region and for maximizing he hroughpu are no equivalen. Using an example ha is moivaed by one-o-many communicaion in wireless newors Figure 1), we nex demonsrae ha a policy ha maximizes he sabiliy region does no maximize he hroughpu. Example 1: Consider he sysem shown in Figure 1. Le λ = 1/2 ǫ,1/2 ǫ), where ǫ is a small posiive real number. Now, consider a policy ha serves each queue whenever i is non-empy. Thus, if only S 1 S 2, resp.) is non-empy, hen will selec service vecor l 2 l 3, resp.) and achieve a reward of 5 1, resp.). If boh queues are non-empy in a slo, hen will selec l 4 and achieve a reward of 4. Clearly, maximizes he sabiliy region. Now, he service process for S 1 is independen of ha for S 2. Using Lile s law, he fracion of slos in which S 1 S 2, resp.) is non-empy and S 2 S 1. resp.) is empy is 1/4 ǫ 2 1/4 ǫ 2, resp.), and he fracion of slos in which boh queues are non-empy is 1/2 ǫ) 2. Thus, Ω b λ) = 5 + 1)1/4 ǫ 2 ) + 41/2 ǫ) 2 10/4. Now, consider a policy ha serves only S 1 when S 1 is non-empy, and serves only S 2 if S 1 is empy and S 2 is non-empy. Noe ha is sable as λ 1 + λ 2 = 1. Thus, whenever S 1 S 2, resp.) is served, he service vecor is l 2 l 3, resp.) and he reward is 5 1, resp.). Since he queues are sable, S 1 and S 2 are served in 1/2 ǫ fracion of slos each. Thus, Ω = 5 + 1)1/2 ǫ) 3. Thus, Ω b λ) < Ω λ). Noe ha in Example 1 always ransmis he maximum number of paces in each slo and also chooses he se of queues whose sum of queue lenghs is he maximum. Tassiulas e al 35 showed ha a policy ha saisfies he laer propery maximizes he sabiliy region in arbirary consrained queueing newors. Bu, Example 1 shows ha does no maximize he hroughpu. This is because does no consider he reward srucure in deciding he service vecor. So, he policies designed o maximize he sabiliy region of he consrained queueing sysem e.g. see 1, 5, 19, 34, 35) need no maximize he hroughpu. Thus, we need alernae mechanism o design hroughpu opimal policies. Now, we consider wo policies, 1, 2, ha see o maximize he reward in a greedy fashion. 1 serves each queue only when he queue can obain is maximum possible reward, and 2 selecs in each slo he service vecor ha aains he maximum possible reward among all valid service vecors in he slo. Simply pu, 1 maximizes he reward per pace, and 2 greedily maximizes he reward in each slo. We show ha 1 does no sabilize he sysem even when he arrival rae vecor is sabilizable, and 2 does no aain he maximum hroughpu among all sable policies. Example 2: Consider he sysem shown in Figure 1. Le λ = 3/4,1/2). Clearly, λ C and policy in Example 1 sabilizes he sysem. Noe ha 1 will never concurrenly serve boh queues. Hence, he sum of he service raes provided o he wo queues is a mos 1. Thus, since λ 1 + λ 2 > 1, 1 does no sabilize he sysem. Noe ha 1 maximizes he reward per pace while serving queues a raes smaller han heir arrival raes and hereby compromises sabiliy. Example 3: Consider he sysem shown in Figure 1 wih he difference ha r 2 l 3 ) = r 2 l 4 ) = 3. Le λ = 1/4,1/4). Noe ha for he above rewards, 2 selecs he same service vecors as. Thus, 2 sabilizes he sysem. Now, he service process for S 1 is independen of ha for S 2. Using Lile s law, he fracion of slos in which S 1 S 2, resp.) is non-empy and S 2 S 1. resp.) is empy is 1/4)3/4) 1/4)3/4)), resp.), and he fracion of slos in which boh queues are non-empy is 1/4) 2. Thus, Ω b λ) = 5 + 3)3/16) + 61/16) = 15/8. Now, consider described in Example 1. Again, is sable as λ 1 + λ 2 < 1. Thus, whenever S 1 S 2, resp.) is served, he service vecor is l 2 l 3, resp.) and he reward is 5 3, resp.). Since he queues are sable, S 1 and S 2 are served in 1/4 fracion of slos each. Thus, Ω = 5 + 3)1/4) = 2. Thus, 2 does no aain he maximum hroughpu among all sable policies. The limiaion of 2 is ha i myopically bases is decision in a slo solely on he aggregae reward in he slo. Thus, even when i is possible o wai and serve queues in muually disjoin slos and achieve a higher reward per pace, 2 serves he queues in he same slo. The examples demonsrae ha 1) a policy ha maximizes he sabiliy region need no maximize he hroughpu

6 6 2) myopically maximizing he reward in each slo or he reward per pace may no maximize he hroughpu or sabilize he sysem 3) he opimal policy should wai jus long enough so as o achieve he highes possible reward per pace while serving each queue a a rae higher han is arrival rae. IV. OPTIMAL POLICIES In his secion, we propose wo policies and prove ha hey are ǫ-hroughpu opimal for every sabilizable arrival rae vecor λ and ǫ > 0. A. Linear program based opimal policy ) The scheduling policy selecs l i L w.p. w i in every slo. If l i is chosen in slo, hen l ) = l i Q) ), i.e., he h queue ransmis a pace if l i = 1 and Q ) > 0. Recall ha l ) is he indicaor vecor for he se of queues served by in slo. Le selec l i in a slo. Then l ) l i. The inequaliy is sric only when some queues in l i are empy in, and hen, as discussed in Secion II, r l )) r l i ) for each for which Q ) > 0. The probabiliy disribuion w = w 1 w m is compued using he following linear program LP λ,δ). Here, δ is a parameer. LP λ,δ) :- Maximize: U λ,δ) = i w il i r l i ) Subjec o: 1) m i=1 w i = 1 and w i 0 for every i 2) m i=1 w il i = λ + δ for every. Consrain 1) ensures ha w is a valid probabiliy disribuion. When δ > 0, consrain 2) ensures ha each queue is seleced for service a a rae higher han he arrival rae in he queue. Thus, consrain 2) ensures sabiliy. Noe ha w and hence depend on λ and he chosen δ. We indicae his dependence by using he noaions w λ,δ) and λ,δ). Now, alhough LP λ,δ) is well-defined, i need no have any feasible soluion, for arbirary λ R n and δ R. Theorem 1 shows ha for all sabilizable λ and sufficienly small posiive δ, LP λ,δ) is feasible and λ,δ) is ǫ-hroughpu opimal. Noe ha allowing arbirary λ R n and δ R in LP λ,δ) simplifies he proof for Theorem 1. Theorem 1: Le λ be any sabilizable arrival rae vecor. Then, for every ǫ > 0 here exiss a δ such ha for every δ 0, δ), LP λ,δ) is feasible and λ,δ) is ǫ-hroughpu opimal. Furhermore, Ω λ,δ) U λ,δ) δ G Ω max λ) ǫ w.p. 1. 2) We prove Theorem 1 in he appendix. Finally, he sabiliy region can be maximized using arbirary feasible soluions of LP λ,δ) 5. Specifically, if λ,δ) selecs he service vecors as per any probabiliy disribuion ha consiues a feasible soluion of LP λ,δ) for any posiive δ, i sabilizes he sysem provided λ is sabilizable 5. Bu, for aaining he maximum hroughpu among all sable policies, an opimal soluion of LP λ,δ) mus be used. Specifically, for any sabilizable λ and ǫ > 0, λ,δ) is ǫ-opimal for any δ 0,min{δ max λ),ǫ/ n G }, where δ max λ) is he maximum value of δ for which LP λ,δ) has a feasible soluion follows from Theorem 1 and Lemma 5 in he appendix). B. Queue lengh based opimal policy O ) The policy λ,δ) requires he nowledge of λ in order o obain he opimal w λ,δ). The sysem may no however now λ. We now design a policy O ha aains he maximum hroughpu among all sable policies and sabilizes he sysem for any sabilizable λ wihou nowing λ. Recall ha an opimal policy should wai as long as possible o achieve he highes possible reward per pace wihou violaing sysem sabiliy Secion III). Now, λ,δ) uses he nowledge of λ o ensure he above, whereas O ensures he above by using only he value of Q).

7 7 We now describe O. In slo, O selecs he service vecor l O ) such ha { } l O ) = arg max Q ) V G r ) l)) l, 3) l L, l= l Q)) where V is a consan. Noe ha he consrain l = ) l Q) implies ha l = 0 if Q ) = 0 and l {0,1} oherwise. Theorem 2: Le λ be any sabilizable arrival rae vecor. Then, for every V 0, O sabilizes he sysem. Moreover, for every ǫ > 0, here exiss V such ha for every V V, O is ǫ-hroughpu opimal. The above resul implies ha any sable off-line policy ha aes ransmission decisions based on he nowledge of pas, presen and fuure arrivals can no aain hroughpu significanly more han Ω O λ) for every sabilizable λ. This holds even hough O aes ransmission decisions based only on he curren queue lenghs. Now, we describe he inuiion behind his resul. Le W Q, def l) = W 1 Q, l) def = Q l V G r l))l ). 4) Q l. Noe ha inuiively G r l) is he loss of reward of he h queue when service vecor l is used. Thus, in each slo, O selecs he service vecor l ha maximizes he do produc, W Q), l), of l and he difference beween he queue lengh vecor Q) and a scaled loss vecor associaed wih l. Noe ha a policy ) ha selecs he service vecor l ha maximizes he do produc, W 1 Q), l), of l and he queue lengh vecor Q) sabilizes he sysem for every sabilizable λ 1, 18, 35. This is because under he queue lengh process has a negaive drif when n Q ) is sufficienly large for every sabilizable λ. When n Q ) >> V n G, W Q), l) W 1 Q), l) for every l L, and herefore, O and selec similar service vecors. Thus, inuiively, for every sabilizable λ, he queue lengh process under O should also have a negaive drif when n Q ) is sufficienly large. Hence, O also sabilizes he sysem for any sabilizable λ. We have however shown ha all sable policies do no aain equal hroughpu Example 1). So, i is no obvious ha O maximizes he hroughpu among all policies ha sabilize he sysem; we now provide he inuiion behind why his is he case. Noe ha when he queue lenghs are small, high hroughpu can be aained wihou violaing sabiliy by serving he queues only when hey receive high rewards. On he oher hand, sabiliy can be ensured by selecing he queues wih higher queue lenghs and by serving a large number of paces when he queue lenghs are large. We now demonsrae ha O follows boh he above principles. For simpliciy, assume ha V,G,r l) are inegers for all, l L. Now, when Q ) < V, Q ) V G r ) l) 0 only if r l) = G. Then, since O maximizes W Q), l), i will serve he h queue only if he maximum possible reward is achievable. Now, if Q {V,...,2V 1}, hen Q ) V G r ) l) 0 only if r l) G 1. Thus, O will serve he h queue only if he achievable reward is greaer han or equal o G 1. Similarly, if Q ) {G u)v,...,g u + 1)V 1}, hen O will serve he h queue only if he achievable reward is greaer han or equal o u. Summarily, O aains he maximum possible reward for every pace while mainaining sabiliy by dynamically selecing he service vecors based on he queue lenghs. Thus, O aains he maximum hroughpu among all sable policies. Now, we prove Theorem 2 using a combinaion of opimizaion and Lyapunov heories. Neely e al. 27 proposed his proof echnique in a differen conex. Proof: Consider a sabilizable λ. For any policy, Q + 1) = Q ) + Λ ) l ). 5) Now for O, { Q)} 1 is an irreducible, aperiodic and counable Marov chain. Now, consider he Lyapunov funcion f Q)) = Q )) 2. 6)

8 8 Le, M def = na 2 max + 1). From 5) and 6), i follows ha f Q + 1)) f Q)) M + 2Q )Λ ) 2Q )l O Thus, E f Q + 1)) f Q)) Q) M + = M + 2E 2Q )λ E ). 2Q )l O ) Q) 2Q )λ 2E W Q), l O )) Q) Q )l ) V G r ) l ))l ) Q) 2V E G r l O ))l O ) Q). 8) Now, since λ is sabilizable, we can obain small enough posiive δ such ha λ,δ) is ǫ/2-hroughpu opimal Theorem 1). We consider λ,δ) for such a δ. Here, l λ,δ) ) is he service vecor λ,δ) would have used a if i had a queue lengh vecor of Q) a. From definiion of O equaion 3)), for every and, E Thus, from 8), for every δ W Q), l O )) Q) E E W Q), l )) Q). f Q + 1)) f Q)) Q) M + 2Q )λ 2E W Q), l λ,δ) )) Q) 2V E G r l O ))l O ) Q). 9) Now λ,δ) chooses each service vecor l i w.p. w i λ,δ) independen of he queue lenghs, and subsequenly serves only hose queues ha are included in he seleced service vecor and are also non-empy. Thus, if Q ) > 0, E ) Q) m = l i w i λ,δ) = λ + δ. 10) l λ,δ) i=1 In addiion, recall ha if selecs l i, and if he h queue is nonempy i receives a reward of a leas r l i ). Thus, if Q ) > 0, E G r ) l λ,δ) )) l λ,δ) ) Q) m w i λ,δ) G r ) l i ) If Q ) = 0, i=1 = G λ + δ) l i m w i λ,δ)r l i )l i. i=1 E G r ) l λ,δ) )) l λ,δ) ) Q) = 0. 7)

9 9 Now, m i=1 w i λ,δ)r l i )l i G λ + δ) since r l i ) G and m i=1 w i λ,δ)l i = λ + δ. Thus, E G r l λ,δ) )) m G λ + δ) i=1 ) l λ,δ) ) Q) From 4), 10), and 11), i follows ha E W Q), l λ,δ) )) Q) Q ) V G )λ + δ) + V U λ,δ). Hence, from 9) E f Q + 1)) f Q)) Q) M 2δQ ) + 2V G λ + δ) 2V U λ,δ) w i λ,δ)r l i )l i. 11) 2V E G r l O ))l O ) Q). 12) 1) Sabiliy of O : From 12), since 0 r l) G,l 0 for all, l L, i follows ha for every sabilizable λ and every non-negaive V E f Q + 1)) f Q)) Q) M 2δQ ) + 2 V G λ + δ). Le A = { Q : n =0 Q V δ n =0 G λ + δ) + M+1 2δ E f Q + 1)) f Q)) Q) }. Then, { < for all Q), 1 if Q) / A. Thus, since A is finie, by Foser s Theorem Theorem in 20), { Q)} 1 is posiive recurren, and for each queue he expeced queue lengh under is saionary disribuion is finie. Thus, he sysem is sable under O. 2) ǫ-throughpu Opimaliy of O : Taing expecaion on boh sides of 12) wih respec o he saionary disribuion of { Q)} 1, we obain E f Q + 1)) f Q)) M 2δEQ ) + 2V G λ + δ) 2V U λ,δ) 2V E G r l O ))l O ). 13) Now, u=1 l O ) is he number of deparures from queue in 0,) under O. Since he queue lengh process { Q)} 1 under O is a posiive recurren Marov chain, for every 1 u El O ) = lim l O u u v) = λ w.p. 1, 14) v=1

10 and 1 = lim u u Er l O ))l O ) u v=1 r l O v))l O v) w.p. 1 = Ω O λ) from 1)). 15) Moreover, since he expecaions are wih respec o he saionary disribuion of Q), i follows ha From 13), 14), 15) and 16), i follows ha Ω O λ) Ef Q + 1)) = Ef Q)). 16) U λ,δ) δ G M 2V Ω max λ) ǫ 2 M from Theorem 1) 2V 17) Ω max λ) ǫ if V M/ǫ. The resul follows. Finally, we commen on he role of he parameer V in deermining he hroughpu of O. From 17), i can be seen ha if V < M/ǫ = na 2 max + 1)/ǫ, hen no hroughpu guaranee can be provided for O. Noe ha A max deermines he bursiness of he arrival process. Thus, he minimum required value of V is higher for more bursy arrival processes. 10 C. Compuaion ime for and O In he wors case, cardinaliy of L can be 2 n as i may conain all n-dimensional binary vecors. Then, λ,δ) can be compued by solving a linear program wih O2 n ) variables and On) consrains. Thus, he ime and he memory required o compue λ,δ) is O2 n ) in he wors case. Under O, we need o find a l L ha maximizes W Q), l) for every. Since L is O2 n ), he ime required o compue he opimal service vecor in each slo is also O2 n ) unless some addiional srucure on he queueing sysem is assumed. We nex propose wo opimal policies which require polynomial compuaion ime in every slo. V. COMPUTATIONALLY SIMPLE OPTIMAL POLICIES We provide a general framewor for designing compuaionally simple policies for maximizing he hroughpu subjec o aaining sabiliy by considering he noion of inaccurae scheduling Subsecion V-A). We subsequenly uilize his framewor o design wo compuaionally simple policies for maximizing he hroughpu subjec o sabilizing he sysem Subsecions V-B,V-C). Finally, we discuss how hese policies can be implemened using disribued compuaion Subsecion V-D). A. Inaccurae scheduling for maximizing he hroughpu subjec o sabilizing he sysem We firs describe a class of scheduling policies referred o as inaccurae scheduling. Noe ha he noion of inaccurae scheduling has earlier been proposed for designing compuaionally simple policies for maximizing he sabiliy region 23, 32, 34. Our conribuion here is o generalize his noion o aain he goal of maximizing he hroughpu subjec o sabilizing he sysem while using simple compuaions. We consider policies for which he sae { Y ) = Q), l ))} 1 consiues an irreducible, aperiodic and counable Marov chain. This assumpion holds when l ) is compued ieraively based on Q) and l 1). Noe ha hen { Q)} 1 may no be a Marov process. Definiion 8 γ-inaccurae Policy): A policy γ is called γ-inaccurae if in each slo i selecs a service vecor l ) such ha W Q), ) l )) max W Q), l X Y )), 18) l L l= l Q))

11 11 where X Y )) is a random variable ha depends on Y ) i.e., he disribuion or X Y )) is deermined by he curren sysem sae Y )), and if { Y )} 1 has a saionary disribuion hen he expecaion EX Y ) under he saionary disribuion is less han or equal o γ. Any service vecor ha saisfies 18) is called a γ-inaccurae service vecor. Noe ha if γ is large, hen he number of γ-inaccurae service vecors will be large and hence he ime needed o find one such service vecor may be small. We show ha for appropriae choices of V all sable γ-inaccurae policies are ǫ-hroughpu opimal. Theorem 3: Le λ be any sabilizable arrival rae vecor and γ be an arbirary γ-inaccurae policy. Then, for every ǫ > 0 and γ <, here exiss V such ha for every V V, 1) if { Y )} 1 is a posiive recurren Marov chain, hen γ is ǫ-hroughpu opimal, and 2) if EX Y ) Y ) = Y γ for every Y, hen { Y )} 1 is a posiive recurren Marov chain, and γ sabilizes he sysem. Now, we provide he inuiion. For simpliciy of explanaion, we assume ha X Y ) γ for every Y, and hence he condiion in 2) of Theorem 3 holds. We firs explain why γ-inaccurae policies maximize he sabiliy region 34. For large queue lenghs, max l L W Q, l) >> γ, and hence from 18), W Q, l γ ) max l L W Q, l). l= l Q) l= l Q) Thus, O and γ selec similar service vecors when he queue lenghs are large. We have shown ha for every sabilizable λ, O has a negaive drif when he queue lenghs are large. Thus, γ also has a negaive drif for large queue lenghs. Hence, γ sabilizes he sysem whenever λ is sabilizable. Incidenally, oher approximae policies may also maximize he sabiliy region. For example, any policy ha saisfies 18) wih W Q), l )) and W Q), l) replaced by W 1 Q), l )) and W 1 Q), l) respecively, maximize he sabiliy region 23, 32, 34. The ey difference beween only sabilizing he sysem and aaining he maximum possible hroughpu subjec o sabilizing he sysem is ha whereas for he former i is sufficien o appropriaely selec he service vecor when he queue lenghs are large, bu for he laer appropriae selecion of service vecors is required for all values of queue lenghs. Hence, i is no clear ha γ maximizes he hroughpu as well; we now explain why his is in fac somewha couner-inuiive. Noe ha for small queue lenghs max l L W Q, l) may be smaller or comparable l= l Q) wih γ. Then, 18) does no guaranee ha he service vecors seleced by γ and O are similar. Hence, i is no clear ha γ achieves he same hroughpu as O, which aains he maximum hroughpu. We now explain why Theorem 3 holds. We argue ha for proper choice of parameers he queue lenghs and he service vecors under O and γ become similar. Clearly, in he firs slo, boh sysems have he same queue lengh vecor, Q. Now, noe ha for large V, W Q, l1 ) and W Q, l 2 ) significanly differ if l 1 and l 2 are significanly differen. Thus, due o 18), and since W Q, l γ ) max l L W Q, l), l γ l O. Thus, he queue lenghs in l= l Q) he nex slo are also similar in boh sysems. Recursive use of he same argumen shows ha he queue lenghs and he service vecors seleced in each slo are similar in boh sysems. Thus, boh policies aain similar hroughpu. Thus, γ is hroughpu opimal for large V. Nex, we prove Theorem 3. Proof: We assume ha λ is sabilizable. We define he following Lyapunov funcion: f Y )) = Q )) 2. Using analysis similar o ha for obaining 8), E f Y + 1)) f Y )) Y ) M + 2Q )λ 2E W Q), l γ )) Y ) 2V E G r l γ ))l γ ) Y ). 19)

12 12 From 18) and 19), i follows ha E f Y + 1)) f Y )) Y ) M + 2Q )λ + 2E X Y )) Y ) 2E max W Q), l) Y ) l L l= l Q)) 2V E G r l γ ))l γ ) Y ). Using argumens similar o hose in he proof of 12) from 8), we can prove ha E f Y + 1)) f Y )) Y ) M + 2E X Y )) Y ) 2δQ ) +2V G λ + δ) 2V U λ,δ) 2V E G r l γ ))l γ ) Y ), 20) where δ is such ha λ,δ) is ǫ/2-hroughpu opimal. 1) Proof for 1): Le he process { Y )} 1 be a posiive recurren Marov chain. Then his process has a saionary disribuion. Taing expecaion on boh sides of 20) wih respec o his saionary disribuion, we obain E f Y + 1)) f Y )) M + 2E X Y )) 2δEQ ) + 2V G λ + δ) 2V U λ,δ) 2V E G r l γ ))l γ ). 21) Since { Y )} 1 is a posiive recurren Marov chain, 1 u El γ ) = lim l γ u u v) = λ w.p. 1, 22) and 1 = lim u u v=1 Er l γ ))l γ ) u v=1 r l γ v))l γ v) w.p. 1 = Ω γ λ) from 1)). 23) From saionariy, Ef Y + 1)) = Ef Y )). 24)

13 13 From 21), 22), 23) and 24), and since from Definiion 8), E X Y )) γ, i follows ha Ω γ λ) U λ,δ) δ G M + 2γ 2V Ω max λ) ǫ 2 M + 2γ from Theorem 1) 2V Ω max λ) ǫ if V M + 2γ. ǫ The resul follows. 2) Proof for 2): Now, le E X Y )) Y ) = Y γ for all Y. Thus, from 20), for every γ and V 0, E f Y + 1)) f Y )) Y ) M + 2γ 2δQ ) + 2V G λ + δ). Le B = { Y = Q, l) : l L, n =0 Q V n δ =0 G λ + δ) + M+2γ+1 2δ }. Thus, { E f Y + 1)) f Y )) Y for all Y ), ) < 1 if Y ) B. Thus, since B is finie, by Foser s Theorem Theorem in 20), { Y )} 1 is posiive recurren, and he expecaions of he queue lenghs under is saionary disribuion are finie. Hence, γ sabilizes he sysem. The main challenge in compuing γ-inaccurae service vecors is ha W Q, l O ) may no be nown and in mos cases is compuaion is complex. Thus, even he verificaion of wheher a given l is γ-inaccurae may be compuaionally complex. We circumven his challenge by designing a compuaionally simple approach ha obains γ-inaccurae service vecors wihou requiring he nowledge of W Q, l O ). B. Periodic Compuaion of Opimal Schedule We divide he ime axis in inervals of lengh T, i.e., in inervals of he form KT,K + 1)T 1. Le ) l OPT ) def = arg max W Q), l. l L l= l Q)) We consider a policy T ha compues l OPT ) a he beginning of each inerval, i.e., in he slos KT for K 0, and hroughou he inerval serves each seleced queue while i is non-empy. The ime needed o compue T is O2 n /T) in he amorized sense, i.e., sup 1 { u=1 cu)/} is O2n /T) on every sample pah, where cu) is he compuaional complexiy is slo u 17. Thus, if we choose T o be sufficienly large 2 n ), hen T requires O1) compuaion ime in he amorized sense. In he following lemma, we show ha X Y ) γ for all Y where γ = nta max + 2). Lemma 1: Le Q) be he queue lengh vecor under T in. Then, W Q), l T )) max W l L l= l Q)) Q), l ) nta max + 2). Proof: Wihou loss of generaliy le KT,K + 1)T 1 for some K. Now, from 5), Q ) Q KT) + u=kt Λ u) Q KT) + nta max. 25)

14 14 Similarly, from 5), Q ) Now, from 4), 25) and 26), we obain Now, from 3), max W l L l= l Q)) Q KT) u=kt l T u) Q KT) nt. 26) Q), l ) W Q), l T )) Q KT) V G r l OPT )))l OPT ) Q KT) V G r l T )))l T ) +nta max + 1) = W QKT), l OPT )) W QKT), l T )) +nta max + 1) 27) W QKT), l OPT )) W QKT), l OPT KT)). 28) Also, since he service vecor seleced by T changes in he inerval only if some queues empy during he period and hen he change ) is o no serve) hem, l T ) l OPT KT) and if l T ) < lopt KT) hen Q KT) T. Thus, r l T ) r l OPT KT) for all for which Q ) > 0. Hence, W QKT), l T )) W QKT), l OPT KT)) nt. 29) The resul follows from 27), 28) and 29). Bu, Y ) is no a Marov chain. Thus, in spie of Lemma 1, T is no γ accurae. Now, Y T) is an irreducible, aperiodic, Marov chain, and he framewor for γ-inaccurae scheduling can be generalized o such cases. We omi his generalizaion for breviy. Bu, using Lemma 1 and a proof similar o ha for Theorem 3, we can prove ha when λ is sabilizable, T is ǫ-hroughpu opimal for every ǫ > 0. We formally sae his in he following heorem, and prove i in he appendix. Theorem 4: Le λ be any sabilizable arrival rae vecor. Then, for every ǫ > 0, here exiss V such ha for every V V he policy T is ǫ-hroughpu opimal. The main challenge in using T is ha i needs o periodically compue he opimal service vecor. Since he ime required in each such compuaion is exponenial in n, for large n, such compuaions may become infeasible. We nex propose an opimal randomized policy which requires On) compuaion ime in every slo. C. Opimal Randomized Policy R ) We now propose a randomized policy R which has been inspired by a randomized policy proposed by Tassiulas 34. The policy in 34 aains he maximum possible sabiliy region in a consrained queueing newor using linear ime compuaions in each slo. Our conribuion here is o show ha linear-ime compuable randomized policies can also maximize he hroughpu subjec o sabilizing he sysem. We now describe R.) In every slo 0, R generaes a service vecor l) randomly among all service vecors l L such ha l Q) = l as per a disribuion P Q ). In every slo 1, once a random vecor is generaed as above, R obains l R ) ieraively as follows. )) ) l) : W Q), l R 1)) Q) <W Q), l) l R )= ) l R 1)) Q) : oherwise.

15 15 Thus, in any slo, R uses a new service vecor only when i increases he value of W ); oherwise i coninues wih he service vecor used in he previous slo. I is ineresing o observe ha he randomized policy proposed by Tassiulas 34, which maximizes he sabiliy region using linear compuaion ime in each slo, uses a new service vecor only when i increases he value of W 1 ). Noe ha he disribuion PQ ) may depend on he curren queue lengh vecor. We only consider disribuions P Q ) such ha for every Q, P Q l OPT ) µ for some µ > 0. Lemma 2: Le λ be a sabilizable arrival rae vecor. Then { Y ) = Q), l γ ))} 1 is a posiive recurren Marov chain, and R sabilizes he sysem. We prove Lemma 2 in appendix. Now, we show ha R is n µ A max + 2)-inaccurae. Lemma 3: Le Q) be he queue lengh vecor under R in. Then, for any iniial disribuion of Y ), ) E max W Q), l W Q), l R )) n l L µ A max + 2). l= l Q)) Proof: Since P Q l OPT ) µ for every Q, l R ) = l OPT ) infiniely ofen w.p. 1. Le {κ K } K 1 be he slos in which l R ) = l OPT ). Again, since P Q l OPT ) µ for every Q, Eκ K+1 κ K 1/µ. Consider he K for which κ K,κ K+1 1. Lie in Lemma 1, we obain ) max W Q), l W Q), l R )) l L l= l Q)) nκ K+1 κ K )A max + 2). Thus, he resul follows since Eκ K+1 κ K 1/µ. Now, from par 1 of Theorem 3 and Lemmas 2 and 3, i follows ha R is ǫ-hroughpu-opimal for any µ. We formally sae his in he following heorem. Theorem 5: Le λ be a sabilizable arrival rae vecor. Then, for every ǫ > 0, here exiss V such ha for every V V he policy R is ǫ-hroughpu opimal. Now, if µ = 2 n, R can be compued in On) ime in each slo. Each non-empy queue can be seleced w.p. 1/2. If he resuling vecor is no in L, hen no queue is served. D. Disribued Implemenaion of R and T Disribued scheduling can be defined in differen ways. One definiion is o consider a policy as disribued if each node selecs is acion based on is observaion, sae and he informaion i acquires by exchanging messages wih is neighbors. Such policies are hen evaluaed on he basis of heir performance and he frequency and he amoun of message exchange. Anoher definiion is o consider a policy as disribued if each node selecs is acion based on is observaion, sae and he saes and acions of nodes in a cerain neighborhood. We firs describe how R and T can be implemened as per he firs definiion. The ime axis can be divided in periods of lengh T. Each node can broadcas is queue lengh a he beginning of every period. The period lengh T should be seleced so ha he broadcass in a period reach oher nodes in he same period. For execuing T, each node compues he opimal service vecor a he beginning of every period based on he broadcass i receives in he previous period. For execuing R, each node randomly selecs a service vecor a he beginning of each period, and subsequenly chooses beween he service vecors seleced in he curren and previous periods based on he broadcass i receives in he previous period, and finally uses he chosen service vecor hroughou he period. All nodes use he same seeds in he random number generaors and herefore obain he same random selecions. For boh policies, each node s compuaions depend on he queue lenghs of oher nodes in he previous period. Theorems 4 and 5 sill hold. The message exchange complexiy can be made arbirarily small in boh cases by increasing T. Deerminaion of an opimal policy which is disribued as per he second definiion for disribued scheduling remains open. Noe ha he design of such scheduling policies in he precursor problem, ha of maximizing he sabiliy region, is sill no compleely undersood, alhough some illuminaing resuls have been obained recenly

16 16 10, 29, 36. We hope ha he opimaliy resuls in his paper and he recen advances in conex of disribued scheduling will moivae furher exploraion of he above open problem. Finally, Ross e. al. has obained local search based policies, which are liely o be compuaionally simple in pracice, for maximizing he sabiliy region of cerain classes of consrained queueing newors 31. I will be ineresing o deermine wheher he hroughpu can be maximized subjec o sabilizing he sysem using similar local search policies, and how he compuaion ime required by he γ-accurae policies we propose compare wih hose for he resuling local search policies. VI. DISCUSSIONS AND GENERALIZATIONS We now generalize our framewor so as o obain opimal policies when some of he assumpions made in Secion II do no hold. Firs, we have so far assumed ha a pace is discarded only afer i is ransmied. We discuss how our framewor can be generalized o allow a queue o discard some or all paces before ransmiing hem, and examine he advanages and disadvanages of his opion Subsecion VI-A). We nex describe how our framewor can be generalized o accommodae random rewards and random ses of valid service vecors L Subsecion VI-B). Finally, we discuss how L and reward funcions can be chosen so as o aain cerain performance goals in an imporan applicaion domain for his framewor ha of wireless newors Subsecion VI-C). A. Discarding paces before ransmission In Secion II, we have assumed ha each pace is discarded from is queue only afer i is served once. However, in pracice, a pace may be discarded from is queue even before i is served. The availabiliy of his opion enhances he sabiliy region, and is judicious use increases he hroughpu. For example, in Example 1 in Secion III, when λ = 1 ǫ,1 ǫ) where ǫ is a small posiive number, Ωmax λ) 4 Ω b 4). Now, if S 2 can discard paces before serving hem, is sable and aains a hroughpu close o 5. Bu, clearly, indiscriminae use of his opion subsanially reduces he hroughpu. We now show ha appropriae augmenaion of L allows us o design policies ha aain he maximum possible hroughpu in presence of his opion. Le O be he original sysem ha does no allow paces o be discarded before ransmission, and le Ô be he new sysem which allows he above. In Ô a queue is said o be served when a pace is removed from is queue. The service vecors in Ô have 2n 0 1 componens. The firs n componens denoe which queues are being served and he remaining componens denoe wheher he paces from he queues ha are being served are ransmied or discarded before ransmission. We obain he se L of valid service vecors of Ô from he corresponding se L of O as follows. Le l L and le l have i 0 componens where 0 i n. Now, l corresponds o 2 i service vecors in L, and each of hese service vecors a) ransmi paces from he queues l were serving in O and b) discard paces from a cerain possibly empy) subse of queues which l were no serving in O. Noe ha he se of queues l were no serving in O has 2 i subses. Thus, he number of service vecors generaed by l is 2 i. Le l1 be one such service vecor generaed by l. Since l and l1 ransmi paces from he same queues, r l1 ) = r l) for each {1,...,n}. The sabiliy region of Ô is a possibly improper) superse of ha of O. This is because as long as he arrival rae of a queue is less han 1 i can be sabilized in Ô by simply discarding all is paces before ransmission. Thus, he sabiliy region of Ô is a superse of { λ : 0 λ i < 1 i = 1,...,n} and a subse of { λ : 0 λ i 1 i = 1,...,n}. For any λ ha is sabilizable in O, he maximum hroughpu of a sable policy in O is less han or equal o ha of he maximum hroughpu of a sable policy in Ô. This is because every policy in O is a valid policy in Ô, since for each l L here exiss l L ha does no discard paces from any queue before ransmission, and ransmis paces from he same queues which l serves. Noe ha, O, γ, T and R can be defined similar o ha in O; he only difference is ha L mus be subsiued by L. The performance guaranees for hese generalized versions, i.e., Theorems 1 o 5, hold in Ô. are he same as hose for O. However, noe ha higher hroughpu and sabiliy region can be aained in Ô while sacrificing fairness. Specifically, for any λ ha is sabilizable in O, if an ǫ hroughpu opimal policy ˆ in Ô aains a hroughpu which is higher han ha of an ǫ hroughpu opimal policy in O, ˆ discards paces before ransmission from some queues. Thus, ˆ aains higher hroughpu by being unfair o some queues. Also, in communicaion newors, in presence of his opion, some receivers may only receive a small fracion of paces ransmied by he corresponding

17 17 sources, which will in urn preven hem from successfully decoding he ransmied informaion. Thus, his opion is no liely o be widely used refer o Secion VI-C). B. Random rewards and random L We have so far assumed ha he reward received by he h queue in is compleely deermined by he service vecor chosen in. We now allow he rewards o be random variables r.v. s) ha depend on an exernal random componen in addiion o he service vecor Subsecion VI-B). This generalizaion is relevan in conex of wireless newors, where he success of a ransmission is a random even whose probabiliy depends on he fading sae of he channels. Thus, in one-o-many or one-o-one communicaion he reward is a r.v. whose disribuion depends on he service vecor and he channel fading sae beween each sender-receiver pair. We generalize, O, γ, T and R so as o maximize he hroughpu subjec o sabilizing he sysem in presence of random rewards. We firs formally describe he generalizaion. We consider a random process {S)} 1 which in any slo is in sae S i wih probabiliy b i, i = 1,2,...,z, independen of is sae in any oher slo and also independen of he arrival process in any slo. Here, b i > 0 for each i {1,...,z}. The policy nows S) a he beginning of slo. The reward received by he h queue when l is he service vecor and he process S) is in sae S is a random variable, R l,s), whose disribuion depends on l and S. Le r l,s) def = ER l,s) G for every l and S. We assume ha R l,s) = 0 if l = 0, and r l 1,S) r l 2,S) if l 1 l 2 and l > 0. Thus, he hroughpu Ω λ) under a policy and arrival rae vecor λ is Ω λ) = lim inf = lim inf 1 1 u=1 u=1 R l u), Su)) R l u), Su))l u). 30) Finally, when he arrival rae vecor is λ, he maximum hroughpu of any sable policy is Ω max λ). We nex elucidae he above formalisms wih a specific example. Example 4: In Fig.1 assume ha he channel o each receiver is in good bad, resp.) sae w.p ), and each receiver can decode he pace w.p , resp.) when is channel is in good bad, resp.) sae and i is no in he range of any oher sender ha is ransmiing paces. The sae of a channel in a slo is independen of ha in oher slos and also independen of he saes of oher channels in any slo. In each slo, he sysem nows he saes of all channels, bu does no now wheher a receiver can decode he pace is sender ransmis. Thus, he sysem has 64 saes corresponding o differen combinaions of channel saes. Now, if l = l 2 l = l 4, resp.) R 1 l,s) equals he number of receivers in he se {R 1,...,R 5 } {R 3,R 4,R 5 }, resp.) ha can decode he pace S 1 ransmis and if l { l 1, l 3 }), R 1 l,s) = 0. Nex, R 2 l,s) = 1 if l 2 = 1 and R 6 can decode he pace, R 2 l,s) = 0 oherwise. Thus, R 1 l,s),r 2 l,s) are random variables whose disribuions depend on l,s. For example, r 1 l 2,S) = 4.5 if S is such ha he channels o R 1,...,R 5 are in good sae, r 2 l,s) = 0.9 if l 2 = 1 and S is such ha he channel o R 6 is in good sae. Firs, noe ha since L does no depend on S), he sabiliy region of he sysem remains he same. Now, we presen he opimaliy resuls. We firs describe how λ,δ) can be generalized. In any slo in which S) = S z, he generalized ) policy λ,δ) selecs l i w.p. w iz. If l i is seleced in slo, hen he sysem selecs service vecor li Q) i.e., l b λ,δ) ) = l i Q) )). The probabiliy disribuion w z = w 1z w mz for every z = 1,...,Z is compued using he following linear program. LP λ,δ) :- Maximize: Û λ,δ) = Z z=1 b m n z i=1 w izl i r l i,s z ) Subjec o: 1) m i=1 w iz = 1 for every z {1,...,Z} 2) w iz 0 for every i {1,...,m} and z {1,...,Z} 3) Z z=1 b z m i=1 w izl i = λ + δ for every {1,...,n}. Noe ha LP λ,δ) is similar o LP λ,δ); he only difference is ha he disribuion for selecing he service vecors depends on he sae S) of he sysem.

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Lecture 4 Notes (Little s Theorem)

Lecture 4 Notes (Little s Theorem) Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,

More information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

INTRODUCTION TO MACHINE LEARNING 3RD EDITION

INTRODUCTION TO MACHINE LEARNING 3RD EDITION ETHEM ALPAYDIN The MIT Press, 2014 Lecure Slides for INTRODUCTION TO MACHINE LEARNING 3RD EDITION alpaydin@boun.edu.r hp://www.cmpe.boun.edu.r/~ehem/i2ml3e CHAPTER 2: SUPERVISED LEARNING Learning a Class

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 2, May 2013, pp. 209 227 ISSN 0364-765X (prin) ISSN 1526-5471 (online) hp://dx.doi.org/10.1287/moor.1120.0562 2013 INFORMS On Boundedness of Q-Learning Ieraes

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs PROC. IEEE CONFERENCE ON DECISION AND CONTROL, 06 A Primal-Dual Type Algorihm wih he O(/) Convergence Rae for Large Scale Consrained Convex Programs Hao Yu and Michael J. Neely Absrac This paper considers

More information

Basic definitions and relations

Basic definitions and relations Basic definiions and relaions Lecurer: Dmiri A. Molchanov E-mail: molchan@cs.u.fi hp://www.cs.u.fi/kurssi/tlt-2716/ Kendall s noaion for queuing sysems: Arrival processes; Service ime disribuions; Examples.

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

Optimality Conditions for Unconstrained Problems

Optimality Conditions for Unconstrained Problems 62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x

More information

Christos Papadimitriou & Luca Trevisan November 22, 2016

Christos Papadimitriou & Luca Trevisan November 22, 2016 U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream

More information

Resource Allocation in Visible Light Communication Networks NOMA vs. OFDMA Transmission Techniques

Resource Allocation in Visible Light Communication Networks NOMA vs. OFDMA Transmission Techniques Resource Allocaion in Visible Ligh Communicaion Neworks NOMA vs. OFDMA Transmission Techniques Eirini Eleni Tsiropoulou, Iakovos Gialagkolidis, Panagiois Vamvakas, and Symeon Papavassiliou Insiue of Communicaions

More information

Learning a Class from Examples. Training set X. Class C 1. Class C of a family car. Output: Input representation: x 1 : price, x 2 : engine power

Learning a Class from Examples. Training set X. Class C 1. Class C of a family car. Output: Input representation: x 1 : price, x 2 : engine power Alpaydin Chaper, Michell Chaper 7 Alpaydin slides are in urquoise. Ehem Alpaydin, copyrigh: The MIT Press, 010. alpaydin@boun.edu.r hp://www.cmpe.boun.edu.r/ ehem/imle All oher slides are based on Michell.

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

On the Stability Region of Multi-Queue Multi-Server Queueing Systems with Stationary Channel Distribution

On the Stability Region of Multi-Queue Multi-Server Queueing Systems with Stationary Channel Distribution On he Sabiliy Region of Muli-Queue Muli-Server Queueing Sysems wih Saionary Channel Disribuion Hassan Halabian, Ioannis Lambadaris, Chung-Horng Lung Deparmen of Sysems and Compuer Engineering, Carleon

More information

Learning a Class from Examples. Training set X. Class C 1. Class C of a family car. Output: Input representation: x 1 : price, x 2 : engine power

Learning a Class from Examples. Training set X. Class C 1. Class C of a family car. Output: Input representation: x 1 : price, x 2 : engine power Alpaydin Chaper, Michell Chaper 7 Alpaydin slides are in urquoise. Ehem Alpaydin, copyrigh: The MIT Press, 010. alpaydin@boun.edu.r hp://www.cmpe.boun.edu.r/ ehem/imle All oher slides are based on Michell.

More information

Decentralized Stochastic Control with Partial History Sharing: A Common Information Approach

Decentralized Stochastic Control with Partial History Sharing: A Common Information Approach 1 Decenralized Sochasic Conrol wih Parial Hisory Sharing: A Common Informaion Approach Ashuosh Nayyar, Adiya Mahajan and Demoshenis Tenekezis arxiv:1209.1695v1 [cs.sy] 8 Sep 2012 Absrac A general model

More information

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Mean-square Stability Control for Networked Systems with Stochastic Time Delay JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455

More information

Optimal Server Assignment in Multi-Server

Optimal Server Assignment in Multi-Server Opimal Server Assignmen in Muli-Server 1 Queueing Sysems wih Random Conneciviies Hassan Halabian, Suden Member, IEEE, Ioannis Lambadaris, Member, IEEE, arxiv:1112.1178v2 [mah.oc] 21 Jun 2013 Yannis Viniois,

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

Comments on Window-Constrained Scheduling

Comments on Window-Constrained Scheduling Commens on Window-Consrained Scheduling Richard Wes Member, IEEE and Yuing Zhang Absrac This shor repor clarifies he behavior of DWCS wih respec o Theorem 3 in our previously published paper [1], and describes

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

h[n] is the impulse response of the discrete-time system:

h[n] is the impulse response of the discrete-time system: Definiion Examples Properies Memory Inveribiliy Causaliy Sabiliy Time Invariance Lineariy Sysems Fundamenals Overview Definiion of a Sysem x() h() y() x[n] h[n] Sysem: a process in which inpu signals are

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

Lecture 2 October ε-approximation of 2-player zero-sum games

Lecture 2 October ε-approximation of 2-player zero-sum games Opimizaion II Winer 009/10 Lecurer: Khaled Elbassioni Lecure Ocober 19 1 ε-approximaion of -player zero-sum games In his lecure we give a randomized ficiious play algorihm for obaining an approximae soluion

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

Appendix 14.1 The optimal control problem and its solution using

Appendix 14.1 The optimal control problem and its solution using 1 Appendix 14.1 he opimal conrol problem and is soluion using he maximum principle NOE: Many occurrences of f, x, u, and in his file (in equaions or as whole words in ex) are purposefully in bold in order

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Approximation Algorithms for Unique Games via Orthogonal Separators

Approximation Algorithms for Unique Games via Orthogonal Separators Approximaion Algorihms for Unique Games via Orhogonal Separaors Lecure noes by Konsanin Makarychev. Lecure noes are based on he papers [CMM06a, CMM06b, LM4]. Unique Games In hese lecure noes, we define

More information

Expert Advice for Amateurs

Expert Advice for Amateurs Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he

More information

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011 Mainenance Models Prof Rober C Leachman IEOR 3, Mehods of Manufacuring Improvemen Spring, Inroducion The mainenance of complex equipmen ofen accouns for a large porion of he coss associaed wih ha equipmen

More information

PROOF FOR A CASE WHERE DISCOUNTING ADVANCES THE DOOMSDAY. T. C. Koopmans

PROOF FOR A CASE WHERE DISCOUNTING ADVANCES THE DOOMSDAY. T. C. Koopmans PROOF FOR A CASE WHERE DISCOUNTING ADVANCES THE DOOMSDAY T. C. Koopmans January 1974 WP-74-6 Working Papers are no inended for disribuion ouside of IIASA, and are solely for discussion and informaion purposes.

More information

Solutions for Assignment 2

Solutions for Assignment 2 Faculy of rs and Science Universiy of Torono CSC 358 - Inroducion o Compuer Neworks, Winer 218 Soluions for ssignmen 2 Quesion 1 (2 Poins): Go-ack n RQ In his quesion, we review how Go-ack n RQ can be

More information

Martingales Stopping Time Processes

Martingales Stopping Time Processes IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765. Volume 11, Issue 1 Ver. II (Jan - Feb. 2015), PP 59-64 www.iosrjournals.org Maringales Sopping Time Processes I. Fulaan Deparmen

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

Families with no matchings of size s

Families with no matchings of size s Families wih no machings of size s Peer Franl Andrey Kupavsii Absrac Le 2, s 2 be posiive inegers. Le be an n-elemen se, n s. Subses of 2 are called families. If F ( ), hen i is called - uniform. Wha is

More information

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem. Noes, M. Krause.. Problem Se 9: Exercise on FTPL Same model as in paper and lecure, only ha one-period govenmen bonds are replaced by consols, which are bonds ha pay one dollar forever. I has curren marke

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems Essenial Microeconomics -- 6.5: OPIMAL CONROL Consider he following class of opimizaion problems Max{ U( k, x) + U+ ( k+ ) k+ k F( k, x)}. { x, k+ } = In he language of conrol heory, he vecor k is he vecor

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

3.1 More on model selection

3.1 More on model selection 3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1 SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

CMU-Q Lecture 3: Search algorithms: Informed. Teacher: Gianni A. Di Caro

CMU-Q Lecture 3: Search algorithms: Informed. Teacher: Gianni A. Di Caro CMU-Q 5-38 Lecure 3: Search algorihms: Informed Teacher: Gianni A. Di Caro UNINFORMED VS. INFORMED SEARCH Sraegy How desirable is o be in a cerain inermediae sae for he sake of (effecively) reaching a

More information

Homework 4 (Stats 620, Winter 2017) Due Tuesday Feb 14, in class Questions are derived from problems in Stochastic Processes by S. Ross.

Homework 4 (Stats 620, Winter 2017) Due Tuesday Feb 14, in class Questions are derived from problems in Stochastic Processes by S. Ross. Homework 4 (Sas 62, Winer 217) Due Tuesday Feb 14, in class Quesions are derived from problems in Sochasic Processes by S. Ross. 1. Le A() and Y () denoe respecively he age and excess a. Find: (a) P{Y

More information

Presentation Overview

Presentation Overview Acion Refinemen in Reinforcemen Learning by Probabiliy Smoohing By Thomas G. Dieerich & Didac Busques Speaer: Kai Xu Presenaion Overview Bacground The Probabiliy Smoohing Mehod Experimenal Sudy of Acion

More information

Games Against Nature

Games Against Nature Advanced Course in Machine Learning Spring 2010 Games Agains Naure Handous are joinly prepared by Shie Mannor and Shai Shalev-Shwarz In he previous lecures we alked abou expers in differen seups and analyzed

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

BU Macro BU Macro Fall 2008, Lecture 4

BU Macro BU Macro Fall 2008, Lecture 4 Dynamic Programming BU Macro 2008 Lecure 4 1 Ouline 1. Cerainy opimizaion problem used o illusrae: a. Resricions on exogenous variables b. Value funcion c. Policy funcion d. The Bellman equaion and an

More information

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

More information

Chapter 7: Solving Trig Equations

Chapter 7: Solving Trig Equations Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions

More information

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013 Mahemaical Theory and Modeling ISSN -580 (Paper) ISSN 5-05 (Online) Vol, No., 0 www.iise.org The ffec of Inverse Transformaion on he Uni Mean and Consan Variance Assumpions of a Muliplicaive rror Model

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive

More information

An Excursion into Set Theory using a Constructivist Approach

An Excursion into Set Theory using a Constructivist Approach An Excursion ino Se Theory using a Consrucivis Approach Miderm Repor Nihil Pail under supervision of Ksenija Simic Fall 2005 Absrac Consrucive logic is an alernaive o he heory of classical logic ha draws

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,

More information

Intelligent Packet Dropping for Optimal Energy-Delay Tradeoffs in Wireless Downlinks

Intelligent Packet Dropping for Optimal Energy-Delay Tradeoffs in Wireless Downlinks IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 54, NO. 3, PP. 565-579, MARCH 2009. 1 Inelligen Packe Dropping for Opimal Energy-Delay Tradeoffs in Wireless Downlinks Michael J. Neely Universiy of Souhern

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3

More information

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990),

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990), SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F Trench SIAM J Marix Anal Appl 11 (1990), 601-611 Absrac Le T n = ( i j ) n i,j=1 (n 3) be a real symmeric

More information

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL: Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

A Shooting Method for A Node Generation Algorithm

A Shooting Method for A Node Generation Algorithm A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan

More information

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model

More information

Orientation. Connections between network coding and stochastic network theory. Outline. Bruce Hajek. Multicast with lost packets

Orientation. Connections between network coding and stochastic network theory. Outline. Bruce Hajek. Multicast with lost packets Connecions beween nework coding and sochasic nework heory Bruce Hajek Orienaion On Thursday, Ralf Koeer discussed nework coding: coding wihin he nework Absrac: Randomly generaed coded informaion blocks

More information

4 Sequences of measurable functions

4 Sequences of measurable functions 4 Sequences of measurable funcions 1. Le (Ω, A, µ) be a measure space (complee, afer a possible applicaion of he compleion heorem). In his chaper we invesigae relaions beween various (nonequivalen) convergences

More information

Transform Techniques. Moment Generating Function

Transform Techniques. Moment Generating Function Transform Techniques A convenien way of finding he momens of a random variable is he momen generaing funcion (MGF). Oher ransform echniques are characerisic funcion, z-ransform, and Laplace ransform. Momen

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

For example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l,

For example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l, Comb Filers The simple filers discussed so far are characeried eiher by a single passband and/or a single sopband There are applicaions where filers wih muliple passbands and sopbands are required The

More information

Convergence of the Neumann series in higher norms

Convergence of the Neumann series in higher norms Convergence of he Neumann series in higher norms Charles L. Epsein Deparmen of Mahemaics, Universiy of Pennsylvania Version 1.0 Augus 1, 003 Absrac Naural condiions on an operaor A are given so ha he Neumann

More information

DIFFERENTIAL GEOMETRY HW 5

DIFFERENTIAL GEOMETRY HW 5 DIFFERENTIAL GEOMETRY HW 5 CLAY SHONKWILER 3. Le M be a complee Riemannian manifold wih non-posiive secional curvaure. Prove ha d exp p v w w, for all p M, all v T p M and all w T v T p M. Proof. Le γ

More information

Examples of Dynamic Programming Problems

Examples of Dynamic Programming Problems M.I.T. 5.450-Fall 00 Sloan School of Managemen Professor Leonid Kogan Examples of Dynamic Programming Problems Problem A given quaniy X of a single resource is o be allocaed opimally among N producion

More information

Intelligent Packet Dropping for Optimal Energy-Delay Tradeoffs in Wireless Downlinks

Intelligent Packet Dropping for Optimal Energy-Delay Tradeoffs in Wireless Downlinks PROC. OF 4TH INT. SYMPOSIUM ON MODELING AND OPTIMIZATION IN MOBILE, AD HOC, AND WIRELESS NETWORKS (WIOPT), APRIL 2006 1 Inelligen Packe Dropping for Opimal Energy-Delay Tradeoffs in Wireless Downlinks

More information

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning

More information

Lecture Notes 3: Quantitative Analysis in DSGE Models: New Keynesian Model

Lecture Notes 3: Quantitative Analysis in DSGE Models: New Keynesian Model Lecure Noes 3: Quaniaive Analysis in DSGE Models: New Keynesian Model Zhiwei Xu, Email: xuzhiwei@sju.edu.cn The moneary policy plays lile role in he basic moneary model wihou price sickiness. We now urn

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS

More information

Logic in computer science

Logic in computer science Logic in compuer science Logic plays an imporan role in compuer science Logic is ofen called he calculus of compuer science Logic plays a similar role in compuer science o ha played by calculus in he physical

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Complexity of Inference in Topic Models

Complexity of Inference in Topic Models Complexiy of Inference in Topic Models David Sonag, Daniel M. Roy Massachuses Insiue of Technology {dsonag,droy}@csail.mi.edu Absrac We consider he compuaional complexiy of finding he MAP assignmen of

More information

Stationary Distribution. Design and Analysis of Algorithms Andrei Bulatov

Stationary Distribution. Design and Analysis of Algorithms Andrei Bulatov Saionary Disribuion Design and Analysis of Algorihms Andrei Bulaov Algorihms Markov Chains 34-2 Classificaion of Saes k By P we denoe he (i,j)-enry of i, j Sae is accessible from sae if 0 for some k 0

More information

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive. LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis

More information

Unit Root Time Series. Univariate random walk

Unit Root Time Series. Univariate random walk Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he

More information

Hamilton Jacobi equations

Hamilton Jacobi equations Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion

More information

Lower and Upper Bounds on FIFO Buffer Management in QoS Switches

Lower and Upper Bounds on FIFO Buffer Management in QoS Switches Lower and Upper Bounds on FIFO Buffer Managemen in QoS Swiches Mahias Engler Mahias Wesermann Deparmen of Compuer Science RWTH Aachen 52056 Aachen, Germany {engler,marsu}@cs.rwh-aachen.de Absrac We consider

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS Exam: ECON4325 Moneary Policy Dae of exam: Tuesday, May 24, 206 Grades are given: June 4, 206 Time for exam: 2.30 p.m. 5.30 p.m. The problem se covers 5 pages

More information