Theano矩阵逆

Question

import numpy as np
from numpy.linalg import inv
import theano.tensor as T
from theano.compile import function
from theano.tensor.nlinalg import matrix_inverse, det
from theano import shared
from theano import config

A = shared(np.matrix('1 0 0; 1 2 0; 0 0 3', dtype=config.floatX))
print(A.get_value())
invA = matrix_inverse(A)
print(invA)

print(invA)语句输出"MatrixInverse.0“。矩阵是有逆的。有人能帮我打印出矩阵逆吗？提前谢谢。

Answer 1

您缺少一个计算符号矩阵逆的Theano函数。

A = T.dmatrix('A')
invA = matrix_inverse(A)
f = theano.function([A], invA)   
print(f(np.matrix('1 0 0; 1 2 0; 0 0 3', dtype=config.floatX)))
#[[ 1.          0.          0.        ]
#[-0.5         0.5         0.        ]
#[ 0.          0.          0.33333333]]