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中的正则表达式的用法
Apr 09 Python
详解Python中的相对导入和绝对导入
Jan 06 Python
Python队列的定义与使用方法示例
Jun 24 Python
Python多进程入门、分布式进程数据共享实例详解
Jun 03 Python
python操作excel让工作自动化
Aug 09 Python
PyQt5+Pycharm安装和配置图文教程详解
Mar 24 Python
基于打开pycharm有带图片md文件卡死问题的解决
Apr 24 Python
tensorflow dataset.shuffle、dataset.batch、dataset.repeat顺序区别详解
Jun 03 Python
Python基于Webhook实现github自动化部署
Nov 28 Python
python不同版本的_new_不同点总结
Dec 09 Python
Python实现简单的猜单词
Jun 15 Python
浅谈Python3中datetime不同时区转换介绍与踩坑
Aug 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
PHP 5.3新特性命名空间规则解析及高级功能
2010/03/11 PHP
PHP函数之error_reporting(E_ALL ^ E_NOTICE)详细说明
2011/07/01 PHP
PHP类与对象中的private访问控制的疑问
2012/11/01 PHP
php中防止SQL注入的最佳解决方法
2013/04/25 PHP
PHP判断图片格式的七种方法小结
2013/06/03 PHP
PHP中使用localhost连接Mysql不成功的解决方法
2014/08/20 PHP
laravel 5 实现模板主题功能
2015/03/02 PHP
为Plesk PHP7启用Oracle OCI8扩展方法总结
2019/03/29 PHP
JavaScript 闭包在封装函数时的简单分析
2009/11/28 Javascript
jQuery.extend 函数详解
2012/02/03 Javascript
JavaScript中统计Textarea字数并提示还能输入的字符
2014/06/10 Javascript
jQuery Masonry瀑布流插件使用详解
2014/11/17 Javascript
angularJs关于指令的一些冷门属性详解
2016/10/24 Javascript
vue cli使用绝对路径引用图片问题的解决
2017/12/06 Javascript
React 高阶组件入门介绍
2018/01/11 Javascript
解决linux下node.js全局模块找不到的问题
2018/05/15 Javascript
Vue.js项目中管理每个页面的头部标签的两种方法
2018/06/25 Javascript
JavaScript实现点击自制菜单效果
2021/02/02 Javascript
JavaScript WeakMap使用详解
2021/02/05 Javascript
[01:09]DOTA2次级职业联赛 - 99战队宣传片
2014/12/01 DOTA
pymssql数据库操作MSSQL2005实例分析
2015/05/25 Python
Python脚本简单实现打开默认浏览器登录人人和打开QQ的方法
2016/04/12 Python
python编写朴素贝叶斯用于文本分类
2017/12/21 Python
python2 与python3的print区别小结
2018/01/16 Python
浅谈pandas筛选出表中满足另一个表所有条件的数据方法
2019/02/08 Python
Python面向对象总结及类与正则表达式详解
2019/04/18 Python
用python写测试数据文件过程解析
2019/09/25 Python
python多进程 主进程和子进程间共享和不共享全局变量实例
2020/04/25 Python
html5 canvas合成海报所遇问题及解决方案总结
2017/08/03 HTML / CSS
教师读书活动总结
2014/05/07 职场文书
四风自我剖析材料思想汇报
2014/10/01 职场文书
图书借阅制度范本
2015/08/06 职场文书
祝福语集锦:朋友新店开业祝福语
2019/12/10 职场文书
Centos系统通过Docker安装并搭建MongoDB数据库
2022/04/12 MongoDB
Python Django / Flask如何使用Elasticsearch
2022/04/19 Python
go goth封装第三方认证库示例详解
2022/08/14 Golang