首页 > 解决方案 > 如何在python中找到对数(非线性)模型的参数?

问题描述

所以我需要提出模型参数,以便能够使用模型进行预测(系数等)。对不起,如果我的问题听起来很无聊,但我正在学习 :) 使用的整体模型应该是这样的:y=p1*np.log(x)+p2。所以问题是 p1 & p2 是否由 popt[0] & popt 1表示?log(x) - 只是日期数(数据集从 2010-08-16 运行到 2021-04-14)?y -price 应根据 p1 & p2 和 x 进行预测。

import pandas as pd
import numpy as np
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt

df = pd.read_csv('C:/Data_crypto/Log regression_bands/BCHAIN-MKPRU.csv')
df = df[ df['Value'] > 0]
df = df.iloc[::-1]
df['Date'] = pd.to_datetime(df['Date'])
df

在此处输入图像描述

def func(x, p1,p2):
    return p1*np.log(x)+p2
ydata = np.log(df['Value'])
xdata = np.array([x + 1 for x in range(len(df))])
popt, pcov = curve_fit(func, xdata, ydata)
print(popt)

print(pcov)

popt[0]

popt[1]

xdata

## y axis will be loged and x axis wont be

plt.semilogy(df["Date"], df['Value'])
plt.show()

## plot fitted data and expinatiate data

fittedydata = func(xdata, popt[0], popt[1])


popt[0]

plt.semilogy(df['Date'], df['Value'])
plt.plot(df['Date'], np.exp(fittedydata))
plt.ylim(bottom = 1)

plt.show()

在此处输入图像描述

标签: pythonregressionbitcoinnon-linear

解决方案


推荐阅读