Opportunistic Scheduling with Reliability Guarantees in Cognitive Radio Networks

Size: px
Start display at page:

Download "Opportunistic Scheduling with Reliability Guarantees in Cognitive Radio Networks"

Transcription

1 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL 008 Opporunisic Scheduling wih Reliabiliy Guaranees in Cogniive Radio Neworks Rahul Urgaonkar, ichael J. Neely Universiy of Souhern California, Los Angeles, CA hp://www-scf.usc.edu/ urgaonka Absrac We develop opporunisic scheduling policies for cogniive radio neworks ha maximize he hroughpu uiliy of he secondary (unlicensed) users subjec o maximum collision consrains wih he primary (licensed) users. We consider a cogniive nework wih saic primary users and poenially mobile secondary users. We use he echnique of Lyapunov Opimizaion o design an online flow conrol, scheduling and resource allocaion algorihm ha mees he desired objecives and provides explici performance guaranees. Index Terms Cogniive Radio, Queueing Analysis, Resource Allocaion, Lyapunov Opimizaion I. INTRODUCTION Cogniive radio neworks have recenly emerged as a promising echnique o improve he uilizaion of he exising radio specrum. The key enabler is he cogniive radio [ ha can dynamically adjus is operaing poins over a wide range depending on specrum availabiliy. The main idea behind a cogniive nework is for he unlicensed users o exploi he spaially and/or emporally under-uilized specrum by ransmiing opporunisically. However, a basic requiremen is o ensure ha he exising licensed users are no adversely affeced by such ransmissions. In his paper, we develop opporunisic scheduling schemes ha maximize he hroughpu uiliy of he secondary (or unlicensed) users subjec o maximum collision consrains wih he primary (or licensed) users in a cogniive radio nework. A survey on he echnology, design issues and recen work in cogniive radio neworks is provided in [, [3. The problem of opimal specrum assignmen o secondary users in saic neworks is reaed in [4 [7 where i is assumed ha scheduling is aware of primary user ransmissions. Scheduling he secondary users under parial channel sae informaion is considered in [8, [9 which uses a probabilisic maximum collision consrain wih he primary users. In his paper, we use he echniques of adapive queueing and Lyapunov Opimizaion o design an online flow conrol, scheduling and resource allocaion algorihm for a cogniive nework ha maximizes he hroughpu uiliy of he secondary users subjec o a maximum rae of collisions wih he primary users. This algorihm operaes wihou knowing he mobiliy paern of he secondary users and provides explici performance bounds. Lyapunov Opimizaion echniques were This maerial is suppored in par by one or more of he following: he DARPA IT-ANET program gran W9NF , he Naional Science Foundaion gran OCE Fig.. 3 Channel 4 4 Channel 4 3 y Channel Channel 3 5 x Primary User Secondary User Access Poin Example cogniive nework showing primary and secondary users perhaps firs applied o wireless neworks in he landmark paper [0, where Lyapunov drif is used o develop a join opimal rouing and scheduling algorihm. This mehod has since been exended o rea problems of join sabiliy and uiliy opimizaion in general sochasic neworks in [ [4 and wireless mesh neworks in [5. The analysis presened in his paper applies he sochasic nework opimizaion framework of [4. The main conribuions of his work are described below. This paper: Develops hroughpu opimal conrol policies for cogniive neworks wih general inerference and mobiliy models. Inroduces he noion of collision queues ha are used o provide srong reliabiliy bounds in erms of he wors case number of collisions suffered by a primary user in any ime inerval. In paricular, he collision queue mehod here is adaped from he virual power queue echnique of [3. However, he collision queues developed here are designed o ensure reliabiliy consrains, raher han average power consrains. Differen from [3, his requires he inpus o he virual queues o be random collision variables ha can be evaluaed only afer packe ransmission has aken place. Develops easier o implemen consan facor approximaions o he opimal resource allocaion problem.

