One-Sided Adaptation for Infinite-Horizon Linear Quadratic N-person Nonzero-sum Dynamic Games

Size: px
Start display at page:

Download "One-Sided Adaptation for Infinite-Horizon Linear Quadratic N-person Nonzero-sum Dynamic Games"

Transcription

1 Proceedngs of the European Control Conference 7 Kos, Greece, July -5, 7 TuB45 One-Sded Adaptaton for Infnte-Horzon Lnear Quadratc N-person Nonzero-sum Dynamc Games Xaohuan Tan and Jose B Cruz, Jr, Lfe Fellow, IEEE Abstract In ths paper, we consder a class of nfntehorzon dscrete-tme lnear quadratc N-person games, n whch one of the players lac the complete nformaton of the game Wth the assumptons on perfect state nformaton pattern and steady state feedbac strateges, we convert the orgnal game problem nto a multvarable adaptve control problem by tang use of the concept of fcttous play and the scheme of adaptve control For the proposed adjustment procedure, we prove that each element of the estmates converges to ts correspondng true value under the condton of persstent exctaton I INTRODUCTION Most research n game theory has been devoted to the class of games wth complete nformaton, n whch all decson maers have the nowledge of the mathematcal structures as well as the exact values of the parameters of games In ths paper, however, we are more nterested n the class of games wth ncomplete nformaton For ths class of games, Harsany [], [], [3] ntroduced the noton of types and proposed to formulate the problem based on Bayes rules Although mathematcally sold, how to obtan the nformaton of the probablty densty dstrbuton over compettors types remans a bg challenge n applcatons As a wdely used model of learnng, fcttous play [4] s another method to deal wth the games wth uncertantes, n whch decson maers behave as f they thn that they are facng a set of statonary but unnown strateges dstrbuton of compettors In a dynamc case, the decson maers may change ther belefs about the compettors strateges after each round of play, whch leaves whether or not play converges a ey queston If t does not, then t not plausble that players wll mantan that assumpton about compettors strateges Fudenberg and Kreps [5] examned some suffcent condtons under whch fcttous play converges Adaptve control [6], as a mature branch of control theores, provdes decson maers the ablty of self-adjustment or self-modfcaton n accordance wth changng envronment or structure The convergence analyss of adaptve control also gves us lots of nsght on how to deal wth the convergence problem concerned n fcttous play In ths paper, we propose to ncorporate the concept of fcttous play wth the technque of adaptve control to solve a class of game problems wth ncomplete nformaton To ths end, we frst consder a class of nfnte-horzon dscrete-tme lnear quadratc N-person game and complete The authors are wth the Department of Electrcal and Computer Engneerng, The Oho State Unversty, Columbus, OH 43, USA emal: tan4@osuedu; jbcruz@eeeorg) nformaton, gvng out the state feedbac Nash equlbrum solutons under perfect state nformaton pattern n Secton II For such games wth one player who lacs the complete nformaton of the game, we propose an adaptaton scheme for the player n Secton III Persstent exctaton and the analyss of parameter convergence are also gven Smulatons n Secton V numercally verfy our theoretcal results and the conclusons are drawn n the last secton II INFINITE-HORIZON DISCRETE-TIME LQ N -PERSON NONZERO-SUM DYNAMIC GAME Let N {,,,N} be the player s set and K be the maxmum number of stages n a dscrete-tme game such that K {,,,K} Denote the strategy sets by {Γ, N, K} and the decson control) varables by {u U, N, K} wth no chance moves allowed Defne the system dynamcs and the cost functonal)s J x KQ Kx K + x + Ax + NB u ) K x Q x + j Nu j Rj uj ) where x R n, u U, A R n n, B Rm R n m, Q Rn n, R j and R j Rm m >, Q > for all,j N, j, K For any gven control strateges {γ Γ, N, K}, the decsons of the players are completely determned by the relatons u γη ) 3) where η N, K) denotes the nformaton set of P N) [7] Under the perfect state nformaton pattern γ {x }, the feedbac Nash equlbrum can be obtaned by followng the steps of dynamc programmng It s gven n [7] that the game defned by ) and ) admts a unque state feedbac Nash equlbrum strateges γ x ) N, K) u : arg mn J γx ) L x 4) u U where L N, K) are matrces of approprate dmensons, satsfyng the set of lnear matrx equatons L R + B S+B ) B S+ A + B j L j 5) j N,j ISBN:

