【正文】 問題二的求解 附表 21 a=[49271 53292 58729 63508 66362 71494 73655 73861 78067 83529 93853 105930 109381 104372 107077 98738 89455 81649 73484 67759 65225] b=[ ] c=zeros(1,21)。 X0=data。 X1(i)=X1(i1)+X0(i)。 end B=zeros(T1,2)。 elseif j1 B(i,j)=1。 end HCS=inv(B’*B)*B’*Y’ H=HCS’。 end for i=21:T+N K(i20)=XR1(i)XR1(i1)。 運行結果 data = +07 * Columns 1 through 8 Columns 9 through 16 Columns 17 through 21 HCS = +07 * ans = +08 * 求 HCS=[a,b] 參數和預測 20202020 年死亡人數的程序 function GM=huise(data,N) T=length(data)。 for i=2:T X1(1)=X0(1)。 end for i=1:T1 M(i)=(*(X1(i)+X1(i+1)))。 for i=1:T1 for j=1:2 if j2 B(i,j)=M(i)。 end end end for i=2:T Y(i1)=X0(i)。 for i=2:T+N XR1(i)=(X0(1)H(2)/H(1))*exp(1*H(1)*(i1))+H(2)/H(1)。 end GM=K。209。202。202。190。177。180。212。180。206。177。206。188。208。 n=6。170。198。175。189。249。196。238。253。 end yhat1 m2=length(yhat1)。 end yhat2 plot(1:21,y223,’*’) a21=2*yhat1(end)yhat2(end) b21=2*(yhat1(end)yhat2(end))/(n1) y2020=a21+b21 y2020=a21+2*b21 y2020=a21+3*b21 y2020=a21+4*b21 y2020=a21+5*b21 運行結果 yhat1 = +008 * yhat2 = +008 * a21 = +008 b21 = +007 y2020 = +008 y2020 = +008 y2020 = +008 y2020 = +008 y2020 = +008 clc,clear load %176。212。188。253。221。163。230。218。191。196。190。196。 214。 m1=length(y221)。 %n 206。210。182。198。190。181。207。202。 for i=1:m1n+1 yhat1(i)=sum(y221(i:i+n1))/n。 for i=1:m2n+1 yhat2(i)=sum(yhat1(i:i+n1))/n。209。202。202。190。177。180。212。180。206。177。206。188。208。 n=6。170。198。175。189。249。196。238。253。 end yhat1 m2=length(yhat1)。 end yhat2 plot(1:21,y222,’*’) a21=2*yhat1(end)yhat2(end) b21=2*(yhat1(end)yhat2(end))/(n1) y2020=a21+b21 y2020=a21+2*b21 y2020=a21+3*b21 y2020=a21+4*b21 y2020=a21+5*b21 運行結果 yhat1 = yhat2 = a21 = b21 = y2020 = y2020 = y2020 = y2020 = y2020 = 附表 41 用灰色系統(tǒng)模型預測 20202020事故起數與直接經濟損失的程序 data = +05 * Columns 1 through 8 Columns 9 through 14 HCS = +04 * ans = +05 * Untitled10 data = +05 * Columns 1 through 8 Columns 9 through 14 HCS = +04 * ans = +05 * function GM=huise(data,N) T=length(data)。 for i=2:T X1(1)=X0(1)。 end for i=1:T1 M(i)=(*(X1(i)+X1(i+1)))。 for i=1:T1 for j=1:2 if j2 B(i,j)=M(i)。 end end end for i=2:T Y(i1)=X0(i)。 for i=2:T+N XR1(i)=(X0(1)H(2)/H(1))*exp(1*H(1)*(i1))+H(2)/H(1)。 end GM=K。 n=6。 end yhat1 m2=length(yhat1)。 end yhat2 plot(1:14,y3,’*’) a21=2*yhat1(end)yhat2(end) b21=2*(yhat1(end)yhat2(end))/(n1) y2020=a21+b21 y2020=a21+2*b21 y2020=a21+3*b21 y2020=a21+4*b21 y2020=a21+5*b21 y2020=a21+6*b21 y2020=a21+7*b21 運行結果 yhat1 = +005 * yhat2 = +005 * a21 = +005 b21 = +004 y2020 = +005 y2020 = +005 y2020 = +005 y2020 = +006 y2020 = +006 y2020 = +006 y2020 = +006 clc,clear load m1=length(y1)。 for i=1:m1n+1 yhat1(i)=sum(y1(i:i+n1))/n。 for i=1:m2n+1 yhat2(i)=sum(yhat1(i:i+n1