解四元线性方程组

更新时间:2023-04-02 12:57:09

问题阐述

问题:我需要用Python解这些方程。

a + 3b + 2c + 2d = 1
2a + b + c + 2d = 0
3a + b + 2c + d = 1
2a + c + 3d = 0

这样我就可以得到a、b、c和d的值。有没有办法以分数形式显示它们?

我的代码:

import numpy as np
A = np.array([[1,3,2,2],[2,1,1,2]])
B = np.array([1,0,1,0])
X2 = np.linalg.solve(A,B)

错误:

LinAlgError: Last 2 dimensions of the array must be square

精准答案

您没有将问题的最后两个方程式添加到A矩阵:

import numpy as np
A = np.array([[1,3,2,2],[2,1,1,2],[3,1,2,1],[2,0,1,3]])
B = np.array([1,0,1,0])
X2 = np.linalg.solve(A,B)

给予:

array([-0.27272727, -0.18181818,  1.09090909, -0.18181818])

这应该是可行的。