2 TuB45 and S N, K) are obtaned recursvely from S A + j N +Q + j N B j L j S + A + j N B j L j L j Rj Lj 6) wth SK Q K The game defned above becomes a game wth nfnte horzon when K We drop the tme ndces from the weghtng matrces n ) and let S N) denote the lmt of S K) as K for fxed assumng that ths lmt exsts and s ndependent of Lewse, let L N) be the lmt of L K) as K Then these lmtng matrces necessarly satsfy the followng two algebrac matrx equatons L R + B S B ) B S F 7) S Q + L R L + F + B L ) S F + B L ) 8) where F, Q N) are defned by F A + B j L j 9) j N,j Q Q + j N,j L j R j L j ) The correspondng steady state feedbac strategy for P N) s then u : arg mn J γ x ) L x,,, ) u U Remar So far there s no exstence proof for an nfntehorzon games n general Obtanng the condtons on the parameters of the game that wll guarantee the exstence of a soluton set to 7) and 8) s qute a challenge and beyond the scope of ths research Hereby we wll assume that the soluton set exsts In the nfnte-horzon case, t s natural to assume that the problem s well-defned, n the sense that there exsts at least one set of control sequences that renders t a fnte cost Condtons whch ensure ths are stablzablty of the matrx par F,B ) N) and detectablty of the par F,C ) N) where C s a matrx such that C C Q The latter condton s also referred to detectablty of the par F, Q ) N) When players strateges are fxed at γ x ) L x and requrng that the resultng nfntehorzon game problem s well-defned, S s a soluton to a dscrete-tme algebrac Rccat equaton, as concluded n the followng lemma Lemma For N, f the par A, B ) s stablzable, A, Q ) s detectable and B B s nvertble, then S s the unque postve sem-defnte soluton to the dscrete-tme algebrac Rccat equaton DARE) S Q + F S F F S B R +B S B ) B S F Proof: Tang expanson of 8) and pluggng 7) nto t, we obtan S Q + L R + B S B )L + F S F + F S B L + L B S F Q + F S F + F S B R + B S B ) R + B S B )R + B S B ) B S F F S B R + B S B ) B S F F S B R + B S B ) B S F Q + F S F F S B R + B S B ) B S F ) It can be found n any textboo on lnear control systems that the unque nonnegatve-defnte soluton to ) exsts when the par F,B ) s stablzable and the par F, Q ) s detectable We wll show next that when B B s nvertble, the stablzablty on A,B ) mples the stablzablty on F,B ), and the detectablty on A, Q ) mples the detectablty on F, Q ) When the par A,B ) s stablzable, there exsts some L such that A + B L s Hurwtz When B B s nvertble, we may construct L L B B ) B B L B j L j 3) j N,j such that F + B L s also Hurwtz, snce F + B L A + B L 4) e, F,B ) s also stablzable Smlarly, the detectablty of A, Q ) means that there exsts some M such that A + C M s Hurwtz, where C C Q Snce Q >, C >, such that there also exsts M M C ) C M B j L j ) 5) satsfyng j N,j F + C M A + C M 6) whch mples that F, Q ) s detectable Ths completes the proof of the lemma III ONE-SIDED ADAPTATION FOR N -PERSON GAMES In many practcal cases such as strategc bddng problems n the deregulated electrcty maret, although all players commt to Nash equlbrum strateges [8], not all players are able to access the complete nformaton of the game, especally the weghtng parameters of compettors cost functons as these are prvate and essental nformaton The lac of complete nformaton of a game maes t mpossble for the decson maers to determne ther strateges In ths case, one may tae use of the concept of fcttous play that each player starts wth belevng that the compettors are playng stable strateges and then each player chooses the best response to the belefs The process s repeated 734

