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操作日期和时间的方法
Mar 11 Python
Python编程实现双链表,栈,队列及二叉树的方法示例
Nov 01 Python
儿童python练习实例
May 27 Python
python selenium 获取标签的属性值、内容、状态方法
Jun 22 Python
python实现自动解数独小程序
Jan 21 Python
python使用pygame模块实现坦克大战游戏
Mar 25 Python
Python pandas.DataFrame调整列顺序及修改index名的方法
Jun 21 Python
Pycharm使用远程linux服务器conda/python环境在本地运行的方法(图解))
Dec 09 Python
python如何将两张图片生成为全景图片
Mar 05 Python
详解Python中的编码问题(encoding与decode、str与bytes)
Sep 30 Python
python 实现图片批量压缩的示例
Dec 18 Python
详解使用python爬取抖音app视频(appium可以操控手机)
Jan 26 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关于IE下的iframe跨域导致session丢失问题解决方法
2013/10/10 PHP
微信公众平台接口开发入门示例
2014/12/24 PHP
php使用gettimeofday函数返回当前时间并存放在关联数组里
2015/03/19 PHP
php使用ob_flush不能每隔一秒输出原理分析
2015/06/02 PHP
PHP中类的继承和用法实例分析
2016/05/24 PHP
php中10个不同等级压缩优化图片操作示例
2016/11/14 PHP
幻宇的层模拟窗口效果-提供演示和下载
2007/01/20 Javascript
上传的js验证(图片/文件的扩展名)
2013/04/25 Javascript
JS设置获取cookies的方法
2014/01/26 Javascript
深入理解JavaScript系列(40):设计模式之组合模式详解
2015/03/04 Javascript
jquery动态增加删减表格行特效
2015/11/20 Javascript
浅谈json取值(对象和数组)
2016/06/24 Javascript
JavaScript实现Java中Map容器的方法
2016/10/09 Javascript
javascript按顺序加载运行js方法
2017/12/01 Javascript
Vue cli构建及项目打包以及出现的问题解决
2018/08/27 Javascript
vue中选中多个选项并且改变选中的样式的实例代码
2020/09/16 Javascript
JS实现纸牌发牌动画
2021/01/19 Javascript
[02:28]DOTA2 2015国际邀请赛中国区预选赛首日现场百态
2015/05/26 DOTA
python实现将excel文件转化成CSV格式
2018/03/22 Python
Python装饰器原理与用法分析
2018/04/30 Python
基于Django与ajax之间的json传输方法
2018/05/29 Python
python 接收处理外带的参数方法
2018/12/03 Python
Python进阶之@property动态属性的实现
2019/04/01 Python
python selenium 执行完毕关闭chromedriver进程示例
2019/11/15 Python
将tensorflow模型打包成PB文件及PB文件读取方式
2020/01/23 Python
Django多层嵌套ManyToMany字段ORM操作详解
2020/05/19 Python
Python 如何实现访问者模式
2020/07/28 Python
澳大利亚在线家具、灯饰和家居装饰店:LivingStyles
2018/11/20 全球购物
捷克移动配件网上商店:ProMobily.cz
2019/03/15 全球购物
美国亚马逊旗下男装网站:East Dane(支持中文)
2019/09/25 全球购物
施工安全协议书
2013/12/11 职场文书
总裁助理岗位职责
2014/02/17 职场文书
八一建军节感言
2014/02/28 职场文书
2014年优秀班主任工作总结
2014/12/16 职场文书
Mysql分库分表之后主键处理的几种方法
2022/02/15 MySQL
vue3.0 数字翻牌组件的使用方法详解
2022/04/20 Vue.js