The Sylvester equation is a matrix equation of the form $AX-XB=C,$ where $A,B,C$ are given matrices of dimension $m\times m,n\times n$ and $m\times n$ and $X$ is an unknown matrix of dimension $m\times n.$ It is a well known fact that the equation has an unique solution if and only if the matrices $A$ and $B$ have disjoint spectrum. If they do not have disjoint spectrum, then the result in general depends on $C.$

While determining perturbation of eigenvalues in certain context I was naturally drawn to the the problem of determining the minimum, $min_{X}||AX-XB-C||,$ where $||.||$ is the Frobenius norm. Clearly, if the spectrum of $A$ and $B$ is disjoint then there is a choice of $X$ for which the norm is zero. Otherwise, we need to resort to certain optimization techniques. One approach could be to vectorize the matrices using Kronecker products and determine the minimum of a linear system.

The problem is: "What is the choice of $X$ for which the norm $||AX-XB-C||$ attains minimum (if it exists) when the spectra of $A$ and $B$ are not disjoint?"

I have not found any literature on discussion about similar problems. I would be very thankful for any references or suggestions in this direction.

Is the question asking to find the minimum of the norm as a function of the matrices $A$, $B$ and $C$ in case the spectra of $A$ and $B$ are disjoint?
–
Igor KhavkineDec 26 '12 at 19:57

the minimization version of the problem is a least-squares problem whose solution can be essentially characterized using the SVD of $I \otimes A - B^T\otimes I$....how helpful that is, I'm not sure.
–
SuvritDec 26 '12 at 20:12

Oops, I meant "not disjoint", since that is the non-trivial case.
–
Igor KhavkineDec 27 '12 at 0:16

@Igor Khavkine Yes, my interest is in the case when the matrices have overlapping spectra. I have edited the question to reflect this.
–
UdayDec 27 '12 at 5:52

This equation may or may not have a unique solution, but the unique least-norm vector $x$ that solves it is given by $x=H^+c$, where $H^+$ is the Moore-Penrose pseudoinverse of $H$. Thus, clearly, the optimization problem has a solution, which answers the question as asked.

Of course, due to the numerical expense of computing the above solution (via the truncated SVD of $H$) might be too high. The OP might be interested in consulting the numerical analysis literature on how to deal with such a case.

@Suvrit Thank you for the answer. Your answer settles the fact that there exists a solution. But, my original interest, which I unfortunately did not make it clear in my post, is to get some quantitative information on the bounds of the solution w.r.t the eigenvalues of $A$ and $B$.
–
UdayDec 27 '12 at 16:35