Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Formulating Linear Constraints on a Matrix

  1. Jun 1, 2010 #1
    Hi everyone. I have this problem which I am trying to formulate. Basically, I have the following linear constraints:

    p_{11} = 2
    p_{22} = 5
    2p_{13} =2
    2p_{23} = 0

    And these are for the symmetric matrix

    \mathbf{P} =
    \left( \begin{array}{ccc}
    p_{11} & p_{12} & p_{13} \\
    p_{12} & p_{22} & p_{23} \\
    p_{13} & p_{23} & p_{33} \end{array} \right)

    I would like to formulate a way to represent the linear constraints and [tex]\mathbf{P}[/tex] as a matrix at the same time.

    I can do this using [tex]\mathbf{P}[/tex] or a vector of the entries of [tex]\mathbf{P}[/tex]. The linear constraints are easy if I use a vector ([tex]\mathbf{Ap}=\mathbf{b}[/tex], but then I don't know how to represent [tex]\mathbf{P}[/tex] as a matrix from the vector! And if I leave [tex]\mathbf{P}[/tex] as a matrix, all the constraints are easy to formulate except [tex]p_{33}+2p_{12}=-1[/tex]. Can anyone help me figure this out?

    If anyone's curious, I'm trying to solve for [tex]\mathbf{P}[/tex] over the cone of PSD matrices using SDP. But I am entirely new to SDP and I'm scratching my head formulating this problem. I feel stupid right now :'(
  2. jcsd
  3. Jun 1, 2010 #2
    I answered my own question... and the answer has to do with algebraic dependence in the structure of the problem, and how SDP problems are formulated. I just thought I was missing something obvious ... hence my convoluted search.
  4. Jun 1, 2010 #3


    User Avatar
    Science Advisor

    Well, I'm glad you answered it. It looks simple to me: I you let [itex]p_{12}= p[/itex] then [itex]p_{33}= -1- 2p[/itex] and all other values are essentially given:
    [tex]\begin{pmatrix}2 & p & 1 \\ p & 5 & 0 \\1 & 0 & -1-2p\end{pmatrix}[/tex]
  5. Jun 1, 2010 #4
    That's exactly what came out in the end. I'm not sure why I didn't see it at first... I always seem to manage to do things in the most roundabout way.
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook