Time-Inconsistent Mean-Field Stochastic Linear-Quadratic Optimal Control

Size: px
Start display at page:

Download "Time-Inconsistent Mean-Field Stochastic Linear-Quadratic Optimal Control"

Transcription

1 Time-Inconsistent Mean-Fied Stochastic Linear-Quadratic Optima Contro I Yuan-Hua 1 1. Department of Mathematics Schoo of Sciences Tianjin Poytechnic University Tianjin P. R. China E-mai: yhni@amss.ac.cn Abstract: This paper is concerned with a ind of time-consistent soution of the time-inconsistent mean-fied stochastic inearquadratic optima contro. Different from standard stochastic inear-quadratic probems both the system matrices and the weighting matrices are dependent on initia time and the conditiona expectations of the contro and state enter quadraticay into the cost functiona. Such features wi ruin the Beman optimaity principe and resut in the time-inconsistency of the optima contro. Due to the dynamica nature of contro probems a ind of open-oop time-consistent equiibrium contro is thoroughy investigated in this paper. It is shown that the existence of an open-oop time-consistent equiibrium contro for a fixed initia pair is equivaent to the sovabiity of a set of forward-bacward stochastic difference equations with stationary conditions and convexity conditions. By decouping the forward-bacward stochastic difference equations the existence of the open-oop e- quiibrium contro for a the initia pairs is shown to be equivaent to the sovabiity of a set of couped constraint generaized difference Riccati equations and a set of couped constraint inear difference equations. Key Words: Time-inconsistency mean-fied contro stochastic inear-quadratic optima contro 1 Introduction The Beman optimaity principe is an essentia property in optima contro theory which provides the theoretica foundation of the dynamic programming approach. From the Beman optimaity principe an optima contro for an initia pair is aso optima aong the optima trajectory. Such a phenomenon is referred as the time-consistency of the optima contro which ensures that one needs ony to sove an optima contro probem for a given initia pair and the obtained optima contro is aso optima aong the optima trajectory. However in reaity the time-consistency fais quite often. For instance when the initia time or initia state enters into the system dynamics or cost functiona expicity or even more the conditiona expectations of the state or contro enters nonineary into the cost functiona the corresponding probems are time-inconsistent. See exampes in 13] and 5] about the hyperboic discounting and quasi-geometric discounting. The probem with noninear term of the conditiona expectation in cost functiona is caed as mean-fied stochastic optima contro. In this case the smoothing property of the conditiona expectation wi not be in effect to ensure the time-consistency of the optima contro. A particuar exampe of this case is the nown mean-variance utiity 3] and 5]. Due to the dynamica nature of the contro probems it is reasonabe to hande the time-inconsistency in a dynamic manner. Instead of seeing an optima contro some ind of equiibrium contros are concerned with. This is mainy motivated by practica appications such as in mathematica finance and economics and has recenty attracted considerabe interest and efforts. The mathematica formuation of the time-inconsistency was first reported by 22] and a quaitative anaysis was given in 21]. Foowing 22] the wors 10] 13] 14] and 19] are for discrete dynamic system- Yuan-Hua i is supported by the ationa atura Science Foundation of China the China Postdoctora Science Foundation 2014M T s or simpe ordinary differentia equations ODEs. Subsequenty 6] and 7] studied the non-exponentia discounting probems both for simpe ODEs and stochastic differentia equations and introduced the notion of time-consistent contro. 5] discussed the probems of genera Marovian time-inconsistent stochastic optima contro. 24] and 25] addressed the deterministic continuous-time inear-quadratic LQ optima contro by an essentiay cooperative game approach. 27] considered the stochastic LQ probem of meanfied case which is caed the cosed-oop formuation there. Different from 24] and 25] 12] studied another ind of time-consistent equiibrium contro which is an infinitesimay open-oop optima contro. Roughy speaing the controer given by 12] can commit to the equiibrium contro in an infinitesima manner. 27] thoroughy investigated both the open-oop and the cosed-oop time-consistent soutions for genera mean-fied stochastic LQ probems. It is shown 27] that the existence of the open-oop equiibrium contro and the cosed-oop equiibrium strategy is ensured via the sovabiity of certain systems of Riccati-type differentia equations. However a these existing resuts focus on the definite LQ probems. In this paper we sha investigate a time-inconsistent mean-fied stochastic LQ optima contro probem whose system dynamics and cost functiona are aso dependent on the initia time. It is worth noting that the definiteness constraints are not posed on the state and contro weighting matrices. Here we adopt a dynamic manner to attac the time-inconsistency and intend seeing a ind of open-oop equiibrium contro of Probem LQ. After giving the definition of open-oop equiibrium pair its existence for a fixed initia pair is shown to be equivaent to the sovabiity of a set of forward-bacward stochastic difference equations FBS Es with stationary conditions and convexity conditions. If for a fixed initia pair Probem LQ admits an openoop equiibrium pair then a set of constrained inear difference equations LDEs is sovabe and the open-oop equiibrium contro admits a cosed-oop representation. Here the cosed-oop representation is a inear feedbac of the cur-

