Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

Let $A$ be a $2n\times 2n$ real matrix with the following structure \begin{equation} A = \left(\begin{matrix} 0 & -I \\ K & S \end{matrix}\right) \end{equation} with all sub-matrices of size $n\times n$: $I$ is the identity matrix, $K$ is symmetric positive definite and $S$ is diagonal but singular. I am interested in the (numerical) solution of the continuous time Lyapunov equation for the $2n\times 2n$ matrix R: \begin{equation} R A^\text{T} + A R = \left(\begin{matrix} 0 & 0 \\ 0 & \Gamma \end{matrix}\right) \end{equation} where $\Gamma$ is a diagonal and singular $n\times n$ real matrix.

However, I only need a few elements of $R$. More specifically, writing

\begin{equation} R = \left(\begin{matrix} X & C \\ C^\text{T} & V \end{matrix}\right) \end{equation} all I really want are the diagonal entries of $V$ (also $n\times n$).

Is there anyway I could reduce this problem to some other (probably Sylvester) equation for $V$ or, better yet, only it's diagonal? I don't really know how to approach this problem.

share|improve this question
    
what is $C$?what kind of matrix –  dato datuashvili Aug 9 '12 at 12:45
    
The solution of the Lyapunov equation, $R$, will be symmetric and positive definite (it is actually a covariance matrix). What I can prove is that $C$ is anti-symmetric ($C^\text{T} = -C$) –  Gabriel Landi Aug 9 '12 at 13:00
    
put instead of $C^T$ ;$-C$ if you get identity ,then it should solve problem –  dato datuashvili Aug 9 '12 at 13:14
    
@dato, Sorry I don't understand. –  Gabriel Landi Aug 9 '12 at 13:22

1 Answer 1

up vote 2 down vote accepted

Just doing the matrix multiplications gives four equations: $$\begin{align} C^T - C &=0 \\ XK + CS &= V \\ KX - SC &= V \\ KC - CK + VS + SV &= \Gamma \end{align}$$ the last equation of course gives $VS + SV = \Gamma - KC + CK$, so if we know what $C$ is we get a Sylvester equation for $V$. If we substitute the LHS of the second and third equation in place for $V$ in the fourth equation we get: $$K(XS + C) + (SX - C)K=\Gamma$$ and assuming $X$ is symmetric, we get $(XS + C)^T = SX - C$, so setting $B = XS + C$, the equation turns into $$KB + B^TK = \Gamma$$ which we want to solve for $B$. Now, this is not really a standard Sylvester equations since the matrix we seek is transposed, but there are methods to solve these kinds of equations too.

Once $B$ has been found we get $$C = \frac{1}{2}(B-B^T)$$ Since $B$ is the sum of one antisymmetric ($C$) and one symmetric matrix ($XS$), the antisymmetric part of $B$ has to be $C$.

Thus, one can now calculate the RHS in $$VS + SV = \Gamma - KC + CK$$ and then solve for $V$ as a standard Sylvester equation.

So, using this method you just solve two $n \times n$ Sylvester-ish equations instead of one $2n \times 2n$ equation. There might be something smarter lurking in here but I have not found it yet.

share|improve this answer
    
Wow, this is great. Thank you! I'll start testing this idea right away and post some results as soon as they come out. –  Gabriel Landi Aug 14 '12 at 0:27
    
Actually, I have been looking at the web and haven't yet found a method for solving the eq. for B. Any ideas on where to look? –  Gabriel Landi Aug 14 '12 at 15:03
    
Try this article: H. W. Braden, The Equations $A^TX \pm X^TA = B$, SIAM J. Matrix Anal. Appl. Vol. 20, No. 2, pp. 295-302. –  Calle Aug 14 '12 at 15:38

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.