Robust Control of Linear Quantum Systems

Size: px
Start display at page:

Download "Robust Control of Linear Quantum Systems"

Transcription

1 B. Ross Barmish Workshop Robust Control of Linear Quantum Systems Ian R. Petersen School of ITEE, University of New South the Australian Defence Force Academy, Based on joint work with Matthew James, A.J. Shaiju, Hendra Nurdin, and Aline Maalouf.

2 B. Ross Barmish Workshop Introduction Recent developments in quantum and nano technology have provided a great impetus for research in quantum feedback control systems. In particular, it is being realized that robustness is critical in quantum feedback control systems, as it is in non-quantum feedback control systems. We present an overview of some recent results in the area of robust feedback control of linear quantum systems. We consider a class of quantum systems described by linear Heisenberg dynamics driven by quantum Gaussian noise processes, and controlled by a linear feedback controller which is also a quantum system.

3 B. Ross Barmish Workshop The most common area in which such systems arise is in the area of quantum optics; e.g., Photo courtesy of Elanor Huntington (UNSW@ADFA)

4 B. Ross Barmish Workshop Analysis of Linear Quantum Systems In the analysis of linear quantum systems, we consider noncommutative stochastic models of the form dx(t) = Ax(t)dt + [ B G ][ dw(t) T dv(t) T ] T ; dz(t) = Cx(t)dt + [ D H ][ dw(t) T dv(t) T ] T (1) where x(t) = [ x 1 (t)... x n (t) ] T is a vector of self-adjoint possibly non-commutative system variables. Also, A, B, G, C D are real matrices. The vector quantity w describes the input signals and is assumed to admit the decomposition dw(t) = β w (t)dt + d w(t) (2) where w(t) is the noise part of w(t) and β w (t) is a self adjoint process.

5 B. Ross Barmish Workshop In this quantum system, the input channel has two components, dw = β w dt + d w which represents disturbance signals, and dv, which represents additional noise sources. The process β w (t) serves to represent variables of other systems which may be passed to the system. This system may represent the closed loop quantum classical system in a quantum H control system.

6 B. Ross Barmish Workshop Definition. The above quantum stochastic system is said to be Strictly Bounded Real with disturbance attenuation g if there exists a positive operator valued quadratic form V (x) = x T Xx (where X is a real positive definite symmetric matrix) and constants λ > 0, ǫ > 0 such that V (x(t)) + t V (x(0)) + λt t > 0, 0 β T z β z (g 2 ǫ)β T wβ w + ǫx T x ds for all Gaussian states ρ for the initial variables x(0). Here we use the shorthand notation for quantum expectation over all initial variables and noises.

7 B. Ross Barmish Workshop Theorem. [Quantum Strict Bounded Real Lemma] The following statements are equivalent (i) The above quantum stochastic system is strictly bounded real with disturbance attenuation g. (ii) A is a stable matrix and C(sI A) 1 B + D < g. (iii) g 2 I D T D > 0 and there exists a positive definite matrix X > 0 : A T X + XA + C T C + ( XB + C T D) (g 2 I D T D) 1 (B T X + D T C) < 0. (iv) g 2 I D T D > 0 and the algebraic Riccati equation A T X + XA + C T C + (XB + C T D) (g 2 I D T D) 1 (B T X + D T C) = 0 has a stabilizing solution X 0. Furthermore, if these statements hold then X < X.

8 B. Ross Barmish Workshop H of linear quantum systems In the quantum H control problem we consider a plant which is described by noncommutative stochastic models of the form dx(t) = Ax(t)dt + B 0 dv(t) + B 1 dw(t) + B 2 du(t); x(0) = x; dz(t) = C 1 x(t)dt + D 12 du(t); dy(t) = C 2 x(t)dt + D 20 dv(t) + D 21 dw(t) We assume that the control input is of the form du(t) = β u (t)dt + dũ(t) where ũ(t) is the noise part of u(t).

9 B. Ross Barmish Workshop Assumptions 1. D T 12D 12 = E 1 > D 21 D21 T = E 2 > 0. [ ] A jωi B2 3. The matrix is full rank for all ω 0. C 1 D 12 [ ] A jωi B1 4. The matrix is full rank for all ω 0. C 2 D 21

