迭代法与高斯消元法.pptVIP

迭代法与高斯消元法 计算地球化学 主要内容 1. 一般迭代法 2. 牛顿迭代法 3. 任意函数求根 4. 高斯消元法 1. 一般迭代法 例1. 用迭代法求方程 定理. 牛顿迭代法的基本思想: 牛顿法的误差估计: 牛顿迭代法 例2. 用牛顿法求方程 例2的计算程序 3. 任意函数求根 Subroutine ZREAL(F, ERRABS, ERRREL, EPS, ETA, NROOT, ITMAX, XGUESS, X, INFO) Real, external :: F Real ERRABS Real ERRREL Real EPS Real ETA Integer NROOT Integer ITMAX Real XGUESS(NROOT) Real X(NROOT) Integer INFO(NROOT) 例 3 例 4：求复数根 Subroutine ZANLY (F, ERRABS, ERRREL, NKNOWN, NNEW, NGUESS, ZINIT, ITMAX, Z, INFO) 4. 高斯消元法 矩阵的加、减、乘法 如：real a(m,n), b(m,n), c(m,n) c = a + b c = a – b 如：real a(m,n), b(n,L), c(m,L) c = matmul(a,b) 高斯消元法的子程序 例 5：高斯消元法 解法2：逆矩阵法 解法3：lin_sol_gen法 作 业 1. 用牛顿迭代法编程计算方程 sin(x)+cos(x) = 0的一个近似实数根，误差不超过0.00001。 2. 用高斯消元法编程解方程组： x+y+2z+w=4 2x+5y+8z-w=15 3x+6y+9z+2w=18 4x+5y+6z+7w = 30 X(I)=B(I)-T 70 CONTINUE 100 FORMAT(1X,FAIL) JS(N)=N DO 150 K=N,1,-1 IF(JS(K).NE. K)THEN T=X(K) X(K)=X(JS(K)) X(JS(K))=T ENDIF 150 CONTINUE RETURN END Subroutine AGAUS program main Implicit none double precision :: A(3,3) = (/ 1,2,3, 1,5,6, 2,8,9 /) double precision :: B(3) = (/ 4, 15, 18 /), C(3) ! 求解 x+y+2z=4 ! 2x+5y+8z=15 ! 3x+6y+9z=18 Call AGAUS(A,B,3,C) write(*,*) C stop end program program main use IMSL real :: A(3,3) = (/ 1,2,3, 1,5,6, 2,8,9 /) real :: B(3) = (/ 4, 15, 18 /), C(3) ! 求解 x+y+2z=4 ! 2x+5y+8z=15 ! 3x+6y+9z=18 C = A .ix. B ! invert(A) * B write(*,*) C stop end program program main use IMSL implicit none real :: A(3,3) = (/ 1,3,2, 1,2,1, 2,1,3 /) real :: B(3,1) = (/ 4,6,6 /) real :: X(3,1) call lin_sol_gen(A,B,X) ! A*X=B,解X write(*,(3F5.2)) X stop end program 注意：接受的参数必须是二维数组，故B和X声明成二维数组。 * * 在隔根区 按递推公式 则 ? 即为原方程的根 . ① ①称为迭代格式 , 初值 . 否则称为发散 .