python用fsolve、leastsq对非线性方程组求解


Posted in Python onDecember 15, 2018

背景:

实现用python的optimize库的fsolve对非线性方程组进行求解。可以看到这一个问题实际上还是一个优化问题,也可以用之前拟合函数的leastsq求解。下面用这两个方法进行对比:

代码:

from scipy.optimize import fsolve,leastsq
from math import sin,cos
 
def f(x):
 x0 = float(x[0])
 x1 = float(x[1])
 x2 = float(x[2])
 return [
 5*x1+3,
 4*x0*x0 - 2*sin(x1*x2),
 x1*x2-1.5
 ]
 
x0 = [1,1,1]
result = fsolve(f,x0)
 
print("===================")
print()
print("求解函数名称:",fsolve.__name__)
print("解:",result)
print("各向量值:",f(result))
#拟合函数来求解
h = leastsq(f,x0)
 
print("===================")
print()
print("求解函数名称:",leastsq.__name__)
print("解:",h[0])
print("各向量的值:",f(h[0]))

结果:

===================

求解函数名称: fsolve
解: [-0.70622057 -0.6        -2.5       ]
各向量值: [0.0, -9.126033262418787e-14, 5.329070518200751e-15]
===================

求解函数名称: leastsq
解: [-0.70622057 -0.6        -2.5       ]
各向量的值: [0.0, -2.220446049250313e-16, 0.0]

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
paramiko模块安装和使用(远程登录服务器)
Jan 27 Python
用pickle存储Python的原生对象方法
Apr 28 Python
详解Python map函数及Python map()函数的用法
Nov 16 Python
利用Python如何制作好玩的GIF动图详解
Jul 11 Python
在python中pandas的series合并方法
Nov 12 Python
python元组和字典的内建函数实例详解
Oct 22 Python
Python数据分析pandas模块用法实例详解
Nov 20 Python
Python tensorflow实现mnist手写数字识别示例【非卷积与卷积实现】
Dec 19 Python
Python使用GitPython操作Git版本库的方法
Feb 29 Python
利用Python裁切tiff图像且读取tiff,shp文件的实例
Mar 10 Python
python爬虫scrapy图书分类实例讲解
Nov 23 Python
Pytorch中Softmax和LogSoftmax的使用详解
Jun 05 Python
python实现一组典型数据格式转换
Dec 15 #Python
python判断计算机是否有网络连接的实例
Dec 15 #Python
Django model反向关联名称的方法
Dec 15 #Python
django orm 通过related_name反向查询的方法
Dec 15 #Python
关于Django ForeignKey 反向查询中filter和_set的效率对比详解
Dec 15 #Python
django 外键model的互相读取方法
Dec 15 #Python
Django之Mode的外键自关联和引用未定义的Model方法
Dec 15 #Python
You might like
论建造顺序的重要性
2020/03/04 星际争霸
十天学会php(2)
2006/10/09 PHP
PHP 日期时间函数的高级应用技巧
2009/10/10 PHP
php array_keys 返回数组的键名
2016/10/25 PHP
thinkphp5+layui实现的分页样式示例
2019/10/08 PHP
基于jquery实现的一个选择中国大学的弹框 (数据、步骤、代码)
2012/07/26 Javascript
整理的比较全的event对像在ie与firefox浏览器中的区别
2013/11/25 Javascript
JavaScript中的函数的两种定义方式和函数变量赋值
2014/05/12 Javascript
jquery 页眉单行信息滚动显示实现思路及代码
2014/06/26 Javascript
javascript操作符"!~"详解
2015/02/10 Javascript
javascript中Math.random()使用详解
2015/04/15 Javascript
Bootstarp风格的toggle效果分享
2016/02/23 Javascript
详解AngularJS 模态对话框
2016/04/07 Javascript
JS中取二维数组中最大值的方法汇总
2016/04/17 Javascript
JavaScript中如何使用cookie实现记住密码功能及cookie相关函数介绍
2016/11/10 Javascript
js中DOM三级列表(代码分享)
2017/03/20 Javascript
ES6 Object属性新的写法实例小结
2019/06/25 Javascript
介绍Python的Django框架中的静态资源管理器django-pipeline
2015/04/25 Python
Python中在脚本中引用其他文件函数的实现方法
2016/06/23 Python
python+pandas分析nginx日志的实例
2018/04/28 Python
Python 绘制酷炫的三维图步骤详解
2019/07/12 Python
Python实现使用dir获取类的方法列表
2019/12/24 Python
东南亚地区最大的购物网站Lazada新加坡站点:Lazada.sg
2016/07/17 全球购物
Oakley官网:运动太阳镜、雪镜和服装
2016/09/30 全球购物
加拿大时尚床上用品零售商:QE Home | Quilts Etc
2018/01/22 全球购物
ASOS亚洲:ASOS Asia
2018/03/04 全球购物
C++:memset ,memcpy和strcpy的根本区别
2013/04/27 面试题
数据库的约束含义
2012/09/09 面试题
领班岗位职责范文
2014/02/06 职场文书
商场消防演习方案
2014/02/12 职场文书
高中微机老师自我鉴定
2014/02/16 职场文书
中药专业自荐信范文
2014/03/18 职场文书
申论倡议书范文
2014/05/13 职场文书
公司员工离职证明书
2014/10/04 职场文书
vue3使用vue-router的完整步骤记录
2021/06/20 Vue.js
mysql联合索引的使用规则
2021/06/23 MySQL