10 B. Ross Barmish Workshop Riccati Equations Our quantum H result is stated in terms of the following pair of algebraic Riccati equations: (A B 2 E1 1 DT 12C 1 ) T X + X(A B 2 E1 1 DT 12C 1 ) +X(B 1 B1 T g 2 B 2 E1 1 B 2)X +g 2 C1 T (I D 12E1 1 DT 12 )C 1 = 0; (A B 1 D T 21E 1 2 C 2)Y + Y (A B 1 D T 21E 1 2 C 2) +Y (g 2 C T 1 C 1 C T 2 E 1 2 C 2)Y +B 1 (I D T 21E 1 2 D 21)B T 1 = 0.

11 B. Ross Barmish Workshop Controller Our quantum H control is defined by state space matrices constructed from the Riccati equations as follows: A K = A + B 2 C K B K C 2 + (B 1 B K D 21 )B T 1 X; B K = (I Y X) 1 (Y C T 2 + B 1 D T 21)E 1 2 ; C K = E 1 1 (g2 B T 2 X + D T 12C 1 ).

12 B. Ross Barmish Workshop Theorem. Necessity. If there exists a controller such that the resulting closed loop system is strictly bounded real with disturbance attenuation g, then the Riccati equations will have stabilizing solutions X 0 and Y 0 such that the matrix XY has a spectral radius strictly less than one. Sufficiency. Suppose the Riccati equations have stabilizing solutions X 0 and Y 0 such that the matrix XY has a spectral radius strictly less than one. If the controller is such that the matrices A K, B K, C K are as defined in terms of the Riccati solutions as above, then the resulting closed loop system will be strictly bounded real with disturbance attenuation g.

13 B. Ross Barmish Workshop Physical Realizability of Linear Quantum Systems In the quantum H control problem, the controller can be a classical system constructed from digital or analog electronics, a quantum system constructed from cavities, beamsplitters and phase shifters, or a combination of both. If the controller is to be a purely quantum system (coherent quantum feedback control), the question arises as to whether a given controller transfer function can be physically realized as a quantum system.

14 B. Ross Barmish Workshop Example Quantum H control system consisting of an optical cavity plant and an optical cavity controller: Photo courtesy of Hideo Mabuchi, Stanford.

15 B. Ross Barmish Workshop Physical Realizability of Quantum System Models To consider the issue of whether a state space quantum system model (controller) is physically realizable, it is most convenient to consider complex quantum systems of the form da = F 0 [ a a ] db = H 0 [ a a ] dt + G 0 [ dv dv ], dt + K 0 [ dv dv ]. We assume without loss of generality that the number of inputs equals the number of outputs.

16 B. Ross Barmish Workshop Also, we write F 0 = [ ] F 1 F 2 C n a 2n a, G 0 = [ ] G 1 G 2 C n a 2n v, H 0 = [ ] H 1 H 2 C n v 2n a, K 0 = [ ] K 1 K 2 C n v 2n v, [ ] [ ] F1 F F = 2 G1 G F2 F1, G = 2 G 2 G, 1 [ ] [ ] H1 H H = 2 K1 K, K = 2. H 2 H 1 K 2 K 1

17 B. Ross Barmish Workshop This complex linear quantum system is equivalent to a real linear quantum system of the form considered previously via the substitutions: x k = a k + a k, x n a +k = i(a k a k ). y l = b l + b l, y nv +l = i(b l b l), 1 l n v, w j = v j + v j, w nv +j = i(v j v j ), 1 j n v.

18 B. Ross Barmish Workshop This leads to the following system with self-adjoint variables and real coefficient matrices: dx = Ax dt + B dw, dy = Cx dt + D dw, where x = [x 1 x 2 x 2na ], y = [y 1 y 2 y 2nv ], A = 1 2 ΦFΦ, B = 1 2 ΦGΦ, C = 1 2 ΦHΦ, D = 1 2 ΦKΦ, [ ] [ ] I I I 0 Φ =, J =. ii ii 0 I

