双指数函数

待拟合曲线为 y(x) = bepx + ceqx

import matplotlib.pyplot as plt

x = ([0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5,

0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0])

y = ([0.33, 0.26, 0.18, 0.16, 0.12, 0.09, 0.08, 0.07, 0.06, 0.06,

0.06, 0.07, 0.09, 0.1, 0.15, 0.19, 0.25, 0.36, 0.47, 0.68])

plt.scatter(x, y)

plt.show()

拟合

import numpy as np

import matplotlib.pyplot as plt

from scipy.optimize import curve_fit

def double_exp(x, b, c, p, q):

x = np.array(x)

return b*np.exp(p*x) + c*np.exp(q*x)

x = ([0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5,

0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0])

y = ([0.33, 0.26, 0.18, 0.16, 0.12, 0.09, 0.08, 0.07, 0.06, 0.06,

0.06, 0.07, 0.09, 0.1, 0.15, 0.19, 0.25, 0.36, 0.47, 0.68])

popt, pcov = curve_fit(double_exp, x, y, [1, 1, 1, 1])

print(popt)

b = popt[0]

c = popt[1]

p = popt[2]

q = popt[3]

y_fit = double_exp(x, b, c, p, q)

plt.scatter(x, y)

plt.plot(x, y_fit, color='red', linewidth=1.0)

plt.show()

numpy 库,实现列表转矩阵,得以进行数学运算。matplotlib.pyplot 库,绘制图像。scipy.optimize 库,curve_fit() 函数,使用非线性最小二乘法拟合曲线。curve_fit()popt,拟合结果,在这里指b, c, p, q 的值。povc,该拟合结果对应的协方差。

拟合结果参数原函数拟合结果误差b0.00110.00110c0.42000.42420.42%p6.39986.49881.55%q-5.1551-5.21641.19%

误差可以满意。

经过测试,如果将初始参数设置为原函数参数(保留 4 位小数),拟合得到的结果并未发生变化。

经过测试,拟合使用的三种方法,"trf","lm" 和 "dogbox" 对该函数拟合结果影响微乎其微。

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持我们。

python函数拟合编程_Python应用实现双指数函数及拟合代码实例相关推荐

  1. python提取pdf文字_python基于pdfminer库提取pdf文字代码实例

    安装pdfminer 库 windows 下安装pdfminer3k pip install pdfminer3k Liunx 下安装pdfminer pip install pdfminer 代码 ...

  2. python监听剪贴板_Python监听剪切板实现方法代码实例

    第一种 import win32clipboard import time #速度快 容易出错 class niubi(): def lihai(self): while True: #jiantin ...

  3. python 大智慧365 数据格式_Python爬取365好书中小说代码实例

    需要转载的小伙伴转载后请注明转载的地址 需要用到的库 from bs4 import BeautifulSoup import requests import time 获取小说的每个章节的名称和章节 ...

  4. python函数参数列表_python函数列表

    广告关闭 腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! python函数函数是python里组织代码的最小单元,python函数包含以下 ...

  5. python使用turtle库绘制一个100长度的十字架_使用python图形模块turtle库绘制樱花、玫瑰、圣诞树代码实例...

    使用python图形模块turtle库绘制樱花.玫瑰.圣诞树代码实例 发布时间:2020-08-22 15:57:19 来源:脚本之家 阅读:76 作者:Soul fragments 今天为大家介绍几 ...

  6. python函数做菜单_PYTHON图形化操作界面的编程七__创建菜单

    PYTHON图形化操作界面的编程七__创建菜单 十八.创建菜单 1.水平菜单的创建 创建菜单需要多条语句,所以这里通过实例来说明水平菜单的创建方法: 下面的语句可以在窗口中添加水平菜单,其中前四行语句 ...

  7. python函数五要素_Python安装及关键要素

    一.Python定义 shell编程: 控制语言:胶水语言 框架:web应用开发 二.Python性能优化工具 Psyco: python语言是一个扩展模块,可以即时对程序代码进行专业的算法优化,可以 ...

  8. python函数五要素_python之基础篇(二)

    防伪码:忘情公子著 一.面向过程与面向对象 面向过程: 以指令为中心,由指令去处理数据 只要考虑如何组织代码去解决问题 面向对象: 以数据为中心,所有的处理代码都围绕数据展开 要考虑如何设计数据结构组 ...

  9. python函数手册 下载_python函数手册

    广告关闭 腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! 这个时候,为了避免反复编写相同的代码,我们可以使用一个函数来对某段代码块进行封装 ...

最新文章

  1. 淘淘商城项目mysql服务器_SpringMVC+Spring+Mybatis+Mysql+Maven+Svn[ 淘淘商城项目环境搭建 ]...
  2. 数据加密辅助类(默认编码UTF-8)
  3. php 字符串 字典序序排序,C++ 怎么实现字典序排序法,自然排序
  4. mso-number-format
  5. 【Kafka】Kafka No serviceName defined in either JAAS or Kafka config
  6. 今天听阿里P8讲完SpringBoot,瞬间感觉自己白学了……
  7. python实现logistic增长模型
  8. Ubuntu中使用dnw工具:没有找到/dev/secbulk0
  9. Mac 应用开发--Cocoa运用程序显示或隐藏MainNemu
  10. c语言算法单循环球队比赛安排,单循环赛赛程安排算法的研究.doc
  11. 高等数学张宇18讲 第三讲 一元函数微分学的概念与计算
  12. 盘点新加坡旅游宣传片中不可错过的景点
  13. python icon生成小工具
  14. latex 去掉(不显示)空白页的页码与页眉
  15. 智能制造,从smart到intelligent
  16. linux里的.mht是什么文件,mht是什么文件?打开mht文件的方法
  17. 泰勒展开式求sinx近似值
  18. linux工作区目录,Linux(以RHEL7为例)下添加工作区的方法|| The Way To Add Workspace On Linux...
  19. vim 安装YouCompleteMe 自动补全插件
  20. User32.dll详细介绍

热门文章

  1. 华强北耳机为啥老是有人翻车?
  2. 9、Flutter 实现 生成二维码
  3. 织梦DEDE网站后台如何上传附件
  4. android使用menu需要重写的方法,Android – 正确使用invalidateOptionsMenu()
  5. centos7 编译安装nginx 设置自启动服务 支持https
  6. r语言rank降序_R语言rank函数详细解析
  7. 数据库实验四 用户权限管理
  8. 什么是框架?spring mvc框架
  9. curl 怎么在xp下使用_Http Post 快速使用
  10. struts2文件上传,下载