Uplink capacity of self-organizing clustered orthogonal CDMA networks in flat fading channels

Size: px
Start display at page:

Download "Uplink capacity of self-organizing clustered orthogonal CDMA networks in flat fading channels"

Transcription

1 Auhor manuscrip, published in "ITW-29, Ialy 29)" Uplink capaciy of self-organizing clusered orhogonal CDMA neworks in fla fading channels Romain Couille ST-Ericsson, Supelec 55 Roue des Lucioles 656 Sophia Anipolis, France Mérouane Debbah Alcael-Lucen Chair, Supelec Plaeau de Moulon, 3 rue Jolio-Curie 992 Gif sur Yvee, France merouane.debbah@supelec.fr hal-4453, version - 3 Jan 2 Absrac In his paper, we derive a deerminisic equivalen of he Shannon-ransform of cerain ype of large uniary random marices. This approximaion is exploied o evaluae he uplink channel capaciy of clusered orhogonal CDMA nework. When non-uniform power allocaion among he users of each cluser is allowed, we derive an explici ieraive waerfilling algorihm which, upon convergence, achieves he muliuser decoding capaciy. In paricular, we show ha, in a selforganizing clusered orhogonal CDMA nework, each cluser can opimize is power allocaion policy independenly of he oher clusers a he expense of a small feedback overhead. Simulaions corroborae he heoreical derivaions. I. INTRODUCTION Previous releases of he universal mobile elecommunicaions sysem UMTS) aimed o ensure qualiy of service for all users in he cells. The recen developmens of 3G proocols however push more and more owards he maximizaion of he user sum raes in he cell. 3G sysems are based on code division muliple access CDMA) which use random i.i.d. codes in he uplink and incur inerference among he users. Orhogonal codes which reduce drasically he inerference) are unforunaely no used as i would require some imporan overhead signaling o synchronize he users. One way of dealing wih his issue, which is par of he recen cogniive radio incenive [2], proposes o allow users o self-organize in small clusers e.g. here CDMA-based clusers) or coaliions [], [2] o opimize heir uplink communicaions o a cenral base saion [3], [4]. Due o he very shor disance beween he users wihin he cluser, non-inerfering ransmissions can be performed hrough he use of orhogonal codes per cluser. Those neworks are referred o as self-organizing neworks, which we shall consider in he following. Among he pas conribuions in he muli-cell CDMA case, a performance sudy of muli-cell orhogonal CDMA neworks in frequency selecive channels was carried ou in [3]. The auhors sugges ha frequency seleciviy has a sronger negaive impac on CDMA performance han he pah loss; herefore muliple small CDMA cells migh urn ou more advanageous han a large single CDMA cell. The same conclusions were given in a similar uplink random CDMA sudy in [4]. In [5], he performance of single-cell CDMA in fla-fading channels wih random codes was analyzed for differen ypes of decoders when power opimizaion is performed a he ransmier. Mos of hese sudies derive resuls under he assumpion ha he number of users per cell is large, since simpler derivaions are possible using he heory of large random marices [6]. In his paper we derive resuls for he uplink capaciy of a clusered CDMA self-organizing nework when he CDMA code marices are large and isomeric. Those resuls are based on ools from random marix heory [6] and free probabiliy [7], [8]. We will in paricular derive wo heorems which relae respecively he η-ransform and he Shannon-ransform see definiions hereafer) of such a channel model o deerminisic approximaions independen of he code marix realizaion. Anoher deerminisic equivalen of he η-ransform has been considered recenly in he frequency selecive case by Peacock e al. [] using a differen approach, referred o as incremenal marix expansions. However, he form of heir deerminisic approximaions does no allow for a rivial calculus of he Shannon-ransform, which is provided here. Up o he auhors undersanding, boh our and Peacock s resuls are differen, hough hey are shown o end o he same asympoic expression. Our resul on he Shannon-ransform allows us firs o derive explici expressions for he uplink clusered orhogonal CDMA capaciy and more imporanly o produce an ieraive waer-filling algorihm for he cell users o achieve nework capaciy under per-cluser sum power consrain where he power wihin a cluser is seen as a long erm exchangeable uiliy, see [2] for more jusificaions). We will also propose an algorihm for each cluser o auomaically reconfigure heir power allocaions, independenly of he oher clusers, when a user connecs or disconnecs o he cell. The remainder of his paper unfolds as follows: in Secion II, he fundamenal mahemaical resuls of his paper are provided. In Secion III, he clusered CDMA self-organizing nework model is inroduced. In Secion IV, he capaciy achieving power allocaion algorihm is presened. In Secion V, numerical simulaions are carried ou. Finally, in Secion VI, we draw our conclusions. Noaion: Capial boldface characers denoe marices I N by isomeric, we mean a non-necessarily square marix wih orhogonal columns.

