4.1 行列式

import numpy as np

# 代码4-1
# 方法一：使用det函数求行列式
arr = np.array([[1, 1],[30, 20]])      # 创建分母的二维数组
arr1 = np.array([[80, 1],[2050, 20]])  # 创建分子的二维数组
arr2 = np.array([[1, 80],[30, 2050]])  # 创建分子的二维数组
# 求解行列式
D = np.linalg.det(arr)
D1 = np.linalg.det(arr1)
D2 = np.linalg.det(arr2)
print('方程组的解x1为：', D1 / D)
print('方程组的解x2为：', D2 / D)

方程组的解x1为： 44.99999999999996
方程组的解x2为： 35.000000000000014

# 方法二：使用solve函数解线性方程组
D = np.array([[1, 1],[30, 20]])  # 创建系数行列式
arr = np.array([80, 2050])
x = np.linalg.solve(D, arr)
print('方程组的解为：', x)

方程组的解为： [45. 35.]

# 代码4-2
arr = np.array([[4, 6, 8],[4, 6, 9],[5, 6, 8]])  # 创建三阶行列式
print('行列式的解为：', np.linalg.det(arr))

行列式的解为： 5.999999999999996

# 代码4-3
# 方法一：使用d