3 TuB45 wth updated observatons In the sequel, we adopt the dea and propose to tae use of avalable adaptve control technque to estmate other decson maers strateges before mang each player s own decson We wll frst propose an adjustment procedure n Secton III-A, then dscuss the parameter convergence problem n Secton III-B A Adaptaton Mechansm Desgn Consder an nfnte-horzon dscrete-tme LQ N-person nonzero-sum game defned n Secton II Suppose all players commt to play the steady state feedbac Nash equlbrum strateges u γ x ) L x N,,,), such that the system dynamcs s x + Ax + B u + B j L j x 7) j N,j where for N and,, we assume that x R n, A R n n, B R n m, u U Rm, and L R m n Also suppose P N) s the only decson maer who s unable to access the complete nformaton of the game, and thus has to estmate all other players control strateges before determne hs/her own controls Defne a subset of players by N {,,,+,,N} Let ˆL j j N,,,) be the estmate of L j at the th tme step, and let ˆx be the estmate of x by P Wthout loss of generosty, we let The state estmate calculated by P s then ˆx + Ax + B u + B j ˆLj x 8) j N Denote L j ˆL j Lj, the dfference between the system dynamcs 7) and the state estmate 8) s ˆx + x + [ L B B N] x 9) Let m θ N m and suppose m θ to be such that m θ n to ensure that the column vectors of the n m θ matrx [B B N ] are lnear ndependent, e, there s no redundant control varables If the m θ m θ matrx [B B N ] [B B N ] s nvertble, then we defne the estmaton error vector e + R m θ as [B e + B N] [ B B N]) [ B B N] ˆx+ x + ) ) Also defne the unnown parameter matrx Θ + R n m θ as Θ + [ˆL ] ˆLN ) and ntroduce the regressor vector φ + z [x + ] x R n, where z s a delay operator, and the nomnal L N vector Θ + [L L N ] R n m θ, such that 9) s reduced to e Θ φ, Θ Θ Θ ) whch s ndeed a multvarable lnear parametrc model Let the ntal guess to be Θ, and choose the normalzed gradent adaptve algorthm for updatng Θ to be Θ + Θ Γφ e m ) 3) where < Γ Γ < I n s a gan matrx, and m κ + φ φ wth κ > beng a desgn parameter Lemma The adaptve law 3) guarantees that lm tr[ Θ Θ ] exsts and s fnte Proof: Introduce the postve defnte functon V Θ ) tr[ Θ Γ Θ ], where tr[ ] s the trace operator for a matrx The tme ncrement of V along 3) s V Θ + ) V Θ ) [ tr m ) [ tr m ) e Θ φ e + e φ Θ e φ Γφ e m ) φ Γφ ) ] m ) e Let λ max [Γ], ) be the maxmum egenvalue of Γ whch satsfes the condton < Γ Γ < I n By lnear algebra theory, we have the nequalty < φ Γφ λ max [Γ]φ φ Wth the desgn parameter κ >, we have < φ Γφ m ) φ Γφ κ + φ φ Denote α λ max [Γ] >, then λ max[γ]φ φ φ φ e e V Θ + ) V Θ ) αtr[ ] m ) λ max [Γ] )] 4) whch mples that lm V Θ ) lm tr[ Θ ΓΘ ] exsts and s fnte, snce we have shown that V Θ ) s nonncreasng and bounded from below V Θ ) > ) Furthermore, matrx Γ > can be decomposed nto Γ DΛD, where Λ s the dagonal matrx consstng of the egenvalue of Γ, and D s the orthogonal matrx consstng of correspondng egenvectors, such that V Θ ) tr[ Θ DΛD Θ ] tr[λ Θ Θ ] The exstence of lm V Θ ) thus mples the exstence of lm tr[ Θ Θ ] snce < λ[γ ] < Ths completes the proof of Lemma 735

