先把方程变型成如下形式：
aX+Xb=-c
这个是Sylvester方程,其中的a,b,c你自己转换吧,因为有矩阵逆的符号,我不好表示.（注意c前面有个负号）
然后你用MATLAB自带的工具箱中的lyap()函数求解.命令输入如下：
X=lyap(a,b,c)