2
$\begingroup$

Let $A\in\mathbb{R}^{n\times n}$ be a stable matrix (i.e., the eigenvalues of $A$ have negative real parts). Consider the following optimization problem in $X \in \mathbb{R}^{n \times n}$

$$\begin{array}{ll} \text{maximize} & \mbox{tr} \left(P X \right)\\ \text{subject to} & \mbox{tr}(X) = 1\\ & X \succeq 0\end{array}$$

where $X \succeq 0$ denotes that $X$ is positive semidefinite, and $P$ is the solution of the following Lyapunov equation $$ AP+PA^\top =-X. $$

Question: Does the above-formulated problem admit a closed-form solution?


An equivalent formulation. We can rewrite the above optimization problem in the following equivalent way. Consider the optimization problem in $P \in \mathbb{R}^{n \times n}$, $P\succeq 0$,

$$\begin{array}{ll} \text{minimize} & \mbox{tr} \left(PAP \right)\\ \text{subject to} & \mbox{tr}(PA) = -1/2\\ & AP+PA^\top \preceq 0\\ & P \succeq 0\end{array}$$

$\endgroup$
3
  • $\begingroup$ Yes, thanks. I've rewritten the problem in order to acknowledge this. $\endgroup$
    – Ludwig
    Commented Nov 22, 2017 at 10:23
  • $\begingroup$ @RodrigodeAzevedo: I added an equivalent formulation of the problem at the end of the question. $\endgroup$
    – Ludwig
    Commented Nov 22, 2017 at 12:50
  • $\begingroup$ @RodrigodeAzevedo: Absolutely right! I've corrected these typos. $\endgroup$
    – Ludwig
    Commented Nov 22, 2017 at 19:19

1 Answer 1

1
$\begingroup$

This is not really an answer, but an observation too long for a comment. The OP may already know this ...

From the solution of Lyapunov equation $P = \int_{0}^{\infty} e^{At}Xe^{A^{\top}t} {\mathrm{d}}t$, we have ${\rm{tr}}(PX) = \int_{0}^{\infty} {\rm{tr}}(YX){\mathrm{d}}t$, where $Y := e^{At}Xe^{A^{\top}t} \succeq 0$.

One can then get a bound by using the fact that for positive semi-definite $Y,X$, we have $\lambda_{\min}(Y){\rm{tr}}(X) \leq {\rm{tr}}(YX) \leq \lambda_{\max}(Y){\rm{tr}}(X)$. Using this inequality with ${\rm{tr}}(X)=1$ gives an upper bound:

$$ {\rm{tr}}(PX) \leq \int_{0}^{\infty} \parallel e^{At} X^{1/2}\parallel_{2}^{2} {\rm{d}}t \leq \lambda_{\max}(X)\int_{0}^{\infty} \parallel e^{At}\parallel_{2}^{2} {\rm{d}}t$$

where the last step used sub-multiplicative property of 2-norm. Not sure how tight is this.

$\endgroup$

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .