2017-09-22 44 views
0

A是N乘N矩陣。 I是(N-2)乘以(N-2)的單位矩陣。 B是另一個N乘N矩陣,其被定義爲有沒有可以解決以下特徵值的Fortran解算器?

B=[I 0 0; 
     0 0 0; 
     0 0 0] 

x是一個有N個元素的數組。如何解決以下形式的本徵值 X = C X,其中c是本徵值, 通過使用特徵值求解器?

+1

另外谷歌爲[廣義特徵值問題](https://en.wikipedia.org/wiki/Eigendecomposition_of_a_matrix#Generalized_eigenvalue_problem) –

回答

1

您可以查看Lapack庫,它提供了廣義矩陣的特徵值問題的解決方案。根據您的數據類型和矩陣類型,您將需要使用不同的子程序。 在這方面看看here。也可以看看here作爲用作矩陣類型的術語。

最後,我之前寫了this module來給出一些Lapack功能的實現示例,其中包括特徵值問題。在那裏你可以找到一個單值的廣義矩陣(sgeev)。