python - 使用 sympy 的行列式

Question

以下代码根据 x、y、z 生成一个 3x3 矩阵。我想生成矩阵的行列式。却做不到。

import numpy as np
import sympy as sp
from sympy import *
from sympy.matrices import Matrix
x,y,z =sp.symbols('x,y,z')
H1=np.array([[x,y,z]])
H2=np.array([[0.0630,0.0314,-0.0001],[0.0314,96.1659,-0.0001],[-0.0001,-0.0001,0.0001]])
H3=H1.T
H=H1*H2*H3
print H

要找到上述矩阵的行列式，我使用以下命令。

H.det()

但是显示错误

AttributeError: 'numpy.ndarray' object has no attribute 'det' 

Answer 1

首先需要将你的numpy n维数组转换为sympy矩阵，然后进行符号行列式的计算。 @hildensia 所说的是行不通的，因为 H 是一个 numpy 对象，它不适用于符号实体。

>>> M = sp.Matrix(H)
>>> M
Matrix([
[ 0.063*x**2,   0.0314*x*y, -0.0001*x*z],
[ 0.0314*x*y, 96.1659*y**2, -0.0001*y*z],
[-0.0001*x*z,  -0.0001*y*z, 0.0001*z**2]])
>>> M.det()
0.000604784913*x**2*y**2*z**2