最近几天在复习线性代数，想用Python求解线性方程组，最开始想用SciPy，但是SciPy求解方程组好像要求系数矩阵为方阵，经过几番搜索，在知乎上找到了一点线索，原知乎问题请点我。

根据知乎的答案实践之后发现此方法甚好，不仅可以求解系数矩阵为方阵的线性方程组，还可以求解含自由变量的线性方程组，所以记录在此。

方程组有一组数值解

原方程组如下：

2x+3y+z=4
4x+2y+3z =17
7x+y-z=1

求解方程组的代码Python代码如下：

from sympy import *

a = Matrix([[2,3,1],[4,2,3], [7,1,-1]])
b = Matrix([[4],[17],[1]])
x = symarray('x', 3)
solve(a*x-b)

执行以上代码得到

{x_1: -1, x_2: 5, x_0: 1}

求解存在自由变量的方程组

代码如下：

from sympy import *

a = Matrix([[3,0,-1,0],[8,0,0,-2], [0,2,-2,-1]])
x = symarray('x', 4)
solve(a*x)

得到：

{x_1: 5x_3/4, x_2: 3x_3/4, x_0: x_3/4}

更多请参照 SymPy Doc.