一、问题背景
以秦皇岛为例，对煤炭的影响因素分析，了解煤炭价格的变化，并提出合理的政策建议。
二、问题分析
2.1问题一的分析
问题一要求通过量化分析的方法，给出不超过10种的影响煤炭价格的主要因素，同时以秦皇岛港动力煤价格为例，给出从2019年5月1日至2020年4月30日，影响秦皇岛港动力煤价格的主要因素的排序。
2.2问题二的分析
问题二要求结合秦皇岛港动力煤价格的历史数据和影响煤炭价格的主要因素，预测未来31天、35周、36个月的煤炭价格。
2.3问题三的分析
问题三要求综合考虑未来各种情况引起的煤炭价格影响因素在结构性和重要性方面的变化并给出预测结果。
2.4问题四的分析
问题四要求在保障我国未来煤炭市场的平稳发展的前提下，结合问题三的模型，向政府部门提供相关的政策建议。
三、模型假设
1.假设国家对煤炭价格的干预呈现周期性；
四、模型建立与求解
4.1指标的选取
本文选取供给因素、需求因素、政策因素、综合因素作为影响煤炭价格的一级指标。
如图示二级指标的选取
图1 指标
4.1.1排序
将指标进行从排序结果如下表所示：
表1 指标优化度排序
指标 排名 指标
优化度
煤炭库存量 1 需求量
6 政策 2 气候变化
7 煤炭成本 3 原煤产量
8 国际煤炭价格 4 下游煤炭需求
9 排放约束 5
铁路运力
10
4.2价格预测1
利用灰色预测可以得到响应方程为：
(1)
(0)0.0012(1)[(1)415608.33]e 415608.33k x k x ∧+=+-
得到结果 4.2.1价格预测2
多元方程为：
0181,(,1,1)8,(,1,1)t t t t t t t y a a x a x ∧
∧
∧
∧
-+-+=++
+
表2预测结果 周 预测价格 周 预测价格 周
预测价格 5月6日 659.09 7月29日 663.13 10月21日 667.17 5月13日 659.42 8月5日 663.46 10月28日 667.51 5月20日 659.76
8月12日 663.80
11月4日 667.84
6月3日660.438月26日664.4811月18日668.52
6月10日660.779月2日664.8111月25日668.86
6月17日661.119月9日665.1512月2日669.19
6月24日661.449月16日665.4912月9日669.53
7月1日661.789月23日665.8212月16日669.87
7月8日662.129月30日666.1612月23日670.20
7月15日662.4510月7日666.5012月30日670.54
7月22日662.7910月14日666.83
4.3预测3
36月结果如表所示
表3 未预测结果3
月预测价格月预测价格月预测价格2020年5月640.002021年5月642.642022年5月645.55 2020年6月640.212021年6月642.872022年6月645.80 2020年7月640.422021年7月643.102022年7月646.06 2020年8月640.642021年8月643.342022年8月646.32 2020年9月640.852021年9月643.582022年9月646.58 2020年10月641.072021年10月643.822022年10月646.84 2020年11月641.292021年11月644.062022年11月647.10 2020年12月641.512021年12月644.302022年12月647.37 2021年1月641.732022年1月644.552023年1月647.64 2021年2月641.952022年2月644.792023年2月647.91 2021年3月642.182022年3月645.042023年3月648.18 2021年4月642.412022年4月645.292023年4月648.45
4.4价格预测4
不妨考虑易量化的指标，如：煤炭库存量能源消耗、排放约束因素对其进行分析。
5.9建议
需要对煤炭价格的波动提出相关建议，本文从国家层面、市场层面和企业层面并结合问题三中的概率组合预测模型提出建议。
六、模型检验
考虑到实际应用层面的因素，对模型进行检验其结果均符合实际要求。
七、模型评价
7.1.模型优点
思想新颖，贴近实际较为合适。
7.2模型缺点
主观性太强。
附录
Shujv=MyData(1,1);
for j=1:1:3
n=1;
for i=1:1:size(M,1)
for k=0:1:MyData(i,1)-MyData(i+1,1)
if k~=0
n=n+1;
end
Shujv(n,j+1)=MyData(i,j+1)+M(i,j).*k;
end
end
Shujv(:,1)=linspace(MyData(1,1),MyData(end,1),size(Shujv,1)) c1=Ydata(1:40,3);
x1=linspace(1,length(a1),length(a1));。