2 hal-4453, version - 3 Jan 2 is he N N ideniy marix). Hermiian ranspose is denoed ) H. The operaor dex represens he deerminan of he marix X. The funcion δx) is he Kronecker Dirac funcion such ha δx) = if x = and equals oherwise. The firs derivaive of he funcion f is denoed f. II. MATHEMATICAL PRELIMINARIES Our main resuls rely on ools of random marix heory analysis, he η-ransform and he Shannon-ransform, which we inroduce in he following. Definiion η-ransform): Le µ be a probabiliy densiy funcion. The η-ransform ηx) of µ is defined as ηx) = µd) ) + x Noe ha he η-ransform is linked o he beer known Sieljes ransform Sx) via he relaion S /x) = xηx). Theorem : Le K, N be posiive inegers, {W i } i=,...,k be K independen N N Haar disribued complex random marices and {D i } i=,...,k be K diagonal N N non-negaive marices. Denoe µ i he empirical disribuion of he enries of D i. Then, for N large, K fixed, and for some x, he η-ransform ηx) of he empirical eigenvalue disribuion of B N = K W kd k Wk H saisfies approximaely ηx) = + x β k x)) 2) where he funcions β k x), k {,...,K} saisfy he K fixed-poin equaions β k x) = xηx)β k x) + xηx) µ id) 3) Proof: The proof relies mosly on he asympoic freeness [6] of he marices W k D k Wk H when N grows o infiniy. In his case, he R-ransform Rx) of B N is asympoically equal o he sum of he R-ransform R k x) of W k D k Wk H. This R- ransform is defined wih respec o he η-ransform as 2 R xη) = x η ) 4) η R + ) = xr + 5) x Denoe η k and R k he η- and R-ransform of W k D k Wk H. Then, from ) and 5), we have xr k + = µ k d) R k + x which leads o R k x) = x R k x) + x µ kd) 6) and, in paricular, defining β k x) = R k xη), we have β k x) = xηβ k + xη µ kd) 7) 2 for readabiliy, we will use he shorcu noaion f = fx) when no confusion is possible. which is exacly 3). Now, since R = K R k, using 4), his gives ) ηx) = + x R k xη)) = + x β k which complees he proof. The η-ransform is of no direc pracical use in our curren analysis hough i migh be essenial in he analysis of MMSE decoders for CDMA downlink, see e.g. []). To esablish he capaciy of uplink ransmissions in CDMA neworks, we need a deerminisic approximaion of he so-called Shannonransform, defined as follows. Definiion 2 Shannon-ransform): Le µ be some probabiliy densiy funcion. We define he Shannon-ransform Vx) of µ as Vx) = log + x)µd) 9) The Shannon-ransform is inimaely linked o he capaciy of a muli-dimensional channel whose marix empirical eigenvalue densiy funcion is µ. Noe also ha he Shannon-ransform is an inegral of /x ηx)). From his observaion, equivalenly o Theorem, we derive in he following heorem a deerminisic approximaion for Vx). Theorem 2: Le B N be an N N marix as defined in he condiions of Theorem. We have he following asympoic approximaion of he Shannon-ransform Vx) of B N as ) Vx) =log + x β k x) 8) + log + xηx)[ β k x)]) µ k d) ) where β k x) and ηx) are given by 2) and 3). Proof: The Shannon-ransform saisfies [6] Vx) = x ηu))du u we herefore seek an inegral form for he η-ransform. Noe firs ha x η) = K x + x β k ) ) = β k η ) Also noe from 3) ha xηβ k = and herefore ha = xηβ k xηβ k + xη µ kd) 2) xηβ k + xη µ kd) 3) Now, for any k, he derivaive along x of C k = log xηβ k + xη)µ k d) is C k = [ ηβk xη β k xηβ k ] + [η + xη ] µ k d) 4) xηβ k + xη

