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 相关文章推荐
python3使用urllib示例取googletranslate(谷歌翻译)
Jan 23 Python
python中的五种异常处理机制介绍
Sep 02 Python
python版微信跳一跳游戏辅助
Jan 11 Python
对python中array.sum(axis=?)的用法介绍
Jun 28 Python
详解python异步编程之asyncio(百万并发)
Jul 07 Python
代码详解django中数据库设置
Jan 28 Python
使用django的ORM框架按月统计近一年内的数据方法
Jul 18 Python
SELENIUM自动化模拟键盘快捷键操作实现解析
Oct 28 Python
基于Python中isfile函数和isdir函数使用详解
Nov 29 Python
Python实现屏幕录制功能的代码
Mar 02 Python
python正则表达式 匹配反斜杠的操作方法
Aug 07 Python
python实现移动木板小游戏
Oct 09 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
php正则表达式(regar expression)
2011/09/10 PHP
PHP file_exists问题杂谈
2012/05/07 PHP
PHP中的替代语法简介
2014/08/22 PHP
php防止sql注入的方法详解
2017/02/20 PHP
非阻塞动态加载javascript广告实现代码
2010/11/17 Javascript
深入理解JavaScript系列(16) 闭包(Closures)
2012/04/12 Javascript
杨氏矩阵查找的JS代码
2013/03/21 Javascript
使用apply方法实现javascript中的对象继承
2013/12/16 Javascript
javascript面向对象之定义成员方法实例分析
2015/01/13 Javascript
jQuery中wrapInner()方法用法实例
2015/01/16 Javascript
jQuery实现的经典滑动门效果
2015/09/22 Javascript
纯js实现瀑布流布局及ajax动态新增数据
2016/04/07 Javascript
微信小程序 Record API详解及实例代码
2016/09/30 Javascript
Angular使用ng-messages与PHP进行表单数据验证
2016/12/28 Javascript
使用openSpeDiv方法实现Ecshop登录弹窗框效果
2017/03/13 Javascript
Angular2+国际化方案(ngx-translate)的示例代码
2017/08/23 Javascript
详解基于 Nuxt 的 Vue.js 服务端渲染实践
2017/10/24 Javascript
webpack打包js文件及部署的实现方法
2017/12/18 Javascript
vue2单元测试环境搭建
2018/05/24 Javascript
Angular 2使用路由自定义弹出组件toast操作示例
2019/05/10 Javascript
Vue 3.0 全家桶抢先体验
2020/04/28 Javascript
[41:21]夜魇凡尔赛茶话会 第三期02:看图识人
2021/03/11 DOTA
python解析中国天气网的天气数据
2014/03/21 Python
将字典转换为DataFrame并进行频次统计的方法
2018/04/08 Python
python3+PyQt5使用数据库窗口视图
2018/04/24 Python
浅析python中的迭代与迭代对象
2018/10/08 Python
python3实现zabbix告警推送钉钉的示例
2019/02/20 Python
Python转换时间的图文方法
2019/07/01 Python
python实现几种归一化方法(Normalization Method)
2019/07/31 Python
详解centos7+django+python3+mysql+阿里云部署项目全流程
2019/11/15 Python
css3新单位vw、vh的使用教程
2018/03/23 HTML / CSS
德国原装品牌香水、化妆品和手表网站:BRASTY.DE
2016/10/16 全球购物
澳大利亚购买太阳镜和眼镜网站:Glamoureyes
2020/09/22 全球购物
网络教育自我鉴定
2013/11/01 职场文书
住宅质量保证书
2014/04/29 职场文书
实现AJAX异步调用和局部刷新的基本步骤
2022/03/17 Javascript