python 行列式_行列式计算with python
问题源于同济大学《线性代数》(第五版)，先看理论：
问题是这样的：
教材上的算术求解方法如下：
接下来我们用python编程求解，代码如下：
import numpy as np
#定义矩阵
D=np.array([[2,1,-5,1],[1,-3,0,-6],[0,2,-1,2],[1,4,-7,6]])
#输出矩阵D行列式运算的值
d=np.linalg.det(D)
print "D=",d
D1=np.array([[8,1,-5,1],[9,-3,0,-6],[-5,2,-1,2],[0,4,-7,6]])
#输出矩阵D1行列式运算的值
d1=np.linalg.det(D1)
print "D1=",d1
D2=np.array([[2,8,-5,1],[1,9,0,-6],[0,-5,-1,2],[1,0,-7,6]])
#输出矩阵D2行列式运算的值
d2=np.linalg.det(D2)
print "D2=",d2
D3=np.array([[2,1,8,1],[1,-3,9,-6],[0,2,-5,2],[1,4,0,6]])
#输出矩阵D3行列式运算的值
d3=np.linalg.det(D3)
print "D2=",d3
D4=np.array([[2,1,-5,8],[1,-3,0,9],[0,2,-1,-5],[1,4,-7,0]])
#输出矩阵D4行列式运算的值
d4=np.linalg.det(D4)
print "D2=",d4
#求解x1,x2,x3,x4
x1=d1/d
x2=d2/d
x3=d3/d
x4=d4/d
print "x1=",x1
print "x2=",x2
print "x3=",x3
print "x4=",x4
问题源于同济大学《线性代数》(第五版)，先看理论： 问题是这样的： 教材上的算术求解方法如下： 接下来我们用python编程求解，代码如下： import numpy as np #定义矩阵 D=np.array([[2,1,-5,1],[1,-3,0,-6],[0,2,-1,2],[1,4,-7,6]]) #输出矩阵D行列式运算的值 d=np.linalg.det(D) print "D=",d D1=np.array([[8,1,-5,1],[9,-3,0,-6],[-5,2,-1,2],[0,4,-7,6]]) #输出矩阵D1行列式运算的值 d1=np.linalg.det(D1) print "D1=",d1 D2=np.array([[2,8,-5,1],[1,9,0,-6],[0,-5,-1,2],[1,0,-7,6]]) #输出矩阵D2行列式运算的值 d2=np.linalg.det(D2) print "D2=",d2 D3=np.array([[2,1,8,1],[1,-3,9,-6],[0,2,-5,2],[1,4,0,6]]) #输出矩阵D3行列式运算的值 d3=np.linalg.det(D3) print "D2=",d3 D4=np.array([[2,1,-5,8],[1,-3,0,9],[0,2,-1,-5],[1,4,-7,0]]) #输出矩阵D4行列式运算的值 d4=np.linalg.det(D4) print "D2=",d4 #求解x1,x2,x3,x4 x1=d1/d x2=d2/d x3=d3/d x4=d4/d print "x1=",x1 print "x2=",x2 print "x3=",x3 print "x4=",x4