2 rent vaue of equiibrium state whose gains are computed via the soutions of the constraint LDEs 11 and a set of generaized difference Riccati equations GDREs 14. For any initia pair Probem LQ admitting an open-oop equiibrium pair is shown to be equivaent to that the constraint LDEs 15 and a set of constraint GDREs 16 are sovabe. Different from is of constrained GDREs. Interestingy if sovabe the set of GDREs 16 does not have symmetric structure i.e. its soution is not symmetric. In 16] a simpified version of Probem LQ is considered where the system dynamics and the cost functiona do not contain mean-fied terms. Hence this paper is a continuation of 16]. If the system dynamics and cost functiona are independent of the initia time the corresponding LQ probem is a dynamic version of that considered in 17] where a static version is studied with the conditiona expectation operator being repaced by the expectation operator. For detais on mean-fied stochastic optima contro and reated meanfied games we refer to for exampe 8] 17] 18] 26] 4] 11] 15] and the references therein. The rest of the paper is organized as foows. Section 2 presents the resuts mentioned above and Section 3 gives some concuding remars. 2 Open-oop time-consistent equiibrium contro Consider the foowing controed stochastic difference e- quation S E X+1 t A t X t + ĀtE t X t + B t u + B t E t u + C t X t + C t E t X t 1 + D t u + D t E t u w Xt t x T t t T where T {0... 1} T t {t 1} and A t Āt C t C t R n n B t B t D t D t R n m are deterministic matrices; {X t T t } X t and {u T t } u with T t being {0... } are the state process and contro process respectivey. Here the initia time t is parameterized in the matrices and state to emphasize that the matrices and state may change according to the initia time t. Simiar notations wi be used throughout the paper. In 1 E t is the conditiona mathematica expectation E F t ] with respect to F t {x 0 w 0 1 t 1} and F 1 is understood as { Ω}. The noise {w T} is assumed to be a martingae difference sequence defined on a probabiity space Ω F P with E +1 w +1 2 ] 1 T. 2 The cost functiona associated with system 1 is Jt x; u t E t X t T Q t X t + E t X t T Qt E t X t + u T R t u + E t u T Rt E t u ] + E t X t T G t X t + E t X t T Ḡ t E t X t 3 where Q t Q t R t R t T t G t Ḡt are deterministic symmetric matrices of appropriate dimensions. In 1 x is in L 2 F t; Rn which is a set of random variabes such that any ξ L 2 F t; Rn is F t -measurabe and E ξ 2 <. Let L 2 F T t; H be a set of H-vaued processes such that for any ν {ν T t } L 2 F T t; H ν is F -measurabe and t E ν 2 <. Then we pose the foowing optima contro probem. Probem LQ. Concerned with 13 and the initia pair t x find a u L 2 F T t; R m such that Jt x; u inf Jt x; u. 4 u L 2 F Tt;Rm Instead of soving Probem LQ for the static precommitted optima contro we adopt the concept of dynamic equiibrium contro which is optima in an infinitesima sense and consistent with the dynamica nature of Probem LQ. Definition 2.1 Given t T and x L 2 F t; Rn a statecontro pair X u with u L 2 F T t; R m is caed an open-oop equiibrium pair of Probem LQ for the initia pair t x if X t x and J X ; u T J X ; u u T+1 5 hods for any u L 2 F ; Rm and T t. Here u T and u T+1 are the restrictions of u on T and T +1 respectivey. Furthermore such a u is caed an openoop equiibrium contro for the initia pair t x. For a u L 2 F T t; R m the requirement that u is F 1 - measurabe is parae to the standard statement on the admissibe contros of continuous-time stochastic optima contro; see 9] 23] for detais. In other words u is determined from { x 0 w 0... w 1 } ony irrespective of how the state process deveops. From this and the standard arguments about the open-oop contro 2] u L 2 F T t; R m can be viewed as an open-oop contro. Hence we ca u an open-oop equiibrium contro. Furthermore noting that u T u u T+1 the contro u u T+1 on the right hand of the inequaity of 5 differs from u T ony at stage. Therefore 5 is viewed as a oca optimaity condition. Simiary to 16] we can show that {u t... u } is the ash equiibrium of a muti-person game with hierarchica structure. By its definition an open-oop equiibrium contro u is timeconsistent in the sense that for any T t u T is an open-oop equiibrium contro for the initia pair X. In other words u is time-consistent aong the trajectory of the equiibrium state X. By a formua of the difference of cost functionas we can derive the foowing necessary and sufficient condition to the existence of the open-oop equiibrium pair for a given initia pair whose proof is omitted due to space imitations. Theorem 2.1 Given t T and x L 2 F t; Rn the foowing statements are equivaent. i There exists an open-oop equiibrium pair of Probem LQ for the initia pair t x. ii There exists a u L 2 F T t; R m such that for any T t the foowing FBS E admits a soution

