Mathematica中的酉对角化

mathematica中是否有一个内置函数或已生成的代码，用于计算给定矩阵的酉对角化

也就是说，如果矩阵是：A，那么我想找到p和D，这样：

---

A=p*dp，其中PP*=I，*是复共轭转置。

这样的矩阵通常不存在。只有正规矩阵（用它们自己的转置/伴随交换）可以用正交变换对角化。对于其他的，你能得到的最好结果是D是上三角形。参见Smith范式和实Smith范式。特征分解算法只是产生这种标准形式。

那么对于标准矩阵，是否已经存在一个代码来计算其酉对角化？正如我所说，这正是所有从Francis QR算法派生的算法所做的。我们必须付出额外的努力来累积所有迭代的旋转，这通常不是在算法的基变量中完成的。