19 B. Ross Barmish Workshop Our notion of physical realizability is defined in terms of the preservation of certain quantum commutation relations and the equivalence of a quantum linear system to a quantum harmonic oscillator. An alternative algebraic condition is given in the following theorem. Theorem. The above complex linear quantum system is physically realizable if and only if there exists a hermitian matrix Ω = Ω such that FΩ + ΩF GJG = 0, G = ΩH J, K = I.

20 B. Ross Barmish Workshop Our main result on physical realizability relates this property to the system theory notion of a (J, J)-unitary transfer function. Consider the complex transfer function matrix corresponding to the above complex linear quantum system Ψ(s) = H(sI F) 1 G + K. Definition. A system with transfer function matrix Ψ(s) is said to be is said to be dual (J,J)-unitary if for s C : Re(s) 0. Ψ(s)JΨ ( s ) = J,

21 B. Ross Barmish Workshop Theorem. The above complex linear quantum system is physically realizable if and only if its transfer function matrix is dual (J,J)-unitary and K = I. This theorem shows the equivalence between the quantum mechanical notion of physical realizability and the frequency domain systems theory notion of dual (J, J)-unitary.

Robust Linear Quantum Systems Theory

Robust Linear Quantum Systems Theory Workshop On Uncertain Dynamical Systems Udine 2011 1 Robust Linear Quantum Systems Theory Ian R. Petersen School of Engineering and Information Technology, University of New South Wales @ the Australian

More information

Quantum Linear Systems Theory

Quantum Linear Systems Theory RMIT 2011 1 Quantum Linear Systems Theory Ian R. Petersen School of Engineering and Information Technology, University of New South Wales @ the Australian Defence Force Academy RMIT 2011 2 Acknowledgments

More information

Control and Robustness for Quantum Linear Systems

Control and Robustness for Quantum Linear Systems CCC 2013 1 Control and Robustness for Quantum Linear Systems Ian R. Petersen School of Engineering and Information Technology, UNSW Canberra CCC 2013 2 Introduction Developments in quantum technology and

More information

Feedback Control of Quantum Systems. Matt James ANU

Feedback Control of Quantum Systems. Matt James ANU Feedback Control of Quantum Systems Matt James ANU Dedicated to Slava Belavkin who pioneered the way Outline What is quantum control? A Little History Optimal Control Using Quantum Langevin Equations H-Infinity

More information

Chapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control

Chapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design

More information

Quantum Linear Systems Theory

Quantum Linear Systems Theory Proceedings of the 19th International Symposium on Mathematical Theory of Networks and Systems MTNS 010 5 9 July, 010 Budapest, Hungary Quantum Linear Systems Theory Ian R Petersen Abstract This paper

More information

arxiv: v1 [quant-ph] 18 Sep 2015

arxiv: v1 [quant-ph] 18 Sep 2015 Tangential Interpolatory Projection for Model Reduction of Linear Quantum Stochastic Systems Onvaree Techakesari and Hendra I. Nurdin September 20, 2018 arxiv:1509.05534v1 [quant-ph] 18 Sep 2015 Abstract

More information

Semidefinite Programming Duality and Linear Time-invariant Systems

Semidefinite Programming Duality and Linear Time-invariant Systems Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,

More information

Introduction to numerical simulations for Stochastic ODEs

Introduction to numerical simulations for Stochastic ODEs Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical

More information

arxiv: v1 [quant-ph] 16 Mar 2016

arxiv: v1 [quant-ph] 16 Mar 2016 Quantum Linear Systems heory Ian R. Petersen arxiv:163.495v1 quant-ph 16 Mar 16 Abstract his paper surveys some recent results on the theory of quantum linear systems and presents them within a unified

More information

COHERENT CONTROL VIA QUANTUM FEEDBACK NETWORKS

COHERENT CONTROL VIA QUANTUM FEEDBACK NETWORKS COHERENT CONTROL VIA QUANTUM FEEDBACK NETWORKS Kavli Institute for Theoretical Physics, Santa Barbara, 2013 John Gough Quantum Structures, Information and Control, Aberystwyth Papers J.G. M.R. James, Commun.

More information

The norms can also be characterized in terms of Riccati inequalities.

