如何找到具有一个常数的矩阵的逆矩阵使用数组？

如何解决如何找到具有一个常数的矩阵的逆矩阵使用数组？

我正在尝试找到它的逆，但是numpy似乎只为数字找到逆。我有alpha不等于零。

我该怎么做？

from sympy import *
import numpy as np 
import scipy.linalg 

A = np.array([[1,2,3],[2,[3,4,5]])
B = np.array([[1,[8,9,5]])

alpha = symbols('a')
K = A - alpha*B
print(np.linalg.inv(K))

解决方法

您是正确的，NumPy inv函数需要数字输入。关于SymPy矩阵对象的basic methods的SymPy文档说：“要找到矩阵的逆，只需将其提高到-1的幂即可。”

所以...

from sympy import *
import numpy as np 
import scipy.linalg 

A = np.array([[1,2,3],[2,[3,4,5]])
B = np.array([[1,[8,9,5]])

alpha = symbols('a')
K = A - alpha*B

M = Matrix(K)
print("M:")
print(repr(M))
print()
print("inv(M):")
print(repr(M**-1))

输出：

M:
Matrix([
[  1 - a,3 - 3*a],[2 - 8*a,2 - 2*a,[3 - 3*a,4 - 9*a,5 - 5*a]])

inv(M):
Matrix([
[                      (17*a - 2)/(181*a**2 - 33*a + 2),-(27*a - 2)/(181*a**2 - 33*a + 2),6*a/(181*a**2 - 33*a + 2)],[                      (31*a - 1)/(181*a**2 - 33*a + 2),(4*a - 4)/(181*a**2 - 33*a + 2),-(21*a - 3)/(181*a**2 - 33*a + 2)],[-(66*a**2 - 38*a + 2)/(181*a**3 - 214*a**2 + 35*a - 2),(9*a + 2)/(181*a**2 - 33*a + 2),-(2*a**2 + 12*a - 2)/(181*a**3 - 214*a**2 + 35*a - 2)]])

正如@OscarBenjamin在评论中指出的那样，您还可以使用inv方法，即M.inv()。