
where
, are independent, identically distributed random variables with
and
is independent of the process
.
The matrices
.

https://github.com/mkhajenejad/Mohammad-Khajenejad/commit/a34713575cd8ae9831cb5b7eb759d0fd45a8c37f
The optimal
returns an upper bound on the
gain of the system. .
It is straightforward to apply scaling method [Boyd, sec.6.3.4] to obtain component-wise results.