多语言展示
当前在线:641今日阅读:103今日分享:49

python 线性代数:[7]解多元一次方程

用python的numpy包中的linalg.solve()方法解多元一次方程,如果你对矩阵解方程非常熟悉,那么现在只是学习一下这个函数就好了。如果你不是很熟悉用矩阵解方程,你需要看一下线性代数的课本。首先看一下我们要解的方程,将这个方程格式调整好,按照x-y-z-常数项的顺序排列将未知数的系数写下来,排列成一个矩阵a,如下常数项构成一个一维数组(向量)使用linalg.solve方法解方程,参数a指的是系数矩阵,参数b指的是常数项矩阵我们得到的解对不对呢?使用点乘的方法可以验证一下,系数乘以未知数可以得到常数项今天用到的所有代码如下>>> a=[[1,2,1],[2,-1,3],[3,1,2]]>>> a=np.array(a)>>> aarray([[ 1, 2, 1], [ 2, -1, 3], [ 3, 1, 2]])>>>>>>>>> b=[7,7,18]>>> b=np.array(b)>>> barray([ 7, 7, 18])>>>>>>>>>>>> x=np.solve(a,b)Traceback (most recent call last): File '', line 1, in x=np.solve(a,b)AttributeError: 'module' object has no attribute 'solve'>>> x=np.linalg.solve(a,b)>>> xarray([ 7., 1., -2.])>>>>>>>>>>>>>>> np.dot(a,x)array([ 7., 7., 18.])
推荐信息