LMIs in Control/pages/Multi-Criterion LQG: Difference between revisions

From testwiki
Jump to navigation Jump to search
imported>Rezajamesahmed
 
(No difference)

Latest revision as of 20:45, 12 December 2019

LMIs in Control/pages/Multi-Criterion LQG

The Multi-Criterion Linear Quadratic Gaussian (LQG) linear matrix inequality will allow one to form an optimized controller, similar to that in an LQR framework, for a state space system with gaussian noise based on several different criterions defined in the Q and R matrices, that are optimized as a part of the arbitrary cost function. Just like traditional LQR, the cost matrices must be tuned in much a similar fashion as traditional gains in classical control. In the LQR and LQG framework however, the gains are more intuitive as each correlates directly to a state or an input.

The System

The system is a linear time-invariant system, that can be represented in state space as shown below:

x˙=Ax+Bu+w,y=Cx+v,z=[Q1/200R1/2][xu]

where xRn,yRl,zRm represent the state vector, the measured output vector, and the output vector of interest, respectively, wRp is the disturbance vector, and A,B,C,Q,R are the system matrices of appropriate dimension. To further define: x is Rn and is the state vector, A is Rn*n and is the state matrix, B is Rn*r and is the input matrix, w is Rr and is the exogenous input, C,Q,R is Rm*n and are the output matrices, and y and z are Rm and are the output and the output of interest, respectively.


Q0 and R>0, and the system is controllable and observable.

The Data

The matrices A,B,C,Q,R,W,V and the noise signals w,v.

The Optimization Problem

In the Linear Quadratic Gaussian (LQG) control problem, the goal is to minimize a quadratic cost function while the plant has random initial conditions and suffers white noise disturbance on the input and measurement.

There are multiple outputs of interest for this problem. They are defined by

z=[Q1/200R1/2][xu],Qi0,Ri>0,i=0,...,p.

For each of these outputs of interest, we associate a cost function:

JLQGi=limtEzi(t)Tzi(t),i=0,...,p.

Additionally, the matrices XLQG and YLQG must be found as the solutions to the following Riccati equations:

ATXLQG+XLQGA=XLQGBR1BTXLQG+Q=0AYLQG+YLQGATYLQGCTV1CYLQG+W=0

The optimization problem is to minimize JLQG0 over u subject to the measurability condition and the constraints JLQGi<γi,i=0,...,p.. This optimization problem can be formulated as:

maxtrace(XLQGU+QYLQG)i=1pγiτi,

over τ1,...,τp, with:

Q=Q0+i=1pτiQi,R=R0+i=1pτiRi.

The LMI: Multi-Criterion LQG

max:trace(XU+(Q0+i=1pτiQi)YLQG)i=1pγiτi,

over X,τ1,...τp, subject to the following constraints:

X>0,τ10,...,τp0,ATX+XAXB(R0+i=1pτiRi)1BTX+Q0+i=1pτiQi0.

Conclusion:

The result of this LMI is the solution to the aforementioned Ricatti equations:

ATXLQG+XLQGA=XLQGBR1BTXLQG+Q=0AYLQG+YLQGATYLQGCTV1CYLQG+W=0

Implementation

This implementation requires Yalmip and Sedumi.

https://github.com/rezajamesahmed/LMImatlabcode/blob/master/multicriterionquadraticproblems.m

  1. Inverse Problem of Optimal Control

Return to Main Page:

Template:BookCat