网站建设中网站图片如何修改,扁平化设计风格的网站模板,an网站建设,网站建设公司平台咨询电话一、问题的描述
问题及要求
1、搜集各个级别世界女子举重比赛的实际数据。分别建立女子举重比赛总成绩的线性模型、幂函数模型、幂函数改进模型#xff0c;并最终建立总冠军评选模型。
应用以上模型对最近举行的一届奥运会女子举重比赛总成绩进行排名#xff0c;并对模型及…
一、问题的描述
问题及要求
1、搜集各个级别世界女子举重比赛的实际数据。分别建立女子举重比赛总成绩的线性模型、幂函数模型、幂函数改进模型并最终建立总冠军评选模型。
应用以上模型对最近举行的一届奥运会女子举重比赛总成绩进行排名并对模型及结果进行必要的分析。
2、请对以上模型作进一步的改进或提出更好的模型并应用模型进行排名。
3、同样对最近举行的一届奥运会所有男女举重冠军行进评选选出1名最优秀的举重运动员给出你的定量评选方法并合理论述。
4、由以上研究过程及结果结论的启发你认为这个模型可以推广应用到生活或学习中的哪些情况下或者应该对生活或学习中哪些事物现有的处理方式需要做出改进请给出你的分析。 说明完成一篇建模论文。具体资料及提交要求见说明文件
二、搜集资料和数据
1.级别的划分
以下括号内容均为网上的查阅资料
男子举重原有8个级别是56、62、69、77、85、94、105公斤级和105公斤以上级于2018年由国际举联调整新的10个级别是55非奥、61、67、73、81、89非奥、96、102非奥、109公斤级和109公斤以上级共10个级别。女子举重原来7个级别是48、53、58、63、69、75公斤级和75公斤级以上级新的10个级别也调整为45非奥、49、55、59、64、71非奥、76、81非奥、87公斤级和87公斤以上级。男女6个非奥级别只在奥运会以外的赛事举行。
在本题目中我们按照奥运会的规定来进行求解规定为8个级别在下面做的时候用的是10个级别
男运动员61、67、73、81、96、109、109以上
女运动员49、55、59、64、76、87、87以上
2.搜集数据
数据来源女子举重世界纪录_国家体育总局 (sport.gov.cn)旧纪录不可用这是一个注意点一定要注意要用2018以后的数据新标准之下的记录 图1.旧标准的记录
可用数据网站需要科学上网https://zh.wikipedia.org/zh-hans/%E8%88%89%E9%87%8D%E4%B8%96%E7%95%8C%E7%B4%80%E9%8C%84%E5%88%97%E8%A1%A8#%E5%A5%B3%E5%AD%90 在这个网站里面我们需要选取当前记录这一栏目千万千万千万记住不要用旧纪录这一栏因为在2018年之后男子和女子的分类的级别已经改变。 图2.男子新记录 图3.女子新纪录
三、第一问
1.分别绘制抓举、挺举、总和下的散点图如下 图4.体重和举重重量的散点趋势图 分析图形可能为线性关系、非线性关系指数、对数等。
散点图代码如下python
#引入相应的库
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns#读取数据
datapd.read_excel(data.xlsx,sheet_name女子)
#选取抓举、挺举、总和
data_1data[data[class]抓举].sort_values(weight)
data_2data[data[class]挺举].sort_values(weight)
data_3data[data[class]总和].sort_values(weight)
#筛选出自变量和因变量
x1data_1[weight].values
y1data_1[record].values
x2data_2[weight].values
y2data_2[record].values
x3data_3[weight].values
y3data_3[record].values#解决中文乱码和符号显示的问题
plt.rcParams[font.sans-serif][SimHei]
plt.rcParams[font.sans-serif][Microsoft YaHei]
#绘制散点图
plt.figure(figsize(10,8),dpi300)
plt.title(体重和举重重量散点趋势图,fontsize16)
plt.xlabel(体重/kg,fontsize14)
plt.ylabel(举重重量/kg,fontsize14)
plt.scatter(x1,y1,colorred,label抓举)
plt.scatter(x2,y2,colorblue,label挺举)
plt.scatter(x3,y3,colorgreen,label总和)
plt.grid(True)
plt.legend()
plt.savefig(1.png)
2.线性拟合
拟合图如下 计算结果 从上面的结果中我们可以看出当多项式的最高次数大于4次的时后虽然模型的R的平方较好但是不符合实际情况属于过拟合现象上面的参数RMSE和MAE来计算误差的话无法排除离群点的影响所以采用了MAPE来计算误差的话这样就排除离群点对模型的影响。这个多项式拟合可以用到第二问当中去作为模型的改进。
2.幂函数拟合
相关的理论支持 上图中的参数是举例子
拟合效果图如下 模型的评估 上面的参数RMSE和MAE来计算误差的话无法排除离群点的影响所以采用了MAPE来计算误差的话这样就排除离群点对模型的影响。
#幂函数进行拟合#引入相应的库
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit# 自定义的拟合函数
def func2(x,k):return k*np.power(x,2/3)# 开始拟合
popt, pcov curve_fit(func2, x3, y3)# 提取拟合结果
k popt[0]#幂函数的系数
print(幂函数的系数为{:}\n.format(k))
print(y3的预测结果为)
y_predictfunc2(x3,k)#预测的结果值
print(y_predict)x_new np.linspace(0, 100, 1000) # 生成新的x轴数据
#绘制拟合曲线效果
plt.figure(figsize(10,8),dpi300)
plt.title(幂函数的拟合效果图,fontsize16)
plt.xlabel(体重/kg,fontsize14)
plt.ylabel(运动员所举的重量,fontsize14)
plt.plot(x3,y3,ro,label真实值)
plt.plot(x_new,func2(x_new,k),colorblue,label预测值)
plt.grid(True)
plt.legend(locupper left)
plt.savefig(3.png)#误差分析
rmse_poly np.sqrt(np.mean((y3 - func2(x3,k))**2))
mae_poly np.mean(np.abs(y3 - func2(x3,k)))
mape_polynp.mean(np.abs(y3-func2(x3,k))/func2(x3,k))# 打印RMSE和MAE
print(幂函数拟合的RMSE和MAE分别为{:.2f}{:.2f}.format(rmse_poly, mae_poly))
print(幂函数拟合的MAPE分别为{:.2f}.format(mape_poly))
y_poly poly(x_new)
# 计算R平方
r2_poly r2_score(y3, func2(x3,k))
# 打印R平方
print(幂函数拟合的R平方为{:.2f}.format(r2_poly))
print(\n)
3.改进幂函数的拟合
相关的理论支持 上图中的参数是举例子 改进幂函数的拟合效果图 误差分析 上面的参数RMSE和MAE来计算误差的话无法排除离群点的影响所以采用了MAPE来计算误差的话这样就排除离群点对模型的影响。
#改进幂函数
#幂函数进行拟合#引入相应的库
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit# 自定义的拟合函数
def func2(x,k,w0,l):return k*np.power(x-w0,l)# 开始拟合
popt, pcov curve_fit(func2, x3, y3)# 提取拟合结果
k popt[0]#改进幂函数的系数k
w0popt[1]#改进幂函数的w0
lpopt[2]#改进幂函数的l
print(改进幂函数的系数为{:}\n.format(k))
print(改进幂函数的w0为{:}\n.format(w0))
print(改进幂函数的l为{:}\n.format(l))
print(y3的预测结果为)
y_predictfunc2(x3,k,w0,l)#预测的结果值
print(y_predict)x_new np.linspace(0, 100, 1000) # 生成新的x轴数据
#绘制拟合曲线效果
plt.figure(figsize(10,8),dpi300)
plt.title(改进幂函数的拟合效果图,fontsize16)
plt.xlabel(体重/kg,fontsize14)
plt.ylabel(运动员所举的重量,fontsize14)
plt.plot(x3,y3,ro,label真实值)
plt.plot(x_new,func2(x_new,k,w0,l),colorblue,label预测值)
plt.grid(True)
plt.legend(locupper left)
plt.savefig(4.png)#误差分析
rmse_poly np.sqrt(np.mean((y3 - func2(x3,k,w0,l))**2))
mae_poly np.mean(np.abs(y3 - func2(x3,k,w0,l)))
mape_polynp.mean(np.abs(y3-func2(x3,k,w0,l))/func2(x3,k,w0,l))# 打印RMSE和MAE
print(改进幂函数拟合的RMSE和MAE分别为{:.2f}{:.2f}.format(rmse_poly, mae_poly))
print(改进幂函数拟合的MAPE分别为{:.2f}.format(mape_poly))
y_poly poly(x_new)
# 计算R平方
r2_poly r2_score(y3, func2(x3,k,w0,l))
# 打印R平方
print(改进幂函数拟合的R平方为{:.2f}.format(r2_poly))
print(\n) 综上所述线性模型一次、幂函数模型、改进幂函数模型的效果大致相同因此可以选择其中一个模型来进行下面的操作为保持与二题的连贯性建议选择线性模型其次是改进幂函数模型最后为幂函数模型。
4.标准统一化
相关理论基础 本文选取中间级别71公斤的作为标准使所以级别的数据转换到71公斤的标准下。
下面是总和成绩记录的所有数据 选用线性模型一次当然也可以选择其它模型改进幂函数等 即y2.52x88.26 将级别6的数据代入即可算出u267.18
本文选取东京奥运会的举重数据进行计算数据网站如下
百度百科-验证 数据处理在当年东京奥运会有的级别中选取总冠军以该级别的最高边界值作为选手的体重。
下表为东京奥运会上各个级别的女子总冠军。 计算得出得结果如下 可视化得分结果如下 从上图可知中国台北的郭婞淳为总冠军
代码
#导库
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
#读取数据
datapd.read_excel(data.xlsx,sheet_name东京奥运会女子举重)
x4data[weight].values
y4data[power].values
#定义标准化函数
coefficients np.polyfit(x3, y3, 1) # 用1次多项式进行拟合
poly np.poly1d(coefficients) # 生成多项式函数
def zfunc(x,y):return 267.18*y/(2.52*x88.26)#求解运动员标准得分该分数越大证明运动员成绩越好
zzfunc(x4,y4)
#将数据读入到data中
data[标准化得分]z#画图
x5data[name].values
y5z
plt.figure(figsize(15,8),dpi300)
plt.title(2021年东京奥运会女子举重统一标准下得分图,fontsize16)
plt.xlabel(运动员姓名,fontsize14)
plt.ylabel(标准得分,fontsize14)
plt.plot(x5,y5,color#FF8066,markero,label标准得分)
plt.grid(True)
plt.legend(locupper left)
plt.savefig(5.png)
print(得分结果为)
print(z)
三、第二问多项式拟合 from scipy.interpolate import interp1d
import numpy as np
import matplotlib.pyplot as plt
from sklearn.metrics import r2_scoreplt.rcParams[font.sans-serif][SimHei] #用来正常显示中文标签
plt.rcParams[axes.unicode_minus] False #用来正常显示负号
plt.figure(figsize(20,15),dpi300)
plt.title(总成绩拟合图,fontsize16)
plt.xlabel(体重,fontsize14)
plt.ylabel(举重重量,fontsize14)
for i in range(1,9):x_new np.linspace(40, 100, 1000) # 生成新的x轴数据coefficients np.polyfit(x3, y3, i) # 用5次多项式进行拟合poly np.poly1d(coefficients) # 生成多项式函数# 绘图plt.subplot(2,4,i)plt.subplots_adjust(leftNone, bottom0.1, rightNone, topNone,wspaceNone, hspaceNone)plt.grid(True)plt.scatter(x3,y3,markero,colorred, label原始数据)plt.plot(x_new, poly(x_new), labelstr(i)次多项式拟合)plt.legend(locupper left)
plt.savefig(2.png)
# 计算拟合后的数值
for i in range(1,9):coefficients np.polyfit(x3, y3, i) # 用5次多项式进行拟合poly np.poly1d(coefficients) # 生成多项式函数y_poly poly(x_new)# # 计算RMSE和MAE# rmse_cubic np.sqrt(np.mean((y - f_cubic(x))**2))# mae_cubic np.mean(np.abs(y - f_cubic(x)))rmse_poly np.sqrt(np.mean((y3 - poly(x3))**2))mae_poly np.mean(np.abs(y3 - poly(x3)))mape_polynp.mean(np.abs(y3-poly(x3))/poly(x3))# 打印RMSE和MAEprint(str(i)次多项式拟合的RMSE和MAE分别为{:.2f}{:.2f}.format(rmse_poly, mae_poly))print(str(i)次多项式拟合的MAPE分别为{:.2f}.format(mape_poly))print(poly)y_poly poly(x_new)# 计算R平方r2_poly r2_score(y3, poly(x3))# 打印R平方print(str(i)次多项式拟合的R平方为{:.2f}.format(r2_poly))print(\n) 从上面的结果中我们可以看出当多项式的最高次数大于4次的时后虽然模型的R的平方较好但是不符合实际情况属于过拟合现象上面的参数RMSE和MAE来计算误差的话无法排除离群点的影响所以采用了MAPE来计算误差的话这样就排除离群点对模型的影响。这个多项式拟合可以用到第二问当中去作为模型的改进因此我们进行的优化为多项式优化。
可能大家会有疑问为什么要用多项式拟合而不去用其它函数呢在这里就是用到泰勒公式的简单应用用多项式近似地表示函数的公式称为泰勒公式并且根据余项表达式的不同而有不同的形式。 得名于英国数学家布鲁克·泰勒他在1712年的一封信里首次叙述了这个公式。嘿嘿这个大家大一的高数应该都是学过滴
关于相关数学理论请看这篇文章其实可以忽略想了解的可以看下记住结论就行不用难为自己
泰勒公式简单应用多项式近似表示任意函数 - 知乎 (zhihu.com)
基于4次多项式的改进模型 从图中可以清晰的看出在4次多项式的条件下候志惠的标准的分最高。
为什么和前面的不一样呢?大家可能有疑问但是在下面可以做一个实验当多项式的最高次大于等于1时模型所评选的女子冠军均为候志惠这说明多项式的拟合精度更高更加具有稳定性和可信性。 在这就不再一一列举了从上面可以观察出随着最高次数的增加模型逐渐趋于稳定。
代码
#导库
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt#读取数据
datapd.read_excel(data.xlsx,sheet_name东京奥运会女子举重)
x4data[weight].values
y4data[power].values#改进的模型4次多项式拟合
coefficients np.polyfit(x3, y3, 5) # 只需改动这里面参数即可拟合多种多项式
poly np.poly1d(coefficients) # 生成多项式函数
print(poly)
upoly(71)
print(u)#定义改进后的标准得分函数
def zfunction(x,y):return (270.5482219490648*y)/poly(x)
z1zfunction(x4,y4)#画图
x5data[name].values
y5z1
plt.figure(figsize(15,8),dpi300)
plt.title(5次改进02021年东京奥运会女子举重统一标准下得分图,fontsize16)
plt.xlabel(运动员姓名,fontsize14)
plt.ylabel(标准得分,fontsize14)
plt.plot(x5,y5,colorb,markero,label标准得分)
plt.grid(True)
plt.legend(locupper left)
plt.savefig(6.png)
print(得分结果为)
print(z1) 至此解决重要的两问。三问、四问都比较容易解决。
四、第三问选取男女冠军 这一问就比较简单啦在这里我就给大家说一种思路 首先按照上面的1、2问的思路对男子数据进行建模选出男子冠军然后再和我们选出的女子冠军进行比较最后选出总冠军关于比较的方法。 我查阅了相关资料 ‘’科学界对于男性与女性在青春期形成的生理差异已经普遍达成共识。在哈伯德已经经历过的男性青春期期间男性的肌肉量会得到增加。 体育科学家塔克Ross Tucker称男性青春期出现的一系列生理变化会在男性与女性之间带来显著的功能优势他称在游泳与自行车运动上这种差异在10-12%之间而在举重等涉及上肢肌肉力量的项目上可能达到“30-40%”这个结论我用往年的数据检验了基本上是符合的因此可以用了。 在模型假设部分假设同体重下男子的上肢力量要比女子高35%这个比例位于30%-40%就行。这句话一定要带上 即是选出的女子总冠军的标准化成绩乘上10.35再去和男子冠军的标准化成绩进行比较。根据我提供的数据男子标准到73kg、女子标准到71kg由于男子的标准女子大为了弥补这个差值我们可以让女子和男子的上肢力量设置的大一些也就是比0.35大一些这样就弥补了这一缺陷啦。 代码上面都有了直接套用就行。
四、第四问语文建模外加一些机理分析 这一问就更加简单了 这一问可以从摔跤、拳击、赛艇按体重分等级的体育项目来考虑 体测的男女的标准转化也是一个思路 至此全部更新完全本人水平有限如有错误可以提出欢迎交流。
