Блог

Следующий код инвертирует матрицу A в sympy.

 import sympy
from sympy import MatrixSymbol

A = MatrixSymbol('A', 3, 3)    
A.I
 

Я хотел бы знать, как получить псевдоинверсию A. Это связано с тем, что я выполняю вычисления с использованием универсальных матриц сингулярных квадратов A и продуктов типа A * A. Появляется pinv, и я не хочу, чтобы sympy упрощал их до матрицы идентичности, что происходит, если я использую обратную вместо псевдоинверсии.

Обратите внимание, что в sympy есть функция pinv (), однако, насколько я могу судить, она работает только для матрицы (т. Е. Матриц с указанными элементами), а не для MatrixSymbol (т. Е. универсальных матриц).

Любая помощь будет признательна!