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求导数的方法
May 09 Python
Python ValueError: invalid literal for int() with base 10 实用解决方法
Jun 21 Python
Windows中安装使用Virtualenv来创建独立Python环境
May 31 Python
python3解析库pyquery的深入讲解
Jun 26 Python
Python 16进制与中文相互转换的实现方法
Jul 09 Python
python 除法保留两位小数点的方法
Jul 16 Python
详解Pandas之容易让人混淆的行选择和列选择
Jul 10 Python
Python Tornado之跨域请求与Options请求方式
Mar 28 Python
浅析python 通⽤爬⾍和聚焦爬⾍
Sep 28 Python
python使用ctypes库调用DLL动态链接库
Oct 22 Python
深入浅析Django MTV模式
Sep 04 Python
PyCharm 配置SSH和SFTP连接远程服务器
May 11 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基础之运算符的使用方法
2013/04/28 PHP
10个php函数实用却不常见
2015/10/13 PHP
PHP版本升级到7.x后wordpress的一些修改及wordpress技巧
2015/12/25 PHP
微信开发之网页授权获取用户信息(二)
2016/01/08 PHP
Laravel 5.1 on SAE环境开发教程【附项目demo源码】
2016/10/09 PHP
PHP实现的简单sha1加密功能示例
2017/08/27 PHP
laravel实现前后台路由分离的方法
2019/10/13 PHP
浅谈tudou土豆网首页图片延迟加载的效果
2010/06/23 Javascript
javascript学习笔记(十八) 获得页面中的元素代码
2012/06/20 Javascript
jquery动态调整div大小使其宽度始终为浏览器宽度
2014/06/06 Javascript
jquery图片滚动放大代码分享(2)
2015/08/28 Javascript
js图片轮播手动切换效果
2015/11/10 Javascript
javascript每日必学之继承
2016/02/23 Javascript
JS上传组件FileUpload自定义模板的使用方法
2016/05/10 Javascript
JavaScript实现复制或剪切内容到剪贴板功能的方法
2016/05/23 Javascript
JS优化与惰性载入函数实例分析
2017/04/06 Javascript
Vue如何引入远程JS文件
2017/04/20 Javascript
微信小程序 request接口的封装实例代码
2017/04/26 Javascript
详解前后端分离之VueJS前端
2017/05/24 Javascript
angular动态删除ng-repaeat添加的dom节点的方法
2017/07/20 Javascript
canvas绘制爱心的几种方法总结(推荐)
2017/10/31 Javascript
vue超时计算的组件实例代码
2018/07/09 Javascript
vue组件数据传递、父子组件数据获取,slot,router路由功能示例
2019/03/19 Javascript
react PropTypes校验传递的值操作示例
2020/04/28 Javascript
Python实现类的创建与使用方法示例
2017/07/25 Python
Python GUI Tkinter简单实现个性签名设计
2018/06/19 Python
Django开发中的日志输出的方法
2018/07/02 Python
Python设计模式之享元模式原理与用法实例分析
2019/01/11 Python
如何利用Pyecharts可视化微信好友
2019/07/04 Python
python pptx复制指定页的ppt教程
2020/02/14 Python
解决pytorch 保存模型遇到的问题
2021/03/03 Python
OPPO手机官方商城:中国手机市场出货量第一品牌
2017/10/18 全球购物
英国和世界各地鲜花速递专家:Arena Flowers
2018/02/10 全球购物
党员学习新党章思想汇报
2014/10/25 职场文书
会议简讯范文
2015/07/20 职场文书
python 命令行传参方法总结
2021/05/25 Python