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 相关文章推荐
python实现通过shelve修改对象实例
Sep 26 Python
用PyQt进行Python图形界面的程序的开发的入门指引
Apr 14 Python
python实现红包裂变算法
Feb 16 Python
python数据结构学习之实现线性表的顺序
Sep 28 Python
将python运行结果保存至本地文件中的示例讲解
Jul 11 Python
在python中用url_for构造URL的方法
Jul 25 Python
python根据多个文件名批量查找文件
Aug 13 Python
python3 BeautifulSoup模块使用字典的方法抓取a标签内的数据示例
Nov 28 Python
Python pandas库中的isnull()详解
Dec 26 Python
Pytorch 计算误判率,计算准确率,计算召回率的例子
Jan 18 Python
利用python汇总统计多张Excel
Sep 22 Python
解决tensorflow模型压缩的问题_踩坑无数,总算搞定
Mar 02 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
多文件上传的例子
2006/10/09 PHP
让PHP更快的提供文件下载的代码
2012/06/13 PHP
使用PHP和JavaScript判断请求是否来自微信内浏览器
2015/08/18 PHP
php foreach如何跳出两层循环(详解)
2016/11/05 PHP
phpcms的分类名称和类别名称的调用
2017/01/05 PHP
Laravel框架实现model层的增删改查(CURD)操作示例
2018/05/12 PHP
ThinkPHP框架下微信支付功能总结踩坑笔记
2019/04/10 PHP
一个简单的JavaScript数据缓存系统实现代码
2010/10/24 Javascript
EXTJS记事本 当CompositeField遇上RowEditor
2011/07/31 Javascript
js不能跳转到上一页面的问题解决方法
2013/03/01 Javascript
js浏览器本地存储store.js介绍及应用
2014/05/13 Javascript
js实现字符串转日期格式的方法
2015/05/20 Javascript
在JavaScript中call()与apply()区别
2016/01/22 Javascript
JS 组件系列之Bootstrap Table 冻结列功能IE浏览器兼容性问题解决方案
2017/06/30 Javascript
深入理解ES6中let和闭包
2018/02/22 Javascript
CKeditor富文本编辑器使用技巧之添加自定义插件的方法
2019/06/14 Javascript
微信小程序利用Canvas绘制图片和竖排文字详解
2019/06/25 Javascript
vue-video-player实现实时视频播放方式(监控设备-rtmp流)
2020/08/10 Javascript
Python中异常重试的解决方案详解
2017/05/05 Python
浅谈Python实现Apriori算法介绍
2017/12/20 Python
python 日期操作类代码
2018/05/05 Python
python中将正则过滤的内容输出写入到文件中的实例
2018/10/21 Python
Django基础知识 web框架的本质详解
2019/07/18 Python
selenium+Chrome滑动验证码破解二(某某网站)
2019/12/17 Python
pytorch 模拟关系拟合——回归实例
2020/01/14 Python
Opencv常见图像格式Data Type及代码实例
2020/11/02 Python
Python3.9.0 a1安装pygame出错解决全过程(小结)
2021/02/02 Python
Python对excel的基本操作方法
2021/02/18 Python
美国男士内衣品牌:Tommy John
2017/12/22 全球购物
倩碧澳大利亚官网:Clinique澳大利亚
2019/07/22 全球购物
几道数据库的面试题或笔试题
2014/05/31 面试题
求职信模板标准格式范文
2014/02/23 职场文书
小学生新年寄语
2014/04/03 职场文书
2015年妇委会工作总结
2015/05/22 职场文书
Java无向树分析 实现最小高度树
2022/04/09 Javascript
在MySQL中你成功的避开了所有索引
2022/04/20 MySQL