3 X tx Z tx X tx +1 A X tx Z tx X tx Z tx T + B u + C X tx + D u + ĀE X tx + B E u + C E X tx + D E u w A T E Z tx +1 + ĀT E Z tx +1 + C T E Z tx +1 w + C T E Z tx +1 w + Q X tx X G X tx with the stationary condition 0 R + R u + Q E X tx + ḠE X tx + B + B T E Z tx +1 + D + D T E Z tx +1 w 7 and the convexity condition 0 E u T R + R u ] + + E Y ū + E Y ū E Y ū T Q Y ū ] T Q E Y ū T G Y ū 6 ] + E Y ū T Ḡ E Y ū.8 In the above Y ū and X are given respectivey by Y ū +1 A Y ū + ĀE Y ū + C Y ū + C E Y ū w T +1 Y ū +1 B + B ū + D + D ū w Y ū 0 and X ĀX +1 A + + C + C X X t x T t. + B + B u + D + D u ] w In this case X u given in ii is an open-oop equiibrium pair. Though the means of the state and contro appear in both the system dynamics and the cost functiona the stationary condition 7 does not contain the means. This is different from the case of static mean-fied LQ optima contro 17]. Here the static mean-fied LQ optima contro is a variant of Probem LQ where the conditiona expectation operator E t is repaced by the expectation operator E and the matrices in the system and cost functiona are independent of the initia time. Another interesting phenomenon is that the system 9 for the equiibrium pair doesn t contain the mean-fied terms too. To simpify the notations et A A + Ā B B + B C C + C D D + D Q Q + Q R R + R G G + Ḡ T t T. 9 Reca the pseudo-inverse of a matrix. By 20] for a given matrix M R n m there exists a unique matrix in R m n denoted by M such that { MM M M M MM M MM T MM M M T M 10 M. This M is caed the Moore-Penrose inverse of M. The foowing emma is from 1]. Lemma 2.1 Let matrices L M and be given with appropriate size. Then LXM has a soution X if and ony if LL MM. Moreover the soution of LXM can be expressed as X L M + Y L LY MM where Y is a matrix with appropriate size. The foowing theorem is concerned with the necessary condition on the existence of the open-oop equiibrium pair for a fixed initia pair. Theorem 2.2 Let Probem LQ for the initia pair t x admit an open-oop equiibrium pair. Then the foowing set of LDEs P Q + A T P +1A + C T P +1C P Q + A T P +1A + C T P +1C P G P G T R + B T P +1B + D T P +1D 0 T t 11 is sovabe in the sense that the foowing inequaities hod { R + B T P +1B + D T P +1D 0 12 T t. Furthermore we have Z tx P X t T + T X E X t + P E X tx E X + T E X and an open-oop equiibrium contro is given by u W H X T t. 13 Here {T T T t } {T T T t } {W T t } and {H T t } are given respectivey by T A T T +1A + C T T +1C A T P +1B + A T T +1B + C T P +1D + C T T +1D W H T A T T +1A + C T T +1C A T P +1B + A T T +1B + C T P +1D + C T T +1D W H T 0 T 0 T T t 14

4 and W R + B T P+1 + T +1 B + D T P+1 + T +1 D H B T P+1 + T +1 A + D T P+1 + T +1 C. Proof. Denote the righthand side of 8 by Ĵ 0; ū. We then have Ĵ 0; ū E Y ū T Q Y ū + ū T R ū + E Y ū + E Y ū T Ḡ E Y ū E Y ū T Q Y ū + Y ū +1 T P +1 Y ū +1 Y ū + E Y ū +1 T P+1 E Y ū + ū T R ū +1 + E Y ū T G Y ū ] T Q E Y ū ] + E Y ū T Q E Y ū T P Y ū +1 E Y ū E E Y ū T Q + A T P +1 A + C T P +1 C P E Y ū + Y ū E Y ū T Q + A T P +1 A + C T P +1 C P Y ū E Y ū ] T P E Y ū ] + ū T R + B T P +1 B + D T P +1 D ū ū T R + B T P +1 B + D T P +1 D ū. By this and 8 we have the sovabiity of 11. Let u be an open-oop equiibrium contro. Then from Theorem 2.1 for any T t the FBS E 6 admits a soution and 7 hods. As Z tx we have by 7 G X tx + ḠE X tx 0 R u + B T G E X ] + g + DG T E X w. Substituting X and by Lemma 2.1 we have u W H X Ψ X where W R + B T G B + D T G D H B T G A + D T G C. Simiary we have Z 2tx P 2 + T 2 X From 7 we have for 2 + P 2 + T 2 E 2 X. 0 R 2 2 u 2 + B 2 2 T P 2 + T 2 E 2 X ] + D 2 2 T P 2 + T 2 E 2 X w 2. Substituting X 2 and by Lemma 2.1 we have u 2 W 2 H 2X 2 Ψ 2X 2 where W 2 R B T 2 2 P 2 + T 2 B D 2 2 T P 2 + T 2 D 2 2 H 2 B 2 2 T P 2 + T 2 A D T 2 2 P 2 + T 2 C 2 2. Bacwardy repeating above procedure we can get 14 and 15. This competes the proof. Theorem 2.3 For any t T and any x L 2 F t; Rn Probem LQ for the initia pair t x admits an open-oop equiibrium pair if and ony if P Q + A T P +1A + C T P +1C P Q + A T P +1A + C T P +1C P G P G T R + B T P +1B + D T P +1D 0 T and the set of GDREs T A T T +1A + C T T +1C A T P +1B + A T T +1B + C T P +1D + C T T +1D W H T A T T +1A + C T T +1C A T P +1B + A T T +1B + C T P +1D + C T T +1D W H T 0 T 0 T W W H H 0 T are sovabe where W R + B T P+1 + T +1 B + D T P+1 + T +1 D H B T P+1 + T +1 A + D T P+1 + T +1 C

5 In this case an open-oop equiibrium contro for the initia pair t x is given in u W H X T t. 18 Proof. Sufficiency. Introduce a dynamics X +1 A B W H X + A B W H and a contro X t x T t ũ W H X w 19 X T t. 20 Then by reversing the first part of the proof of Theorem 2.2 we can show that for any T t the foowing FBS E admits an adapted soution X tx +1 A Xtx tx X tx T + B ũ + C Xtx + D ũ A T E with property and tx 0 R ũ + C T E + Q Xtx + ĀE Xtx + B E ũ + C E Xtx + D E ũ tx +1 + ĀT E tx X tx P + T X tx tx +1 w +1 w + C T E + Q E Xtx X G Xtx tx +1 w + ḠE Xtx E Xtx + P E Xtx E X + T E X + BE T tx +1 + DT tx E +1 w. Furthermore by 12 and 15 we have 8. From Theorem 2.1 Probem LQ for the initia pair t x admits an openoop equiibrium pair and X ũ is an equiibrium pair. ecessity. By Theorem 2.2 we need ony to prove W W H H 0 T. 21 Consider Probem LQ for the initia pair 1 x with x L 2 F 1; Rn. By the proof of Theorem 2.2 we have 0 W u x + H X x. 22 Let e i be a R n -vaued vector with the i-th entry being 1 and other entries 0. Then we have 0 W u e 1... u en + H e1... e n. oting that e 1... e n is the identity matrix and by Lemma 2.1 we have W W H H 0. Considering Probem LQ for the initia pair 2 x with x L 2 F 2; Rn we can simiary prove W 2 W 2 H 2 H 2 0. Continuing above procedure we then achieve the concusion. ote that P P T T +1 are symmetric. If Q Q R R are seected such that Q Q + Q R R + R 0 T T then 15 is sovabe. Furthermore Θ {P P T T π } are used to express Z tx. {P P } is then caed the symmetric part of Θ and {T T } are viewed as the nonsymmetric part. 3 Concusion In this paper the open-oop time-consistent equiibrium contro is investigated for a ind of mean-fied stochastic LQ probem. ecessary sufficient conditions are presented for the cases with a fixed initia pair and a the initia pairs. Furthermore the GDREs and LDEs are introduced to characterize the open-oop equiibrium contro. For future research the cosed-oop time-consistent soution shoud be studied. References 1] M. Ait Rami X. Chen and X. Y. Zhou. Discrete-time indefinite LQ contro with state and contro dependent noises. Journa of Goba Optimization 2002 vo ] K. J. Astrom. Introduction to Stochastic Contro Theory. A- cademic Press Inc. ew Yor Mathematics in Science and Engineering 1970 vo.70. 3] S. Basa and G. Chabaauri. Dynamic mean-variance asset aocation. Rev. Financia Stud vo ] A. Bensoussan J. Frehse and P. Yam. Mean fied games and mean fied type contro theory. Spriger Briefs in Mathematics ] T. Bjor A. Murgoci. A genera theory of Marovian time inconsisiten stochastic contro probem. Woring paper. 6] I. Eeand A. Lazra. Being serious about non-commitment: subgame perfect equiibrium in continuous time. Preprint U- niversity of University British Coumbia ] I. Eeand T. Privu. Investment and consumption without commitment. Preprint University of University British Coumbia ] R.J. Eiott X. Li and Y.H. i. Discrete time mean-fied s- tochastic inear-quadratic optima contro probems. Automatica 2013 vo.49 no ] W.H. Feming H.M. Soner. Controed Marov Processes and Viscosity Soutions Second Edition. Springer Science+Business Media Inc ] S. M. Godman. Consistent pan. Review of Economic Studies 1980 vo ] M. Huang P.E. Caines and R.P. Mahame. Largepopuation cost-couped LQG probems with nonuniform a- gents: individua-mass behavior and decentraized ε-ash equiibria. IEEE Transactions on Automatic Contro 2007 vo ] Y. Hu H. Jin and X.Y. Zhou. Time-inconsistent stochastic inear-quadratic contro. SIAM Journa on Contro and Optimization 2012 vo ] P. Kruse A. A. Smith. Consumption and savings decisions with quasi-georetric discounting. Econometrica 2003 vo.71 no

6 14] D. Laibson. Goden eggs and hyperboic discounting. The Quartery Journa of Economics 1997 vo ] J.M. Lasry and P.L. Lions. Mean-fied games Japanese Journa of Mathematics 2007 vo ] X. Li Y.H. i and J.F. Zhang. On time-consistent soution to time-inconsistent inear-quadratic optima contro of discretetime stochastic systems. Submitted to IEEE Transactions on Automatic Contro ] Y.H. i J.F. Zhang and X. Li. Indefinite mean-fied stochastic inear-quadratic optima contro. IEEE Transactions on Automatic Contro 2015 vo.60 no ] Y.H. i X. Li and J.F. Zhang. Indefinite mean-fied stochastic inear-quadratic optima contro: from finite horizon to infinite horizon. IEEE Transactions on Automatic Contro DOI: /TAC ] I. Paacios-Huerta. Time-inconsistent preferences in Adam Smith and Davis Hume. History of Poitica Economy 2003 Vo ] R. Penrose. A generaized inverse of matrices. Mathematica Proceedings of the Cambridge Phiosophica Society 1955 vo ] A. Smith. The theory of mora sentiments. First Edition 1759; Reprint Oxford University Press ] R.H. Strotz. Myopia and inconsistency in dynamic utiity maximization. Review of Economic Studies 1955 vo ] J.M. Yong and X.Y. Zhou. Stochastic contros: Hamitonian systems and HJB equations. Springer ew Yor ] J.M. Yong. A deterministic inear quadratic time-inconsitent optima contro probem. Mathematica Contro and Reated Rieds 2011 vo.1 no ] J.M. Yong. Deterministic time-inconsistent optima contro probems an essentiay cooperative approach. Acta App. Math. Sinica 2012 vo ] J.M. Yong. A inear-quadratic optima contro probem for mean-fied stochastic differentia equations. SIAM Journa on Contro Optimization 2013 vo ] J.M. Yong. Linear-quadratic optima contro probems for mean-fied stochastic differentia equations time-consistent soutions. Eectronicay pubished by Transactions of the American Mathematica Society DOI: /tran/

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC (January 8, 2003) A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC DAMIAN CLANCY, University of Liverpoo PHILIP K. POLLETT, University of Queensand Abstract

