Fortran을 이용한 직선회귀, 다항식 회귀법

소개글

본 보고서는 Fortran으로 직성 회귀, 다항식 회귀법을 직접 인코딩 해보고 그 결과를 정리하였다.

목차

Question 5.1, 직선회귀법
Question 5.3, 다항식 회귀법

본문내용

Question 5.1, 직선회귀법
C LEAST SQUARE METHOD
DIMENSION X(99), F(99)
READ(5,*) N
READ(5,*)(X(I), F(I), I = 1, N)
XS = 0.0
YS = 0.0
XXS = 0.0
XYS = 0.0
DO 10 I = 1, N
XS = XS + X(I)
YS = YS + F(I)
XXS = XXS + X(I) * X(I)
XYS = XYS + X(I) * F(I)
10 CONTINUE
B = (N * XYS - XS * YS) / (N * XXS - XS * XS)
A = YS / N - B *XS / N
WRITE(6, 200) A, B
200 FORMAT(// , 5X, ` F(X) = `, F10.4, ` + `, F10.4, ` X`/)
STOP
END



Question 5.3, 다항식 회귀법
POLYNOMIAL REGRESSION METHOD
COMMON M, N, X(50), F(50), A(10), AM(11, 12), C(50, 11)
READ (5,*) N
READ(5,*) M
READ(5,*) (X(I), F(I), I = 1, N)
CALL LEAST
WRITE(6,200) (I-1,A(I) ,I = 1, M+1)
200 FORMAT(2X,`A(`,I2,`) = `,F10.6)
WRITE(6,210)
210 FORMAT(6X, `X`,8X, `F(X)`, 9X, `ERROR`)
DO 10 I = 1, N
FX = A(M+1)
DO 20 J = M,1,-1
FX = FX * X(I) + A(J)
20 CONTINUE
ERROR = FX - F(I)
WRITE(6, 220) X(I),FX,ERROR
220 FORMAT(2X, 2F10.6,E15.6/)

참고 자료

없음