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单元测试unittest实例详解
May 11 Python
使用Python3制作TCP端口扫描器
Apr 17 Python
Python对ElasticSearch获取数据及操作
Apr 24 Python
python提取log文件内容并画出图表
Jul 08 Python
django之使用celery-把耗时程序放到celery里面执行的方法
Jul 12 Python
Python多线程模块Threading用法示例小结
Nov 09 Python
python-web根据元素属性进行定位的方法
Dec 13 Python
Python如何使用内置库matplotlib绘制折线图
Feb 24 Python
Django REST framwork的权限验证实例
Apr 02 Python
Python smtp邮件发送模块用法教程
Jun 15 Python
Python requests模块安装及使用教程图解
Jun 30 Python
python 算法题——快乐数的多种解法
May 27 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优化及注意事项
2016/07/04 PHP
php获取ip及网址的简单方法(必看)
2017/04/01 PHP
PHP日期和时间函数的使用示例详解
2020/08/06 PHP
Javascript-Mozilla和IE中的一个函数直接量的问题
2007/01/09 Javascript
JS实现程序暂停与继续功能代码解读
2013/10/10 Javascript
ECMAScript 5中的属性描述符详解
2015/03/02 Javascript
jquery处理checkbox(复选框)是否被选中实例代码
2017/06/12 jQuery
通过nodejs 服务器读取HTML文件渲染到页面的方法
2018/05/17 NodeJs
JS实现把一个页面层数据传递到另一个页面的两种方式
2018/08/13 Javascript
node读写Excel操作实例分析
2019/11/06 Javascript
python基础教程之数字处理(math)模块详解
2014/03/25 Python
python对字典进行排序实例
2014/09/25 Python
Python实现的文本简单可逆加密算法示例
2017/05/18 Python
CentOS 7下Python 2.7升级至Python3.6.1的实战教程
2017/07/06 Python
Python 多线程Threading初学教程
2017/08/22 Python
python中numpy的矩阵、多维数组的用法
2018/02/05 Python
Python对切片命名的实现方法
2018/10/16 Python
Python中三元表达式的几种写法介绍
2019/03/04 Python
一文了解Python并发编程的工程实现方法
2019/05/31 Python
使用Python OpenCV为CNN增加图像样本的实现
2019/06/10 Python
python 搜索大文件的实例代码
2019/07/08 Python
python文件操作seek()偏移量,读取指正到指定位置操作
2020/07/05 Python
记一次django内存异常排查及解决方法
2020/08/07 Python
python 实现图片批量压缩的示例
2020/12/18 Python
HTML5画渐变背景图片并自动下载实现步骤
2013/11/18 HTML / CSS
英国在线电子和小工具商店:TecoBuy
2018/10/06 全球购物
开水果连锁店创业计划书
2013/12/29 职场文书
课改先进个人汇报材料
2014/01/26 职场文书
财务人员求职自荐书范文
2014/02/10 职场文书
电脑售后服务承诺书
2014/03/27 职场文书
保护环境倡议书300字
2014/05/19 职场文书
卖车协议书范本4篇
2014/10/01 职场文书
地道战观后感
2015/06/04 职场文书
创业计划书之烤红薯
2019/09/26 职场文书
golang 实现菜单树的生成方式
2021/04/28 Golang
python面向对象版学生信息管理系统
2021/06/24 Python