More information

Asynchronous Control for Coupled Markov Decision Systems

Asynchronous Control for Coupled Markov Decision Systems INFORMATION THEORY WORKSHOP (ITW) 22 Asynchronous Contro for Couped Marov Decision Systems Michae J. Neey University of Southern Caifornia Abstract This paper considers optima contro for a coection of

More information

Theory of Generalized k-difference Operator and Its Application in Number Theory

Theory of Generalized k-difference Operator and Its Application in Number Theory Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication

More information

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

Mat 1501 lecture notes, penultimate installment

Mat 1501 lecture notes, penultimate installment Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the

More information

Disturbance decoupling by measurement feedback

Disturbance decoupling by measurement feedback Preprints of the 19th Word Congress The Internationa Federation of Automatic Contro Disturbance decouping by measurement feedback Arvo Kadmäe, Üe Kotta Institute of Cybernetics at TUT, Akadeemia tee 21,

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

The Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation

The Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation The Symmetric Antipersymmetric Soutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B 2 + + A X B C Its Optima Approximation Ying Zhang Member IAENG Abstract A matrix A (a ij) R n n is said to be symmetric

More information

Reliability: Theory & Applications No.3, September 2006

Reliability: Theory & Applications No.3, September 2006 REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate

More information

Global Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations

Global Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations Goba Optimaity Principes for Poynomia Optimization Probems over Box or Bivaent Constraints by Separabe Poynomia Approximations V. Jeyakumar, G. Li and S. Srisatkunarajah Revised Version II: December 23,

More information

King Fahd University of Petroleum & Minerals

King Fahd University of Petroleum & Minerals King Fahd University of Petroeum & Mineras DEPARTMENT OF MATHEMATICAL SCIENCES Technica Report Series TR 369 December 6 Genera decay of soutions of a viscoeastic equation Saim A. Messaoudi DHAHRAN 3161

More information

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled. imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,

More information

Improving the Reliability of a Series-Parallel System Using Modified Weibull Distribution

Improving the Reliability of a Series-Parallel System Using Modified Weibull Distribution Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution

More information

BP neural network-based sports performance prediction model applied research

BP neural network-based sports performance prediction model applied research Avaiabe onine www.jocpr.com Journa of Chemica and Pharmaceutica Research, 204, 6(7:93-936 Research Artice ISSN : 0975-7384 CODEN(USA : JCPRC5 BP neura networ-based sports performance prediction mode appied

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

Analysis of the problem of intervention control in the economy on the basis of solving the problem of tuning

Analysis of the problem of intervention control in the economy on the basis of solving the problem of tuning Anaysis of the probem of intervention contro in the economy on the basis of soving the probem of tuning arxiv:1811.10993v1 q-fin.gn] 9 Nov 2018 1 st eter Shnurov Nationa Research University Higher Schoo

More information

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation

More information

Intuitionistic Fuzzy Optimization Technique for Nash Equilibrium Solution of Multi-objective Bi-Matrix Games

Intuitionistic Fuzzy Optimization Technique for Nash Equilibrium Solution of Multi-objective Bi-Matrix Games Journa of Uncertain Systems Vo.5, No.4, pp.27-285, 20 Onine at: www.jus.org.u Intuitionistic Fuzzy Optimization Technique for Nash Equiibrium Soution of Muti-objective Bi-Matri Games Prasun Kumar Naya,,

More information

Combining reaction kinetics to the multi-phase Gibbs energy calculation

Combining reaction kinetics to the multi-phase Gibbs energy calculation 7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation

More information

Effect of transport ratio on source term in determination of surface emission coefficient

Effect of transport ratio on source term in determination of surface emission coefficient Internationa Journa of heoretica & Appied Sciences, (): 74-78(9) ISSN : 975-78 Effect of transport ratio on source term in determination of surface emission coefficient Sanjeev Kumar and Apna Mishra epartment

More information

Mode in Output Participation Factors for Linear Systems

Mode in Output Participation Factors for Linear Systems 2010 American ontro onference Marriott Waterfront, Batimore, MD, USA June 30-Juy 02, 2010 WeB05.5 Mode in Output Participation Factors for Linear Systems Li Sheng, yad H. Abed, Munther A. Hassouneh, Huizhong

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t. Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes

More information

EXPONENTIAL DECAY OF SOLUTIONS TO A VISCOELASTIC EQUATION WITH NONLINEAR LOCALIZED DAMPING

EXPONENTIAL DECAY OF SOLUTIONS TO A VISCOELASTIC EQUATION WITH NONLINEAR LOCALIZED DAMPING Eectronic Journa of Differentia Equations, Vo. 24(24), No. 88, pp. 1 1. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (ogin: ftp) EXPONENTIAL DECAY

More information

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations) Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =

More information

General Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping

General Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping Journa of Mathematica Research with Appications Jan.,, Vo. 3, No., pp. 53 6 DOI:.377/j.issn:95-65...7 Http://jmre.dut.edu.cn Genera Decay of Soutions in a Viscoeastic Equation with Noninear Locaized Damping

More information

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications Pricing Mutipe Products with the Mutinomia Logit and Nested Logit Modes: Concavity and Impications Hongmin Li Woonghee Tim Huh WP Carey Schoo of Business Arizona State University Tempe Arizona 85287 USA

More information

Homogeneity properties of subadditive functions

Homogeneity properties of subadditive functions Annaes Mathematicae et Informaticae 32 2005 pp. 89 20. Homogeneity properties of subadditive functions Pá Burai and Árpád Száz Institute of Mathematics, University of Debrecen e-mai: buraip@math.kte.hu

More information

Convergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems

Convergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,

More information

Conditions for Saddle-Point Equilibria in Output-Feedback MPC with MHE

Conditions for Saddle-Point Equilibria in Output-Feedback MPC with MHE Conditions for Sadde-Point Equiibria in Output-Feedback MPC with MHE David A. Copp and João P. Hespanha Abstract A new method for soving output-feedback mode predictive contro (MPC) and moving horizon

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Indirect Optimal Control of Dynamical Systems

Indirect Optimal Control of Dynamical Systems Computationa Mathematics and Mathematica Physics, Vo. 44, No. 3, 24, pp. 48 439. Transated from Zhurna Vychisite noi Matematiki i Matematicheskoi Fiziki, Vo. 44, No. 3, 24, pp. 444 466. Origina Russian

More information

<C 2 2. λ 2 l. λ 1 l 1 < C 1

<C 2 2. λ 2 l. λ 1 l 1 < C 1 Teecommunication Network Contro and Management (EE E694) Prof. A. A. Lazar Notes for the ecture of 7/Feb/95 by Huayan Wang (this document was ast LaT E X-ed on May 9,995) Queueing Primer for Muticass Optima

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

Distributed average consensus: Beyond the realm of linearity

Distributed average consensus: Beyond the realm of linearity Distributed average consensus: Beyond the ream of inearity Usman A. Khan, Soummya Kar, and José M. F. Moura Department of Eectrica and Computer Engineering Carnegie Meon University 5 Forbes Ave, Pittsburgh,

More information

Another Class of Admissible Perturbations of Special Expressions

Another Class of Admissible Perturbations of Special Expressions Int. Journa of Math. Anaysis, Vo. 8, 014, no. 1, 1-8 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.014.31187 Another Cass of Admissibe Perturbations of Specia Expressions Jerico B. Bacani

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

An Extension of Almost Sure Central Limit Theorem for Order Statistics

An Extension of Almost Sure Central Limit Theorem for Order Statistics An Extension of Amost Sure Centra Limit Theorem for Order Statistics T. Bin, P. Zuoxiang & S. Nadarajah First version: 6 December 2007 Research Report No. 9, 2007, Probabiity Statistics Group Schoo of

More information

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

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

More information

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1 Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rues 1 R.J. Marks II, S. Oh, P. Arabshahi Λ, T.P. Caude, J.J. Choi, B.G. Song Λ Λ Dept. of Eectrica Engineering Boeing Computer Services University

More information

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES

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

More information

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0 Bocconi University PhD in Economics - Microeconomics I Prof M Messner Probem Set 4 - Soution Probem : If an individua has an endowment instead of a monetary income his weath depends on price eves In particuar,

More information

Discrete Bernoulli s Formula and its Applications Arising from Generalized Difference Operator

Discrete Bernoulli s Formula and its Applications Arising from Generalized Difference Operator Int. Journa of Math. Anaysis, Vo. 7, 2013, no. 5, 229-240 Discrete Bernoui s Formua and its Appications Arising from Generaized Difference Operator G. Britto Antony Xavier 1 Department of Mathematics,

More information

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5]. PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

Appendix for Stochastic Gradient Monomial Gamma Sampler

Appendix for Stochastic Gradient Monomial Gamma Sampler 3 4 5 6 7 8 9 3 4 5 6 7 8 9 3 4 5 6 7 8 9 3 3 3 33 34 35 36 37 38 39 4 4 4 43 44 45 46 47 48 49 5 5 5 53 54 Appendix for Stochastic Gradient Monomia Gamma Samper A The Main Theorem We provide the foowing

More information

Differential equations associated with higher-order Bernoulli numbers of the second kind

Differential equations associated with higher-order Bernoulli numbers of the second kind Goba Journa of Pure and Appied Mathematics. ISS 0973-768 Voume 2, umber 3 (206), pp. 2503 25 Research India Pubications http://www.ripubication.com/gjpam.htm Differentia equations associated with higher-order

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),

More information

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,

More information

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier

More information

Explicit overall risk minimization transductive bound

Explicit overall risk minimization transductive bound 1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,

More information

Coupling of LWR and phase transition models at boundary

Coupling of LWR and phase transition models at boundary Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of

More information

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport

More information

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE Progress In Eectromagnetics Research, PIER 30, 73 84, 001 CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE W. Lin and Z. Yu University of

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

Homework 5 Solutions

Homework 5 Solutions Stat 310B/Math 230B Theory of Probabiity Homework 5 Soutions Andrea Montanari Due on 2/19/2014 Exercise [5.3.20] 1. We caim that n 2 [ E[h F n ] = 2 n i=1 A i,n h(u)du ] I Ai,n (t). (1) Indeed, integrabiity

More information

