python如何解方程的三种方法

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了python如何解方程的三种方法，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含2315字，纯文字阅读大概需要4分钟。

内容图文

相关推荐：《python视频》

Numpy求解方程组

x + 2y = 3
4x ＋ 5y = 6

当然我们可以手动写出解析解，然后写一个函数来求解，这实际上只是用 Python 来单纯做“数值计算”. 但实际上，numpy.linalg.solve 可以直接求解线性方程组.

一般地，我们设解线性方程组形如 Ax=b，其中 A 是系数矩阵，b 是一维(n 维也可以，这个下面会提到)，x 是未知变量. 再拿上面地最简单的二元一次方程组为例，我们用 numpy.linalg.solve 可以这样写：

In [1]: import numpy as np
  ...: A = np.mat('1,2; 4,5')    # 构造系数矩阵 A
  ...: b = np.mat('3,6').T       # 构造转置矩阵 b （这里必须为列向量）
  ...: r = np.linalg.solve(A,b)  # 调用 solve 函数求解
  ...: print r
  ...:
Out[1]: [[-1.]
      [ 2.]]

那么前面提到的“ n 维”情形是什么呢？实际上就是同时求解多组形式相同的二元一次方程组，例如我们想同时求解这样两组：

x + 2y = 3
4x ＋ 5y = 6

和

x + 2y = 7
4x ＋ 5y = 8

就可以这样写：

In [2]: import numpy as np
  ...: A = np.mat('1,2; 4,5')          # 构造系数矩阵 A
  ...: b = np.array([[3,6], [7,8]]).T  # 构造转置矩阵 b （这里必须为列向量），
  ...: 注意这里用的是 array
  ...: r = np.linalg.solve(A,b)        # 调用 solve 函数求解
  ...: print r
  ...:
Out[2]: [[-1.         -6.33333333]
      [ 2.          6.66666667]]

SciPy 求解非线性方程组

一般来说，我们只需要用到 func 和 x0 就够了. func 是自己构造的函数，也就是需要求解的方程组的左端（右端为 0），而 x0 则是给定的初值.

我们来看一个具体的例子，求解：

x + 2y + 3z - 6 = 0
5 * (x ** 2) + 6 * (y ** 2) + 7 * (z ** 2) - 18 = 0
9 * (x ** 3) + 10 * (y ** 3) + 11 * (z ** 3) - 30 = 0

就可以这么写：

In [3]: from scipy.optimize import fsolve
  ...:
  ...: def func(i):
  ...:     x, y, z = i[0], i[1], i[2]
  ...:     return [
  ...:             x + 2 * y + 3 * z - 6,
  ...:             5 * (x ** 2) + 6 * (y ** 2) + 7 * (z ** 2) - 18,
  ...:             9 * (x ** 3) + 10 * (y ** 3) + 11 * (z ** 3) - 30
  ...:            ]
  ...:
  ...: r = fsolve(func,[0, 0, 0])
  ...: print r
  ...:
Out[3]: [ 1.00000001  0.99999998  1.00000001]

当然，SciPy 也可以用来求解线性方程组，这是因为 scipy.optimize.fsolve 本质上是最小二乘法来逼近真实结果.

SymPy 求解方程组

例如求解一个：

x + 2 * (x ** 2) + 3 * (x ** 3) - 6 = 0

直接就是：

In [4]: from sympy import *
  ...: x = symbols('x')
  ...: solve(x + 2 * (x ** 2) + 3 * (x ** 3) - 6, x)
Out[4]: [1, -5/6 - sqrt(47)*I/6, -5/6 + sqrt(47)*I/6]

以上就是python如何解方程的三种方法的详细内容，更多请关注Gxl网其它相关文章！

内容总结

以上是互联网集市为您收集整理的python如何解方程的三种方法全部内容，希望文章能够帮你解决python如何解方程的三种方法所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端