Python应用实现双指数函数及拟合代码实例
作者:kurrrr 时间:2023-04-14 18:17:14
双指数函数
待拟合曲线为 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" 对该函数拟合结果影响微乎其微。
来源:https://www.cnblogs.com/kurrrr/p/13127011.html
标签:Python,应用,双指数,函数,拟合
0
投稿
猜你喜欢
Python3中小括号()、中括号[]、花括号{}的区别详解
2022-05-31 09:20:36
基于Python实现简单的定时器详解
2021-04-12 01:21:28
Bootstrap Paginator分页插件使用方法详解
2023-09-05 23:08:42
Python中表达式x += y和x = x+y 的区别详解
2022-02-24 04:11:19
Golang自定义结构体转map的操作
2024-05-08 10:21:39
MySQL中SQL Mode的查看与设置详解
2024-01-26 12:14:08
Golang 中的json.Marshal问题总结(推荐)
2024-02-16 04:20:22
Mysql Explain 详解
2010-12-03 16:09:00
一文彻底理解JS回调函数
2024-04-10 10:52:10
利用python画一颗心的方法示例
2021-06-04 01:14:33
90行Python代码开发个人云盘应用
2021-12-17 13:44:12
python 利用zmail库发送邮件
2023-02-25 11:06:03
python2.7 mayavi 安装图文教程(推荐)
2022-05-12 05:39:13
通过python顺序修改文件名字的方法
2023-05-20 08:05:29
Asp WinHttp.WinHttpRequest.5.1 对象使用详解
2012-05-02 10:15:27
Python使用matplotlib和pandas实现的画图操作【经典示例】
2023-03-24 16:07:16
Golang使用channel实现一个优雅退出功能
2024-04-28 09:14:28
python 解决print数组/矩阵无法完整输出的问题
2021-05-13 00:45:15
如何尽快释放掉Connection对象建立的连接?
2009-12-16 18:38:00
解析SQL server与asp 互操作的时间处理
2009-02-05 16:13:00