2 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL 008 II. NETWORK ODEL We consider a cogniive radio nework consising of primary users and N secondary users as shown in Fig.. Each primary user has a unique licensed channel and hese are orhogonal in frequency and/or space. Thus, he primary users can send daa over heir own licensed channels o heir respecive access poins simulaneously. The secondary users do no have any such channels and opporunisically ry o send heir daa o he access poins by uilizing idle primary channels. Such opporuniies are called specrum holes. We consider a ime-sloed model. The primary users are assumed o be saic. However, he secondary users could be mobile so ha he se of channels hey can access can change over ime. In a imeslo, a secondary user can access a subse of he primary channels poenially depending on is curren locaion. This informaion is concisely represened by an N binary channel accessibiliy marix H() = h nm () N where: if sec. user n can access channel m in slo h nm () = 0 else For example, he channel accessibiliy marix for he example nework in Fig. is given by: H() = We assume ha he mobiliy process for he secondary users is such ha he resuling H() process is arkovian and has a well defined seady sae disribuion. Le S() = (S (), S (),..., S ()) represen he curren primary user occupancy sae of he channels. Here, S i () 0, (for i,,..., ) wih he inerpreaion ha S i () = 0 if channel i is occupied by primary user i in imeslo and S i () = if primary user i is idle in imeslo. We assume ha exacly packe can be ransmied over any channel in a imeslo. A secondary user can aemp ransmission over a mos channel subjec o he consrains in H(). This ransmission is successful only when he channel is no being used by is primary user or any oher secondary user. If a secondary user ransmis on a channel which is busy, here is a collision and boh packes are los. To capure he inerference ha a secondary user ransmission may cause on oher channels, for all n,,..., N, m,,...,, we define I nm as he se of channels ha secondary user n inerferes wih when i uses channel m. We include m in he se I nm. We furher define he following indicaor variables (o be used laer): Inm k if k Inm k,,..., = 0 else Clearly, Inm m = for all m, n. Under his inerference model, he following wo condiions are necessary for a We assume ha muli-user deecion/inerference cancelaion is no available so ha if he secondary user aemps o ransmi is own daa when some oher user is also ransmiing, here is enough inerference a he access poin and no daa is successfully received. ransmission by secondary user n on channel m in slo o be successful: ) S m () = ) For all oher secondary users i ransmiing on a channel j,,...,, we have m / I ij (where i,,..., N \ n) This inerference model is general enough o capure scenarios in which he channels may no be orhogonal wih respec o he secondary user ransmissions alhough hey are orhogonal for he primary user ransmissions. In mos pracical siuaions, he cardinaliy of he inerference ses I nm would be small. An imporan special case is when he channels are indeed orhogonal for all secondary user ransmissions, so ha I nm = m n. We assume ha he ses I nm are fixed. In a more general model, hese ses could also change over ime (possibly depending on he secondary user locaion). Our formulaion can be exended o rea hose cases. We assume ha he primary user channel occupancy sae process S() evolves according o a finie sae ergodic arkov chain on he space 0,. I is furher assumed o be independen of he secondary user mobiliy process H(). The channel sae informaion available o he secondary users is described by a probabiliy vecor P() = (P (), P (),..., P ()) where P i () is he probabiliy ha primary user i is idle in imeslo and is defined as P m () = E S m () S( ) m,,...,. We assume ha his informaion is obained eiher by sensing he channels, or hrough a knowledge of he raffic saisics of he primary users, or a combinaion of hese. This models siuaions where he exac channel sae may no be available o he secondary users (e.g. due o limiaions in carrier sensing). These probabiliies capure he inheren sensing measuremen errors associaed wih any primary ransmission deecion algorihm. Inuiively, he closer P() is o S(), he smaller he chances of collisions. Each primary user m receives exogenous daa a a rae ν m packe/slo and can olerae a maximum ime average rae of collisions given by ρ m, where ρ m is a known consan. A. Queueing Dynamics and Conrol Decisions Each secondary user n,,..., N receives daa according o an i.i.d arrival process A n () ha has rae λ n packes/slo. We assume ha he maximum number of arrivals o any secondary user n is upper bounded by a consan value A max every imeslo. This daa arrives a he ranspor layer and flow conrol decisions on how many packes o admi o he nework layer are aken by each secondary user. We assume ha here are no ranspor layer buffers and add/drop decisions are aken immediaely. Le U n () be he backlog in he nework layer queue of secondary user n a he beginning of imeslo. Le R n () be he number of new packes admied ino his queue in slo. Define µ nm () as he number of aemped packe ransmissions when a conrol acion allocaes channel m o In addiion, predicion based echniques could be also be used o ge his informaion.

