《gm(1,1)程序编写及数据文件》由会员分享，可在线阅读，更多相关《gm(1,1)程序编写及数据文件（2页珍藏版）》请在金锄头文库上搜索。
1、 integer N,IINTEGER,PARAMETER : N_MAX=20REAL U0(N_MAX),U1(N_MAX),BU1(N_MAX),C1(N_MAX),C2(N_MAX),y(N_MAX)REAL B11,B12,B21,B22,bB11,bB12,bB21,bB22,a,b,M,d1,d2,d3,d4open(unit=10,file=out.txt)open(unit=20,file=date.txt)WRITE (*,*) PLEASE INPUT YOUR NUMBER: N=WRITE (10,*) PLEASE INPUT YOUR NUMBER: N=!REA
2、D (*,*) NREAD (20,*) NWRITE (10,*) NWRITE (*,*) PLEASE INPUT YOUR INITAL NUMBER :WRITE (10,*) PLEASE INPUT YOUR INITAL NUMBER :DO I=1,n ! READ (*,*) U0(I)READ (20,*) U0(I)END DO u1(1)=u0(1)DO I=2,NU1(I)=U0(i)+U1(I-1)END DODO I=2,N BU1(I)=-(U1(i-1)+U1(i)/2END DO b11=0.0b12=0.0b21=0.0b22=0.0DO I=2,N B
3、11=BU1(I)*2+B11B12=BU1(I)+B12b21=BU1(I)+B21B22=1.0+b22END DO M=B11*B22-B12*B21bB11=b22bB12=-B12bB21=-B21bB22=B11DO I=2,N C1(I)=bB11*BU1(I)+bB12C2(I)=bB21*BU1(I)+bB22END DOa=0.0b=0.0DO I=2,N a=c1(i)*u0(i)/M+ab=c2(i)*u0(i)/M+bEND DO d1=(b11*bb11+b12*bb21)/md2=(b11*bb12+b12*bb22)/md3=(b21*bb11+b22*bb21
4、)/md4=(b21*bb12+b22*bb22)/mwrite (*,*) a,bwrite (*,*) d1,d2,d3,d4write (*,*) bb11,bb12,bb21,bb22write (*,*) b11,b12,b21,b22,mwrite (10,*) youhuahou:u1write (10,*) u0(1)DO I=2,N+5 y(i)=(1-exp(a)*(u0(1)-(b/a)*exp(-a*(I-1)write (*,*) y(I)write (10,*) y(I)END DO end 以下另保存一个新的 txt 格式，命名 date10 0.66 1.15 1.57 2.19 3.04 3.99 5.42 7.48 9.54 12.9
