【正文】 ),xx]。[b,bint,r,rint,stats]=regress(yy,xx1,)rcoplot(r,rint)二、問(wèn)題二程序function [U,xx0,XXX,E,e,X0_ave,S1,E_ave,S2,C,P]=GM11(X0)%（1）由原始數(shù)據(jù)列計(jì)算一次累加序列X1[m,n]=size(X0)。X1(1)=X0(1)。for i=2:n X1(i)=X1(i1)+X0(i)。end%（2）建立矩陣B，yB=ones(n1,2)。for j=1:n1 B(j)=((X1(j)+X1(j+1)))/2。 y(j,1)=X0(j+1)。end%(3)計(jì)算C=（BTB）1CC=inv(B39。*B)%(4)由U=((BTB)1)BTy,估值a,u.%得時(shí)間響應(yīng)方程X1(k+1)=(X1(1)u/a)*exp(a*k)+u/a。U=CC*B39。*y。%輸出a,u的值a=U(1,1)。u=U(2,1)。disp(39。a,u值：39。)%(5)計(jì)算擬合值xxxx(1)=X1(1)。for i=2:n xx(i)=(X1(1)u/a)*exp(a*(i1))+u/a。end%還原計(jì)算得到模型計(jì)算值xx0（i）；xx0(1)=xx(1)。for k=2:n xx0(n+2k)=xx(n+2k)xx(n+2k1)。end%精度檢驗(yàn)：%【1】殘差檢驗(yàn)%E為殘差，e為相對(duì)殘差E(1)=0。e(1)=0。for i=2:n E(i)=X0(i)xx0(i)。 e(i)=(X0(i)xx0(i))/X0(i)。enddisp(E)disp(e)%【2】后驗(yàn)差檢驗(yàn)均值，方差，殘差均值，殘差方差，后驗(yàn)差比值，小誤差概率X0_ave=sum(X0)/n。for i=1:n X0_cha(i)=(X0(i)sum(X0)/n)^2。 E_cha(i)=(E(i)sum(E)/(n1))^2。endS1=sqrt(sum(X0_cha)/n)。E_ave=sum(E)/(n1)。S2=sqrt(sum(E_cha)/(n1))C=S2/S1。P=0。for i=1:n dd(i)=abs(E(i)sum(E)/(n1))。endfor i=1:nif dd(i)*S1 P=P+1。endendP=P/(n1)。if Pamp。amp。C disp(39。預(yù)測(cè)精度等級(jí)：好39。) elseif Pamp。amp。C disp(39。預(yù)測(cè)精度等級(jí)：合格39。) elseif Pamp。amp。C disp(39。預(yù)測(cè)精度等級(jí)：勉強(qiáng)39。) elseif Pamp。amp。C disp(39。預(yù)測(cè)精度等級(jí)：不合格39。)endXXX=X1(1)。return:clear。clcx=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\39。,39。sheet139。,39。c15:c4639。)39。[U,xx0,XXX,E,e,X0_ave,S1,E_ave,S2,C,P]=GM11(x)。n=64。for i=1:n xx(i)=(XXXU(2,1)/U(1,1))*exp(U(1,1)*(i1))+U(2,1)/U(1,1)。end%還原計(jì)算得到模型計(jì)算值xx0（i）；xx0(1)=XXX。for k=2:n xx0(n+2k)=xx(n+2k)xx(n+2k1)。end3預(yù)測(cè)出生率x=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\39。,39。sheet139。,39。b2:b3239。)39。[U,xx0,XXX,E,e,X0_ave,S1,E_ave,S2,C,P,NNN]=GM11(x)。n=63。for i=1:n xx(i)=(XXXU(2,1)/U(1,1))*exp(U(1,1)*(i1))+U(2,1)/U(1,1)。end%還原計(jì)算得到模型計(jì)算值xx0（i）；xx0(1)=XXX。for k=2:n xx0(n+2k)=xx(n+2k)xx(n+2k1)。endx=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\39。,39。sheet139。,39。d2:d3239。)39。[U,xx0,XXX,E,e,X0_ave,S1,E_ave,S2,C,P,NNN]=GM11(x)。n=63。for i=1:n xx(i)=(XXXU(2,1)/U(1,1))*exp(U(1,1)*(i1))+U(2,1)/U(1,1)。end%還原計(jì)算得到模型計(jì)算值xx0（i）；xx0(1)=XXX。for k=2:n xx0(n+2k)=xx(n+2k)xx(n+2k1)。endx=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\39。,39。sheet139。,39。b60:j6039。)。[U,xx0,XXX,E,e,X0_ave,S1,E_ave,S2,C,P,NNN]=GM11(x)。n=63。for i=1:n xx(i)=(XXXU(2,1)/U(1,1))*exp(U(1,1)*(i1))+U(2,1)/U(1,1)。end%還原計(jì)算得到模型計(jì)算值xx0（i）；xx0(1)=XXX。for k=2:n xx0(n+2k)=xx(n+2k)xx(n+2k1)。endx=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\39。,39。sheet139。,39。 b64:j6439。)。[U,xx0,XXX,E,e,X0_ave,S1,E_ave,S2,C,P,NNN]=GM11(x)。n=63。for i=1:n xx(i)=(XXXU(2,1)/U(1,1))*exp(U(1,1)*(i1))+U(2,1)/U(1,1)。end%還原計(jì)算得到模型計(jì)算值xx0（i）；xx0(1)=XXX。for k=2:n xx0(n+2k)=xx(n+2k)xx(n+2k1)。end：傅里葉擬合程序clear。clcy1=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\39。,39。sheet239。,39。b1:b3139。)。x1=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\39。,39。sheet239。,39。c1:c3139。)。p1=cftool(x1,y1,8) for x=1:63 f(x) = a0 + a1*cos(x*w) + b1*sin(x*w) + ... a2*cos(2*x*w) + b2*sin(2*x*w) + a3*cos(3*x*w) + b3*sin(3*x*w)+... a4*cos(4*x*w) + b4*sin(4*x*w) + a5*cos(5*x*w) + b5*sin(5*x*w)+... a6*cos(6*x*w) + b6*sin(6*x*w) + a7*cos(7*x*w) + b7*sin(7*x*w)+... a8*cos(8*x*w) + b8*sin(8*x*w)。 end三、問(wèn)題三程序clear。clccc1=xlsread(39。D:\0000第一次模擬訓(xùn)練\第二問(wèn)數(shù)據(jù)整理結(jié)果\39。,39。自然增長(zhǎng)率39。,39。K5:K3439。)。y1=xlsread(39。D:\0000第一次模擬訓(xùn)練\第二問(wèn)數(shù)據(jù)整理結(jié)果\39。,39。自然增長(zhǎng)率39。,39。O5:O3439。)。y2=xlsread(39。D:\0000第一次模擬訓(xùn)練\第二問(wèn)數(shù)據(jù)整理結(jié)果\39。,39。自然增長(zhǎng)率39。,39。Q5:Q3439。)。[m,n]=size(cc1)。alpha=ones(m,n)*5。cc2=cc1+alpha。pp1=polyfit(cc1,y1,1)。pp2=polyfit(cc1,y2,2)。yy1=polyval(pp1,cc1)。yy2=polyval(pp2,cc1)。yy3=polyval(pp1,cc2)。yy4=polyval(pp2,cc2)。plot(cc1,yy1,39。r39。,cc1,yy2,39。r39。,cc2,yy3,39。b39。,cc2,yy4,39。b39。)。四、問(wèn)題四程序x1=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\固定資產(chǎn)投資39。,39。sheet239。,39。c2:o1639。)。x2=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\固定資產(chǎn)投資39。,39。sheet239。,39。c20:o2739。)。y1=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\固定資產(chǎn)投資39。,39。sheet239。,39。b2:b1639。)。[m,n]=size(x1)。%計(jì)算樣本的相關(guān)系數(shù)矩陣Rfor i=1:m for j=1:n xx1(i,j)=x1(i,j)/x1(1,j)。 % xx為對(duì)原數(shù)據(jù)矩陣無(wú)量綱化后的矩陣 endend[mm,nn]=size(x2)。for i=1:mm for j=1:nn xx2(i,j)=x2(i,j)/x1(1,j)。 % xx為對(duì)原數(shù)據(jù)矩陣無(wú)量綱化后的矩陣 endend [m,n]=size(y1)。%計(jì)算樣本的相關(guān)系數(shù)矩陣Rfor i=1:m for j=1:n yy1(i,j)=y1(i,j)/y1(1,j)。 % xx為對(duì)原數(shù)據(jù)無(wú)量綱化后的矩陣 endendxxx1=[ones(15,1),xx1]。[b,bint,r,rint,stats]=regress(yy1,xxx1,)xxx2=[ones(8,1),xx2]。yy2=xxx2*b。yy3=xxx1*b。38 / 38