3 hal-4453, version - 3 Jan 2 Fig.. λ λ 2 λ 22 λ 23 λ 3 λ 32 Self-organizing clusered CDMA nework Recalling now 3) and 3), his yields C k = ηβ k xη β k xηβ k ) +η+xη ) β k = xηβ k 5) We also have log + x β k )) = η β k + x ηβ k 6) Adding his las expression o k C k, we end up wih he desired K β kη. Verifying ha V) =, we finally obain ). A. Model III. SYSTEM MODEL Consider a nework composed of K clusers 3 and a single base saion, each cluser being composed of a maximum of N users each. The users of cluser k {,...,K} are assigned orhogonal spreading codes forming he columns of a uniary N N Haar disribued marix W k. Each cluser operaes independenly regarding he allocaion of spreading codes, and herefore he W k marices are considered independen. Depending on is posiion in cluser k, user n {,...,N} induces a pah-loss λ kn in he uplink ransmission o he base saion 4. The pah-losses in cluser k are gahered in he diagonal marix Λ k = diagλ k,...,λ kn ). Moreover we assume ha he ransmission channel is frequency fla, i.e. he propagaion channel fade is consan over he ransmission bandwidh. User n in cluser k ransmis wih power p kn, he collecion of which is gahered in he marices P k = diagp k,..., p kn ). We addiionally assume he base saion experiences a background noise disribued as CN, σ 2 ). This siuaion is depiced in Figure. B. Uplink capaciy Denoing D k = Λ k P k, he uplink per-ime sample) capaciy C reads ) Cσ 2 ) = N log de I N + σ 2 W k D k Wk H 7) 3 we assume ha he clusers are formed in a self-organized manner based on he disance, see [2] 4 when here are less han N users in a specific cluser, we can freely wrie λ kn = for he virual users. Denoing µ k ) = N N n= δ λ knp kn ), from Theorem 2, when N is large, we have approximaely 5 ) Cσ 2 ) = log + K σ 2 β k σ 2 ) + log + σ 2 ησ 2 )[ β k σ 2 )] ) µ k d) 8) wih ηx) and β k x) saisfying he implici Equaion 3). We herefore have an expression of he uplink capaciy which only depends on he pah losses λ kn, k {,..., K}, n {,..., N}. This suggess ha he choice of he spreading codes W k s has asympoically no impac on he achievable uplink sum rae. We will now seek a power allocaion sraegy which maximizes he capaciy under cluser power consrains. IV. SELF-ORGANIZED POWER ALLOCATION In he following, based on relaive fairness among clusers and on individual mobile erminal capabiliies, cluser k is allocaed a oal allowed ransmi P k ; herefore N n= p kn = P k. A. Opimal power allocaion To maximize he sysem capaciy wih respec o he KN scalars p kn, we shall firs show he following lemma, Lemma : The power allocaion policy p kn = p kn opimizing he deerminisic capaciy approximaion 8) saisfies, for all k, n, p kn α = k σ2 η βk ) + λ kn η 9) where η, β k are he respecive values of η and β k when C achieves is maximum, and α k is such ha k p kn = P k. Proof: To prove his lemma, we follow Prop. 4 in []. We need o ensure successively ha, denoing Cσ 2 ) = V p,..., p KN, β,..., β K, η), 6 i) V/ β k = V/ η =. ii) C is sricly concave in he p kn. Par ii) is in fac obvious once i) is shown. Indeed, in his case, he successive derivaives of C along p kn equal hose of V along p kn plus oher erms he derivaives of V along η and β k ) which are null. Therefore, since V is he sum of sricly concave funcions of he p kn, C is sricly concave. To show i), we simply differeniae and observe ha V η = σ 2 β k σ 2 ηβ k σ 2 + ησ 2 µ kd) = 2) k V = ησ 2 ησ 2 + β k ηβ k σ 2 + ησ 2 µ kd) = 2) boh equaliies semming from he remarks 3) and 3). 5 noe ha he capaciy is he Shannon-ransform of σ 2 and no of σ 2. 6 V can be seen as a funcion of N + )K + independen variables.