The norms can also be characterized in terms of Riccati inequalities. 9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements

More information

arxiv: v1 [quant-ph] 13 Nov 2016

arxiv: v1 [quant-ph] 13 Nov 2016 Mixed LQG and H Coherent Feedback Control for Linear Quantum Systems arxiv:1611.04140v1 [quant-ph] 13 Nov 016 Lei Cui, Zhiyuan Dong, Guofeng Zhang and Heung Wing Joseph Lee Department of Applied Mathematics,

More information

6.241 Dynamic Systems and Control

6.241 Dynamic Systems and Control 6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May

More information

Linear Filtering of general Gaussian processes

Linear Filtering of general Gaussian processes Linear Filtering of general Gaussian processes Vít Kubelka Advisor: prof. RNDr. Bohdan Maslowski, DrSc. Robust 2018 Department of Probability and Statistics Faculty of Mathematics and Physics Charles University

More information

Schrodinger Bridges classical and quantum evolution

Schrodinger Bridges classical and quantum evolution .. Schrodinger Bridges classical and quantum evolution Tryphon Georgiou University of Minnesota! Joint work with Michele Pavon and Yongxin Chen! LCCC Workshop on Dynamic and Control in Networks Lund, October

More information

arxiv: v1 [quant-ph] 20 Jun 2017

arxiv: v1 [quant-ph] 20 Jun 2017 A Direct Couling Coherent Quantum Observer for an Oscillatory Quantum Plant Ian R Petersen arxiv:76648v quant-h Jun 7 Abstract A direct couling coherent observer is constructed for a linear quantum lant

More information

EL2520 Control Theory and Practice

EL2520 Control Theory and Practice EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller

More information

Robust Control 9 Design Summary & Examples

Robust Control 9 Design Summary & Examples Robust Control 9 Design Summary & Examples Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/6/003 Outline he H Problem Solution via γ-iteration Robust stability via

More information

Schrodinger Bridges classical and quantum evolution

Schrodinger Bridges classical and quantum evolution Schrodinger Bridges classical and quantum evolution Tryphon Georgiou University of Minnesota! Joint work with Michele Pavon and Yongxin Chen! Newton Institute Cambridge, August 11, 2014 http://arxiv.org/abs/1405.6650

More information

Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013

Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,

More information

Modern Optimal Control

Modern Optimal Control Modern Optimal Control Matthew M. Peet Arizona State University Lecture 22: H 2, LQG and LGR Conclusion To solve the H -optimal state-feedback problem, we solve min γ such that γ,x 1,Y 1,A n,b n,c n,d

More information

LQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin

LQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system

More information

6.241 Dynamic Systems and Control

6.241 Dynamic Systems and Control 6.41 Dynamic Systems and Control Lecture 5: H Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 11, 011 E. Frazzoli (MIT) Lecture 5: H Synthesis May 11, 011

More information

Introduction to Quantum Linear Systems