4 TuB45 B Persstent Exctaton and Parameter Convergence We observe from the estmaton error equaton ) that when e, t s suffcent but not necessary that Θ To mae Θ Θ Θ small when e s small, the regressor φ should be suffcently rch, whch s crucal for ensurng parameter convergence wth the gradent algorthm We thus requre that condton of persstent exctaton [9] s satsfed for the regressor φ, upon whch we prove next n Theorem that each element n the matrx ˆL j j N,,, ) wll converge to ts correspondng true value exponentally Note that we let n the prevous secton for clarfcaton purpose Here we wll leave t as a varable snce any player who s dong estmaton may adopt the same adjustment mechansm as 3) wth approprately defned e and Θ Theorem If the sgnal ψ φ m s persstently exctng over the tme sequence set {,,}, the adaptve law 3) ensures that lm ˆL j ) ql L j ) ql exponentally, where ˆL j ) ql and L j ) ql are the q,l) th element of the matrces ˆL j ) ql and L j ) ql, respectvely j N, q, l n) Proof: Consder the adaptve law 3) wth Γ I, whch can be rewrtten as Θ + Θ φ φ m ) Θ I ψ ψ ) Θ 5) where ψ φ m can be vewed as the observaton vector For a postve defnte functon V Θ ) tr[ Θ Θ ], V Θ + ) V Θ ) tr[ Θ I ψ ψ ) I ψ ψ ) Θ ] tr[ Θ Θ ] tr[ Θ ψ ψ + ψ ψ ) ψ ψ ) Θ ] tr[ Θ ψ ψ Θ ] 6) Denote the state transton matrx as Φσ+δ,σ) σ+δ σ I ψ ψ ) such that Θ σ+δ Φσ +δ,σ) Θ σ Then for arbtrary δ >, σ > t we have [ σ+δ ] V Θ σ+δ ) V Θ σ ) tr Θ ψ ψ Θ where σ V Θ σ ) tr[ Θ σ S Θ σ ] 7) σ+δ S Φ,σ)) ψ ψ Φ, σ) 8) σ If ψ s persstently excted over the nterval σ,σ + δ), the matrx S s stll postve defnte because the observablty property of the system s the same Decompose S > nto Γ DΛ S D, such that tr[ Θ σ S Θ σ ] tr[ Θ σ DΛ S D Θσ ] tr[λ S Θ σ Θσ ] Denote the mnmum egenvalue of S as λ mn [S] >, then tr[ Θ σ S Θ σ ] tr[λ S Θ σ Θσ ] λ mn [S]tr[ Θ σ Θ σ ] whch results n the nequalty 7) to be reduced to V Θ σ+δ ) λ mn [S])V Θ σ ) βv Θ σ ) 9) where β,) as V Θ σ ) > for Θ σ Furthermore, note that the d th d N ) dagonal bloc of the m θ m θ matrx Θ Θ, L j j L j N ), s n L j ) l n L j j L L j ) l n L j ) l where L j ) ql s the q,l) th element of the n matrx L j j N, q, l n,,,) Therefore tr[ Θ Θ ] and 9) s equvalent to j N q l j N q l n L j σ+δ ) ql β It follows that m j n L j t +Tδ ) ql β T j N q l whch mples that lm T + j N q l n L j ) ql 3) j N q l j N q l exponentally Ths s equvalent to say that lm + j N q l n L j σ) ql 3) n L j t ) ql 3) n L j t +Tδ ) ql 33) n L j ) ql 34) snce we have nown from Lemma that the lmt exsts and s fnte By the defnton of lmt, 34) means that ε >, K ε such that > K ε, we have j N q l n L j ) ql < ε 35) Snce L j ) ql for all j N, q, l n,,,, the worst case s when all but one element, say L j ) ql, are zeros In ths case, 35) turns out to be L j ) ql < ε 36) whch agan by defnton mples that lm + L j ) ql, e, lm + ˆL j ) ql L j ) ql Ths completes the proof of Theorem To satsfy the condton on PE, we may ntroduce reference sgnals r and u r N,,,) to be traced by the 736