4 hal-4453, version - 3 Jan 2 From ii), C admis a unique maximum and from i), when β k and η have reached his maximum, he power allocaion of he p kn s reduces o maximizing K log + σ 2 β k )+ log + σ 2 η[ β k ] ) µ k d) 22) independenly of η and he β k s. This is, for each k, i reduces o maximizing log + σ 2 η[ β k ] ) µ k d), whose soluion is he waer-filling soluion 9). To achieve he capaciy opimizing power allocaion policy, similarly o [], we propose he following ieraive waerfilling algorihm, A iniializaion, for all k, p kn = P k N, η =, β k =. while he p kn s have no converged do for k {,...,K} do Se η, β k ) as soluion of 2), 3) for n =...,N do ) +, wih αk such ha Se p kn = α k σ2 ηβ k λ kn η n p kn = P k. end while This algorihm canno be proved o converge. However, from he proof of Lemma and following he argumens in [], upon convergence, he algorihm is proved o converge o he opimal soluion. The power allocaion algorihm allows now he base saion o indicae o all users he uplink ransmission powers, provided ha i is aware of all he λ kn s. We will show in he following ha such a cenralized approach is unnecessary o achieve opimal power allocaion; precisely, we will show ha each cluser k can independenly perform power allocaion regardless of he λ jn, j k, provided ha each cluser feedbacks a single real parameer o a neighboring cluser afer he individual opimizaion process. B. Self-organized capaciy opimizaion The main ineres of he sysem model under consideraion lies in he independence of every cluser in erms of code allocaion: each cluser can behave auonomously, wihou he need for he base saion o inervene in he code allocaion policy. Now, in order o achieve he opimal uplink capaciy hroughou he K clusers, user n in cluser k mus be allocaed a ransmi power p kn = p kn. Those p kn s can be compued explicily a he base saion or by any cenral eniy aware of all he λ kn s. However, his requires a large amoun of overhead daa o be ransferred back and forh hrough he nework; his comes a a non-negligible cos and goes agains he philosophy of self-organized neworks. I is acually possible o circumven his issue by performing successive local opimizaions in every cluser and by successively ransmiing a single parameer o he oher clusers: his parameer is he updaed evaluaion of he η- ransform of B N. More explicily, we consider he following algorihm, Capaciy Cσ 2 ) [bis/s/hz] Simulaion, uniform power allocaion Simulaion, Mone Carlo allocaion Theory, uniform power allocaion Theory, opimal power allocaion SNR [db] Fig. 2. Uplink capaciy of a K = 3-clusered orhogonal CDMA, uniform and opimal power allocaion, N = 8. For user n in cluser k, λ kn = α k [n ]/[N ]) db, α = db, α 2 = 2 db, α 3 = 3 db. A iniializaion, in cluser k, se p kn = P k N, η = and β k =. while η has no converged do for cluser k {,..., K} do Compue η as in 2). Se β k as soluion of 3), from he updaed η. for n =...,N do Se p kn = n p kn = P k. α k σ2 ηβ k λ kn η ) +, wih αk such ha Transmi he updaed value η o cluser k + wih convenion K + ). end while This algorihm only requires for he clusers o ransmi he updaed value of η a each sep. Noe ha compuing η requires only o know he updaed value of k β k. This value does no need be addiionally ransmied if we assume ha each cluser sores he las value of η and infers from i he updaed k β k. This way, when a new user regisers or a curren user disconnecs from a cluser, all he ransmied powers can be quickly updaed hroughou he nework a a low feedback cos. Moreover, as will be observed in Secion V, he convergence ime of he algorihm is in general so fas ha a single round of opimizaion over he K clusers is sufficien and herefore he ouer while loop is no needed) o achieve a saisfying uplink rae. In he nex secion, he heoreical resuls discussed so far are confroned o simulaions. V. SIMULATION AND RESULTS In his secion, we verify he heoreical formula derived in Theorem 2, he power allocaion algorihm derived from Lemma and he self-organizing successive power allocaion proposed in Secion IV-B. Throughou his secion, we