Introduction to Quantum Linear Systems The Hong Kong Polytechnic University Guofeng.Zhang@polyu.edu.hk 29 July 2014 Quantum Control Engineering, Cambridge, 21 July - 15 August 2014 Classical linear systems mathematical model (via first-principle

More information

Robust Control 5 Nominal Controller Design Continued

Robust Control 5 Nominal Controller Design Continued Robust Control 5 Nominal Controller Design Continued Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 4/14/2003 Outline he LQR Problem A Generalization to LQR Min-Max

More information

Quantum Dissipative Systems and Feedback Control Design by Interconnection

Quantum Dissipative Systems and Feedback Control Design by Interconnection 1 Quantum Dissipative Systems and Feedback Control Design by Interconnection Matthew R. James, John Gough arxiv:0707.1074v3 [quant-ph] 7 May 2009 Abstract The purpose of this paper is to extend J.C. Willems

More information

Quantum Dissipative Systems and Feedback Control Design by Interconnection

Quantum Dissipative Systems and Feedback Control Design by Interconnection 1 Quantum Dissipative Systems and Feedback Control Design by Interconnection Matthew R. James, Fellow, IEEE John Gough Abstract The purpose of this paper is to extend J.C. Willems theory of dissipative

More information

Unconditionally stable scheme for Riccati equation

Unconditionally stable scheme for Riccati equation ESAIM: Proceedings, Vol. 8, 2, 39-52 Contrôle des systèmes gouvernés par des équations aux dérivées partielles http://www.emath.fr/proc/vol.8/ Unconditionally stable scheme for Riccati equation François

More information

Coherent controllers for optical-feedback cooling of quantum oscillators

Coherent controllers for optical-feedback cooling of quantum oscillators Coherent controllers for optical-feedback cooling of quantum oscillators Ryan Hamerly and Hideo Mabuchi Edward L. Ginzton Laboratory, Stanford University, Stanford, CA 94305 (Dated: June 12, 2012) We study

More information

Control, Stabilization and Numerics for Partial Differential Equations

Control, Stabilization and Numerics for Partial Differential Equations Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua

More information

Linear Quadratic Gausssian Control Design with Loop Transfer Recovery

Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Leonid Freidovich Department of Mathematics Michigan State University MI 48824, USA e-mail:leonid@math.msu.edu http://www.math.msu.edu/

More information

ME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ

ME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and

More information

arxiv: v3 [quant-ph] 30 Oct 2014

arxiv: v3 [quant-ph] 30 Oct 2014 Quantum feedback: theory, experiments, and applications Jing Zhang a,b,c, Yu-xi Liu d,b,c, Re-Bing Wu a,b,c, Kurt Jacobs e,c, Franco Nori c,f a Department of Automation, Tsinghua University, Beijing 100084,

More information

Week 9 Generators, duality, change of measure

Week 9 Generators, duality, change of measure Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov

More information

Ian R. Petersen. University of New South Wales, Canberra

Ian R. Petersen. University of New South Wales, Canberra A Direct Coupling Coherent Quantum Observer 1 Ian R. Petersen University of New South Wales, Canberra A Direct Coupling Coherent Quantum Observer 2 This paper considers the problem of constructing a direct

More information

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors /88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix

More information

Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection

Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute

More information

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition

More information

10 Transfer Matrix Models

10 Transfer Matrix Models MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important

More information

BV functions in a Gelfand triple and the stochastic reflection problem on a convex set

BV functions in a Gelfand triple and the stochastic reflection problem on a convex set BV functions in a Gelfand triple and the stochastic reflection problem on a convex set Xiangchan Zhu Joint work with Prof. Michael Röckner and Rongchan Zhu Xiangchan Zhu ( Joint work with Prof. Michael

More information

The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications

The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department

More information

Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra

Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear

More information

Denis ARZELIER arzelier

Denis ARZELIER   arzelier COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15

More information

Controls Problems for Qualifying Exam - Spring 2014

Controls Problems for Qualifying Exam - Spring 2014 Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function

More information

A Quantum Particle Undergoing Continuous Observation

A Quantum Particle Undergoing Continuous Observation A Quantum Particle Undergoing Continuous Observation V.P. Belavkin and P. Staszewski y December 1988 Published in: Physics Letters A, 140 (1989) No 7,8, pp 359 {362 Abstract A stochastic model for the

More information

The semantics of algebraic quantum mechanics and the role of model theory.

The semantics of algebraic quantum mechanics and the role of model theory. The semantics of algebraic quantum mechanics and the role of model theory. B. Zilber University of Oxford August 6, 2016 B.Zilber, The semantics of the canonical commutation relations arxiv.org/abs/1604.07745

More information

NLControl - package for modelling, analysis and synthesis of nonlinear control systems

NLControl - package for modelling, analysis and synthesis of nonlinear control systems NLControl - package for modelling, analysis and synthesis of nonlinear control systems Maris Tõnso Institute of Cybernetics Tallinn University of Technology Estonia maris@cc.ioc.ee May, 2013 Contents 1

More information

Stability of Stochastic Differential Equations

Stability of Stochastic Differential Equations Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010

More information

Control of elementary quantum flows Luigi Accardi Centro Vito Volterra, Università di Roma Tor Vergata Via di Tor Vergata, snc Roma, Italy

Control of elementary quantum flows Luigi Accardi Centro Vito Volterra, Università di Roma Tor Vergata Via di Tor Vergata, snc Roma, Italy Control of elementary quantum flows Luigi Accardi Centro Vito Volterra, Università di Roma Tor Vergata Via di Tor Vergata, snc 0033 Roma, Italy accardi@volterra.mat.uniroma2.it Andreas Boukas Department

More information

arxiv: v1 [quant-ph] 3 Sep 2008

arxiv: v1 [quant-ph] 3 Sep 2008 Frequency Locking of an Optical Cavity using LQG Integral Control S. Z. Sayed Hassen, M. Heurs, E. H. Huntington, I. R. Petersen University of New South Wales at the Australian Defence Force Academy, School

More information

where r n = dn+1 x(t)

where r n = dn+1 x(t) Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution

More information

Stochastic Calculus. Kevin Sinclair. August 2, 2016

Stochastic Calculus. Kevin Sinclair. August 2, 2016 Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed

More information

CONTROLLABILITY OF QUANTUM SYSTEMS. Sonia G. Schirmer

CONTROLLABILITY OF QUANTUM SYSTEMS. Sonia G. Schirmer CONTROLLABILITY OF QUANTUM SYSTEMS Sonia G. Schirmer Dept of Applied Mathematics + Theoretical Physics and Dept of Engineering, University of Cambridge, Cambridge, CB2 1PZ, United Kingdom Ivan C. H. Pullen

More information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

1. Find the solution of the following uncontrolled linear system. 2 α 1 1 Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +

More information

Linear Quadratic Zero-Sum Two-Person Differential Games

Linear Quadratic Zero-Sum Two-Person Differential Games Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard To cite this version: Pierre Bernhard. Linear Quadratic Zero-Sum Two-Person Differential Games. Encyclopaedia of Systems and Control,

More information

On Stochastic Adaptive Control & its Applications. Bozenna Pasik-Duncan University of Kansas, USA

On Stochastic Adaptive Control & its Applications. Bozenna Pasik-Duncan University of Kansas, USA On Stochastic Adaptive Control & its Applications Bozenna Pasik-Duncan University of Kansas, USA ASEAS Workshop, AFOSR, 23-24 March, 2009 1. Motivation: Work in the 1970's 2. Adaptive Control of Continuous

More information

Solution of Stochastic Optimal Control Problems and Financial Applications

Solution of Stochastic Optimal Control Problems and Financial Applications Journal of Mathematical Extension Vol. 11, No. 4, (2017), 27-44 ISSN: 1735-8299 URL: http://www.ijmex.com Solution of Stochastic Optimal Control Problems and Financial Applications 2 Mat B. Kafash 1 Faculty

More information

Numerical treatment of stochastic delay differential equations

Numerical treatment of stochastic delay differential equations Numerical treatment of stochastic delay differential equations Evelyn Buckwar Department of Mathematics Heriot-Watt University, Edinburgh Evelyn Buckwar Roma, 16.10.2007 Overview An example: machine chatter

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Short Course in Quantum Information Lecture 2

Short Course in Quantum Information Lecture 2 Short Course in Quantum Information Lecture Formal Structure of Quantum Mechanics Course Info All materials downloadable @ website http://info.phys.unm.edu/~deutschgroup/deutschclasses.html Syllabus Lecture

More information

Decentralized LQG Control of Systems with a Broadcast Architecture

Decentralized LQG Control of Systems with a Broadcast Architecture Decentralized LQG Control of Systems with a Broadcast Architecture Laurent Lessard 1,2 IEEE Conference on Decision and Control, pp. 6241 6246, 212 Abstract In this paper, we consider dynamical subsystems

More information

7.3 State Space Averaging!

7.3 State Space Averaging! 7.3 State Space Averaging! A formal method for deriving the small-signal ac equations of a switching converter! Equivalent to the modeling method of the previous sections! Uses the state-space matrix description

More information

Stochastic Differential Equations.

Stochastic Differential Equations. Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)

More information

Here represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.

Here represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities. 19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that

More information

Identification and Estimation of Quantum Linear Input-Output Systems

Identification and Estimation of Quantum Linear Input-Output Systems University of Nottingham arxiv:1712.08531v1 [quant-ph] 22 Dec 2017 Doctoral Thesis Identification and Estimation of Quantum Linear Input-Output Systems Author: Supervisors: Matthew Levitt Mădălin Guță

More information

Lecture 3. Econ August 12

Lecture 3. Econ August 12 Lecture 3 Econ 2001 2015 August 12 Lecture 3 Outline 1 Metric and Metric Spaces 2 Norm and Normed Spaces 3 Sequences and Subsequences 4 Convergence 5 Monotone and Bounded Sequences Announcements: - Friday

More information

1 Linear Algebra Problems

1 Linear Algebra Problems Linear Algebra Problems. Let A be the conjugate transpose of the complex matrix A; i.e., A = A t : A is said to be Hermitian if A = A; real symmetric if A is real and A t = A; skew-hermitian if A = A and

More information

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7) EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and

More information

SEA s workshop- MIT - July 10-14

SEA s workshop- MIT - July 10-14 Matrix-valued Stochastic Processes- Eigenvalues Processes and Free Probability SEA s workshop- MIT - July 10-14 July 13, 2006 Outline Matrix-valued stochastic processes. 1- definition and examples. 2-

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechanics Rajdeep Sensarma sensarma@theory.tifr.res.in Quantum Dynamics Lecture #2 Recap of Last Class Schrodinger and Heisenberg Picture Time Evolution operator/ Propagator : Retarded

More information

Decay rates for partially dissipative hyperbolic systems

Decay rates for partially dissipative hyperbolic systems Outline Decay rates for partially dissipative hyperbolic systems Basque Center for Applied Mathematics Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Numerical Methods

More information

AN INTRODUCTION CONTROL +

AN INTRODUCTION CONTROL + AN INTRODUCTION CONTROL + QUANTUM CIRCUITS AND COHERENT QUANTUM FEEDBACK CONTROL Josh Combes The Center for Quantum Information and Control, University of New Mexico Quantum feedback control Parameter

More information

Zeros and zero dynamics

Zeros and zero dynamics CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)

