第四章-课后练习5-修正指数曲线模型--Python实现

介绍:

在统计预测与决策中,学习了时间序列的趋势外推法,其中修正指数曲线模型的用处广泛。

散点图:

我们通过数据观察,通过python汇制散点图,代码如下:

复制代码
#修正指数曲线预测
import pandas as pd
import matplotlib.pyplot as plt
#导入数据
data = pd.read_excel('P112-5.xlsx',sheet_name='data')
data.head()
#绘制数据图
#显示中文标签
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False

plt.figure(figsize=(10,6))
plt.plot(data, label='data', color='b', marker='o')
plt.title('销售量关于时间的拟合图数据')
plt.xlabel('Year')
plt.xticks(data.index)
#设置y轴范围
plt.ylim(45000, 60000)
plt.show()

修正指数曲线模型:

参数求解:

复制代码
#修正指数曲线函数
def func(t, a, b, c):
    return a+ b*c**t

#拟合
#使用分组法估计参数

#分组
n = 3
#读取前3个数据
data1 = data.iloc[:3]
data2 = data.iloc[3:6]
data3 = data.iloc[6:]

#分组求和
data1_sum = data1['销售量'].sum()
data2_sum = data2['销售量'].sum()
data3_sum = data3['销售量'].sum()

#计算参数
c = ((data3_sum-data2_sum)/(data2_sum-data1_sum))**(1/n)
b = (data2_sum-data1_sum)*((c-1)/(c**n-1)**2)
a = (data1_sum-b*(c**n-1)/(c-1))*(1/n)

#输出参数
print('a = %f, b = %f, c = %f' % (a,b,c))

#输出拟合函数
print('拟合函数为:y = %f + %f * %f^t' % (a,b,c))

data3 = pd.read_excel('P112-5.xlsx',sheet_name='data')
x = data3.销售量

#计算第2022年的销售量
d9 = a + b * c**9
print('2022年的销售量为:%d' % d9)
#绘图
plt.figure(figsize=(10,6))
plt.plot(x, label='data', color='b', marker='o')
#添加点
plt.scatter(9, d9, color='r')
plt.plot(data.index, func(data.index,a,b,c), 'r', label='fit')
#拟合直线过添加的点
plt.scatter(9, d9, color='r')
plt.text(9, d9, '(%d, %d)' % (9, d9), ha='left', va='bottom')
plt.ylim(45000, 60000)
plt.title('销售量关于时间的拟合图')
plt.xlabel('Year')
plt.legend()
plt.show()

求解参数结果:

分组求和法

通过python计算得出:

拟合曲线:

相关推荐
董先生_ad986ad2 小时前
C# 中的 `lock` 关键字本质
开发语言·c#
元亓亓亓2 小时前
Java后端开发day36--源码解析:HashMap
java·开发语言·数据结构
道剑剑非道2 小时前
QT 打包安装程序【windeployqt.exe】报错c000007d原因:Conda巨坑
开发语言·qt·conda
小邓儿◑.◑3 小时前
C++武功秘籍 | 入门知识点
开发语言·c++
码银5 小时前
Java 集合:泛型、Set 集合及其实现类详解
java·开发语言
大G哥5 小时前
PHP标签+注释+html混写+变量
android·开发语言·前端·html·php
傻啦嘿哟5 小时前
HTTP代理基础:网络新手的入门指南
开发语言·php
fish_study_csdn5 小时前
pytest 技术总结
开发语言·python·pytest
咖啡调调。5 小时前
使用Django框架表单
后端·python·django
BO_S__5 小时前
python调用ffmpeg对截取视频片段,可批量处理
python·ffmpeg·音视频