5 hal-4453, version - 3 Jan 2 Capaciy Cσ 2 ) [bis/s/hz] Uniform power allocaion Ieraive waer-filling, round Ieraive waer-filling, opimal Local power allocaion SNR [db] Fig. 3. Uplink capaciy of a K = 3-cluser orhogonal CDMA, ieraive waer-filling power allocaion, N = 8. For user n in cluser k, λ kn = α k [n ]/[N ]) db, α = db, α 2 = 2 db, α 3 = 3 db. consider a K = 3-cluser self-organizing orhogonal CDMA nework wih N = 8 users per cluser. The uplink signal from user n in cluser k experiences a pah loss λ kn = α k [n ]/[N ]) db, wih α = db, α 2 = 2 db, α 3 = 3 db; his allows o have in each cluser a linear range in db) of users wih relaive pah losses varying from db o α k ) db. In Figure 2, we presen simulaion resuls based on he averaging of Mone Carlo simulaions when uniform power allocaion is applied or when he allocaed powers are issued from random Mone Carlo rials bu saisfy he power consrains); for he laer, for each of he simulaions, we define he sysem capaciy as he maximal rae achieved over, random realizaions of vecors {p,...,p KN } ordered in such a way ha p kj < p k,j+. 7 Those are compared agains he soluions given by Theorem 2 wih uniform power allocaion and he ieraive waer-filling algorihm, respecively. We observe a perfec fi, even for his low N = 8 value, beween he simulaed and heoreical resuls in he uniform power allocaion case. As for he opimal power allocaion, we observe an almos perfec fi, which will surely be more accurae if one increases he number of random power allocaions of he Mone Carlo mehod. In Figure 3, we apply he self-organizing power allocaion algorihm under he same consrains as in Figure 3, when, or an infinie number of rounds a leas a sufficien number of rounds o ensure convergence of η) are performed). We observe surprisingly enough ha a single round is enough o already achieve he opimal capaciy limi of he complee algorihm. As a consequence, he algorihm is exremely fas a providing high performance power allocaion among he K clusers. For he sake of comparison, we also propose in Figure 3 he so-called local power allocaion policy applied 7 his allows o significanly reduce he number of required simulaions by discarding obviously wrong soluions. by each cluser as hough i were alone in he nework, which is maximized by he classical waer-filling soluion, i.e. wih powers ) + k σ2 p kn = α 23) λ kn where α k ensures ha n p kn = P k. We observe ha, while in he low SNR region, his local sraegy is reasonable, i is no longer he case in he high SNR regime. Since he cos of he ransmission of he daum η is obviously very low in his regime, he self-organizing power allocaion policy is even more ineresing. VI. CONCLUSION In his paper, we provided wo novel heorems relaing he uplink capaciy of a large clusered orhogonal frequency fla CDMA nework communicaing o a single base saion o deerminisic capaciy approximaions. Those acually very accurae approximaions only depend on he pah losses of he users uplink signals. Moreover, we provided an efficien algorihm for he clusers o perform local power allocaion in order o maximize he sysem hroughpu a he cos of a single daum exchange in he nework. Simulaion resuls show a perfec fi beween Mone Carlo and heoreical resuls even when he sysem dimension is no very large. REFERENCES [] M. J. M. Peacock, I. B. Collings and M. L. Honig, Eigenvalue disribuions of sums and producs of large random marices via incremenal marix expansions, IEEE Trans. on Informaion Theory, vol. 54, no. 5, pp. 223, 28. [2] J. Miola III and G. Q. Maguire Jr, Cogniive radio: making sofware radios more personal, IEEE personal communicaions, vol. 6, no. 4, pp. 3-8, 999. [3] M. Debbah, Downlink CDMA: o cell or no o cell, 2h European Signal Processing Conference, Vienna, Ausria, 24 [4] O. Somekh, B.J. Zaidel and S. Shamai, Specral Efficiency of Join Muliple Cell-Sies Processors for Randomly Spread DS-CDMA, ISIT, Chicago, USA, 24. [5] S. Shamai and S. Verdu, The Impac of Frequency-Fla Fading on he Specral Efficiency of CDMA, IEEE Trans. on Informaion Theory, pp , May 2. [6] A. Tulino, S. Verdú, Random Marix Theory and Wireless Communicaions, Now Publishers, vol., issue, 24 [7] D.V. Voiculescu, K.J. Dykema and A. Nica Free random variables, American Mahemaical Sociey, 992. [8] P. Biane, Free probabiliy for probabiliss, High Magneic Fields: Science and Technology, World Scienific, pp. 55, 23. [9] T.M. Cover, and J.A. Thomas, Elemens of informaion heory, Wiley New York, 99. [] J. Dumon, S. Lasaulce, W. Hachem, Ph. Loubaon and J. Najim, On he Capaciy Achieving Covariance Marix for Rician MIMO Channels: An Asympoic Approach, submied o IEEE Trans. on Informaion Theory, Ocober 27, revised in Ocober 28. [] A. Coso, S. Savazzi, U. Spagnolini, and C. Ibars, Virual MIMO channels in cooperaive muli-hop wireless sensor neworks, in Proc. Conf. on Informaion Sciences and Sysems, Mar. 26, pp [2] W. Saad, Z. Han, M. Debbah, and A. Hjørungnes, A disribued merge and spli algorihm for fair cooperaion in wireless neworks, in Proc. In. Conf. on Communicaions, Workshop on Cooperaive Communicaions and Neworking, Beijing, China, 28. [3] N. Jindal, U. Mira, and A. Goldsmih, Capaciy of ad-hoc neworks wih node cooperaion, in Proc. In. Symp. on Informaion Theory, Chicago, IL, USA, Jun. 24, pp. 27. [4] J. Cai and U. Pooch, Allocae fair payoff for cooperaion in wireless ad hoc neworks using shapley value, in Proc. In. Parallel and Disribued Processing Symposium, Sana Fe, NM, USA, Apr. 24, pp

