如何像 Matlab 一样使 LAPACK 特征向量正交?

计算科学 线性代数 特征值 本征系统 拉帕克
2021-12-17 21:34:55

我正在使用 LAPACKzgeev计算高维对称复数矩阵的特征向量(n2000)。我需要这些特征向量来满足

m=1n|mm¯|m¯|m=In

在哪里|m是一个特征向量并且m¯|是相同的特征向量转置(注意:未共轭)和In是大小的单位矩阵n.

我不知道是否所有对称复数矩阵都满足这一点,但这个矩阵具有特定的结构,它总是应该:

n=6

在此处输入图像描述

(生成这些矩阵的算法太复杂了,无法分享,但我认为可以从这两个例子中得到它们的模式)

n=9

在此处输入图像描述

但是我得到的特征向量zgeev不遵循我提到的方程,它们不是正交的。

我已经在 Matlab 中测试了这些矩阵,Matlab 生成的特征向量满足这个条件并且是正交的。(IEmi|m¯j=δij)

如何从 LAPACK 获得相同的特征向量zgeev

0个回答
没有发现任何回复~