More information

Lecture 15: H Control Synthesis

Lecture 15: H Control Synthesis c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control for Linear Systems: Lecture 15 p. 1/14 Lecture 15: H Control Synthesis Example c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control

More information

LECTURE 7. k=1 (, v k)u k. Moreover r

LECTURE 7. k=1 (, v k)u k. Moreover r LECTURE 7 Finite rank operators Definition. T is said to be of rank r (r < ) if dim T(H) = r. The class of operators of rank r is denoted by K r and K := r K r. Theorem 1. T K r iff T K r. Proof. Let T

More information

Return Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems

Return Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems Spectral Properties of Linear- Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2018! Stability margins of single-input/singleoutput (SISO) systems! Characterizations

More information

Regularization by noise in infinite dimensions

Regularization by noise in infinite dimensions Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of

More information

A Comparative Study on Automatic Flight Control for small UAV

A Comparative Study on Automatic Flight Control for small UAV Proceedings of the 5 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'18) Niagara Falls, Canada June 7 9, 18 Paper No. 13 DOI: 1.11159/cdsr18.13 A Comparative Study on Automatic

More information

Stochastic Processes. A stochastic process is a function of two variables:

Stochastic Processes. A stochastic process is a function of two variables: Stochastic Processes Stochastic: from Greek stochastikos, proceeding by guesswork, literally, skillful in aiming. A stochastic process is simply a collection of random variables labelled by some parameter:

