from scipy.optimize import curve_fit
import matplotlib.pyplot as plt
import numpy as np
def func_exp_fb_1(x, a,b,d): # 自定义函数
# return a*np.exp(-b*np.log(x+15)**0.4)+d
# return a*(1/x)+1/(b*np.log(x**0.2+10))+d
# return a*np.exp(-b*x*x*x**0.2)+d +a*(1/x)+1/(b*np.log(x+50))+d
# return a*np.exp(-b*x*x*x**0.2)+d
# return -a*(1/((x+5)**2.2+b+d))
# return a*np.exp(-b*x*x*x**0.2)+d +a*(1/x)+1/(b*np.log(x+50))+d +(-a*(1/((x+5)**2.2+b+d)) + a * (x ** b)-40)/2
# return a * (x ** b)-40
# return (-a*(1/((x+5)**2.2+b+d)) + a * (x ** b)-40)/2
# return a*np.exp(-b*x**0.2) +a*(1/x)+1/(b*np.log(x+1))
# return a*np.exp(-b*x*x*x**0.2)+d +a*(1/x)+1/(b*np.log(x+50))+d -a*(1/((x+5)**2.2+b+d))
# return a*np.exp(-b*x*x*x**0.2)+d +a*(1/x)+1/(b*np.log(x+1))+d
# return a*np.exp(-b*x**0.2) +a*(1/x)+1/(b*np.log(x-9.5))
return a*np.exp(-b*x**0.8)
def tran_num(data):
## 换序号
y = data
x = list(range(len(y)))
return x,y
data = [5,4,3,2,1,0.8] # 你的数据
plt.figure(figsize=(12,6))
x,y = tran_num(data)
x = np.linspace(1,len(y),len(y),dtype=np.int) #这个是从1开始的
popt,pcov = curve_fit(func_exp_fb_1,x,y,maxfev=5000000)
a,b,d = popt
print(a,b,d)
y_pre1 = func_exp_fb_1(x,a,b,d)
x_future = np.linspace(x.shape[0],60,60-x.shape[0]+1) #x的长度设为了60,可以根据自己的情况修改
y_future = func_exp_fb_1(x_future,a,b,d)
plt.plot(x, y, 'ko', label="Original Data")
plt.plot(x, y_pre1, 'r-', label=" Fitting Curve")
plt.plot(x_future, y_future, 'green', label=" prediction Curve")
plt.legend()
## 换序号和时间
plt.title("plot")
plt.xlabel("x")
plt.ylabel("y")
plt.ylim(0,max(data)+1)
plt.show()
python 根据点以及自定义函数预测
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- Mapbox GL JS MapboxGLJS是一个JavaScript库,它使用WebGL呈现交互式地图vect...
- R语言使用plot函数可视化数据、使用type参数自定义设置可视化的类型(数据点和线关系的类型)、设置type参数...
- 最近项目中有需要应用到高德地图的模块,在参考别的app地图相关模块实现过程中,自己产生了一些想法。首先说明要实现的...
- 前言 写这篇文章是因为自己做小程序的时候遇到一个需求,根据登录的角色来展示不同的tabbar,中间遇到了tabba...
- 转自http://www.h3blog.com/article/python-autobuild-image/[h...