2022 APMCM亚太数学建模竞赛 C题全球是否变暖思路及代码实现(持续更新中)

news2024/9/27 23:25:56

2022 APMCM亚太数学建模竞赛 C题全球是否变暖思路及代码实现(持续更新中)

1 题目

全球变暖与否? 加拿大49.6°C的高温为地球北纬50°以上地区创造了新的气温记录，一周内就有数百人死于高温;美国加利福尼亚州死亡谷54.4°C，是地球上有记录以来的最高温度;科威特53.5°C，阳光直射下甚至超过70°C，中东很多国家超过50°C。

图1所示。北美的高温。今年以来，我们看到了大量惊人的气温报告。地球正在燃烧的事实是毋庸置疑的。继这些地区从6月底到7月初的可怕高温之后，意大利再次创下欧洲气温纪录，达到惊人的48.8°C，许多国家宣布进入紧急状态。

全球气候变暖是一种与自然有关的现象。正是由于温室效应的不断积累，导致地球大气系统吸收和排放的能量失衡，能量在地球大气系统中不断积累，导致气温上升，全球气候变暖。

在工业革命之前，二氧化碳(CO2大气中的二氧化碳一直保持在280 ppm左右。2004年3月，大气中的二氧化碳浓度达到了377.7 ppm，这是到那时为止10年平均增幅最大的一次。根据美国国家海洋和大气管理局(NOAA)和斯克里普斯海洋学研究所(SIO)的科学家们的研究，月平均CO2 浓度水平在2022年5月达到峰值421 ppm。经济合作与发展组织(OECD)的一份报告预测CO2 到2050年，二氧化碳浓度将达到685 PPM

APMCM组委会已要求您的团队解决这些关于当前报告和未来全球温度水平预测的主张。他们提供了数据集2022_APMCM_C_Data.csv，其中包含239177条记录，以协助你们的研究。

需求 1.你同意全球温度的说法吗?使用附件中的2022_APMCM_C_Data.csv和你们团队收集的其他数据集来分析全球温度变化。

a)你是否同意2022年3月全球温度的上升导致了比以往任何10年期间观测到的更大的增长?为什么?为什么不? b)基于历史数据，请建立两个或两个以上的数学模型来描述过去并预测未来的全球温度水平。 c)分别使用1(b)中的每个模型来预测2050年和2100年的全球温度。你们的任何一个模型是否与2050年或2100年观测点的全球平均温度将达到20.00°C的预测一致?如果不是在2050年或2100年，你们预测模型中观测点的平均温度什么时候会达到20.00°C?

d)你在1(b)中建立的哪个模型你认为最准确?为什么? 2.影响温度变化的原因是什么? a)利用问题1的结果和附件2022_ APMCM_C_ data .csv和你们团队收集的其他数据集中的数据，建立一个数学模型来分析全球温度、时间和位置之间的关系(如果有)，并解释它们之间的关系或证明它们之间没有关系。

b)请收集相关数据，分析自然灾害(如火山爆发、森林火灾和新冠肺炎)的因素。对全球气温有影响吗?

2 题目分析

2.1 问题一

（a）你是否同意2022年3月全球温度的上升导致了比以往任何10年期间观测到的更大的增长?为什么?为什么不?

思路：第一步，挑选几个城市，选取2022年的一整年的数据，将温度变化以折线图可视化，分析3月是不是存在异常值。

第二步，计算2022年排除3月的平均温度，和加上3月后的平均温度，计算2011年-2021年的平均温度，以折线图可视化。对比每年的增长，与2022年有无3月份数据的增长。

（b）基于历史数据，请建立两个或两个以上的数学模型来描述过去并预测未来的全球温度水平。

从（c）知道，需要以年为跨度，建立时间序列预测模型。通过分析虽然有100个城市，但在题目中，没指定说明哪个观测点。

随机挑一个城市即可。以年为时间单位，统计每年的平均温度，得到一个时间序列。选择以下时间序列预测模型进行预测即可。注意题目需要两个模型。

有哪些时间序列预测模型

自回归 (AR)
移动平均线
自回归移动平均线
自回归积分移动平均线 (ARIMA)
季节性自回归积分移动平均线 (SARIMA)
具有外生回归量的季节性自回归综合移动平均线 (SARIMAX)
具有 ARIMA 误差的回归模型
向量自回归 (VAR)
GARCH 模型
Glostan、Jagannathan 和 Runkle GARCH 模型

（c）利用以上模型，预测未来100年的平均温度，并将这100个数据，可视化，分析什么时候能达到20度。

（d）模型的准确性，用评价指标来评价

模型的评价指标有

（1）MSE（mean squared error），是观测值（observed values）与预测值（predicted values）的误差的平方和的均值，即SSE/n。它是误差的二阶矩，包含估计量的方差（variance）及其偏差（bias），是衡量估计量质量的指标，其公式为：

$MSE(y,\hat{y})=1n∑(yi−\hat{yi})2$

（2）RMSE（root mean squared error），也称作RMSD(root mean square deviation)，是MSE的算数平方根。由于每个误差（each error）对RMSD的影响与误差的平方（squared error）成正比，因此较大的误差会对RMSE影响过大，RMSE对异常值很敏感。其公式为：

$RMSE(y,\hat{y})=1n∑(yi−\hat{yi})2$

（3）MAE（mean absolute error），是时间序列分析中预测误差常用的指标，由于MAE使用的是与被测数据相同的尺度（scale）,因此不能用于比较两个不同尺度的序列。MAE又被称为L1范数损失函数（就是可以做为损失函数），是真实数据与预测数据之差的绝对值的均值。

$MAE(y,\hat{y})=1n∑|yi−\hat{yi}|$

（4）MAPE（mean absolute percentage error），也被称为MAPD（mean absolute percentage deviation），是一种衡量预测方法的预测准确性的指标。MAPE在解释相对误差（relative error）方面非常直观，在评价模型时MAPE通常用作回归（regression）问题的损失函数（loss function）。

$MAPE(y,\hat{y})=1n∑|yi−\hat{yi}|y$