10
$\begingroup$

How does one solve the matrix equation $AX+XB=C$ for $X$? It doesn't seem too difficult. I tried many times but failed.

I'm an adult student... I am now vexed about Gilbert Strang - An Introduction to Linear Algebra. I don't even understand a single word in Wikipedia: Sylvester equation. If you have ever use some nice workable materials or lecture notes? You can generously upload and share the links of the lecture notes and assignments. Different subjects/ topics are welcome, as long as you deem they are nice and workable.

The problem origins from a system of diff equation, using undetermined coefficients (matrix) to find the particular solution. Try $y_p=X\begin{pmatrix} e^{\alpha t} \\ e^{\beta t} \end{pmatrix}$

$\dot{y}+Ay=C\begin{pmatrix} e^{\alpha t} \\ e^{\beta t}\end{pmatrix}$

$\dot{y_p}=X\begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix} \begin{pmatrix}e^{\alpha t} \\ e^{\beta t}\end{pmatrix}$

substitute $\dot{y_p}$ and $y_p$ into the original differential equation..

$X\begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix}+AX=C$

  • 4
    What you have there is a [Sylvester equation](http://en.wikipedia.org/wiki/Sylvester_equation) for which many solution methods are known.2012-08-25
  • 1
    Are in interested in solving the matrix equation to solve the differential equation, or in its own right?2012-08-25
  • 1
    I think this question was asked here on math.SE before, but I can't find the link.2012-08-25
  • 0
    "I don't even understand a single word in Wikipedia" - fine... do you at least know what a Kronecker product is? Or, better yet, what part of the Wiki article do you not understand?2012-08-25
  • 0
    @JenniferDylan: I think the closest thing is this, but it's not exactly a duplicate, since it has no free term: http://math.stackexchange.com/questions/39906/solving-a-matrix-equation-ax-xb-in-a-cas2012-08-25
  • 0
    It seems I'm seriously incapable. It is better to suggest me some materials or books rather than just give a solution to my problem. I’d like to contact J.M. also...2012-08-26

2 Answers 2

3

It seems $AX+XB=[A,I_1]X[I_2,B]^T$, where $I_1$ and $I_2$ are identity matrices and subscript 'T' is transpose of a matrix. Thus, the problem becomes how to solve $[A,I_1]X[I_2,B]^T = C $. Obviously, $X=[A,I_1]^\dagger C {[I_2,B]^T}^\dagger$, where subscript '$\dagger$' means the Pseudo-inverse.

  • 1
    This is incorrect. $AX+XB$ is $[A,I]\begin{bmatrix}X\\ &X\end{bmatrix}\begin{bmatrix}I\\ B\end{bmatrix}$, not $[A,I]X[I,B]^T$. Note that $[A,I]X$ is not a legitimate matrix product, as $[A,I]$ has $2n$ columns but $X$ has $n$ rows.2015-01-16
  • 0
    Yeah, I made a mistake. Thank you for your pointing out.2015-01-18
0

Use the superoperator formalism: $$ AXB=C \mapsto ( {B}^T\otimes A)\text{vec} X = C $$ (see here for a definition of $\text{vec}(X)\;$ and here for more information: Kronecker product)

Rewrite your equation (using $AXB\to (B^T\otimes A) \text{vec}\;X$) to $$ AX+XB=C \to (1\otimes A)\text{vec}\;X + (B^T\otimes 1)\text{vec}\;X=\text{vec}\;C\\ \text{vec}\;X=\left((1\otimes A) + (B^T\otimes 1)\right)^{-1}\text{vec}\;C, $$ assuming that the last inverse exists...