博客
关于我
python数据分析之曲线拟合(3):指数函数拟合
阅读量:172 次
发布时间:2019-02-28

本文共 796 字,大约阅读时间需要 2 分钟。

引入

指数形式的曲线也是工程实践中经常遇到的。比如指数衰减。

处理流程

  1. 获取实验数据x, y
  2. 利用scipy.optimize.curve_fit()进行指数函数拟合。
    curve_fit本质是提供一个目标函数和初值,通过优化算法去搜索出最佳的拟合参数。可以提供一个初值,使得拟合更快更准。
  3. 得到拟合出的系数,进行后续的数据处理。

实例

已知一组类似指数衰减数据,形如: y = a 0 e x / a 1 + a 2 y = a_0e^{x/a_1}+a_2 y=a0ex/a1+a2,需拟合出系数。

import numpy as npimport matplotlib.pyplot as pltimport scipy.optimize as optimizepi = np.pi# 模拟生成一组实验数据x = np.arange(0, 100, 0.2)y = np.exp(-x / 51.3)noise = np.random.uniform(0, 0.1, len(x))y += noisefig, ax = plt.subplots()ax.plot(x, y, 'b--')# 拟合指数曲线def target_func(x, a0, a1, a2):    return a0 * np.exp(-x / a1) + a2a0 = max(y) - min(y)a1 = x[round(len(x) / 2)]a2 = min(y)p0 = [a0, a1, a2]print(p0)para, cov = optimize.curve_fit(target_func, x, y, p0=p0)print(para)y_fit = [target_func(a, *para) for a in x]ax.plot(x, y_fit, 'g')plt.show()

在这里插入图片描述

转载地址:http://tthd.baihongyu.com/

你可能感兴趣的文章
Netstat端口占用情况
查看>>
Netty 4的内存管理:sun.misc.Unsafe
查看>>
Netty channelRegistered\ChannelActive---源码分析
查看>>
Netty NIO transport && OIO transport
查看>>
netty php,netty
查看>>
Netty WebSocket客户端
查看>>
netty 主要组件+黏包半包+rpc框架+源码透析
查看>>
Vue过渡 & 动画---vue工作笔记0014
查看>>
Netty 异步任务调度与异步线程池
查看>>
Netty 的 Handler 链调用机制
查看>>
Netty 编解码器详解
查看>>
Netty 解决TCP粘包/半包使用
查看>>
Netty 调用,效率这么低还用啥?
查看>>
Netty+Protostuff实现单机压测秒级接收35万个对象实践经验分享
查看>>
Netty+SpringBoot+FastDFS+Html5实现聊天App详解(一)
查看>>
netty--helloword程序
查看>>
Netty5.x 和3.x、4.x的区别及注意事项(官方翻译)
查看>>
netty——bytebuf的创建、内存分配与池化、组成、扩容规则、写入读取、内存回收、零拷贝
查看>>
netty——Channl的常用方法、ChannelFuture、CloseFuture
查看>>
netty——Future和Promise的使用 线程间的通信
查看>>