Mathematica中的单一对角化

Question

是否有内置函数或者已经在mathematica中编写了一个计算给定矩阵的单一对角线的代码？

即，如果矩阵是：A，那么我想找到P和D，以便：

A = P * D P，其中PP * = I，*是复共轭转置。

Answer 1

这种矩阵一般不存在。只有正规矩阵（用他们自己的转置/伴随通勤）可以用正交变换对角化。对于其他人来说，你能得到的最好的是D是上三角形。参见史密斯正规形式和真实史密斯正规形式。特征分解算法只产生这种正常形式。