Uplink capacity of self-organizing clustered orthogonal CDMA networks in flat fading channels

Uplink capacity of self-organizing clustered orthogonal CDMA networks in flat fading channels Uplink capacity of self-organizing clustered orthogonal CDMA networks in flat fading channels Romain Couillet, Merouane Debbah To cite this version: Romain Couillet, Merouane Debbah. Uplink capacity of

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

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

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

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

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

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

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

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

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

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013 IMPLICI AND INVERSE FUNCION HEOREMS PAUL SCHRIMPF 1 OCOBER 25, 213 UNIVERSIY OF BRIISH COLUMBIA ECONOMICS 526 We have exensively sudied how o solve sysems of linear equaions. We know how o check wheher

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

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

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

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

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

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

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

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

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

Energy-Efficient Link Adaptation on Parallel Channels

Energy-Efficient Link Adaptation on Parallel Channels Energy-Efficien Link Adapaion on Parallel Channels Chrisian Isheden and Gerhard P. Feweis Technische Universiä Dresden Vodafone Chair Mobile Communicaions Sysems D-0069 Dresden Germany Email: (see hp://www.vodafone-chair.com/

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

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

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

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

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

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

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

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

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

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

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

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

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

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

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

. 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

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

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a

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

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

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

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

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

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

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

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

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.

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

arxiv: v1 [cs.dc] 8 Mar 2012

arxiv: v1 [cs.dc] 8 Mar 2012 Consensus on Moving Neighborhood Model of Peerson Graph arxiv:.9v [cs.dc] 8 Mar Hannah Arend and Jorgensen Jos Absrac In his paper, we sudy he consensus problem of muliple agens on a kind of famous graph,

More information

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

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

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

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

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

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

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

Basic Circuit Elements Professor J R Lucas November 2001

Basic Circuit Elements Professor J R Lucas November 2001 Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

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 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

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

Solutions to the Exam Digital Communications I given on the 11th of June = 111 and g 2. c 2

Solutions to the Exam Digital Communications I given on the 11th of June = 111 and g 2. c 2 Soluions o he Exam Digial Communicaions I given on he 11h of June 2007 Quesion 1 (14p) a) (2p) If X and Y are independen Gaussian variables, hen E [ XY ]=0 always. (Answer wih RUE or FALSE) ANSWER: False.

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

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

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

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1. Karl Henrik Johansson Alberto Speranzon,2 Sandro Zampieri

ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1. Karl Henrik Johansson Alberto Speranzon,2 Sandro Zampieri ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1 Karl Henrik Johansson Albero Speranzon, Sandro Zampieri Deparmen of Signals, Sensors and Sysems Royal Insiue of Technology Osquldas

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

ADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91

ADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91 ddiional Problems 9 n inverse relaionship exiss beween he ime-domain and freuency-domain descripions of a signal. Whenever an operaion is performed on he waveform of a signal in he ime domain, a corresponding

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

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

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

Signal and System (Chapter 3. Continuous-Time Systems)

Signal and System (Chapter 3. Continuous-Time Systems) Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

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

Problem Set #1. i z. the complex propagation constant. For the characteristic impedance:

Problem Set #1. i z. the complex propagation constant. For the characteristic impedance: Problem Se # Problem : a) Using phasor noaion, calculae he volage and curren waves on a ransmission line by solving he wave equaion Assume ha R, L,, G are all non-zero and independen of frequency From

More information

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients mahemaics Aricle A Noe on he Equivalence of Fracional Relaxaion Equaions o Differenial Equaions wih Varying Coefficiens Francesco Mainardi Deparmen of Physics and Asronomy, Universiy of Bologna, and he

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

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

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

Intermediate Macro In-Class Problems

Intermediate Macro In-Class Problems Inermediae Macro In-Class Problems Exploring Romer Model June 14, 016 Today we will explore he mechanisms of he simply Romer model by exploring how economies described by his model would reac o exogenous

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 6: Wiener Process

Lecture 6: Wiener Process Lecure 6: Wiener Process Eric Vanden-Eijnden Chapers 6, 7 and 8 offer a (very) brief inroducion o sochasic analysis. These lecures are based in par on a book projec wih Weinan E. A sandard reference for

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

BBP-type formulas, in general bases, for arctangents of real numbers

BBP-type formulas, in general bases, for arctangents of real numbers Noes on Number Theory and Discree Mahemaics Vol. 19, 13, No. 3, 33 54 BBP-ype formulas, in general bases, for arcangens of real numbers Kunle Adegoke 1 and Olawanle Layeni 2 1 Deparmen of Physics, Obafemi

More information

GMM - Generalized Method of Moments

GMM - Generalized Method of Moments GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................

More information

Performance of Stochastically Intermittent Sensors in Detecting a Target Traveling between Two Areas

Performance of Stochastically Intermittent Sensors in Detecting a Target Traveling between Two Areas American Journal of Operaions Research, 06, 6, 99- Published Online March 06 in SciRes. hp://www.scirp.org/journal/ajor hp://dx.doi.org/0.436/ajor.06.60 Performance of Sochasically Inermien Sensors in

More information

CHERNOFF DISTANCE AND AFFINITY FOR TRUNCATED DISTRIBUTIONS *

CHERNOFF DISTANCE AND AFFINITY FOR TRUNCATED DISTRIBUTIONS * haper 5 HERNOFF DISTANE AND AFFINITY FOR TRUNATED DISTRIBUTIONS * 5. Inroducion In he case of disribuions ha saisfy he regulariy condiions, he ramer- Rao inequaliy holds and he maximum likelihood esimaor

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

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

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018 MATH 5720: Gradien Mehods Hung Phan, UMass Lowell Ocober 4, 208 Descen Direcion Mehods Consider he problem min { f(x) x R n}. The general descen direcions mehod is x k+ = x k + k d k where x k is he curren

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

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