5 TuB45 state and the controls respectvely, such that the cost functon J N) s revsed to J + [x r ) Q x r ) + u j ujr ) R j uj ujr )] 37) j N where r s a nown trajectory by all players Snce the problem s well-defned, we need to restrct the control references u r N) to be B u r r + Ar 38) N such that n steady state ) when x r and u u r, the system can stll be controlled at a fnte cost The Nash equlbrum soluton to the game wth cost functon 37) can stll be obtaned after tang a coordnate transformaton, e, denote y x r and v u ur The state feedbac Nash soluton for N n the new coordnates s v : arg mn J γ y ) L y 39) v U whch s equvalent to u : arg mn J γ x ) L x + u r u L r ) 4) U n the orgnal coordnates, where L s defned n 7) Gven the form of control n 4) and the restrcton 38), the state equaton ) can be expressed as x + A + N L ) x + r + A + N L ) r 4) whch shows that PE condton can be satsfed by choosng r properly At the every step of the game, wth all estmates ˆL j j N ), P N) can calculate hs/her best response strategy ˆL from 7) and 8), whch n ths adaptaton s ˆL R + B ŜB ) B Ŝ ˆ F 4) Ŝ ˆ Q + ˆ F Ŝ ˆ F ˆ F ŜB R + B ŜB ) B Ŝ ˆ F 43) where ˆ F, ˆ Q N,,,) are defned by ˆ F A + B j ˆLj 44) j N,j ˆ Q Q + j N,j IV NUMERICAL SIMULATIONS ˆL j Rj ˆLj 45) Now we consder a three-person game wth vector control, n whch P lacs the complete nformaton of the costs of other players and thus employs the adaptaton mechansm proposed to estmate other players strateges before mang hs/her own decsons The parameters of the state equaton are supposed to be A B B B The weghtng matrces n the cost functonals are 5 6 Q 5 Q Q R 8 R 5 R 3 5 R 5 5 R 5 5 R 3 4 R R 3 3 R By choosng the desgn parameters κ 4, Γ I 4 and the tracng sgnal r [7sn5), 5cos9), 5cos9 ), 7sn35)], we see from Fg and Fg 3 that each element n ˆL and ˆL 3 matrces converges well to ther correspondng true values n L and L 3, whch maes the elements n ˆL also converge to ther true values n L as shown n Fg V CONCLUSIONS In ths paper, we proposed to solve a class of nfntehorzon LQ Nash games wth ncomplete nformaton by tang use of the concept of fcttous play and the scheme of adaptve control We frst consdered a class of dscrete-tme lnear quadratc N-person game wth complete nformaton, gvng out the state feedbac Nash equlbrum solutons under perfect state nformaton pattern We proved that under condtons of stablzablty and detectablty, the matrx S s actually the unque postve sem-defnte soluton to the assocated coupled dscrete-tme algebrac equaton, whch lad out the base on whch a decson maer can determne hs/her best response strateges aganst the estmates to the compettors strateges For such a class of N-person games wth vector-valued decson varables, we assumed that all players commt to play the steady state feedbac equlbrum strateges, n whch an arbtrary player P lacs the complete nformaton of the game Under some condtons, we converted the game 737

6 TuB45 3 L ) and ˆL ) L ) 3 and ˆL ) L ) and ˆL ) 4 L ) 4 and ˆL ) L 3 ) 3 and ˆL 3 ) 3 5 L 3 ) and ˆL 3 ) L 3 ) 4 and ˆL 3 ) L 3 ) and ˆL 3 ) 8 L ) and ˆL ) 4 L ) and ˆL ) 4 5 L 3 ) and ˆL 3 ) 5 L 3 ) and ˆL 3 ) L ) 3 and ˆL ) 3 6 tme step L ) 4 and ˆL ) 4 tme step 5 L 3 ) 3 and ˆL 3 ) tme step 4 L 3 ) 4 and ˆL 3 ) tme step Fg Control gan L dashed) and ts estmates ˆL L ) ql or ˆL ) ql) s the q, l) th element of L or ˆL ) sold), where Fg 3 Control gan L 3 dashed) and ts estmates ˆL 3 L 3 ) ql or ˆL 3 ) ql) s the q, l) th element of L 3 or ˆL 3 ) sold), where 5 5 L ) and ˆL ) L ) 3 and ˆL ) L ) and ˆL ) 8 4 L ) 3 and ˆL ) tme step 5 L ) and ˆL ) 5 L ) 4 and ˆL ) L ) and ˆL ) 5 5 L ) 4 and ˆL ) tme step Fg Control gan L dashed) and ts estmates ˆL sold), where L ) ql or ˆL ) ql) s the q, l) th element of L or ˆL ) problem nto a mult-nput mult-output adaptve control problem, for whch P employs an adjustment procedure usng normalzed gradent algorthm to estmate the compettors strateges at every tme step before mang hs/her own decsons We proved n theory that each element of P s estmates converges to ts correspondng equlbrum value under the condton of persstent exctaton Numercal smulatons of a three-person game verfed our clams REFERENCES [] J C Harsany, Games wth ncomplete nformaton played by Bayesan player, Part I: The basc model, Management Scence, vol 4, no 3, pp 59 8, Nov 967 [], Games wth ncomplete nformaton played by Bayesan player, Part II: Bayesan equlbrum ponts, Management Scence, vol 4, no 5, pp 3 334, Jan 968 [3], Games wth ncomplete nformaton played by Bayesan player, Part III: The basc probablty dstrbuton of the game, Management Scence, vol 4, no 7, pp 486 5, Mar 968 [4] G W Brown, Iteratve soluton of games by fcttous play, n Actvty Analyss of Producton and Allocaton, T C Koopmans, Ed New Yor: Wley, 95, pp [5] D Fudenberg and D Kreps, Lectures on learnng and equlbrum n stratgc-from games, n CORE Lecture Seres Mmeo, 99 [6] K J Åström, Theory and applcaton of adaptve control - a survey, Automatca, vol 9, pp , 983 [7] T Başar and G J Olsder, Dynamc Noncoopertve Game Theory, revsed ed New Yor, NY: Academc Press, 998 [8] R W Ferrero, V C Ramesh, and S M Shahdehpour, Applcaton games wth ncomplete nformaton for prcng electrcty n deregulated power pools, IEEE Trans on Power Systems, vol 3, no, pp 84 89, Feb 998 [9] G Tao, Adaptve Control Desgn and Analyss Hoboen, NY: Wley- Interscence, 3 738

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

