Form linear-quadratic-Gaussian (LQG) regulator

`rlqg = lqgreg(kest,k)`

rlqg = lqgreg(kest,k,controls)

`rlqg = lqgreg(kest,k)`

returns the LQG regulator
`rlqg`

(a state-space model) given the Kalman estimator
`kest`

and the state-feedback gain matrix `k`

. The same
function handles both continuous- and discrete-time cases. Use consistent tools to design
`kest`

and `k`

:

Continuous regulator for continuous plant: use

`lqr`

or`lqry`

and`kalman`

Discrete regulator for discrete plant: use

`dlqr`

or`lqry`

and`kalman`

Discrete regulator for continuous plant: use

`lqrd`

and`kalmd`

In discrete time, `lqgreg`

produces the regulator

$$u\left[n\right]=-K\widehat{x}\left[n|n\right]$$ when

`kest`

is the "current" Kalman estimator$$u\left[n\right]=-K\widehat{x}\left[n|n-1\right]$$ when

`kest`

is the "delayed" Kalman estimator

For more information on Kalman estimators, see the `kalman`

reference page.

`rlqg = lqgreg(kest,k,controls)`

handles estimators that have access to
additional deterministic known plant inputs *u _{d}*. The
index vector

`controls`

then specifies which estimator inputs are the
controls `rlqg`

has
Always use *positive* feedback to connect the LQG regulator to the
plant.

See the example LQG Regulation: Rolling Mill Case Study.

`lqgreg `

forms the linear-quadratic-Gaussian (LQG)
regulator by connecting the Kalman estimator designed with `kalman`

and the
optimal state-feedback gain designed with `lqr`

, `dlqr`

, or
`lqry`

. The LQG regulator minimizes some quadratic cost function that
trades off regulation performance and control effort. This regulator is dynamic and relies on
noisy output measurements to generate the regulating commands.

In continuous time, the LQG regulator generates the commands

$$u=-K\widehat{x}$$

where $$\widehat{x}$$ is the Kalman state estimate. The regulator state-space equations are

$$\begin{array}{l}\dot{\widehat{x}}=[A-LC-(B-LD)K]\widehat{x}+Ly\\ u=-K\widehat{x}\end{array}$$

where *y* is the vector of plant output measurements (see
`kalman`

for background and notation). The following diagram shows this
dynamic regulator in relation to the plant.

In discrete time, you can form the LQG regulator using either the delayed state estimate $$\widehat{x}[n|n-1]$$ of *x*[*n*], based on measurements up to
*y*[*n–1*], or the current state estimate $$\widehat{x}[n|n]$$, based on all available measurements including
*y*[*n*]. While the regulator

$$u\left[n\right]=-K\widehat{x}\left[n|n-1\right]$$

is always well-defined, the *current regulator*

$$u\left[n\right]=-K\widehat{x}\left[n|n\right]$$

is causal only when *I*-*KMD* is invertible (see
`kalman`

for the notation). In addition, practical implementations of the
current regulator should allow for the processing time required to compute
*u*[*n*] after the measurements
*y*[*n*] become available (this amounts to a time delay in
the feedback loop).

For a discrete-time plant with equations:

$$\begin{array}{c}x\left[n+1\right]=Ax\left[n\right]+Bu\left[n\right]+Gw\left[n\right]\\ y\left[n\right]=Cx\left[n\right]+Du\left[n\right]+Hw\left[n\right]+v\left[n\right]\text{\hspace{1em}}\left\{\text{Measurements}\right\}\end{array}$$

connecting the "current" Kalman estimator to the LQR gain is optimal only when $$E\left(w\left[n\right]v{\left[n\right]}^{\prime}\right)=0$$ and *y*[*n*] does not depend on
*w*[*n*] (*H* = 0). If these conditions are not satisfied, compute the optimal LQG controller
using `lqg`

.