西尔维斯特方程（Sylvester equation）AW+WB=C求解

可以求解 AW+WB=C 形式的矩阵方程。
参考：https://zhuanlan.zhihu.com/p/37965266

---

#1 Vec(X) 理解

表示矩阵的按列拼装向量化。如果有

则

在matlab中可以使用 X( : ) 的形式得到。

#2 AXB+CXD=E求解

一种朴实的解法是参考sylverster的朴素解法，两边做向量化, 即表示为：
$$(B^T\otimes A+D^T\otimes C)vec(X)=vec(E)，变成解线性方程$$
期间，Bartels-Stewart算法的求解方式能提高性能。
参：https://zhuanlan.zhihu.com/p/105326895