2012-11-26 24 views
2

我很好奇在MATLAB中用於計算SVD的DGESVD函數。據我可以從Gene H. Golub和Charles F. Van Loan的「矩陣計算」中得知,有兩種可能的雙對角化方案 - Householder雙對角化和R-雙對角化。但是,我無法從正在使用的LAPACK文檔中確定。有人有主意嗎?哪個算法在LAPACK中做DGESVD工具?

回答

3

根據LAPACK User's Guide,通過使用Householder反射的常規DGEBRD完成對比對角線形式的縮減。

2.4.6 「奇異值分解」:

例程xGEBRD表示因子形式作爲初級反射器的製品U1和V1,如5.4節中所述。

5.4「的正交或酉矩陣代表」:

一個真正的正交的或複雜的酉矩陣(通常表示Q)在LAPACK通常表示爲基本反射的產物 - - 也被稱爲基本Householder矩陣