3 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL secondary user n in slo. Noe ha in his model µ nm () 0, m, n. Also, here is a successful ransmission on channel m only when he necessary condiions specified earlier are me. Then he queueing dynamics of he secondary user n is described by: U n ( + ) = max[u n () µ nm ()S m (), 0 + R n () where () µ nm () 0, m, n () µ nm () h nm () m, n (3) 0 µ nm () n (4) µ nm () = j= i= i n N Iij m µ ij () = 0 m, n (5) 0 R n () A n () (6) Here, inequaliy (3) represens he consrain imposed by he channel accessibiliy marix H(). Inequaliy (4) represens he consrain ha a secondary user can be allocaed a mos channel. (5) represens he second necessary condiion for successful ransmission expressed in erms of he Inm k variables. In he special case of orhogonal channels, his simplifies o he consrain ha a channel can be allocaed o a mos secondary user, i.e., N 0 µ nm () m (7) B. Discussion of Nework odel The above nework model considers access poin based neworks wih saic (or locally mobile) licensed and fully mobile unlicensed users. Examples of real neworks ha can be modeled like his include Wi-Fi, cellular and mesh neworks wih boh licensed and unlicensed users. In such neworks, he licensed users may no schedule heir ransmissions and hus send a any ime hey desire. The unlicensed users mus make an effor o opporunisically use he specrum holes wihou inerfering oo much wih he licensed users, and hence need sophisicaed scheduling mechanisms. A axonomy of differen approaches o specrum sharing in cogniive neworks is provided in [3. The nework model used in his paper falls ino he specrum overlay approach o specrum sharing. III. AXIU THROUGHPUT OBJECTIVE Le r n denoe he ime average rae of admied daa for secondary user n, i.e., r n = lim R n (τ) Le r = (r,..., r N ) denoe he vecor of hese ime average raes. We define he following collision variables for each primary user m,..., : if here was a collision wih he primary user C m () = in channel m in slo 0 else Le c m denoe he ime average rae of collision for primary user m, i.e., c m = lim C m (τ) Le θ,..., θ N be a collecion of posiive weighs. Then he conrol objecive is o design a flow conrol and scheduling policy ha yields ime average rae r ha solves he following opimizaion problem: aximize: N θ n r n Subjec o: 0 r n λ n n,..., N c m ρ m m,..., r Λ (8) Here, Λ represens he nework capaciy region for he nework model as described above. I is defined as he se of all inpu rae vecors λ = (λ,..., λ N ) of he secondary users for which a scheduling sraegy exiss ha can suppor λ (wihou flow conrol) subjec o he consrains imposed by he nework, including he c m ρ m consrain. The noion of nework capaciy for general neworks wih ime varying channels and energy consrains is formalized in [, [3, [4 where i is shown o be a funcion of he seady sae nework opology disribuion, channel probabiliies, and ime average ransmission raes. Le r = (r,..., rn ) denoe he opimal soluion o he opimizaion problem defined above. In principle, i can be solved if all sysem parameers, including he se Λ are known in advance. However, in pracice, his region may no be known o he nework conroller (e.g., because he mobiliy paerns of he secondary users are unknown) and he above maximizaion problem mus be done for inpu raes eiher inside or ouside of he capaciy region. Even if all sysem parameers are known, he opimal soluion may be difficul o implemen as i may require cenralized coordinaion among all users. We nex presen an online conrol algorihm ha overcomes all of hese challenges. IV. OPTIAL CONTROL ALGORITH We now presen he Cogniive Nework Conrol Algorihm (CNC), a cross-layer conrol sraegy ha can be shown o achieve he opimal soluion r o he nework opimizaion problem presened earlier. I operaes wihou knowledge of wheher he inpu rae is wihin or ouside of he capaciy region Λ. Furher, i provides deerminisic wors case bounds on he maximum backlog a all imes and he maximum number of collisions wih a primary user in a given ime inerval. These are much sronger han probabilisic performance guaranees.

4 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL Finally, i offers a conrol parameer V 0 ha enables an explici rade-off beween he average hroughpu uiliy and delay. This algorihm is similar in spiri o he backpressure algorihms proposed in [3, [5 for problems of energy opimal neworking in wireless ad-hoc and mesh neworks. The algorihm is decoupled ino wo separae componens. The firs componen performs opimal flow conrol a he ranspor layers and is implemened independenly a each secondary user. The second componen deermines a nework wide resource allocaion every slo and needs o be solved collecively by he secondary users. In addiion o he acual queue backlog U n (), his algorihm uses a se of collision queues X m () for each channel m. These queues are virual in ha hey are mainained purely in sofware. These are used o rack he amoun by which he number of collisions suffered by a primary user m exceeds is ime average collision consrain rae ρ m. We assume ha he secondary users are aware of he X m () value for each channel m ha hey can access a ime. We define he collision queue X m () for channel m as follows: X m ( + ) = max[x m () ρ m, 0 + C m () (9) where C m () is he collision variable for channel m as defined in he previous secion. The above equaion represens he queueing dynamics of a single server sysem wih inpu process C m () and consan service rae ρ m. This sysem is rae sable if and only if he service rae exceeds he inpu rae, i.e., c m = lim C m (τ) ρ m This is precisely he collision consrain in he uiliy opimizaion problem saed earlier. Thus, if our policy sabilizes all collision queues as defined above, he maximum average rae of collisions will mee he required consrain. This echnique of urning ime average consrains ino queueing sabiliy problems was inroduced in [3 where i was used for saisfying average power consrains. A. Cogniive Nework Conrol Algorihm (CNC) Le he flow conrol and resource allocaion decision under he CNC algorihm be Rn CNC () and µ CNC nm () respecively. These are deermined according o he following algorihm, which uses a consan parameer V 0 o affec a performance/delay radeoff: Flow Conrol: A each secondary user n, choose he number of packes o admi Rn CNC () as he soluion o he following problem: inimize: R n ()[U n () V θ n Subjec o: 0 R n () A n () (0) This problem has a simple hreshold-based soluion. In paricular, if he curren queue backlog U n () > V θ n, hen Rn CNC () = 0 and no new packes are admied. Else, if U n () V θ n, hen Rn CNC () = A n () and all new packes are admied. Noe ha his can be solved separaely a each user and does no require knowledge of θ n weighs of oher users. Resource Allocaion: Choose a resource allocaion µ CNC nm () ha solves he following problem: ax: [ µ nm () U n ()P m () X k ()( P k ())Inm k k= Subjec o: consrains (), (3), (4), (5) () Afer observing he oucome of his allocaion a he end of he slo, he virual queues are updaed as in (9) based on he feedback received abou a collision wih a primary user or a successful ransmission. Noe ha only collisions wih a primary user affec (9), collisions beween secondary users do no affec he virual collision queues. The above problem is a generalized aximum Weigh ach ( problem where he weigh for a pair (n, m) is given by U n ()P m () ) k= X k()( P k ())Inm k. This is he difference beween he curren queue backlog U n () weighed by he probabiliy ha primary user m is idle and he weighed sum of all collision queue backlogs for he channels ha user n inerferes wih if i uses channel m. The weigh for a collision queue is he probabiliy ha he corresponding primary user will ransmi. Noe ha if his difference is nonposiive, hen he link (n, m) can be removed from he decision opions, simplifying scheduling. This problem is hard o solve in general, hough consan facor approximaions exis ha are easier o implemen. We discuss hese in Sec.VI. For he case when all channels are orhogonal from he poin of view of secondary users (which means a secondary user ransmission on a channel does no cause inerference o oher channels), I nm = m so ha Inm m =, Inm k = 0 k m. Then he above maximizaion simplifies o he following problem: aximize: [ µ nm () U n ()P m () X m ()( P m ()) Subjec o: consrains (), (3), (4), (7) () The above maximizaion requires solving he aximum Weigh ach (W) problem on an N biparie graph of N secondary users and channels. This problem can be solved in polynomial ime, hough his may require cenralized conrol. We discuss simpler consan facor approximaions in Sec.VI. Also, we consider a cell pariioned nework in he simulaions of Sec.VII for which a full maximum weigh mach can be implemened in a disribued manner. B. Performance Analysis We now characerize he performance of he CN C algorihm. This holds for general secondary user mobiliy processes ha are described by finie sae ergodic arkov Chains. Theorem : (Algorihm Performance) Assume ha all queues are iniialized o 0. Suppose all arrivals A n () are upper bounded so ha A n () A max for all n,. Also suppose he H() and P() processes are arkovian and have a well defined seady sae disribuion. Then, implemening

5 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL he CNC algorihm every slo for any fixed conrol parameer V 0 sabilizes all real and virual queues (hereby saisfying he maximum ime average collision consrains) and yields he following performance bounds: ) The wors case queue backlog for all secondary users n is upper bounded by a finie consan U max for all : U n () U max = V θ max + A max (3) where θ max = max n,...,n θ n ) For all m, such ha P m (), le ɛ > 0 be such ha P m () ɛ. 3 Then, he wors case collision queue backlog for all channels m is upper bounded by a finie consan X max : ( ɛ) X m () X max= U max + (4) ɛ Furher, he wors case number of collisions suffered by any primary user m is no more han ρ m T + X max over any finie inerval of size T. 3) The ime average hroughpu uiliy achieved by he CNC algorihm is wihin B/V of he opimal value: lim inf N θ n E R n (τ) N θ n r n B V (5) where B = B + C U + C X + N log V + log V and where B, C U, C X are consans (defined precisely in (8), (8), (9)). The consans C U and C X are deermined by he sochasics of he mobiliy and channel sae probabiliy processes and i is shown in Appendix A ha hese are O(log V ) when hese processes evolve according o a finie sae ergodic arkov model. Therefore, by par (3) of he heorem, he achieved average hroughpu uiliy is wihin O(log V/V ) of he opimal value. This can be pushed arbirarily close o he opimal value by increasing he conrol parameer V. However, his increases he maximum queue backlog bound U max linearly in V, leading o a uiliy-delay radeoff. The above bounds are quie srong. In paricular, he maximum collisions bound in par () gives deerminisic performance guaranees ha hold for any inerval size. This is quie useful in he conex of cogniive neworks since i implies ha he licensed users are guaraneed o suffer a mos hese many collisions. Probabilisic guaranees (e.g. [8) do no provide such bounds. We nex prove he firs wo pars of Theorem. Proof of par (3) uses he echnique of Sochasic Lyapunov Opimizaion and is provided in he nex secion. Proof of par (): Suppose ha U n () U max for all n,..., N for some ime. This is rue for = 0 as all queues are iniialized o 0. We show ha he same holds for ime +. We have cases. If U n () U max A max, hen from (), we have U n (+) U max (because R n () A max for all ). Else, if U n () > U max A max, hen U n () > V θ n + A max A max = V θ n. Then, he flow conrol par of he algorihm chooses R n () = 0, so ha by (): U n ( + ) U n () U max 3 Such an ɛ exiss for any finie sae ergodic arkov Chain. This proves (3). Proof of par (): Suppose ha X m () X max for all m,..., for some ime. This is rue for = 0 as all queues are iniialized o 0. We show ha he same holds for ime +. Firs suppose P m () =. Then, by definiion, here is no collision wih he primary user in channel m in slo so ha C m () = 0. Then, from (9), we have X m ( + ) X max. Nex, suppose P m () <. We again have cases. If X m () X max, hen from (9), we have X m ( + ) X max (because C m () ( ɛ) for all ). Else, if X m () > X max = U max ɛ, hen X m ()ɛ > U max ( ɛ). This implies X m ()( P m ()) > U max P m () U n ()P m () for all n,..., N. Then, since Inm m =, he resource allocaion par of he algorihm chooses µ nm () = 0 for all n. This would yield C m () = 0 (since no collision akes place wih primary user m), so ha by (9): This proves (4). X m ( + ) X m () X max V. STOCHASTIC LYAPUNOV OPTIIZATION Le Q() = (Q (),..., Q K ()) be a vecor process of queue lenghs for a discree ime sochasic queueing nework wih K queues (possibly including some virual queues like he collision queues defined in he previous subsecion). Le L(Q) be any non-negaive scalar valued funcion of he queue lenghs, called a Lyapunov funcion. Define he Lyapunov drif () as follows: () =E L(Q( + )) L(Q()) Suppose he nework accumulaes rewards every imeslo (where rewards migh correspond o uiliy measures of conrol acions). Assume rewards are real valued and bounded, and le he sochasic process f() represen he reward earned during slo. Le f represen he arge reward. The following resul (a varian of relaed resuls from [3, [4) specifies a drif condiion which ensures ha he ime average of he reward process f() is close o meeing or exceeding f. Theorem : (Delayed Lyapunov Opimizaion wih Rewards) Suppose here exis finie consans V > 0, B > 0, d > 0, and a non-negaive funcion L(Q) such ha E L(Q(d)) < and for every imeslo > d, he Lyapunov drif saisfies: hen we have: () V E f() B V f (6) lim inf E f(τ) f B V Proof : Inequaliy (6) holds for all > d. Summing boh sides over τ d,..., yields: E L(Q()) E L(Q(d)) B( d) V ( d)f + V E f(τ) τ=d

6 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL Rearranging erms, dividing by, and using non-negaiviy of L(Q) yields: ( d)f ( d)b V E L(Q(d) V E f(τ) The resul follows by aking limi as. We now use Theorem o prove par (3) of Theorem. This is done by comparing he Lyapunov drif of he CNC algorihm wih ha of a saionary randomized algorihm ha makes conrol decisions every slo purely as a funcion of he curren channel sae informaion P() and H(). We firs obain an expression for he Lyapunov drif under any conrol policy for our cogniive nework model. A. Lyapunov Drif Le Q() = (U (),..., U N (), X (),..., X ()) represen he collecion of all real and virual queue backlogs in he cogniive nework. We define he following Lyapunov funcion: L(Q()) = [ N U n() + Xm() Using queueing dynamics () and (9), he Lyapunov drif () under any conrol policy (including CN C) can be compued as follows: () B E where N ( U n () E X m ()(ρ m C m ()) ) µ nm ()S m () R n () B = N(A max + ) + ρ m + (7) (8) The collision variable C m () can be expressed in erms of he conrol decisions µ ij () and channel sae S() as follows: C m () = N i= j= µ ij ()Iij m [Ui()>0( S m ()) (9) where [Ui()>0 is an indicaor variable of non-zero queue backlog in secondary user i. This follows by observing ha a collision wih he primary user occurs in channel m if he primary user is busy (i.e. S m () = 0) and if µ ij () = for some secondary user i wih non-zero backlog using channel j ha inerferes wih channel m. We will find i useful o define he following relaed variable: Ĉ m () = N i= j= µ ij ()Iij m ( S m ()) (0) For a given conrol parameer V 0, we subrac he N reward meric V E θ nr n () from boh sides of he drif inequaliy (7) and use he fac ha Ĉ m () C m () o ge he following: N () V E θ n R n () B N ( E U n () ) µ nm ()S m () R n () N E X m ()(ρ m Ĉm()) V E θ n R n () B. Opimal Saionary, Randomized Policy () We now describe he saionary, randomized policy ha chooses conrol acions only as a funcion of P() and H() every slo. We have he following fac: Opimal Saionary, Randomized Policy: For any rae vecor (λ,..., λ N ) (inside or ouside of he nework capaciy region Λ), here exiss a saionary randomized scheduling policy n nm ha chooses feasible allocaions R (), µ () for all n,..., N, m,..., every slo as a funcion of he channel sae informaion P() and H() and yields he following seady sae values: E Rn () = rn () µ n = lim E µ nm (τ)s m (τ) rn (3) ĉ m = lim Ĉ E m () ρ m (4) where r = (r,..., rn ) is he soluion o (8). The above fac can be proven using echniques similar o he ones used in [, [3 for showing he exisence of capaciy achieving saionary, randomized policies ha make conrol decisions independen of queue backlog and is omied for breviy. We now prove an imporan propery of he CN C algorihm. Claim: Given a paricular backlog Q(), he CNC algorihm minimizes he righ side of inequaliy () over all alernae feasible decisions ha can be made on slo, including he saionary, randomized policy. Proof : By changing he order of summaions and using (0), he righ side of () can be expressed in a more convenien form: N B ρ m E X m () + E R n ()(U n () V θ n ) E [ µ nm () U n ()S m () X k ()( S k ())Inm k k= (5) Noe ha E S m () S( ) = P m () m. By wriing he las wo erms on he righ hand side as an ieraed expecaion by condiioning on he queue backlog and S( ), i can be seen ha CNC chooses conrol decisions (0) and () ha minimize hese erms for every possible value of he backlog and S( ), so ha he acual expecaion is also minimized.

7 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL Using his fac, we have: N CNC () B N ( E U n () µ E X m ()(ρ m N V E θ n R n nm () ()S m () R AT ĈST m ()) n ) () (6) x y Primary User Secondary User In Appendix A, we show ha his can be expressed as: N CNC () B N V θ n rn This is in a form ha fis (6). Thus, applying Theorem proves (5). VI. DISTRIBUTED IPLEENTATION Here we discuss consan facor approximaions o he resource allocaion problem () ha are easier o implemen in a disribued nework. We focus on he orhogonal channel case in which a secondary user ransmission on a channel does no cause inerference o oher channels. As noed earlier, in his case, he resource allocaion problem () reduces o a aximum Weigh ach (W) problem on an N biparie graph beween N secondary users and channels. An edge exiss beween nodes n and m of his graph if h nm () =, i.e., if secondary user n can access channel m in slo. The weigh of his edge is given by (U n ()P m () X m ()( P m ())). While he W problem can be solved in polynomial ime in a cenralized way, here we are ineresed in simpler implemenaions. In paricular, we use he idea of Greedy aximal ach Scheduling. Similar maximal scheduling algorihms for consan facor guaranees have been invesigaed in recen works [6 [8, [4. A maximal mach is defined as any se of edges (m, n) ha do no inerfere wih each oher such ha adding any new edge o his se necessarily violaes a maching consrain. A Greedy aximal ach can be achieved as follows: Firs selec he edge (m, n) wih he larges posiive weigh and label i acive. Then selec he edge wih he second larges posiive weigh (breaking ies arbirarily) ha does no conflic wih an acive edge and label i acive. Coninue in he same way, unil no more edges can be added. I is no difficul o see ha his final se of edges labeled acive has he desired maximal propery. A Greedy aximal ach can be compued wih much less overhead as compared o he aximum Weigh ach. I can be shown ha using maximal maches insead of he maximum weigh mach every slo can sill suppor any rae wihin Λ. In paricular, resource allocaions µg mn () chosen according o greedy maximal maches have he following Fig.. Example cell-pariioned nework used in simulaion propery: [ µ G nm () U n ()P m () X m ()( P m ()) [ µ nm() U n ()P m () X m ()( P m ()) where µ nm() is he opimal soluion o (). Using his, we have he following resul: Theorem 3: (Performance Bound for Orhogonal Channels wih Greedy aximal ach Scheduling) The ime average hroughpu uiliy achieved by he CN C algorihm wih Greedy aximal ach Scheduling is wihin B/V of N θ nr n: lim inf N θ n E R n (τ) VII. SIULATIONS N θ n r n B V We simulae he CNC algorihm on an example cogniive nework consising of 9 primary users and 8 secondary users as shown in Fig.. We consider a simple cell-pariioned nework wih one primary user per cell. The primary users are saic and each has is own licensed channel ha can be used by hem simulaneously. A secondary user can only aemp o ransmi on he channel associaed wih he primary user in is curren cell. The secondary users move from one cell o anoher according o a arkovian random walk. In paricular, a he end of every slo, a secondary user decides o say in is curren cell wih probabiliy β, else decides o move o an adjacen cell wih probabiliy β/4 (where β = 0.5 for he simulaions). If here is no feasible adjacen cell (e.g, if he previous cell is a corner cell and he new chosen cell does no exis), hen he user remains in he curren cell. I can be shown ha he resuling H() process forms an irreducible, aperiodic arkov Chain where he seady sae locaion disribuion is uniform over all cells. The channel sae process S m () for each primary user m is governed by an ON/OFF arkov chain wih symmeric

8 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL Toal Average Congesion (log scale) V = V = V = 5 V = 0 V = 00 No Flow Conrol Throughpu (packes/slo) V = V = V = 5 V = 0 V = Inpu Rae (packes/slo) Inpu Rae (packes/slo) Fig. 3. Toal average congesion vs. inpu rae for differen values of V Fig. 4. Achieved hroughpu vs. inpu rae for differen values of V ransiion probabiliies beween he ON and OFF saes given by 0. m. The maximum collision rae ρ m = 0. m. New packes arrive a he secondary users according o independen Bernoulli processes, so ha a single packe arrives i.i.d. wih probabiliy λ every slo. We assume here are no ranspor layer sorage buffers, so ha all packes ha are no immediaely admied o he nework layer are necessarily dropped. Flow conrol is performed according o (0) (wih θ n = n) and resource allocaion decisions are made every slo according o (). In his paricular cell-pariioned nework srucure wih one channel per cell, he maximum weigh mach is same as a greedy maximal mach so ha he opimal resource allocaion decisions can be made in each cell in a disribued way. In Fig. 3 we plo he average oal occupancy (summing all packes in he queues of he secondary users) versus he inpu rae λ. Each daa poin represens a simulaion over 500, 000 imeslos, and he differen curves correspond o values of he flow conrol parameer V,, 5, 0, 00, and he case V = (no flow conrol) is also shown. In his case, he average oal occupancy increases wihou bound as he inpu rae approaches nework capaciy. The verical asympoe which appears roughly a λ = 0.85 packes/slo corresponds o his value. Fig. 4 illusraes he achieved hroughpu versus he raw daa inpu rae λ for various V parameers. The achieved hroughpu is almos idenical o he inpu rae λ for small values of λ, and he hroughpu sauraes a a value ha depends on V, being very close o he 0.85 capaciy level when V is large. Also, i was found ha all real and virual queue backlogs are always bounded by he maximum values given in (3) and (4). In paricular, ɛ = 0. for his nework, so ha X m () X max = 4U max + = 4V + 5. Finally, he maximum average rae of collisions was very close o he arge ρ m = 0.. APPENDIX A DELAYED LYAPUNOV DRIFT Here, we use delayed queue backlogs o express he Lyapunov drif of he CN C algorihm in a form ha fis (6). Throughou, he backlogs U n () and X m () represen he values associaed wih he CN C algorihm. Recall ha Rn () and µ nm () denoe he resource allocaion decisions under he saionary, randomized policy inro- duced in Sec. V. We use he following sample pah inequaliies. Specifically, for all > d, we have for each secondary user queue U n () and for each collision queue X m (): U n ( d) + da max U n () U n ( d) d X m ( d) + d X m () X m ( d) dρ m These follow by noing ha he queue backlog a ime canno be smaller han he queue backlog a ime ( d) minus he maximum possible deparures in duraion ( d, d). Similarly, i canno be larger han he queue backlog a ime ( d) plus he maximum possible arrivals in duraion ( d, d). Using hese in (6) and using E Rn () = rn (from ()), we ge: N CNC () B + C U + C X N ( E U n ( d) µ E X m ( d)(ρ m nm where C U and C X are given by: ()S m () R AT ĈST m ()) V n ) () N θ n rn (7) C U = dn + da maxn (8) C X= d ( + ρ m) (9) Using ieraed expecaions, we have he following: N E U n ( d) µ ()S m () = nm N E U n ( d) E E X m ( E AT d)ĉst m () µ nm ()S m () T ( d) = (30) X m ( d) E Ĉ m () T ( d) (3)

9 PROC. IEEE INFOCO, PHOENIX, AZ, APRIL where T ( d) = (H( d), S( d), Q( d)) represens he composie sysem sae a ime ( d) and includes he opology sae and queue backlogs. By he arkovian propery of he H(), S() (and herefore P()) processes, any funcionals of hese saes converge exponenially fas o heir seady sae values (his is formalized in Appendix B). Since he policy makes conrol decisions only as a funcion of P() and H(), he resuling allocaions are funcionals of hese arkovian processes. Thus, here exis posiive consans α, α and 0 < γ, γ < such ha: E µ nm ()S m () T ( d) µ n α γ d ĈST AT E m () T ( d) ĉm + α γ d where µ n, ĉm are he seady sae values as defined in (3), (4). Using hese, he above can be wrien as: E µ nm ()S m () T ( d) rn α γ d (3) ĈST AT E m () T ( d) ρ m + α γ d (33) Thus, using (3), (33) in (30), (3), inequaliy (7) can be expressed as: N CNC () B + C U + C X N + E U n ( d)α γ d + E X m ( d)α γ d V N θ n rn B + C U + C X + NU max α γ d + X max α γ d V N θ n rn (using (3), (4)) Define d = log(αumax) log(/γ ), d = log(αxmax) log(/γ ). Then choosing d = max(d, d ), we have: N CNC () B + C U + C X + N log V + log V V N θ n rn (34) Since U max and X max are O(V ), we have d O(log V ). APPENDIX B CONVERGENCE OF ARKOV CHAINS Le Z() be a finie sae, discree ime ergodic arkov chain. Le S denoe is sae space and le π i i S be he seady sae probabiliy disribuion. Then, for all inegers d 0, here exis consans α, γ such ha: P rz() = j Z( d) = i π j αγ d (35) where α 0 and 0 < γ <. This implies ha he arkov chain converges o is seady sae probabiliy disribuion exponenially fas (see [9). Le f(z()) be a posiive random funcion of Z() (negaive case can be reaed similarly). Define f = j S π jm j where m j = E f(z()) Z() = j. Then: E f(z()) Z( d) = i = j S E f(z()) Z() = j P rz() = j Z( d) = i j S m j (π j + αγ d ) (using (35)) f + sm max αγ d where m max = max j S m j and s = cards. This shows ha funcionals of he saes of a finie sae ergodic arkov chain converge o heir seady sae value exponenially fas. REFERENCES [ J. iola. Sofware radios: wireless archiecure for he s cenury. John Wiley & Sons Inc, 000. [ I. F. Akyildiz, W-Y Lee,. C. Vuran and S. ohany. NeX generaion/dynamic specrum access/cogniive radio wireless neworks: a survey. Compuer Neworks Journal (Elsevier), vol. 50, pp. 7-59, Sep [3 Q. Zhao and B.. Sadler. A survey of dynamic specrum access: signal processing, neworking, and regulaory policy. IEEE Signal Processing agazine, vol. 4, no. 3, pp , ay, 007. [4 C. Peng, H. Zheng, and B.Y. Zhao. Uilizaion and fairness in specrum assignmen for opporunisic specrum access. AC/Springer ONET, vol., issue 4, Aug [5 L. Cao, H. Zheng. Disribued specrum allocaion via local bargaining. Proc. IEEE SECON, Sep [6 W. Wang, X. Liu and H. Xiao. Exploring opporunisic specrum availabiliy in wireless communicaion neworks. Proc. of IEEE VTC, Sep [7 Y.T. Hou, Y. Shi, and H.D. Sherali. Opimal specrum sharing for mulihop sofware defined radio neworks. Proc. IEEE INFOCO, ay 007. [8 Y. Chen, Q. Zhao, and A. Swami. Join design and separaion principle for opporunisic specrum access. under submission o IEEE Transacions on Informaion Theory [9 Q. Zhao, L. Tong, A. Swami, and Y. Chen. Decenralized cogniive AC for opporunisic specrum access in ad hoc neworks: A PODP framework. IEEE Journal on Seleced Areas in Communicaions (JSAC): Special Issue on Adapive, Specrum Agile and Cogniive Wireles Neworks, vol. 5, no. 3, pp , April, 007. [0 L. Tassiulas and A. Ephremides. Sabiliy properies of consrained queueing sysems and scheduling policies for maximum hroughpu in mulihop radio neworks. IEEE Transacaions on Auomaic Conrol, vol. 37, no., pp , Dec. 99. [. J. Neely, E. odiano, and C. E Rohrs. Dynamic power allocaion and rouing for ime varying wireless neworks. IEEE Journal on Seleced Areas in Communicaions, vol. 3, no., pp , Jan [. J. Neely, E. odiano, and C. Li. Fairness and opimal sochasic conrol for heerogeneous neworks. Proc. IEEE INFOCO, arch 005. [3. J. Neely. Energy opimal conrol for ime varying wireless neworks. IEEE Transacions on Informaion Theory, vol. 5, no. 7, July 006. [4 L. Georgiadis,. J. Neely, L. Tassiulas. Resource allocaion and crosslayer conrol in wireless neworks. Foundaions and Trends in Neworking, Vol., no., pp. -49, 006. [5. J. Neely and R. Urgaonkar. Cross layer adapive conrol for wireless mesh neworks. Ad Hoc Neworks (Elsevier), vol. 5, no. 6, pp , Aug [6 X. Lin and N. B. Shroff. The impac of imperfec scheduling on crosslayer rae conrol in wireless neworks. Proc. IEEE INFOCO, arch 005. [7 P. Chaporkar, K. Kar, and S. Sarkar. Throughpu guaranees hrough maximal scheduling in wireless neworks. Proc. of 43rd Alleron Conference, Sep [8 X. Wu and R. Srikan. Bounds on he capaciy region of mulihop wireless neworks under disribued greedy scheduling. Proc. IEEE INFOCO, April 006. [9 S. Ross. Sochasic Processes. John Wiley & Sons, Inc., New York, 996.

Opportunism, Backpressure, and Stochastic Optimization with the Wireless Broadcast Advantage

Opportunism, Backpressure, and Stochastic Optimization with the Wireless Broadcast Advantage ASILOMAR CONFERENCE ON SIGNALS, SYSTEMS, AND COMPUTERS, PACIFIC GROVE, CA, OCTOBER 2008 1 Opporunism, Backpressure, and Sochasic Opimizaion wih he Wireless Broadcas Advanage Michael J. Neely, Rahul Urgaonkar

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

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

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

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

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

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

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

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

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

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

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

Network Utility Maximization over Partially Observable Markovian Channels

Network Utility Maximization over Partially Observable Markovian Channels WIOPT 2011 1 Nework Uiliy Maximizaion over Parially Observable Markovian Channels Chih-ping Li, Suden Member, IEEE and Michael J. Neely, Senior Member, IEEE Absrac This paper considers maximizing hroughpu

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

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

Stable Scheduling Policies for Maximizing Throughput in Generalized Constrained Queueing Systems

Stable Scheduling Policies for Maximizing Throughput in Generalized Constrained Queueing Systems 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

More information

Delay Reduction via Lagrange Multipliers in Stochastic Network Optimization

Delay Reduction via Lagrange Multipliers in Stochastic Network Optimization IEEE TRANS. ON AUTOMATIC CONTROL, OL. 56, ISSUE 4, PP. 842-857, APRIL 2011 1 Delay Reducion via Lagrange Mulipliers in Sochasic Nework Opimizaion Longbo Huang, Michael J. Neely Absrac In his paper, we

More information

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3 Macroeconomic Theory Ph.D. Qualifying Examinaion Fall 2005 Comprehensive Examinaion UCLA Dep. of Economics You have 4 hours o complee he exam. There are hree pars o he exam. Answer all pars. Each par has

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

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

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

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

Energy Storage Benchmark Problems

Energy Storage Benchmark Problems Energy Sorage Benchmark Problems Daniel F. Salas 1,3, Warren B. Powell 2,3 1 Deparmen of Chemical & Biological Engineering 2 Deparmen of Operaions Research & Financial Engineering 3 Princeon Laboraory

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

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

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

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

RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY

RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY ECO 504 Spring 2006 Chris Sims RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY 1. INTRODUCTION Lagrange muliplier mehods are sandard fare in elemenary calculus courses, and hey play a cenral role in economic

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

Some Ramsey results for the n-cube

Some Ramsey results for the n-cube Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain

More information

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.

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

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

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

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

Optima and Equilibria for Traffic Flow on a Network

Optima and Equilibria for Traffic Flow on a Network Opima and Equilibria for Traffic Flow on a Nework Albero Bressan Deparmen of Mahemaics, Penn Sae Universiy bressan@mah.psu.edu Albero Bressan (Penn Sae) Opima and equilibria for raffic flow 1 / 1 A Traffic

More information

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne

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

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

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

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

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

Appendix to Creating Work Breaks From Available Idleness

Appendix to Creating Work Breaks From Available Idleness Appendix o Creaing Work Breaks From Available Idleness Xu Sun and Ward Whi Deparmen of Indusrial Engineering and Operaions Research, Columbia Universiy, New York, NY, 127; {xs2235,ww24}@columbia.edu Sepember

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

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

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

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

arxiv: v1 [math.oc] 11 Sep 2017

arxiv: v1 [math.oc] 11 Sep 2017 Online Learning in Weakly Coupled Markov Decision Processes: A Convergence ime Sudy Xiaohan Wei, Hao Yu and Michael J. Neely arxiv:709.03465v [mah.oc] Sep 07. Inroducion Absrac: We consider muliple parallel

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

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

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

Optimal Power Cost Management Using Stored Energy in Data Centers

Optimal Power Cost Management Using Stored Energy in Data Centers Opimal Power Cos Managemen Using Sored Energy in Daa Ceners Rahul Urgaonkar, Bhuvan Urgaonkar, Michael J. Neely, Anand Sivasubramaniam Dep. of EE, Dep. of CSE Universiy of Souhern California, The Pennsylvania

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

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

Upper confidence bound based decision making strategies and dynamic spectrum access

Upper confidence bound based decision making strategies and dynamic spectrum access Upper confidence bound based decision making sraegies and dynamic specrum access Wassim Jouini wassim.jouini@supelec.fr Damien Erns Universiy of Liège derns@ulg.ac.be Chrisophe Moy chrisophe.moy@supelec.fr

More information

Longest Common Prefixes

Longest Common Prefixes Longes Common Prefixes The sandard ordering for srings is he lexicographical order. I is induced by an order over he alphabe. We will use he same symbols (,

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

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

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

Control of computer chip semi-conductor wafer fabs

Control of computer chip semi-conductor wafer fabs Conrol of compuer chip semi-conducor wafer fabs Cos: x9 $ Reurn: years Cycle ime: 6 weeks WIP: 6, wafers, 8x6 $ Challenge: Conrol he queues a ~5 work seps Conrol of Manufacuring Sysems:! Sochasic racking

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

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,

More information

Online Convex Optimization Example And Follow-The-Leader

Online Convex Optimization Example And Follow-The-Leader CSE599s, Spring 2014, Online Learning Lecure 2-04/03/2014 Online Convex Opimizaion Example And Follow-The-Leader Lecurer: Brendan McMahan Scribe: Sephen Joe Jonany 1 Review of Online Convex Opimizaion

More information

Block Diagram of a DCS in 411

Block Diagram of a DCS in 411 Informaion source Forma A/D From oher sources Pulse modu. Muliplex Bandpass modu. X M h: channel impulse response m i g i s i Digial inpu Digial oupu iming and synchronizaion Digial baseband/ bandpass

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

Scheduling of Crude Oil Movements at Refinery Front-end

Scheduling of Crude Oil Movements at Refinery Front-end Scheduling of Crude Oil Movemens a Refinery Fron-end Ramkumar Karuppiah and Ignacio Grossmann Carnegie Mellon Universiy ExxonMobil Case Sudy: Dr. Kevin Furman Enerprise-wide Opimizaion Projec March 15,

More information

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j = 1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

. Now define y j = log x j, and solve the iteration.

. Now define y j = log x j, and solve the iteration. Problem 1: (Disribued Resource Allocaion (ALOHA!)) (Adaped from M& U, Problem 5.11) In his problem, we sudy a simple disribued proocol for allocaing agens o shared resources, wherein agens conend for resources

More information

Energy-Aware Wireless Scheduling with Near Optimal Backlog and Convergence Time Tradeoffs

Energy-Aware Wireless Scheduling with Near Optimal Backlog and Convergence Time Tradeoffs IEEE/ACM TRANSACTIONS ON NETWORKING, VOL. 24, NO. 4, AUG. 2016, PP. 2223-2236 1 Energy-Aware Wireless Scheduling wih Near Opimal Backlog and Convergence Time Tradeoffs Michael J. Neely Universiy of Souhern

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

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based

More information

Object tracking: Using HMMs to estimate the geographical location of fish

Object tracking: Using HMMs to estimate the geographical location of fish Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging

More information

The field of mathematics has made tremendous impact on the study of

The field of mathematics has made tremendous impact on the study of A Populaion Firing Rae Model of Reverberaory Aciviy in Neuronal Neworks Zofia Koscielniak Carnegie Mellon Universiy Menor: Dr. G. Bard Ermenrou Universiy of Pisburgh Inroducion: The field of mahemaics

More information

The Optimality of Two Prices: Maximizing Revenue in a Stochastic Network

The Optimality of Two Prices: Maximizing Revenue in a Stochastic Network PROC. OF 45TH ANNUAL ALLERTON CONFERENCE ON COMMUNICATION, CONTROL, AND COMPUTING (INVITED PAPER), SEPT. 2007 1 The Opimaliy of Two Prices: Maximizing Revenue in a Sochasic Newor Longbo Huang, Michael

More information

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1.

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1. Timed Circuis Asynchronous Circui Design Chris J. Myers Lecure 7: Timed Circuis Chaper 7 Previous mehods only use limied knowledge of delays. Very robus sysems, bu exremely conservaive. Large funcional

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

Simulating models with heterogeneous agents

Simulating models with heterogeneous agents Simulaing models wih heerogeneous agens Wouer J. Den Haan London School of Economics c by Wouer J. Den Haan Individual agen Subjec o employmen shocks (ε i, {0, 1}) Incomplee markes only way o save is hrough

More information

Exploiting Wireless Channel State Information for Throughput Maximization

Exploiting Wireless Channel State Information for Throughput Maximization Exploiing Wireless Channel ae Informaion for Throughpu Maximizaion V Tsibonis, L Georgiadis, L Tassiulas Absrac We consider he problem of scheduling packes over a number of channels wih ime varying conneciviy

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

Empirical Process Theory

Empirical Process Theory Empirical Process heory 4.384 ime Series Analysis, Fall 27 Reciaion by Paul Schrimpf Supplemenary o lecures given by Anna Mikusheva Ocober 7, 28 Reciaion 7 Empirical Process heory Le x be a real-valued

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

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

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

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

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE

INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE JAMES ALEXANDER, JONATHAN CUTLER, AND TIM MINK Absrac The enumeraion of independen ses in graphs wih various resricions has been a opic of much ineres

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

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

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

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion

More information

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Inventory Control of Perishable Items in a Two-Echelon Supply Chain Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan

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