Python解方程组 scipy.optimize.fsolve()函数如何求解带有循环求和的方程式

 更新时间:2023年06月15日 09:19:11   作者:qq_41846459  
这篇文章主要介绍了Python解方程组 scipy.optimize.fsolve()函数如何求解带有循环求和的方程式,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下

最近做的课题需要求解方程组,方程组还挺复杂,未知参数比较多

简单粗暴介绍fsolve()函数用法:

假如要解方程组

x*y=6
x+y=5

先变换成右边为0的形式

x*y-6
x+y-5

func函数中直接替换要解的方程就可

import scipy
from scipy.optimize import fsolve
def func(i):
    x, y = i[0], i[1]
    return [  # 这里写要求解的方程组式子,变成等于0的形式
        x*y-6,
        x+y-5 ]
if __name__ == '__main__':
    r = scipy.optimize.fsolve(func, [0, 0])
    print(r)

结果

[2. 3.]

再看看有参数的解法:

如果有参数a,b需要传入

axy-6
x+by-5

import scipy
from scipy.optimize import fsolve
def func(i, params):
    x, y = i[0], i[1]
    a, b = params
    return [  # 这里写要求解的方程组式子,变成等于0的形式
        a*x*y-6,
        x+b*y-5 ]
if __name__ == '__main__':
    a = 2
    b = 0.5
    params = [a, b]
    r = scipy.optimize.fsolve(func, [0, 0], args=params)
    print(r)

结果

[0.32055053 9.35889894]

近似解

更复杂的方程式带有循环求和:

附三角函数写法

import math
math.sin(theta)

带求和符号的方程求解(这里举一个简单的例子好理解一些)

这里用sum来计算和

import scipy
from scipy.optimize import fsolve
def func(i, params):
    x, y = i[0], i[1]
    a, b = params
    list_e1 = [x * a[j][0] + y * b[j][0] - 5 for j in range(2)]
    list_e2 = [x * a[j][1] - y * b[j][1] - 2 for j in range(2)]
    return [  # 这里写要求解的方程组式子,变成等于0的形式
        sum(list_e1),
        sum(list_e2)
    ]
if __name__ == '__main__':
    a = [[3, 4], [3, 4]]
    b = [[1, 3], [4, 1]]
    params = [a, b]
    r = scipy.optimize.fsolve(func, [0, 0], args=params)
    print(r)

结果

[0.9375 0.875 ]

到此这篇关于Python解方程组 scipy.optimize.fsolve()函数如何求解带有循环求和的方程式的文章就介绍到这了,更多相关Python scipy.optimize.fsolve()函数内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!

相关文章

  • 教你pycharm运行Django第一个项目

    教你pycharm运行Django第一个项目

    本文主要介绍了教你pycharm运行Django第一个项目的实现,文中通过图文示例介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2021-08-08
  • 浅谈Django自定义模板标签template_tags的用处

    浅谈Django自定义模板标签template_tags的用处

    这篇文章主要介绍了浅谈Django自定义模板标签template_tags的用处,具有一定借鉴价值,需要的朋友可以参考下。
    2017-12-12
  • Python打印scrapy蜘蛛抓取树结构的方法

    Python打印scrapy蜘蛛抓取树结构的方法

    这篇文章主要介绍了Python打印scrapy蜘蛛抓取树结构的方法,实例分析了打印scrapy蜘蛛抓取树结构的技巧,非常具有实用价值,需要的朋友可以参考下
    2015-04-04
  • 详解Python 中的命名空间、变量和范围

    详解Python 中的命名空间、变量和范围

    Python 是一种动态类型语言,在程序执行期间,变量名可以绑定到不同的值和类型,这篇文章主要介绍了Python 中的命名空间、变量和范围,需要的朋友可以参考下
    2022-09-09
  • python使用multiprocessing的详细方法

    python使用multiprocessing的详细方法

    multiprocessing是Python标准库中的一个模块,用于实现多进程编程,它提供了一种简单而高效的方式来利用多核处理器的能力,通过在多个进程中同时执行任务,加快程序的执行速度和提高系统的吞吐量,这篇文章主要介绍了python使用multiprocessing,需要的朋友可以参考下
    2024-03-03
  • linux中使用Python对图片进行批量命名

    linux中使用Python对图片进行批量命名

    大家好,本篇文章主要讲的是linux中使用Python对图片进行批量命名,感兴趣的同学赶快来看一看吧,对你有帮助的话记得收藏一下,方便下次浏览
    2021-12-12
  • VSCode基础使用与VSCode调试python程序入门的图文教程

    VSCode基础使用与VSCode调试python程序入门的图文教程

    这篇文章主要介绍了VSCode基础使用+VSCode调试python程序入门图文教程,本文通过图文并茂的形式给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下
    2020-03-03
  • python Django的web开发实例(入门)

    python Django的web开发实例(入门)

    这篇文章主要介绍了python Django的web开发实例,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2019-07-07
  • Python中的字符串查找操作方法总结

    Python中的字符串查找操作方法总结

    这里我们来整理一下Python中的字符串查找操作方法总结,除了基本的find()方法外,还会讲解到朴素匹配算法和KMP算法的使用:
    2016-06-06
  • pytorch的batch normalize使用详解

    pytorch的batch normalize使用详解

    今天小编就为大家分享一篇pytorch的batch normalize使用详解,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2020-01-01

最新评论