Competitive Diffusion in Social Networks: Quality or Seeding?

Competitive Diffusion in Social Networks: Quality or Seeding? Competitive Diffusion in Socia Networks: Quaity or Seeding? Arastoo Fazei Amir Ajorou Ai Jadbabaie arxiv:1503.01220v1 [cs.gt] 4 Mar 2015 Abstract In this paper, we study a strategic mode of marketing and

More information

Symbolic models for nonlinear control systems using approximate bisimulation

Symbolic models for nonlinear control systems using approximate bisimulation Symboic modes for noninear contro systems using approximate bisimuation Giordano Poa, Antoine Girard and Pauo Tabuada Abstract Contro systems are usuay modeed by differentia equations describing how physica

More information

Math 124B January 17, 2012

Math 124B January 17, 2012 Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia

More information

On Some Basic Properties of Geometric Real Sequences

On Some Basic Properties of Geometric Real Sequences On Some Basic Properties of eometric Rea Sequences Khirod Boruah Research Schoar, Department of Mathematics, Rajiv andhi University Rono His, Doimukh-791112, Arunacha Pradesh, India Abstract Objective

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

arxiv: v1 [math.ap] 8 Nov 2014

arxiv: v1 [math.ap] 8 Nov 2014 arxiv:1411.116v1 [math.ap] 8 Nov 014 ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS FOR m-component REACTION-DIFFUSION SYSTEMS WITH A TRIDIAGONAL SYMMETRIC TOEPLITZ MATRIX OF DIFFUSION COEFFICIENTS SALEM ABDELMALEK

More information

Appendix for Stochastic Gradient Monomial Gamma Sampler

Appendix for Stochastic Gradient Monomial Gamma Sampler Appendix for Stochastic Gradient Monomia Gamma Samper A The Main Theorem We provide the foowing theorem to characterize the stationary distribution of the stochastic process with SDEs in (3) Theorem 3

More information

ANALOG OF HEAT EQUATION FOR GAUSSIAN MEASURE OF A BALL IN HILBERT SPACE GYULA PAP

ANALOG OF HEAT EQUATION FOR GAUSSIAN MEASURE OF A BALL IN HILBERT SPACE GYULA PAP ANALOG OF HEAT EQUATION FOR GAUSSIAN MEASURE OF A BALL IN HILBERT SPACE GYULA PAP ABSTRACT. If µ is a Gaussian measure on a Hibert space with mean a and covariance operator T, and r is a} fixed positive

More information

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs

More information

arxiv: v1 [math.oc] 21 Sep 2018

arxiv: v1 [math.oc] 21 Sep 2018 Differentia Dynamic Programming for Noninear Dynamic Games Boei Di and Andrew Lampersi arxiv:8090830v mathoc Sep 08 Abstract Dynamic games arise when mutipe agents with differing objectives choose contro

More information

Semidefinite relaxation and Branch-and-Bound Algorithm for LPECs

Semidefinite relaxation and Branch-and-Bound Algorithm for LPECs Semidefinite reaxation and Branch-and-Bound Agorithm for LPECs Marcia H. C. Fampa Universidade Federa do Rio de Janeiro Instituto de Matemática e COPPE. Caixa Posta 68530 Rio de Janeiro RJ 21941-590 Brasi

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

Secure Information Flow Based on Data Flow Analysis

Secure Information Flow Based on Data Flow Analysis SSN 746-7659, Engand, UK Journa of nformation and Computing Science Vo., No. 4, 007, pp. 5-60 Secure nformation Fow Based on Data Fow Anaysis Jianbo Yao Center of nformation and computer, Zunyi Norma Coege,

More information

Tikhonov Regularization for Nonlinear Complementarity Problems with Approximative Data

Tikhonov Regularization for Nonlinear Complementarity Problems with Approximative Data Internationa Mathematica Forum, 5, 2010, no. 56, 2787-2794 Tihonov Reguarization for Noninear Compementarity Probems with Approximative Data Nguyen Buong Vietnamese Academy of Science and Technoogy Institute

More information

Minimum Enclosing Circle of a Set of Fixed Points and a Mobile Point

Minimum Enclosing Circle of a Set of Fixed Points and a Mobile Point Minimum Encosing Circe of a Set of Fixed Points and a Mobie Point Aritra Banik 1, Bhaswar B. Bhattacharya 2, and Sandip Das 1 1 Advanced Computing and Microeectronics Unit, Indian Statistica Institute,

More information

EXISTENCE OF SOLUTIONS FOR ONE-DIMENSIONAL WAVE EQUATIONS WITH NONLOCAL CONDITIONS

EXISTENCE OF SOLUTIONS FOR ONE-DIMENSIONAL WAVE EQUATIONS WITH NONLOCAL CONDITIONS Eectronic Journa of Differentia Equations, Vo. 21(21), No. 76, pp. 1 8. ISSN: 172-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (ogin: ftp) EXISTENCE OF SOLUTIONS

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

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

More information

Dynamic Savings Choices with Disagreements

Dynamic Savings Choices with Disagreements Dynamic Savings Choices with Disagreements Dan Cao Georgetown University Iván Werning MIT November 2017 We study a fexibe dynamic savings game in continuous time, where decision makers rotate in and out

More information

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia

More information

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

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

More information

Summation of p-adic Functional Series in Integer Points

