SVD - 奇异值分解
该程序计算矩阵的奇异值分解 A = U * SIGMA * transpose(V)。
用途
SVD 存储过程的语法如下:
- SVD matrixA,matrixU,matrixS,matrixVT
- 参数
- matrixA
- 输入矩阵 A 的名称。
- matrixU
- 输出矩阵 U 的名称。
- matrixS
- 单列输出矩阵 S 的名称。
- matrixVT
- 输出矩阵 transpose(V) 的名称。
详细信息
使用 "调用 nzm..vec_to_diag('S','SIGMA'); "从单列矩阵 S 创建对角矩阵 SIGMA。
示例
call nzm..shape('1,2,3,4,5,6,7,8,9,0',3,3,'A');
call nzm..svd('A', 'U', 'S', 'VT');
call nzm..vec_to_diag('S','SIGMA');
call nzm..gemm('U', 'SIGMA', 'USIGMA');
call nzm..gemm('USIGMA', 'VT', 'A1');
call nzm..subtract('A', 'A1', 'A0');
call nzm..print('A0');
call nzm..delete_matrix('A');
call nzm..delete_matrix('U');
call nzm..delete_matrix('S');
call nzm..delete_matrix('VT');
call nzm..delete_matrix('SIGMA');
call nzm..delete_matrix('USIGMA');
call nzm..delete_matrix('A0');
call nzm..delete_matrix('A1');
SHAPE
-------
t
(1 row)
SVD
-----
t
(1 row)
VEC_TO_DIAG
-------------
t
(1 row)
GEMM
------
t
(1 row)
GEMM
------
t
(1 row)
SUBTRACT
----------
t
(1 row)
PRINT
---------------------------------------------------------------------------------------------------------------
-- matrix: A0 --
-1.7763568394003e-15, -6.2172489379009e-15, -3.1086244689504e-15
-1.7763568394003e-15, -1.7763568394003e-15, -4.4408920985006e-15
-2.6645352591004e-15, -7.105427357601e-15, -7.105427357601e-15
(1 row)
DELETE_MATRIX
---------------
t
(1 row)
DELETE_MATRIX
---------------
t
(1 row)
DELETE_MATRIX
---------------
t
(1 row)
DELETE_MATRIX
---------------
t
(1 row)
DELETE_MATRIX
---------------
t
(1 row)
DELETE_MATRIX
---------------
t
(1 row)
DELETE_MATRIX
---------------
t
(1 row)
DELETE_MATRIX
---------------
t
(1 row)