Ian R. Petersen. University of New South Wales, Canberra
|
|
- Marilynn Stewart
- 6 years ago
- Views:
Transcription
1 A Direct Coupling Coherent Quantum Observer 1 Ian R. Petersen University of New South Wales, Canberra
2 A Direct Coupling Coherent Quantum Observer 2 This paper considers the problem of constructing a direct coupling quantum observer for a closed linear quantum system. The proposed observer is shown to be able to estimate some but not all of the plant variables in a time averaged sense. A simple example and simulations are included to illustrate the properties of the observer.
3 A Direct Coupling Coherent Quantum Observer 3 In most previous papers on quantum observers the coupling between the plant and the observer is via a field coupling. This leads to an observer structure of the form shown below. quantum plant quantum noise quantum observer quantum noise This enables a one way connection between the quantum plant and the quantum observer. Also, since both the quantum plant and the quantum observer are open quantum systems, they are both subject to quantum noise.
4 A Direct Coupling Coherent Quantum Observer 4 However in a recent paper by Zhang and James, a coherent quantum control problem is considered in which both field coupling and direct coupling is considered between the quantum plant and the quantum controller. In this paper, we explore the construction of a coherent quantum observer in which there is only direct coupling between quantum plant and the quantum observer. Furthermore, both the quantum plant and the quantum observer are assumed to be closed quantum systems which means that they are not subject to quantum noise and are purely deterministic systems. This leads to an observer structure of the form shown below. quantum plant quantum observer It is shown that for the case being considered, a quantum observer can be constructed to estimate some but not all of the system variables of the quantum plant. Also, the observer variables converge to the plant variables in a time averaged sense rather than a quantum expectation sense.
5 A Direct Coupling Coherent Quantum Observer 5 We consider linear non-commutative systems of the form ẋ(t) = Ax(t); x() = x (1) whereais a real matrix inr n n, andx(t) = [x 1 (t)... x n (t)] T is a vector of self-adjoint possibly non-commutative system variables. Herenis assumed to be an even number and n 2 quantum system. is the number of modes in the The initial system variablesx() = x are assumed to satisfy the commutation relations [x j (),x k ()] = 2iΘ jk, j,k = 1,...,n, whereθis a real antisymmetric matrix with componentsθ jk. Here, the matrixθis assumed to be of the formθ = diag(j,j,...,j) wherej denotes the real skew-symmetric2 2matrix J = [ 1 1 ].
6 A Direct Coupling Coherent Quantum Observer 6 This linear quantum system is said to be physically realizable if it ensures the preservation the canonical commutation relations (CCRs): x(t)x(t) T (x(t)x(t) T ) T = 2iΘ for allt. This holds when the system corresponds to a collection of closed quantum harmonic oscillators. Such quantum harmonic oscillators are described by a quadratic Hamiltonian H = 1 2 x()t Rx(), whereris a real symmetric matrix. Theorem 1 The above system is physically realizable if and only if: AΘ+ΘA T =. In this case, the corresponding Hamiltonian matrixris given byr = 1 4 ( ΘA+AT Θ). In addition, for a given Hamiltonian matrixr, the corresponding matrixais given by A = 2ΘR.
7 A Direct Coupling Coherent Quantum Observer 7 Note that the above system cannot be asymptotically stable if it is physically realizable. Since it is not possible for a physically realizable quantum system of the form above to be asymptotically stable, we will need a new notion of convergence for our direct coupled quantum observer.
8 A Direct Coupling Coherent Quantum Observer 8 We first consider general closed linear quantum plants described by non-commutative models of the following form: ẋ p (t) = A p x p (t); x p () = x p ; z p (t) = C p x p (t) wherez p denotes the vector of system variables to be estimated by the observer and A p R n p n p,c p R m p n p. It is assumed that this quantum plant is physically realizable and corresponds to a plant Hamiltonian H p = 1 2 x p() T R p x p () where the symmetric matrixr p is given byr p = 1 4 ( ΘA p +A T pθ). Also, we consider a direct coupled linear quantum observer defined a symmetric matrix R o R n o n o, and matricesr c R n p n o,c o R m p n o.
9 A Direct Coupling Coherent Quantum Observer 9 These matrices define an observer Hamiltonian and a coupling Hamiltonian H o = 1 2 x o() T R o x o (), H c = 1 2 x p() T R c x o ()+ 1 2 x o() T R T c x p (). The matrixc o also defines the vector of estimated variables for the observer as z o (t) = C o x o (t).
10 A Direct Coupling Coherent Quantum Observer 1 The augmented quantum linear system consisting of the quantum plant and the direct coupled quantum observer is then a quantum system of the form described by the total Hamiltonian H a = H p + H c + H o = 1 2 x a() T R a x a () wherex a = [ xp x o ] andr a = [ Rp R c Rc T R o ]. Then, it follows that the augmented quantum linear system is described by the equations [ẋp (t) ẋ o (t) wherea a = 2ΘR a. ] = A a [ xp (t) x o (t) z p (t) = C p x p (t); z o (t) = C o x o (t) ] ; x p () = x p ; x o () = x o ;
11 A Direct Coupling Coherent Quantum Observer 11 Definition 1 The matrices R o R n o n o, R c R n p n o, C o R m p n o define a direct coupled linear quantum observer for the quantum linear plant if the corresponding augmented linear quantum system is such that lim T 1 T T (z p (t) z o (t))dt =.
12 A Direct Coupling Coherent Quantum Observer 12 We now describe the construction of a direct coupled linear quantum observer. We assume that the quantum plant is such thata p =. This corresponds tor p = in the plant Hamiltonian. It follows that the plant system variablesx p (t) will remain fixed if the plant is not coupled to the observer. However, when the plant is coupled to the quantum observer this will no longer be the case. We will show that if the quantum observer is suitably designed, the plant quantity to be estimatedz p (t) will remain fixed and the averaged convergence condition will be satisfied.
13 A Direct Coupling Coherent Quantum Observer 13 We also assume thatm p = n p β = andβ i R 2 1 fori = 1,2,..., n p 2. 2 and the matrixc p is of the formc p = β T where β 1 β 2... np Rn p 2 β n p 2 This assumption means that the plant variables to be estimated include only one quadrature for each mode of the plant. We now suppose that the matricesr o,r c,c o are such thatr c = βα T, α R n o n p 2 and the matrixr o is positive definite. [ ] Θ1 Also, we writeθ = whereθ Θ 1 R n p n p andθ 2 R n o n o. 2 [ ] [ ] βα T 2Θ1 βα Then,R a = αβ T anda R a = 2ΘR a = T o 2Θ 2 αβ T 2Θ 2 R o.
14 A Direct Coupling Coherent Quantum Observer 14 Hence, the augmented system equations describing the combined plant-observer system become ẋ p (t) = 2Θ 1 βα T x o (t); ẋ o (t) = 2Θ 2 αβ T x p (t)+2θ 2 R o x o (t); z p (t) = C p x p (t); z o (t) = C o x o (t). We now use Laplace Transforms to solve these equations. It follows that sx p (s) = 2Θ 1 βα T X o (s)+x p (); sx o (s) = 2Θ 2 αβ T X p (s)+2θ 2 R o X o (s)+x o () and hence, sx o (s) = 4 s Θ 2αβ T Θ 1 βα T X o (s)+ 2 s Θ 2αβ T x p () +2Θ 2 R o X o (s)+x o ().
15 A Direct Coupling Coherent Quantum Observer 15 However, β T Θ 1 β = β T 1 Jβ 1 β T 2 Jβ 2... = β T n p Jβ np sincej is skew-symmetric. Therefore, ( ) X o (s) = (si 2Θ 2 R o ) 1 2 s Θ 2αβ T x p ()+x o ().
16 A Direct Coupling Coherent Quantum Observer 16 Taking the inverse Laplace Transform of this equation, we obtain Also, we obtain x o (t) = e 2Θ 2R o t x o ()+2 = e 2Θ 2R o t x o () t o e 2Θ 2R o (t τ) dτθ 2 αβ T x p () e 2Θ 2R o t ( e 2Θ 2R o t I ) Ro 1 Θ 1 2 Θ 2αβ T x p () = e 2Θ 2R o t ( x o ()+Ro 1 αβ T x p () ) R 1 o αβ T x p (). X p (s) = 4 s 2Θ 1βα T (si 2Θ 2 R o ) 1 Θ 2 αβ T x p () + 2 s Θ 1βα T (si 2Θ 2 R o ) 1 x o () + 1 s x p().
17 A Direct Coupling Coherent Quantum Observer 17 Taking the inverse Laplace Transform of this equation, we obtain t x p (t) = 4Θ 1 βα T e 2Θ 2R o (t τ) τdτθ 2 αβ T x p () o t +2Θ 1 βα T e 2Θ 2R o (t τ) dτx o () +x p () o = 2tΘ 1 βα T R 1 o αβ T x p () +Θ 1 βα T R 2 o Θ 2 αβ T x p () Θ 1 βα T e 2Θ 2R o t R 2 o Θ 2 αβ T x p () +Θ 1 βα T R 1 o Θ 2 x o () Θ 1 βα T e 2Θ 2R o t R 1 o Θ 2 x o () +x p (). We now choose the parameters of the quantum observer so thatc o R 1 o α = I.
18 A Direct Coupling Coherent Quantum Observer 18 It follows that the quantitiesz p (t) = C p x p (t) andz o (t) = C o x o (t) are given by z o (t) = C o e 2Θ 2R o t ( x o ()+R 1 o αβ T x p () ) +z p () and z p (t) = z p () where we have used the fact thatc p Θ 1 β = β T Θ 1 β =. That is, the quantityz p (t) remains constant and is not affected by the coupling to the coherent quantum observer.
19 A Direct Coupling Coherent Quantum Observer 19 Note that the matrixa a will have all purely imaginary eigenvalues. We now verify that the averaged convergence condition is satisfied for this quantum observer. We recall that the quantity 1 2 xt R o x remains constant in time for the linear system: ẋ = 2Θ 2 R o x; x() = x. That is 1 2 x(t)t R o x(t) = 1 2 xt R o x t.
20 A Direct Coupling Coherent Quantum Observer 2 However,x(t) = e 2Θ 2R o t x andr o >. Therefore, it follows that e 2Θ 2R o t x λ max (R o ) λ min (R o ) x for allx andt. Hence, for allt. e 2Θ 2R o t λ max (R o ) λ min (R o )
21 A Direct Coupling Coherent Quantum Observer 21 Now sinceθ 2 andr o are non-singular, T and therefore, it follows that ast. e 2Θ 2R o t dt = 1 2 e2θ 2R o T Ro 1 Θ R 1 o Θ T T e 2Θ 2R o t dt = 1 T 1 2 e2θ 2R o T Ro 1 Θ R 1 o Θ T e2θ 2R o T Ro 1 Θ T 2T R 1 o Θ λ max (R o ) λ min (R o ) R 1 o Θ 1 2 2T R 1 o Θ 1 2
22 A Direct Coupling Coherent Quantum Observer 22 Hence, lim T 1 T T z o (t)dt = z p (). Also, lim T 1 T T z p (t)dt = z p (). Therefore, averaged convergence condition is satisfied. Thus, we have established the following theorem.
23 A Direct Coupling Coherent Quantum Observer 23 Theorem 2 Consider a quantum plant as above wherea p =,C p = β T and β 1 β 2 β =... np Rn p 2. β n p 2 Then the matrices R o >, R c, C o will define direct coupled quantum observer for this quantum plant if R c is of the form R c = Cp T α T where α R n o n p 2 and Co T Ro 1 α = I.
24 A Direct Coupling Coherent Quantum Observer 24 We consider the above result for the single mode case withn p = 2,m p = 1, in which C p = [1]. This means that the variable to be estimated by the quantum observer is the position [ ] qp (t) operator of the quantum plant; i.e.,z p (t) = q p (t) wherex p (t) = p p (t) [ ] [ ] 1 1 By choosingn o = 2,R o = I,C o = [1],β = andα =., the conditions of the above theorem will be satisfied and the observer output variable will be the position operator of the quantum observerq o (t); i.e.,z o (t) = q o (t) where [ ] qo (t) x o (t) =. p o (t) Before the quantum observer is connected to the quantum plant, the quantitiesq p (t) andp p (t) will remain constant since we have assumed thata p =.
25 A Direct Coupling Coherent Quantum Observer 25 Now suppose that the quantum observer is connected to the quantum plant at time t =. Then, the plant position operatorq p (t) will remain constant at its initial value q p (t) = q p () but the plant momentum operatorp p (t) will evolve in an time varying and oscillatory way. In addition, the observer position operatorq o (t) will evolve in an oscillatory way but its time average will converge toq p (). This behavior of the quantum observer is similar to the behavior of quantum measurements. This is not surprising since the behavior of the direct coupled quantum observers considered in this paper and the behavior of quantum measurements are both determined by the quantum commutation relations which are fundamental to the theory of quantum mechanics.
26 A Direct Coupling Coherent Quantum Observer 26 We now present some numerical simulations to illustrate the direct coupled quantum observer. We consider the single mode quantum observer considered above wheren p = 2, m p = 1,n o = 2,A p =,C p = [1],R o = I,C o = [1],β = [ ] 1 α =. As described above, the variable to be estimated by the quantum observer is the position operator of the quantum plant; i.e.,z p (t) = q p (t) wherex p (t) = [ 1 ] and [ qp (t) p p (t) Also, the observer output variable will be the position operator of the quantum observer q o (t); i.e.,z o (t) = q o (t) wherex o (t) = [ qo (t) p o (t) ]. ].
27 A Direct Coupling Coherent Quantum Observer 27 Then the augmented plant-observer system is described by the equations q p (t) ṗ p (t) q o (t) ṗ o (t) = A a q p (t) p p (t) q o (t) p o (t) where A a = [ 2Jβα T 2Jαβ T 2JR o ] =
28 A Direct Coupling Coherent Quantum Observer 28 Then, we can write where Φ(t) = q p (t) p p (t) q o (t) p o (t) = Φ(t) q p () p p () q o () p o () φ 11 (t) φ 12 (t) φ 13 (t) φ 14 (t) φ 21 (t) φ 22 (t) φ 23 (t) φ 24 (t) φ 31 (t) φ 32 (t) φ 33 (t) φ 34 (t) φ 41 (t) φ 42 (t) φ 43 (t) φ 44 (t) Thus, the plant variable to be estimatedq p (t) is given by = ea at. q p (t) = φ 11 (t)q p ()+φ 12 (t)p p ()+φ 13 (t)q o ()+φ 14 (t)p o () and we plot the functionsφ 11,φ 12 (t),φ 13 (t),φ 14 (t) below.
29 A Direct Coupling Coherent Quantum Observer φ 11 (t) φ 12 (t) φ 13 (t).6 φ 14 (t) φ 1i Time From this figure, we can see thatφ 11 (t) 1,φ 12 (t),φ 13 (t),φ 14 (t), andq p (t) will remain constant atq p () for allt.
30 A Direct Coupling Coherent Quantum Observer 3 Also, the other plant variablep p (t) is given by p p (t) = φ 21 (t)q p ()+φ 22 (t)p p ()+φ 23 (t)q o ()+φ 24 (t)p o () and we plot the functionsφ 21,φ 22 (t),φ 23 (t),φ 24 (t) below φ 21 (t) φ 22 (t) φ 2i 2 φ 23 (t) φ 24 (t) Time From this figure, we can see thatp p (t) evolves in a time-varying and oscillatory way when the quantum plant is connected to the quantum observer.
31 A Direct Coupling Coherent Quantum Observer 31 We now consider the output variable of the quantum observerq o (t) which is given by q o (t) = φ 31 (t)q p ()+φ 32 (t)p p ()+φ 33 (t)q o ()+φ 34 (t)p o () and we plot the functionsφ 31,φ 32 (t),φ 33 (t),φ 34 (t) below. 4 3 φ 31 (t) φ 32 (t) φ 33 (t) φ 34 (t) 2 φ 3i Time
32 A Direct Coupling Coherent Quantum Observer 32 To illustrate the time average convergence property of the quantum observer, we now plot the quantities φ ave 31 (T) = 1 T φ ave 32 (T) = 1 T φ ave 33 (T) = 1 T φ ave 34 (T) = 1 T T T T T φ 31 (t)dt φ 32 (t)dt φ 33 (t)dt φ 34 (t)dt below.
33 A Direct Coupling Coherent Quantum Observer av φ 31 (t) av φ (t) 32 av φ (t) 33 av φ 34 (t) 1.5 av φ 3i Time From this figure, we can see that the time average ofq o (t) converges toq p () as t.
34 A Direct Coupling Coherent Quantum Observer 34 To investigate the time average property of the other quantum observer variable, we now plot the quantities φ ave 41 (T) = 1 T φ ave 42 (T) = 1 T φ ave 43 (T) = 1 T φ ave 44 (T) = 1 T T T T T φ 41 (t)dt φ 42 (t)dt φ 43 (t)dt φ 44 (t)dt below.
35 A Direct Coupling Coherent Quantum Observer av φ 41 (t) av φ (t) 42 av φ 43 (t).5 av φ 44 (t) av φ 4i Time
36 A Direct Coupling Coherent Quantum Observer 36 In this paper we have introduced a notion of a direct coupling observer for closed quantum linear systems and given a result which shows how such an observer can be constructed. The main result shows the time average convergence properties of the direct coupling observer. We have also presented an illustrative example along with simulations to investigate the behavior of a direct coupling observer when applied to a simple one mode quantum linear system. Future research in this area might involve extending the class of quantum linear systems for a which a direct coupling observer can be designed and also considering the problem of constructing an observer which is optimal in some sense. Also, future research could investigate the role of direct coupling observers in the design of direct coupling coherent quantum feedback control systems.
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 informationQuantum 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 informationControl 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 informationRobust Control of Linear Quantum Systems
B. Ross Barmish Workshop 2009 1 Robust Control of Linear Quantum Systems Ian R. Petersen School of ITEE, University of New South Wales @ the Australian Defence Force Academy, Based on joint work with Matthew
More informationAdvanced Mechatronics Engineering
Advanced Mechatronics Engineering German University in Cairo 21 December, 2013 Outline Necessary conditions for optimal input Example Linear regulator problem Example Necessary conditions for optimal input
More informationRobust 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 informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationQuantum Theory and Group Representations
Quantum Theory and Group Representations Peter Woit Columbia University LaGuardia Community College, November 1, 2017 Queensborough Community College, November 15, 2017 Peter Woit (Columbia University)
More information1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is
ECE 55, Fall 2007 Problem Set #4 Solution The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ Ax + Bu is x(t) e A(t ) x( ) + e A(t τ) Bu(τ)dτ () This formula is extremely important
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More information1. 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 informationECEEN 5448 Fall 2011 Homework #5 Solutions
ECEEN 5448 Fall 211 Homework #5 Solutions Professor David G. Meyer December 8, 211 1. Consider the 1-dimensional time-varying linear system ẋ t (u x) (a) Find the state-transition matrix, Φ(t, τ). Here
More informationSimple Harmonic Oscillator
Classical harmonic oscillator Linear force acting on a particle (Hooke s law): F =!kx From Newton s law: F = ma = m d x dt =!kx " d x dt + # x = 0, # = k / m Position and momentum solutions oscillate in
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationOn complexified quantum mechanics and space-time
On complexified quantum mechanics and space-time Dorje C. Brody Mathematical Sciences Brunel University, Uxbridge UB8 3PH dorje.brody@brunel.ac.uk Quantum Physics with Non-Hermitian Operators Dresden:
More informationAdvanced Control Theory
State Space Solution and Realization chibum@seoultech.ac.kr Outline State space solution 2 Solution of state-space equations x t = Ax t + Bu t First, recall results for scalar equation: x t = a x t + b
More information6.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 informationStatistical Mechanics Solution Set #1 Instructor: Rigoberto Hernandez MoSE 2100L, , (Dated: September 4, 2014)
CHEM 6481 TT 9:3-1:55 AM Fall 214 Statistical Mechanics Solution Set #1 Instructor: Rigoberto Hernandez MoSE 21L, 894-594, hernandez@gatech.edu (Dated: September 4, 214 1. Answered according to individual
More informationLecture 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 informationMatrix Solutions to Linear Systems of ODEs
Matrix Solutions to Linear Systems of ODEs James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 3, 216 Outline 1 Symmetric Systems of
More informationState Variable Analysis of Linear Dynamical Systems
Chapter 6 State Variable Analysis of Linear Dynamical Systems 6 Preliminaries In state variable approach, a system is represented completely by a set of differential equations that govern the evolution
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [
More information8 Periodic Linear Di erential Equations - Floquet Theory
8 Periodic Linear Di erential Equations - Floquet Theory The general theory of time varying linear di erential equations _x(t) = A(t)x(t) is still amazingly incomplete. Only for certain classes of functions
More informationLecture 19 (Nov. 15, 2017)
Lecture 19 8.31 Quantum Theory I, Fall 017 8 Lecture 19 Nov. 15, 017) 19.1 Rotations Recall that rotations are transformations of the form x i R ij x j using Einstein summation notation), where R is an
More informationCreation and Destruction Operators and Coherent States
Creation and Destruction Operators and Coherent States WKB Method for Ground State Wave Function state harmonic oscillator wave function, We first rewrite the ground < x 0 >= ( π h )1/4 exp( x2 a 2 h )
More informationNumerical Linear Algebra Homework Assignment - Week 2
Numerical Linear Algebra Homework Assignment - Week 2 Đoàn Trần Nguyên Tùng Student ID: 1411352 8th October 2016 Exercise 2.1: Show that if a matrix A is both triangular and unitary, then it is diagonal.
More informationAn Exactly Solvable 3 Body Problem
An Exactly Solvable 3 Body Problem The most famous n-body problem is one where particles interact by an inverse square-law force. However, there is a class of exactly solvable n-body problems in which
More informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
More informationSemi-Classical Theory of Radiative Transitions
Semi-Classical Theory of Radiative Transitions Massimo Ricotti ricotti@astro.umd.edu University of Maryland Semi-Classical Theory of Radiative Transitions p.1/13 Atomic Structure (recap) Time-dependent
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationModal Decomposition and the Time-Domain Response of Linear Systems 1
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING.151 Advanced System Dynamics and Control Modal Decomposition and the Time-Domain Response of Linear Systems 1 In a previous handout
More informationChapter 2 The Group U(1) and its Representations
Chapter 2 The Group U(1) and its Representations The simplest example of a Lie group is the group of rotations of the plane, with elements parametrized by a single number, the angle of rotation θ. It is
More informationENSC327 Communications Systems 19: Random Processes. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 19: Random Processes Jie Liang School of Engineering Science Simon Fraser University 1 Outline Random processes Stationary random processes Autocorrelation of random processes
More informationTopic # Feedback Control Systems
Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the
More informationECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77
1/77 ECEN 605 LINEAR SYSTEMS Lecture 7 Solution of State Equations Solution of State Space Equations Recall from the previous Lecture note, for a system: ẋ(t) = A x(t) + B u(t) y(t) = C x(t) + D u(t),
More informationSample Quantum Chemistry Exam 2 Solutions
Chemistry 46 Fall 7 Dr. Jean M. Standard Name SAMPE EXAM Sample Quantum Chemistry Exam Solutions.) ( points) Answer the following questions by selecting the correct answer from the choices provided. a.)
More informationNIU PHYS 500, Fall 2006 Classical Mechanics Solutions for HW6. Solutions
NIU PHYS 500, Fall 006 Classical Mechanics Solutions for HW6 Assignment: HW6 [40 points] Assigned: 006/11/10 Due: 006/11/17 Solutions P6.1 [4 + 3 + 3 = 10 points] Consider a particle of mass m moving in
More informationSketchy Notes on Lagrangian and Hamiltonian Mechanics
Sketchy Notes on Lagrangian and Hamiltonian Mechanics Robert Jones Generalized Coordinates Suppose we have some physical system, like a free particle, a pendulum suspended from another pendulum, or a field
More informationSemi-Classical Dynamics Using Hagedorn Wavepackets
Semi-Classical Dynamics Using Hagedorn Wavepackets Leila Taghizadeh June 6, 2013 Leila Taghizadeh () Semi-Classical Dynamics Using Hagedorn Wavepackets June 6, 2013 1 / 56 Outline 1 The Schrödinger Equation
More informationME 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 informationFirst-Order Low-Pass Filter
Filters, Cost Functions, and Controller Structures Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 218! Dynamic systems as low-pass filters! Frequency response of dynamic systems!
More informationOPTIMAL CONTROL SYSTEMS
SYSTEMS MIN-MAX Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University OUTLINE MIN-MAX CONTROL Problem Definition HJB Equation Example GAME THEORY Differential Games Isaacs
More information6. Linear Quadratic Regulator Control
EE635 - Control System Theory 6. Linear Quadratic Regulator Control Jitkomut Songsiri algebraic Riccati Equation (ARE) infinite-time LQR (continuous) Hamiltonian matrix gain margin of LQR 6-1 Algebraic
More informationGeneralized Riccati Equations Arising in Stochastic Games
Generalized Riccati Equations Arising in Stochastic Games Michael McAsey a a Department of Mathematics Bradley University Peoria IL 61625 USA mcasey@bradley.edu Libin Mou b b Department of Mathematics
More informationPostulates of Quantum Mechanics
EXERCISES OF QUANTUM MECHANICS LECTURE Departamento de Física Teórica y del Cosmos 018/019 Exercise 1: Stern-Gerlach experiment Postulates of Quantum Mechanics AStern-Gerlach(SG)deviceisabletoseparateparticlesaccordingtotheirspinalonga
More informationComplicated behavior of dynamical systems. Mathematical methods and computer experiments.
Complicated behavior of dynamical systems. Mathematical methods and computer experiments. Kuznetsov N.V. 1, Leonov G.A. 1, and Seledzhi S.M. 1 St.Petersburg State University Universitetsky pr. 28 198504
More informationClassifying and building DLMs
Chapter 3 Classifying and building DLMs In this chapter we will consider a special form of DLM for which the transition matrix F k and the observation matrix H k are constant in time, ie, F k = F and H
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 12: I/O Stability Readings: DDV, Chapters 15, 16 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology March 14, 2011 E. Frazzoli
More informationECE 388 Automatic Control
Controllability and State Feedback Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage:
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationOptimal Control. Quadratic Functions. Single variable quadratic function: Multi-variable quadratic function:
Optimal Control Control design based on pole-placement has non unique solutions Best locations for eigenvalues are sometimes difficult to determine Linear Quadratic LQ) Optimal control minimizes a quadratic
More informationSecond quantization: where quantization and particles come from?
110 Phys460.nb 7 Second quantization: where quantization and particles come from? 7.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system? 7.1.1.Lagrangian Lagrangian
More informationThe 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 informationControl, 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 informationQuantum Mechanics C (130C) Winter 2014 Assignment 7
University of California at San Diego Department of Physics Prof. John McGreevy Quantum Mechanics C (130C) Winter 014 Assignment 7 Posted March 3, 014 Due 11am Thursday March 13, 014 This is the last problem
More informationMATH 1140(M12A) Semester I. Ax =0 (1) x 4. x 2. x 3. By using Gaussian Elimination, we obtain the solution
MATH 4(M2A) Semester I Homogeneous Systems These are systems of linear equations of the form Ax = () whereais the coefficient matrix and x is the vector of unknowns. Example x 2x 2 4x 3 3x 4 = 2x +x 2
More informationSYSTEMTEORI - KALMAN FILTER VS LQ CONTROL
SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL 1. Optimal regulator with noisy measurement Consider the following system: ẋ = Ax + Bu + w, x(0) = x 0 where w(t) is white noise with Ew(t) = 0, and x 0 is a stochastic
More informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More informationGENERAL ARTICLE Realm of Matrices
Realm of Matrices Exponential and Logarithm Functions Debapriya Biswas Debapriya Biswas is an Assistant Professor at the Department of Mathematics, IIT- Kharagpur, West Bengal, India. Her areas of interest
More informationNon-Differentiable Embedding of Lagrangian structures
Non-Differentiable Embedding of Lagrangian structures Isabelle Greff Joint work with J. Cresson Université de Pau et des Pays de l Adour CNAM, Paris, April, 22nd 2010 Position of the problem 1. Example
More informationarzelier
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 informationDifferential equations
Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system
More informationKAM for NLS with harmonic potential
Université de Nantes 3rd Meeting of the GDR Quantum Dynamics MAPMO, Orléans, 2-4 February 2011. (Joint work with Benoît Grébert) Introduction The equation : We consider the nonlinear Schrödinger equation
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationChapter 29. Quantum Chaos
Chapter 29 Quantum Chaos What happens to a Hamiltonian system that for classical mechanics is chaotic when we include a nonzero h? There is no problem in principle to answering this question: given a classical
More informationEE221A Linear System Theory Final Exam
EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,
More informationQuantization of scalar fields
Quantization of scalar fields March 8, 06 We have introduced several distinct types of fields, with actions that give their field equations. These include scalar fields, S α ϕ α ϕ m ϕ d 4 x and complex
More informationMin-Max Output Integral Sliding Mode Control for Multiplant Linear Uncertain Systems
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July -3, 27 FrC.4 Min-Max Output Integral Sliding Mode Control for Multiplant Linear Uncertain
More informationHamiltonian Systems of Negative Curvature are Hyperbolic
Hamiltonian Systems of Negative Curvature are Hyperbolic A. A. Agrachev N. N. Chtcherbakova Abstract The curvature and the reduced curvature are basic differential invariants of the pair: Hamiltonian system,
More informationNewtonian Mechanics. Chapter Classical space-time
Chapter 1 Newtonian Mechanics In these notes classical mechanics will be viewed as a mathematical model for the description of physical systems consisting of a certain (generally finite) number of particles
More informationQUALIFYING EXAMINATION, Part 1. Solutions. Problem 1: Mathematical Methods. x 2. 1 (1 + 2x 2 /3!) ] x 2 1 2x 2 /3
QUALIFYING EXAMINATION, Part 1 Solutions Problem 1: Mathematical Methods (a) Keeping only the lowest power of x needed, we find 1 x 1 sin x = 1 x 1 (x x 3 /6...) = 1 ) 1 (1 x 1 x /3 = 1 [ 1 (1 + x /3!)
More informationController Coefficient Truncation Using Lyapunov Performance Certificate
Controller Coefficient Truncation Using Lyapunov Performance Certificate Joëlle Skaf Stephen Boyd Information Systems Laboratory Electrical Engineering Department Stanford University European Control Conference,
More informationSmall oscillations and normal modes
Chapter 4 Small oscillations and normal modes 4.1 Linear oscillations Discuss a generalization of the harmonic oscillator problem: oscillations of a system of several degrees of freedom near the position
More information( ) in the interaction picture arises only
Physics 606, Quantum Mechanics, Final Exam NAME 1 Atomic transitions due to time-dependent electric field Consider a hydrogen atom which is in its ground state for t < 0 For t > 0 it is subjected to a
More informationDifferentiable Functions
Differentiable Functions Let S R n be open and let f : R n R. We recall that, for x o = (x o 1, x o,, x o n S the partial derivative of f at the point x o with respect to the component x j is defined as
More informationAssignment 8. [η j, η k ] = J jk
Assignment 8 Goldstein 9.8 Prove directly that the transformation is canonical and find a generating function. Q 1 = q 1, P 1 = p 1 p Q = p, P = q 1 q We can establish that the transformation is canonical
More informationECEEN 5448 Fall 2011 Homework #4 Solutions
ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with
More informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
More informationModule 02 CPS Background: Linear Systems Preliminaries
Module 02 CPS Background: Linear Systems Preliminaries Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html August
More informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation
High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y
More informationTypicality paradigm in Quantum Statistical Thermodynamics Barbara Fresch, Giorgio Moro Dipartimento Scienze Chimiche Università di Padova
Typicality paradigm in Quantum Statistical Thermodynamics Barbara Fresch, Giorgio Moro Dipartimento Scienze Chimiche Università di Padova Outline 1) The framework: microcanonical statistics versus the
More informationControl Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017
More informationControls 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 informationIntroduction to Integrability
Introduction to Integrability Problem Sets ETH Zurich, HS16 Prof. N. Beisert, A. Garus c 2016 Niklas Beisert, ETH Zurich This document as well as its parts is protected by copyright. This work is licensed
More informationDirac Equation. Chapter 1
Chapter Dirac Equation This course will be devoted principally to an exposition of the dynamics of Abelian and non-abelian gauge theories. These form the basis of the Standard Model, that is, the theory
More informationLecture 3: Review of Linear Algebra
ECE 83 Fall 2 Statistical Signal Processing instructor: R Nowak Lecture 3: Review of Linear Algebra Very often in this course we will represent signals as vectors and operators (eg, filters, transforms,
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science Quantum Optical Communication
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.453 Quantum Optical Communication Date: Tuesday, September 13, 216 Dirac-notation Quantum Mechanics. Introduction
More informationThe 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 informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More informationLecture 3: Review of Linear Algebra
ECE 83 Fall 2 Statistical Signal Processing instructor: R Nowak, scribe: R Nowak Lecture 3: Review of Linear Algebra Very often in this course we will represent signals as vectors and operators (eg, filters,
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationQuantum Dynamics. March 10, 2017
Quantum Dynamics March 0, 07 As in classical mechanics, time is a parameter in quantum mechanics. It is distinct from space in the sense that, while we have Hermitian operators, X, for position and therefore
More informationforms Christopher Engström November 14, 2014 MAA704: Matrix factorization and canonical forms Matrix properties Matrix factorization Canonical forms
Christopher Engström November 14, 2014 Hermitian LU QR echelon Contents of todays lecture Some interesting / useful / important of matrices Hermitian LU QR echelon Rewriting a as a product of several matrices.
More informationON COST MATRICES WITH TWO AND THREE DISTINCT VALUES OF HAMILTONIAN PATHS AND CYCLES
ON COST MATRICES WITH TWO AND THREE DISTINCT VALUES OF HAMILTONIAN PATHS AND CYCLES SANTOSH N. KABADI AND ABRAHAM P. PUNNEN Abstract. Polynomially testable characterization of cost matrices associated
More informationPhysics 581, Quantum Optics II Problem Set #4 Due: Tuesday November 1, 2016
Physics 581, Quantum Optics II Problem Set #4 Due: Tuesday November 1, 2016 Problem 3: The EPR state (30 points) The Einstein-Podolsky-Rosen (EPR) paradox is based around a thought experiment of measurements
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationImplications of Time-Reversal Symmetry in Quantum Mechanics
Physics 215 Winter 2018 Implications of Time-Reversal Symmetry in Quantum Mechanics 1. The time reversal operator is antiunitary In quantum mechanics, the time reversal operator Θ acting on a state produces
More informationUnitary Dynamics and Quantum Circuits
qitd323 Unitary Dynamics and Quantum Circuits Robert B. Griffiths Version of 20 January 2014 Contents 1 Unitary Dynamics 1 1.1 Time development operator T.................................... 1 1.2 Particular
More information