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的异常处理
Jun 19 Python
Python 实现 贪吃蛇大作战 代码分享
Sep 07 Python
python append、extend与insert的区别
Oct 13 Python
Python实现文件信息进行合并实例代码
Jan 17 Python
Python lambda函数基本用法实例分析
Mar 16 Python
Python使用pyautogui模块实现自动化鼠标和键盘操作示例
Sep 04 Python
python重试装饰器的简单实现方法
Jan 31 Python
python 绘制拟合曲线并加指定点标识的实现
Jul 10 Python
解决python3 安装不了PIL的问题
Aug 16 Python
浅谈Python访问MySQL的正确姿势
Jan 07 Python
IDLE下Python文件编辑和运行操作
Apr 25 Python
Python使用Chrome插件实现爬虫过程图解
Jun 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
CI框架在CLI下执行占用内存过大问题的解决方法
2014/06/17 PHP
PHP实现数据分页显示的简单实例
2016/05/26 PHP
php的4种常用运行方式详解
2016/12/22 PHP
jQuery的一些注意
2006/12/06 Javascript
分享别人写的一个小型js框架
2007/08/13 Javascript
javascript 子窗体父窗体相互传值方法
2010/05/31 Javascript
有关javascript的性能优化 (repaint和reflow)
2013/04/12 Javascript
div拖拽插件——JQ.MoveBox.js(自制JQ插件)
2013/05/17 Javascript
Js实现当前点击a标签变色突出显示其他a标签回复原色
2013/11/27 Javascript
动态创建script在IE中缓存js文件时导致编码的解决方法
2014/05/04 Javascript
使用Plupload实现直接上传附件至七牛云存储
2014/12/26 Javascript
JavaScript模块化之使用requireJS按需加载
2017/04/12 Javascript
jQuery树插件zTree使用方法详解
2017/05/02 jQuery
微信小程序上滑加载下拉刷新(onscrollLower)分批加载数据(一)
2017/05/11 Javascript
json 带斜杠时如何解析的实现
2019/08/12 Javascript
Vue 根据条件判断van-tab的显示方式
2020/08/03 Javascript
微信小程序将页面按钮悬浮固定在底部的实现代码
2020/10/29 Javascript
Python模拟登陆淘宝并统计淘宝消费情况的代码实例分享
2016/07/04 Python
python email smtplib模块发送邮件代码实例
2018/04/26 Python
对python实现合并两个排序链表的方法详解
2019/01/23 Python
python实现AES和RSA加解密的方法
2019/03/28 Python
Python TestCase中的断言方法介绍
2019/05/02 Python
Flask配置Cors跨域的实现
2019/07/12 Python
PyCharm无法识别PyQt5的2种解决方法,ModuleNotFoundError: No module named 'pyqt5'
2020/02/17 Python
Python如何使用turtle库绘制图形
2020/02/26 Python
Django 解决model 反向引用中的related_name问题
2020/05/19 Python
Python偏函数Partial function使用方法实例详解
2020/06/17 Python
CSS3的transition和animation的用法实例介绍
2014/08/20 HTML / CSS
工商治理实习生的自我评价
2014/01/15 职场文书
我有一个梦想演讲稿
2014/05/05 职场文书
党员教师个人对照检查材料范文
2014/09/25 职场文书
亮剑精神观后感
2015/06/05 职场文书
贫困证明怎么写
2015/06/16 职场文书
如何在CocosCreator里画个炫酷的雷达图
2021/04/16 Javascript
JavaScript的Set数据结构详解
2022/02/18 Javascript
Python字符串的转义字符
2022/04/07 Python