More information

Convergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics

Convergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics Convergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics Meng Xu Department of Mathematics University of Wyoming February 20, 2010 Outline 1 Nonlinear Filtering Stochastic Vortex

More information

Lyapunov-based control of quantum systems

Lyapunov-based control of quantum systems Lyapunov-based control of quantum systems Symeon Grivopoulos Bassam Bamieh Department of Mechanical and Environmental Engineering University of California, Santa Barbara, CA 936-57 symeon,bamieh@engineering.ucsb.edu

More information

Linear Algebra: Matrix Eigenvalue Problems

Linear Algebra: Matrix Eigenvalue Problems CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given

More information

Optimal Polynomial Control for Discrete-Time Systems

Optimal Polynomial Control for Discrete-Time Systems 1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should

More information

Linear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters

Linear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters Linear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 204 Emo Todorov (UW) AMATH/CSE 579, Winter

More information

Analog Signals and Systems and their properties

Analog Signals and Systems and their properties Analog Signals and Systems and their properties Main Course Objective: Recall course objectives Understand the fundamentals of systems/signals interaction (know how systems can transform or filter signals)

More information

Numerical methods for solving stochastic differential equations

Numerical methods for solving stochastic differential equations Mathematical Communications 4(1999), 251-256 251 Numerical methods for solving stochastic differential equations Rózsa Horváth Bokor Abstract. This paper provides an introduction to stochastic calculus

