成品油定價(jià)體系-資料下載頁(yè)
2025-05-23 18:30本頁(yè)面
【正文】 教育出版社 2003年 【2】孫祥 徐流美 吳清《》 清華大學(xué)出版社 2005年【3】，2002【4】：華中科技大學(xué)出版社，2002.【5】紀(jì)利群 中國(guó)石油大學(xué)報(bào) 第35卷4期 2011【6】【7】附錄：MATLAB源代碼%一次函數(shù)的回歸分析x1=[1:23]39。y=[ ]39。x=[ones(23,1),x1]。[b,bint,r,rint,stats]=regress(y,x)。b,bint,stats,rcoplot(r,rint )%二次函數(shù)的回歸分析x=[1:24]。y=[ ]。A=polyfit(x,y,2)z=polyval(A,x)plot(x,y,39。k+39。,x,z,39。r39。)%GM（1,1）實(shí)現(xiàn)y=input(39。請(qǐng)輸入數(shù)據(jù) 39。)。n=length(y)。yy=ones(n,1)。yy(1)=y(1)。for i=2:nyy(i)=yy(i1)+y(i)。endB=ones(n1,2)。for i=1:(n1) B(i,1)=(yy(i)+yy(i+1))/2。 B(i,2)=1。endBT=B39。for j=1:n1 YN(j)=y(j+1)。endYN=YN39。A=inv(BT*B)*BT*YN。a=A(1)。u=A(2)。t=u/a。t_test=input(39。請(qǐng)輸入需要預(yù)測(cè)個(gè)數(shù)：39。)。i=1:t_test+n。yys(i+1)=(y(1)t).*exp(a.*i)+t。yys(1)=y(1)。for j=n+t_test:1:2 ys(j)=yys(j)yys(j1)。endx=1:n。xs=2:n+t_test。yn=ys(2:n+t_test)。plot(x,y,39。^r39。,xs,yn,39。*b39。)。det=0。for i=2:n det=det+abs(yn(i)y(i))。enddet=det/(n1)。disp([39。預(yù)測(cè)值為： 39。,num2str(ys(n+1:n+t_test))])。然后能將數(shù)據(jù)[ ]輸入程序，預(yù)測(cè)數(shù)目為6個(gè)%層次分析法a=[ 1 6 5 2 1 7 1/21/6 1 1/5 1/5 1/7 1 1/71/5 5 1 1/3 1/4 3 1/41/2 5 3 1 1/2 3 1/21 7 4 2 1 7 11/7 1 1/3 1/3 1/7 1 1/72 7 4 2 1 7 1]。[x,y]=eig(a)。eigenvalue=diag(y)。lamda=eigenvalue(1)。ci1=(lamda7)/6。cr1=ci1/w1=x(:,1)/sum(x(:,1))b1=[ 1 2 2 3 5 1/2 1 1 2 3 1/2 1 1 2 4 1/3 1/2 1/2 1 2 1/5 1/3 1/4 1/2 1]。[x,y]=eig(b1)。eigenvalue=diag(y)。lamda=eigenvalue(1)。ci21=(lamda5)/4。cr21=ci21/w21=x(:,1)/sum(x(:,1))b2=[ 1 2 4 5 7 1/2 1 1/2 2 4 1/4 2 1 3 5 1/5 1/2 1/3 1 2 1/7 1/4 1/5 1/2 1]。[x,y]=eig(b2)。eigenvalue=diag(y)。lamda=eigenvalue(1)。ci22=(lamda5)/4。cr22=ci22/w22=x(:,1)/sum(x(:,1))b3=[ 1 4 5 6 8 1/4 1 2 4 5 1/5 1/2 1 2 4 1/6 1/4 1/2 1 3 1/8 1/5 1/4 1/3 1]。[x,y]=eig(b3)。eigenvalue=diag(y)。lamda=eigenvalue(1)。ci23=(lamda5)/4。cr23=ci23/w23=x(:,1)/sum(x(:,1))b4=[ 1 3 4 5 7 1/3 1 2 3 5 1/4 1/2 1 2 4 1/5 1/3 1/2 1 2 1/7 1/5 1/4 1/2 1]。[x,y]=eig(b4)。eigenvalue=diag(y)。lamda=eigenvalue(1)。ci24=(lamda5)/4。cr24=ci24/w24=x(:,1)/sum(x(:,1))b5=[ 1 2 3 5 61/2 1 2 4 5 1/3 1/2 1 2 41/5 1/4 1/2 1 21/6 1/5 1/4 1/2 1]。[x,y]=eig(b5)。eigenvalue=diag(y)。lamda=eigenvalue(2)。ci25=(lamda5)/4。cr25=ci25/w25=x(:,2)/sum(x(:,2))b6=[ 1 4 4 5 7 1/4 1 3 4 6 1/4 1/3 1 3 5 1/5 1/4 1/3 1 3 1/7 1/6 1/5 1/3 1]。[x,y]=eig(b6)。eigenvalue=diag(y)。lamda=eigenvalue(1)。ci26=(lamda5)/4。cr26=ci26/w26=x(:,1)/sum(x(:,1))b7=[ 1 2 3 4 6 1/2 1 2 4 6 1/3 1/2 1 2 4 1/4 1/4 1/2 1 3 1/6 1/5 1/4 1/3 1]。[x,y]=eig(b7)。eigenvalue=diag(y)。lamda=eigenvalue(1)。ci27=(lamda5)/4。cr26=ci26/w27=x(:,1)/sum(x(:,1))w_sum=[w21,w22,w23,w24,w25,w26,w27]*w1ci=[ci21,ci22,ci23,ci24,ci25,ci26,ci27]。cr=ci*w1/sum(*w1)圖41程序x=1:24。y=[ ]。plot(x,y,39。+39。)圖42程序x1=[1:24]39。y=[ ]39。x=[ones(24,1),x1]。[b,bint,r,rint,stats]=regress(y,x)。b,bint,stats,rcoplot(r,rint )圖49程序x=2003:2011。y=[ ]。plot(x,y,39。LineWidth39。,3)hold ony1=[ ]。plot(x,y1,39。r39。,39。LineWidth39。,3)legend(39。國(guó)外原油價(jià)格39。,39。國(guó)內(nèi)汽油價(jià)格39。)xlabel(39。年份39。)ylabel(39。油的價(jià)格39。)grid on圖411 程序x=2006:2011。y=[11759 13786 15781 17175 19109 20714]。plot(x,y,39。bo39。,39。LineWidth39。,2)hold ony1=[4950 5280 6030 6730 7420 8280]。plot(x,y1,39。rd39。,39。LineWidth39。,2)legend(39。人均收入39。,39。油價(jià)39。)xlabel(39。年份39。)ylabel(39。油的價(jià)格與收入39。)grid on來自中國(guó)統(tǒng)計(jì)年鑒的石油消費(fèi)表72 能源消費(fèi)總量及構(gòu)成 年 份能源消費(fèi)總量占能源消費(fèi)總量的比重 (%)(萬(wàn)噸標(biāo)準(zhǔn)煤)煤 炭石 油天然氣水電、核電、風(fēng)電 197857144198060275198576682199098703 19911037831992109170199311599319941227371995131176 19961351921997135909199813618419991405692000145531 20011504062002159431200318379220042134562005235997 20062586762007280508200829144820093066472010324939
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1