线性方程组Python求解

如下所示,可以用scipy模块和numpy模块求解线性方程组


from scipy import linalg
import numpy as np
A = np.array([[1, 1, 7], [2, 3, 5], [4, 2, 6]])  # A代表系数矩阵
b = np.array([2, 3, 4])  # b代表常数列
x = linalg.solve(A, b)
print(x)

最后的结果如下:

[0.6  0.35 0.15]

但是,为啥我模仿了一下,就不可以了,
代码如下:

from scipy import linalg
import numpy as np
A=np.array([[1/201,1/202],[1/202,1/203]])
b=np.array([1/103,1/104])
x=linalg.solve(A,b)
print(x)

提示的报错如下:

img

是数据类型的问题吗?谢谢大家!

矩阵A是奇异矩阵就等价于A的行列式等于0,如果你用的是python2的话,改成:

A=np.array([[1.0/201.0,1.0/202.0],[1.0/202.0,1.0/203.0]])
b=np.array([1.0/103.0,1.0/104.0])