北京科技大学matlab大作业Word下载.docx
《北京科技大学matlab大作业Word下载.docx》由会员分享,可在线阅读,更多相关《北京科技大学matlab大作业Word下载.docx(5页珍藏版)》请在冰豆网上搜索。
200449.16
200559.84
200669.65
200793.51
200889.40
200968.58
预计2012-2015年世界原油价格。
建立模型
世界油价的走势一定可以用一个函数来拟合它,因为我们已经知道最近几年的油价走势,和2010,2011年的油价情况,所以,我们可以用拟合的函数图象与2010,2011年的原油价进行比对,找出最相近的拟合函数,对他进行分析找出其未来2012-2015年所对应的原油价格,并对他们的图形进行合理分析,看出偏差的原因。
模型求解
首先我们要求出在用不同拟合曲线中最接近2010,与2011年的曲线,然后取出它单独作为分析对象。
(因世界油价取决因素过多,所以我们只能在短时间范围内进行拟合,所以我只取到2015年的世界原油价格)
用matlab求数值拟合函数图象
x=2003:
1:
2009;
y=[40.33,49.16,59.84,69.65,93.51,89.40,95.18];
m2=polyfit(x,y,2);
m3=polyfit(x,y,3);
m5=polyfit(x,y,5);
m9=polyfit(x,y,9);
disp('
二次拟合'
),f2=poly2str(m2,'
x'
)
三次拟合'
),f3=poly2str(m3,'
五次拟合'
),f5=poly2str(m5,'
九次拟合'
),f9=poly2str(m9,'
x1=2003:
2015;
y2=polyval(m2,x1);
y3=polyval(m3,x1);
y5=polyval(m5,x1);
y9=polyval(m9,x1);
wch2=abs(y-polyval(m2,x))./y
wch3=abs(y-polyval(m3,x))./y
wch5=abs(y-polyval(m5,x))./y
wch9=abs(y-polyval(m9,x))./y
pjwch2=mean(wch2)
pjwch3=mean(wch3)
pjwch5=mean(wch5)
pjwch9=mean(wch9)
plot(x,y,'
rp'
x1,y2,'
--'
x1,y3,'
k-'
x1,y5,x1,y9)
legend('
,'
由图可以得出应选择二次拟合图象,再把二次拟合图象单独分析
m=polyfit(x,y,3);
),f=poly2str(m,'
wch3=abs(y-polyval(m,x))./y
得到的函数m输出为
m=-1.2683.*x.^3+7630.5197.*x.^2-15302171.1633.*x+10228931022.1912
由此可拟合估计2012-2015的值
m=-1.2683.*x.^3+7630.5197.*x.^2-15302171.1633.*x+10228931022.1912;
xi=2012:
mi1=interp1(x,m,xi,'
*nearest'
);
mi2=interp1(x,m,xi,'
*linear'
mi3=interp1(x,m,xi,'
*spline'
mi4=interp1(x,m,xi,'
*cubic'
plot(x,m,'
ro'
xi,mi1,'
xi,mi2,'
-'
xi,mi3,'
k.-'
xi,mi4,'
m:
'
原始数据'
'
最近点插值'
线性插值'
样条插值'
立方插值'
油价已成为现今社会人们非常关注的一个话题,我们对油价的估计也许并非很准确,但也能说明一些问题,如下面专家所分析的那样。
综合外电7月29日报道,美国财政部负责国际事务的副部长麦考密克29日称,有关投机因素和美元币值在原油价格上涨中所起作用的争论一直在持续;
供需不平衡是导致原油价格上涨的根本原因。
麦考密克在一份事先准备好的讲稿中表示,由于需求强劲增长而产量增长缓慢,原油价格必须上涨才能令有限的原油得到分配。
麦考密克称,美元币值下跌对原油价格的影响相对较小,而且没有证据表明原油市场存在操纵现象。
麦考密克虽然没有完全否认这些因素的影响,但强调供需因素是导致原油价格高涨的根本原因。