北京科技大学matlab大作业.pdf
分析未来五年世界汽油价格问题机械 1011实验目的实验目的通过本次实验,目的在于提高学生根据实际问题建立函数拟合与插值方面的问题的运用巩固我们的有关知识,在 21 世纪这个对能源问题而苦恼的世纪中,让我们认识到现实,为资源节约尽一份力。实验问题根据以往油价,推测未来五年的油价走势,最终估计按此发展下去,在2012-2015 年的原油价格会达到多少,我们应如何应对这些情况。日期国际原油价(美元/桶)200340.33200449.16200559.84200669.65200793.51200889.40200968.58预计 2012-2015 年世界原油价格。建立模型建立模型世界油价的走势一定可以用一个函数来拟合它,因为我们已经知道最近几年的油价走势,和 2010,2011年的油价情况,所以,我们可以用拟合的函数图象与2010,2011 年的原油价进行比对,找出最相近的拟合函数,对他进行分析找出其未来 2012-2015 年所对应的原油价格,并对他们的图形进行合理分析,看出偏差的原因。模型求解模型求解首先我们要求出在用不同拟合曲首先我们要求出在用不同拟合曲线中最接近线中最接近 20102010,与与 20112011 年的曲线,年的曲线,然后取出它单独作为分析对象。然后取出它单独作为分析对象。(因(因世界油价取决因素过多,所以我们只世界油价取决因素过多,所以我们只能在短时间范围内进行拟合,所以我能在短时间范围内进行拟合,所以我只取到只取到 20152015 年的世界原油价格)年的世界原油价格)用用 matlabmatlab 求数值拟合函数图象求数值拟合函数图象x=2003:1:2009;x=2003:1:2009;y=40.33,49.16,59.84,69.65,93.51,89.40,y=40.33,49.16,59.84,69.65,93.51,89.40,95.18;95.18;m2=polyfit(x,y,2);m2=polyfit(x,y,2);m3=polyfit(x,y,3);m3=polyfit(x,y,3);m5=polyfit(x,y,5);m5=polyfit(x,y,5);m9=polyfit(x,y,9);m9=polyfit(x,y,9);disp(disp(二次拟合二次拟合),f2=poly2str(m2,x),f2=poly2str(m2,x)disp(disp(三次拟合三次拟合),f3=poly2str(m3,x),f3=poly2str(m3,x)disp(disp(五次拟合五次拟合),f5=poly2str(m5,x),f5=poly2str(m5,x)disp(disp(九次拟合九次拟合),f9=poly2str(m9,x),f9=poly2str(m9,x)x1=2003:1:2015;x1=2003:1:2015;y2=polyval(m2,x1);y2=polyval(m2,x1);y3=polyval(m3,x1);y3=polyval(m3,x1);y5=polyval(m5,x1);y5=polyval(m5,x1);y9=polyval(m9,x1);y9=polyval(m9,x1);wch2=abs(y-polyval(m2,x)./ywch2=abs(y-polyval(m2,x)./ywch3=abs(y-polyval(m3,x)./ywch3=abs(y-polyval(m3,x)./ywch5=abs(y-polyval(m5,x)./ywch5=abs(y-polyval(m5,x)./ywch9=abs(y-polyval(m9,x)./ywch9=abs(y-polyval(m9,x)./ypjwch2=mean(wch2)pjwch2=mean(wch2)pjwch3=mean(wch3)pjwch3=mean(wch3)pjwch5=mean(wch5)pjwch5=mean(wch5)pjwch9=mean(wch9)pjwch9=mean(wch9)plot(x,y,rp,x1,y2,-,x1,y3,k-,x1,y5,x1,y9)plot(x,y,rp,x1,y2,-,x1,y3,k-,x1,y5,x1,y9)legend(legend(二次拟合二次拟合,三次拟合三次拟合,五次拟合五次拟合,九次拟合九次拟合)由图可以得出应选择二次拟合图象,再把二次拟合图象单独分析x=2003:1:2009;y=40.33,49.16,59.8469.65,93.51,89.40,95.18;m=polyfit(x,y,3);disp(三次拟合),f=poly2str(m,x)wch3=abs(y-polyval(m,x)./y,得到的函数 m 输出为m=-1.2683.*x.3+7630.5197.*x.2-15302171.1633.*x+10228931022.1912由此可拟合估计 2012-2015 的值x=2003:2015;m=-1.2683.*x.3+7630.5197.*x.2-15302171.1633.*x+10228931022.1912;xi=2012:1:2015;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:)legend(原始数据,最近点插值,线性插值,样条插值,立方插值)油价已成为现今社会人们非常关注的一个话题,我们对油价的估计也许并非很准确,但也能说明一些问题,如下面专家所分析的那样。综合外电 7 月 29 日报道,美国财政部负责国际事务的副部长麦考密克 29 日称,有关投机因素和美元币值在原油价格上涨中所起作用的争论一直在持续;供需不平衡是导致原油价格上涨的根本原因。麦考密克在一份事先准备好的讲稿中表示,由于需求强劲增长而产量增长缓慢,原油价格必须上涨才能令有限的原油得到分配。麦考密克称,美元币值下跌对原油价格的影响相对较小,而且没有证据表明原油市场存在操纵现象。麦考密克虽然没有完全否认这些因素的影响,但强调供需因素是导致原油价格高涨的根本原因。