即它可以看成非独立变量，不必列入回归方程。 (5分) 回归方程可写为：7=b+x+bx3+bx (2) 将表中所列实验数据代入多元线性回归子程序（其中调用高斯消去法子程序解正规方 程组)，通过计算确定方程中的各个系数和复相关系数。 (5分) 2.程序框图(30分) 开始 输入：数据点数N,三种气体的物质的量分数x,2, 粘度7的数值(8分) 调用多元线性回归子程序计算回归方程的系数b0.b1,b2.b3及复相关系数R (其中调用高斯消去法子程序解正规方程组)(15分)N) 输出：b0,b1,b2,b3及R(7分) 结束 3.源程序(30分) (15分) COMMO A10.1i),X50,10us=oL OPEN(6.FILE-2005TB2.TXT.STATUS-UNKNOWN) 的mt 2义可 Y.7X.Yealc.7X:DIFF/IX.66(H-) D030=1,M 30 YCALC+A(J.MI)X(IJ) 40 DIFF WEHEGO XLJU-LMI.YCALC.DIFF5 即它可以看成非独立变量，不必列入回归方程。 （5 分） 回归方程可写为： 0 11 2 2 33 η = b bx bx bx ++ + （2） 将表中所列实验数据代入多元线性回归子程序（其中调用高斯消去法子程序解正规方 程组），通过计算确定方程中的各个系数和复相关系数。 （5 分） 2.程序框图（30 分） 3.源程序（30 分） program main （15 分） DOUBLE PRECISION XV,S,A,X,B0,QT,U,QE,SE,R,F dimension XV(11),S(11) COMMON /AX/A(10,11),X(50,10) OPEN(5,FILE='2005TB1.TXT',STATUS='OLD') OPEN(6,FILE='2005TB2.TXT',STATUS='UNKNOWN') READ(5,*) M,N M1=M+1 READ(5,*) ((X(I,J),J=1,M1),I=1,N) CALL MLR(M,N,XV,S,B0,QT,U,QE,SE,R,F) WRITE(6,20) M,N 20 FORMAT(/4X,'INPUT DATA:','Xn1,'1X,'Xn2,....',1X,'Xnm,Yn(m=1,' & ,i2,',n=1,',i2,')' / 7X,'X1',1X,'X2',1X,'X3',1X,'Y=η'& /1X,66(1H-)/6X,'X1',8X,'X2',8X,'X3',8X, & 'Y',7X,'Ycalc',7X,'DIFF'/1X,66(1H-)) DO 40 I=1,N YCALC=B0 DO 30 J=1,M 30 YCALC=YCALC+A(J,M1)*X(I,J) DIFF=YCALC-X(I,M1) 40 WRITE(6,50) (X(I,J),J=1,M1),YCALC,DIFF 开始 输入：数据点数 N, 三种气体的物质的量分数 x1，x2，x3 粘度η的数值（8 分） 调用多元线性回归子程序计算回归方程的系数 b0,b1,b2,b3 及复相关系数 R （其中调用高斯消去法子程序解正规方程组） （15 分）N） 输出：b0,b1,b2,b3 及 R （7 分） 结束