Python实现的拟合二元一次函数功能示例【基于scipy模块】


Posted in Python onMay 15, 2018

本文实例讲述了Python实现的拟合二元一次函数功能。分享给大家供大家参考,具体如下:

背景:

使用scipy拟合一元二次函数。

参考:

HYRY Studio-《用Python做科学计算

代码:

# -*- coding:utf-8 -*-
#! python3
import numpy as np
from scipy.optimize import leastsq
import pylab as pl
def func(x,p):
  """
  数组拟合函数
  """
  A,k,theta = p
  return A*(x-k)**2+theta
def residuals(p,y,x):
  """
  残差
  """
  return y-func(x,p)
x = np.linspace(0,2,100)
A,k,theta = 10.,1,2. #真实数据参数
y0 = func(x,[A,k,theta]) #真实数据
y1 = y0 + 2 * np.random.randn(len(x)) #加入噪声序列
p0 = [7.,0.2,1.]
plsq = leastsq(residuals,p0,args = (y1,x))
print("真实参数:",[A,k,theta])
print("拟合参数:",plsq[0]) #试验数据拟合后的参数
pl.plot(x,y0,label = "real")
pl.plot(x,y1,label = "real+noise")
pl.plot(x,func(x,plsq[0]),label = "fitting")
pl.legend()
pl.show()

结果:

(貌似这里的求解方法用了智能算法,因为每次的结果都有细小差异。具体资料没见到,以后有精力再找)

真实参数: [10.0, 1, 2.0]
拟合参数: [ 10.83391995   0.98950039   1.63356065]

Python实现的拟合二元一次函数功能示例【基于scipy模块】

Python 相关文章推荐
Python enumerate函数功能与用法示例
Mar 01 Python
python flask安装和命令详解
Apr 02 Python
pandas数据筛选和csv操作的实现方法
Jul 02 Python
Python (Win)readline和tab补全的安装方法
Aug 27 Python
python 队列基本定义与使用方法【初始化、赋值、判断等】
Oct 24 Python
Python使用psutil获取进程信息的例子
Dec 17 Python
Python urlopen()和urlretrieve()用法解析
Jan 07 Python
python实现扫雷小游戏
Apr 24 Python
总结Pyinstaller的坑及终极解决方法(小结)
Sep 21 Python
python代数式括号有效性检验示例代码
Oct 04 Python
Python的scikit-image模块实例讲解
Dec 30 Python
Python代码风格与编程习惯重要吗?
Jun 03 Python
Python实现的多项式拟合功能示例【基于matplotlib】
May 15 #Python
浅析PHP与Python进行数据交互
May 15 #Python
Python使用min、max函数查找二维数据矩阵中最小、最大值的方法
May 15 #Python
django自带的server 让外网主机访问方法
May 14 #Python
Python中利用xpath解析HTML的方法
May 14 #Python
利用pyinstaller将py文件打包为exe的方法
May 14 #Python
django的登录注册系统的示例代码
May 14 #Python
You might like
用定制的PHP应用程序来获取Web服务器的状态信息
2006/10/09 PHP
PHP源代码数组统计count分析
2011/08/02 PHP
php中的注释、变量、数组、常量、函数应用介绍
2012/11/16 PHP
ubuntu10.04配置 nginx+php-fpm模式的详解
2013/06/03 PHP
PHP图片处理之图片旋转和图片翻转实例
2014/11/19 PHP
Laravel重写用户登录简单示例
2016/10/08 PHP
PHP+JS实现的实时搜索提示功能
2018/03/13 PHP
ThinkPHP中获取指定日期后工作日的具体日期方法
2018/10/14 PHP
PHP如何防止用户重复提交表单
2020/12/09 PHP
通过正则格式化url查询字符串实现代码
2012/12/28 Javascript
JavaScript中的noscript元素属性位置及作用介绍
2013/04/11 Javascript
javascript创建和存储cookie示例
2014/01/07 Javascript
谷歌地图打不开的解决办法
2014/08/07 Javascript
jQuery过滤选择器用法分析
2015/02/10 Javascript
如何用js实现鼠标向上滚动时浮动导航
2016/07/18 Javascript
几句话带你理解JS中的this、闭包、原型链
2016/09/26 Javascript
js指定步长实现单方向匀速运动
2017/07/17 Javascript
详解基于webpack&gettext的前端多语言方案
2019/01/29 Javascript
JS中的算法与数据结构之常见排序(Sort)算法详解
2019/08/16 Javascript
Vue使用NProgress进度条的方法
2019/09/21 Javascript
vue项目引入ts步骤(小结)
2019/10/31 Javascript
Python中列表和元组的相关语句和方法讲解
2015/08/20 Python
详解python中requirements.txt的一切
2017/03/03 Python
Python中第三方库Requests库的高级用法详解
2017/03/12 Python
Python Flask-web表单使用详解
2017/11/18 Python
Flask web开发处理POST请求实现(登录案例)
2018/07/26 Python
NLTK 3.2.4 环境搭建教程
2018/09/19 Python
Liu Jo西班牙官网:意大利服装品牌
2019/09/11 全球购物
中专毕业生的自我鉴定
2013/12/01 职场文书
预备党员政审材料
2014/02/04 职场文书
上课玩手机检讨书
2014/02/08 职场文书
认识深刻的检讨书
2014/02/16 职场文书
大龄毕业生求职别忘职业规划
2014/03/11 职场文书
公司总经理岗位职责
2014/03/15 职场文书
刑事上诉状范文
2015/05/22 职场文书
SQL Server查询某个字段在哪些表中存在
2022/03/03 SQL Server