The Second Anti-Mathima on Game Theory

The Second Anti-Mathima on Game Theory The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player

More information

Google PageRank with Stochastic Matrix

Google PageRank with Stochastic Matrix Google PageRank wth Stochastc Matrx Md. Sharq, Puranjt Sanyal, Samk Mtra (M.Sc. Applcatons of Mathematcs) Dscrete Tme Markov Chan Let S be a countable set (usually S s a subset of Z or Z d or R or R d

More information

6. Stochastic processes (2)

6. Stochastic processes (2) Contents Markov processes Brth-death processes Lect6.ppt S-38.45 - Introducton to Teletraffc Theory Sprng 5 Markov process Consder a contnuous-tme and dscrete-state stochastc process X(t) wth state space

More information

6. Stochastic processes (2)

6. Stochastic processes (2) 6. Stochastc processes () Lect6.ppt S-38.45 - Introducton to Teletraffc Theory Sprng 5 6. Stochastc processes () Contents Markov processes Brth-death processes 6. Stochastc processes () Markov process

More information

Singular Value Decomposition: Theory and Applications

Singular Value Decomposition: Theory and Applications Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

First day August 1, Problems and Solutions

First day August 1, Problems and Solutions FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

1 GSW Iterative Techniques for y = Ax

1 GSW Iterative Techniques for y = Ax 1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples

More information

Dynamic Systems on Graphs

Dynamic Systems on Graphs Prepared by F.L. Lews Updated: Saturday, February 06, 200 Dynamc Systems on Graphs Control Graphs and Consensus A network s a set of nodes that collaborates to acheve what each cannot acheve alone. A network,

More information

The lower and upper bounds on Perron root of nonnegative irreducible matrices

The lower and upper bounds on Perron root of nonnegative irreducible matrices Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College

More information

ECE559VV Project Report

ECE559VV Project Report ECE559VV Project Report (Supplementary Notes Loc Xuan Bu I. MAX SUM-RATE SCHEDULING: THE UPLINK CASE We have seen (n the presentaton that, for downlnk (broadcast channels, the strategy maxmzng the sum-rate

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE III LECTURE - 2 EXPERIMENTAL DESIGN MODELS Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 2 We consder the models

More information

Economics 101. Lecture 4 - Equilibrium and Efficiency

Economics 101. Lecture 4 - Equilibrium and Efficiency Economcs 0 Lecture 4 - Equlbrum and Effcency Intro As dscussed n the prevous lecture, we wll now move from an envronment where we looed at consumers mang decsons n solaton to analyzng economes full of

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton

More information

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010 Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

CS286r Assign One. Answer Key

CS286r Assign One. Answer Key CS286r Assgn One Answer Key 1 Game theory 1.1 1.1.1 Let off-equlbrum strateges also be that people contnue to play n Nash equlbrum. Devatng from any Nash equlbrum s a weakly domnated strategy. That s,

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD Matrx Approxmaton va Samplng, Subspace Embeddng Lecturer: Anup Rao Scrbe: Rashth Sharma, Peng Zhang 0/01/016 1 Solvng Lnear Systems Usng SVD Two applcatons of SVD have been covered so far. Today we loo

More information

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

More information

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Hidden Markov Models & The Multivariate Gaussian (10/26/04) CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models

More information

Lecture 10: May 6, 2013

Lecture 10: May 6, 2013 TTIC/CMSC 31150 Mathematcal Toolkt Sprng 013 Madhur Tulsan Lecture 10: May 6, 013 Scrbe: Wenje Luo In today s lecture, we manly talked about random walk on graphs and ntroduce the concept of graph expander,

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Feb 14: Spatial analysis of data fields

Feb 14: Spatial analysis of data fields Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

COS 521: Advanced Algorithms Game Theory and Linear Programming

COS 521: Advanced Algorithms Game Theory and Linear Programming COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Neuro-Adaptive Design - I:

Neuro-Adaptive Design - I: Lecture 36 Neuro-Adaptve Desgn - I: A Robustfyng ool for Dynamc Inverson Desgn Dr. Radhakant Padh Asst. Professor Dept. of Aerospace Engneerng Indan Insttute of Scence - Bangalore Motvaton Perfect system

More information

Deriving the X-Z Identity from Auxiliary Space Method

Deriving the X-Z Identity from Auxiliary Space Method Dervng the X-Z Identty from Auxlary Space Method Long Chen Department of Mathematcs, Unversty of Calforna at Irvne, Irvne, CA 92697 chenlong@math.uc.edu 1 Iteratve Methods In ths paper we dscuss teratve

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Appendix B. Criterion of Riemann-Stieltjes Integrability

Appendix B. Criterion of Riemann-Stieltjes Integrability Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

Lecture 4: Constant Time SVD Approximation

Lecture 4: Constant Time SVD Approximation Spectral Algorthms and Representatons eb. 17, Mar. 3 and 8, 005 Lecture 4: Constant Tme SVD Approxmaton Lecturer: Santosh Vempala Scrbe: Jangzhuo Chen Ths topc conssts of three lectures 0/17, 03/03, 03/08),

More information

Infinitely Split Nash Equilibrium Problems in Repeated Games

Infinitely Split Nash Equilibrium Problems in Repeated Games Infntely Splt ash Equlbrum Problems n Repeated Games Jnlu L Department of Mathematcs Shawnee State Unversty Portsmouth, Oho 4566 USA Abstract In ths paper, we ntroduce the concept of nfntely splt ash equlbrum

More information

for Linear Systems With Strictly Diagonally Dominant Matrix

for Linear Systems With Strictly Diagonally Dominant Matrix MATHEMATICS OF COMPUTATION, VOLUME 35, NUMBER 152 OCTOBER 1980, PAGES 1269-1273 On an Accelerated Overrelaxaton Iteratve Method for Lnear Systems Wth Strctly Dagonally Domnant Matrx By M. Madalena Martns*

More information

Lecture 17 : Stochastic Processes II

Lecture 17 : Stochastic Processes II : Stochastc Processes II 1 Contnuous-tme stochastc process So far we have studed dscrete-tme stochastc processes. We studed the concept of Makov chans and martngales, tme seres analyss, and regresson analyss

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

Random Walks on Digraphs

Random Walks on Digraphs Random Walks on Dgraphs J. J. P. Veerman October 23, 27 Introducton Let V = {, n} be a vertex set and S a non-negatve row-stochastc matrx (.e. rows sum to ). V and S defne a dgraph G = G(V, S) and a drected

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013 ISSN: 2277-375 Constructon of Trend Free Run Orders for Orthogonal rrays Usng Codes bstract: Sometmes when the expermental runs are carred out n a tme order sequence, the response can depend on the run

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

Feature Selection: Part 1

Feature Selection: Part 1 CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?

More information

Entropy of Markov Information Sources and Capacity of Discrete Input Constrained Channels (from Immink, Coding Techniques for Digital Recorders)

Entropy of Markov Information Sources and Capacity of Discrete Input Constrained Channels (from Immink, Coding Techniques for Digital Recorders) Entropy of Marov Informaton Sources and Capacty of Dscrete Input Constraned Channels (from Immn, Codng Technques for Dgtal Recorders). Entropy of Marov Chans We have already ntroduced the noton of entropy

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

Stability and Stabilization for Discrete Systems with Time-varying Delays Based on the Average Dwell-time Method

Stability and Stabilization for Discrete Systems with Time-varying Delays Based on the Average Dwell-time Method Proceedngs of the 29 IEEE Internatonal Conference on Systems, an, and Cybernetcs San Antono, TX, USA - October 29 Stablty and Stablzaton for Dscrete Systems wth Tme-varyng Delays Based on the Average Dwell-tme

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

THE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION

THE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION THE WEIGHTED WEAK TYPE INEQUALITY FO THE STONG MAXIMAL FUNCTION THEMIS MITSIS Abstract. We prove the natural Fefferman-Sten weak type nequalty for the strong maxmal functon n the plane, under the assumpton

More information

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M CIS56: achne Learnng Lecture 3 (Sept 6, 003) Preparaton help: Xaoyng Huang Lnear Regresson Lnear regresson can be represented by a functonal form: f(; θ) = θ 0 0 +θ + + θ = θ = 0 ote: 0 s a dummy attrbute

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Digital Signal Processing

Digital Signal Processing Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

The Prncpal Component Transform The Prncpal Component Transform s also called Karhunen-Loeve Transform (KLT, Hotellng Transform, oregenvector Transfor

The Prncpal Component Transform The Prncpal Component Transform s also called Karhunen-Loeve Transform (KLT, Hotellng Transform, oregenvector Transfor Prncpal Component Transform Multvarate Random Sgnals A real tme sgnal x(t can be consdered as a random process and ts samples x m (m =0; ;N, 1 a random vector: The mean vector of X s X =[x0; ;x N,1] T

More information

Continuous Time Markov Chain

Continuous Time Markov Chain Contnuous Tme Markov Chan Hu Jn Department of Electroncs and Communcaton Engneerng Hanyang Unversty ERICA Campus Contents Contnuous tme Markov Chan (CTMC) Propertes of sojourn tme Relatons Transton probablty

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Communication Complexity 16:198: February Lecture 4. x ij y ij

Communication Complexity 16:198: February Lecture 4. x ij y ij Communcaton Complexty 16:198:671 09 February 2010 Lecture 4 Lecturer: Troy Lee Scrbe: Rajat Mttal 1 Homework problem : Trbes We wll solve the thrd queston n the homework. The goal s to show that the nondetermnstc

More information

Tracking with Kalman Filter

Tracking with Kalman Filter Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Math 217 Fall 2013 Homework 2 Solutions

Math 217 Fall 2013 Homework 2 Solutions Math 17 Fall 013 Homework Solutons Due Thursday Sept. 6, 013 5pm Ths homework conssts of 6 problems of 5 ponts each. The total s 30. You need to fully justfy your answer prove that your functon ndeed has

More information

Convergence of random processes

Convergence of random processes DS-GA 12 Lecture notes 6 Fall 216 Convergence of random processes 1 Introducton In these notes we study convergence of dscrete random processes. Ths allows to characterze phenomena such as the law of large

More information

k t+1 + c t A t k t, t=0

k t+1 + c t A t k t, t=0 Macro II (UC3M, MA/PhD Econ) Professor: Matthas Kredler Fnal Exam 6 May 208 You have 50 mnutes to complete the exam There are 80 ponts n total The exam has 4 pages If somethng n the queston s unclear,

More information

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper Games of Threats Elon Kohlberg Abraham Neyman Workng Paper 18-023 Games of Threats Elon Kohlberg Harvard Busness School Abraham Neyman The Hebrew Unversty of Jerusalem Workng Paper 18-023 Copyrght 2017

More information

THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY. William A. Pearlman. References: S. Arimoto - IEEE Trans. Inform. Thy., Jan.

THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY. William A. Pearlman. References: S. Arimoto - IEEE Trans. Inform. Thy., Jan. THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY Wllam A. Pearlman 2002 References: S. Armoto - IEEE Trans. Inform. Thy., Jan. 1972 R. Blahut - IEEE Trans. Inform. Thy., July 1972 Recall

More information

Randomness and Computation

Randomness and Computation Randomness and Computaton or, Randomzed Algorthms Mary Cryan School of Informatcs Unversty of Ednburgh RC 208/9) Lecture 0 slde Balls n Bns m balls, n bns, and balls thrown unformly at random nto bns usually

More information

ECEN 667 Power System Stability Lecture 21: Modal Analysis

ECEN 667 Power System Stability Lecture 21: Modal Analysis ECEN 667 Power System Stablty Lecture 21: Modal Analyss Prof. Tom Overbye Dept. of Electrcal and Computer Engneerng Texas A&M Unversty, overbye@tamu.edu 1 Announcements Read Chapter 8 Homework 7 s posted;

More information