Summation of p-adic Functional Series in Integer Points Fiomat 31:5 (2017), 1339 1347 DOI 10.2298/FIL1705339D Pubished by Facuty of Sciences and Mathematics, University of Niš, Serbia Avaiabe at: http://www.pmf.ni.ac.rs/fiomat Summation of p-adic Functiona

More information

Many-Help-One Problem for Gaussian Sources with a Tree Structure on their Correlation

Many-Help-One Problem for Gaussian Sources with a Tree Structure on their Correlation Many-Hep-One Probem for Gaussian Sources with a Tree Structure on their Correation Yasutada Oohama arxiv:090988v csit 4 Jan 009 Abstract In this paper we consider the separate coding probem for + correated

More information

Stability Analysis for a Class of Non-Weakly Reversible Chemical Reaction Networks

Stability Analysis for a Class of Non-Weakly Reversible Chemical Reaction Networks Specia Issue on SICE Annua Conference 217 SICE Journa of Contro, Measurement, and System Integration, Vo. 11, No. 3, pp. 138 145, May 218 Stabiity Anaysis for a Cass of Non-eay Reversibe Chemica Reaction

More information

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997 Stochastic utomata etworks (S) - Modeing and Evauation Pauo Fernandes rigitte Pateau 2 May 29, 997 Institut ationa Poytechnique de Grenobe { IPG Ecoe ationae Superieure d'informatique et de Mathematiques

More information

arxiv: v1 [cs.sy] 31 May 2014

arxiv: v1 [cs.sy] 31 May 2014 Cooperative Contro of Linear Muti-Agent Systems via Distributed Output Reguation and Transient Synchronization Georg S. Seyboth, Wei Ren, Fran Agöwer arxiv:46.85v cs.sy 3 May 24 Abstract A wide range of

More information

Gaussian Curvature in a p-orbital, Hydrogen-like Atoms

Gaussian Curvature in a p-orbital, Hydrogen-like Atoms Advanced Studies in Theoretica Physics Vo. 9, 015, no. 6, 81-85 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/astp.015.5115 Gaussian Curvature in a p-orbita, Hydrogen-ike Atoms Sandro-Jose Berrio-Guzman

More information

17 Lecture 17: Recombination and Dark Matter Production

17 Lecture 17: Recombination and Dark Matter Production PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was

More information

ON THE POSITIVITY OF SOLUTIONS OF SYSTEMS OF STOCHASTIC PDES

ON THE POSITIVITY OF SOLUTIONS OF SYSTEMS OF STOCHASTIC PDES ON THE POSITIVITY OF SOLUTIONS OF SYSTEMS OF STOCHASTIC PDES JACKY CRESSON 1,2, MESSOUD EFENDIEV 3, AND STEFANIE SONNER 3,4 On the occasion of the 75 th birthday of Prof. Dr. Dr.h.c. Wofgang L. Wendand

More information

#A48 INTEGERS 12 (2012) ON A COMBINATORIAL CONJECTURE OF TU AND DENG

#A48 INTEGERS 12 (2012) ON A COMBINATORIAL CONJECTURE OF TU AND DENG #A48 INTEGERS 12 (2012) ON A COMBINATORIAL CONJECTURE OF TU AND DENG Guixin Deng Schoo of Mathematica Sciences, Guangxi Teachers Education University, Nanning, P.R.China dengguixin@ive.com Pingzhi Yuan

More information

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES

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

More information

Input-to-state stability for a class of Lurie systems

Input-to-state stability for a class of Lurie systems Automatica 38 (2002) 945 949 www.esevier.com/ocate/automatica Brief Paper Input-to-state stabiity for a cass of Lurie systems Murat Arcak a;, Andrew Tee b a Department of Eectrica, Computer and Systems

More information

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

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

More information

Approximated MLC shape matrix decomposition with interleaf collision constraint

Approximated MLC shape matrix decomposition with interleaf collision constraint Approximated MLC shape matrix decomposition with intereaf coision constraint Thomas Kainowski Antje Kiese Abstract Shape matrix decomposition is a subprobem in radiation therapy panning. A given fuence

More information

Demand in Leisure Markets

Demand in Leisure Markets Demand in Leisure Markets An Empirica Anaysis of Time Aocation Shomi Pariat Ph.D Candidate Eitan Bergas Schoo of Economics Te Aviv University Motivation Leisure activities matter 35% of waking time 9.7%

More information

Quantum Electrodynamical Basis for Wave. Propagation through Photonic Crystal

Quantum Electrodynamical Basis for Wave. Propagation through Photonic Crystal Adv. Studies Theor. Phys., Vo. 6, 01, no. 3, 19-133 Quantum Eectrodynamica Basis for Wave Propagation through Photonic Crysta 1 N. Chandrasekar and Har Narayan Upadhyay Schoo of Eectrica and Eectronics

More information

Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay

Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay Throughput Optima Scheduing for Wireess Downinks with Reconfiguration Deay Vineeth Baa Sukumaran vineethbs@gmai.com Department of Avionics Indian Institute of Space Science and Technoogy. Abstract We consider

More information

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,

More information

Proceedings of the 2012 Winter Simulation Conference C. Laroque, J. Himmelspach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds.

Proceedings of the 2012 Winter Simulation Conference C. Laroque, J. Himmelspach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds. Proceedings of the 2012 Winter Simuation Conference C. aroque, J. Himmespach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds. TIGHT BOUNDS FOR AMERICAN OPTIONS VIA MUTIEVE MONTE CARO Denis Beomestny Duisburg-Essen

More information