More information

arzelier

arzelier COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.1 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS STABILITY ANALYSIS Didier HENRION www.laas.fr/ henrion henrion@laas.fr Denis ARZELIER www.laas.fr/

More information

OPTIMAL CONTROL AND ESTIMATION

OPTIMAL CONTROL AND ESTIMATION OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION

More information

Perturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop

Perturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael

More information

Multivariate Generalized Ornstein-Uhlenbeck Processes

Multivariate Generalized Ornstein-Uhlenbeck Processes Multivariate Generalized Ornstein-Uhlenbeck Processes Anita Behme TU München Alexander Lindner TU Braunschweig 7th International Conference on Lévy Processes: Theory and Applications Wroclaw, July 15 19,

More information

Modern Optimal Control

Modern Optimal Control Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints

More information

Suppose that we have a specific single stage dynamic system governed by the following equation:

Suppose that we have a specific single stage dynamic system governed by the following equation: Dynamic Optimisation Discrete Dynamic Systems A single stage example Suppose that we have a specific single stage dynamic system governed by the following equation: x 1 = ax 0 + bu 0, x 0 = x i (1) where

More information

Definition of a Stochastic Process

Definition of a Stochastic Process Definition of a Stochastic Process Balu Santhanam Dept. of E.C.E., University of New Mexico Fax: 505 277 8298 bsanthan@unm.edu August 26, 2018 Balu Santhanam (UNM) August 26, 2018 1 / 20 Overview 1 Stochastic

More information

A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES

A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES STEFAN TAPPE Abstract. In a work of van Gaans (25a) stochastic integrals are regarded as L 2 -curves. In Filipović and Tappe (28) we have shown the connection

More information

On the linear quadratic optimization problems and associated Riccati equations for systems modeled by Ito linear differential equations

On the linear quadratic optimization problems and associated Riccati equations for systems modeled by Ito linear differential equations Innovativity in Modeling and Analytics Journal of Research vol. 1, 216, pp.13-33, http://imajor.info/ c imajor ISSN XXXX-XXXX On the linear quadratic optimization problems and associated Riccati equations

More information

Risk-Sensitive and Robust Mean Field Games

Risk-Sensitive and Robust Mean Field Games Risk-Sensitive and Robust Mean Field Games Tamer Başar Coordinated Science Laboratory Department of Electrical and Computer Engineering University of Illinois at Urbana-Champaign Urbana, IL - 6181 IPAM

More information

Control Design of a Distributed Parameter. Fixed-Bed Reactor

Control Design of a Distributed Parameter. Fixed-Bed Reactor Workshop on Control of Distributed Parameter Systems p. 1/1 Control Design of a Distributed Parameter Fixed-Bed Reactor Ilyasse Aksikas and J. Fraser Forbes Department of Chemical and Material Engineering

More information

Lifted approach to ILC/Repetitive Control

Lifted approach to ILC/Repetitive Control Lifted approach to ILC/Repetitive Control Okko H. Bosgra Maarten Steinbuch TUD Delft Centre for Systems and Control TU/e Control System Technology Dutch Institute of Systems and Control DISC winter semester

More information

Networked Sensing, Estimation and Control Systems

Networked Sensing, Estimation and Control Systems Networked Sensing, Estimation and Control Systems Vijay Gupta University of Notre Dame Richard M. Murray California Institute of echnology Ling Shi Hong Kong University of Science and echnology Bruno Sinopoli

More information