我正在编写一种基于局部惯性矩处理3D图像的算法。
我有一个3x3的实对称矩阵,需要从中找到特征值。我已经找到了各种用于矩阵对角化的通用算法,但是我不知道是否存在针对该矩阵的3个特征的解析表达式。
精通数学的人会知道吗?
编辑
作为记录,这是我自己在问题上发现的。正如Matthias Odisio所说,一旦有了3x3矩阵,就无法简单地分析一个简单的表达式。
但是,我已经找到了3x3 Hermitian矩阵这种特殊情况的专用纸,其中比较了各种数值专用方法:
http://arxiv.org/abs/physics/0610206
这是本文的C和Fortran代码:
http://www.mpi-hd.mpg.de/personalhomes/globes/3x3/index.html