南开大学计算物理讲义3.pptVIP

* 4. Secant method 割线, 正切 * Convergence on the root using the secant method * Divergence using the secant method * * =1.618 =0.618 * 例：用割线法求方程 在x = 1.5附近的根，使绝对误差精确到10-4。 按上式计算得： 解：取初值x0 = 1.5，x1 = 1.4，按公式 得迭代格式 x2 = 1.33522；x3 = 1.32541； x4 = 1.32472； x5 = 1.32472 取x* ? 1.32472 * f(x)=exp(x)*ln(x)-x2 f(1)=-1 f(2)=1 x?[1,2] A=1.0 B=2.0 DX=(B-A)/10.0 X0=(A+B)/2.0 WRITE(*,*) ISTEP,X0,DX STOP END Call secant(DL,X0,DX,ISTEP) FUNCTION F(X) F=exp(x)*ALOG(x)-X*X RETURN END DL=1.0E-06 * DO WHILE(ABS(DX).GT.DL) D= F(X1)- F(X0) X2=X1-(X1-X0)* F(X1)/D X0=X1 X1=X2 DX=X1-X0 X1=X0+DX ISTEP=0 ISTEP=ISTEP+1 ENDDO SUBROUTINE secant(DL,X0,DX,ISTEP) RETURN END * 5. Direct iteration （直接迭代法） 迭代法是求方程近似根的一个重要方法，也是计算方法中的一种基本方法，它的算法简单，应用范围广，故学好这一方法对以后的学习也很有帮助。 * 简单迭代法 从而构成序列 简单迭代法的基本思想是：将方程f (x) = 0化为一个等价的方程 给定一个初值x0，由上式可算得x1 = ? (x0)，再将x1代入上式的右端又可得x2 =? (x1)，…。我们称{xk}为迭代序列，而称 x=? (x) 中的? (x)为迭代函数，上式为迭代格式。如果? (x)连续，迭代序列{xk}收敛于x*，则x*就是方程 x=? (x) 的解。 * 方程f (x) = 0 等价于方程 * 迭代法的几何意义 交点的横坐标 y=x ? í ì = = T = ) ( ) ( x y x y x x j j y=?(x) 方程f (x) = 0 * 求以下方程在[0,1]内的根的近似值 现将原方程改为同解方程 由此得迭代格式 取初始值x0 = 1，可逐次算得 x1 = 0.4771, x2 = 0.3939 … x6 = 0.3758, x7 =0.3758 因为x6 和x7 已趋于一致，所以取x7 =0.3758 为原方程在[0, 1]内的一个根的近似值。 * 一个方程的迭代格式并不是唯一的，且迭代也不总是收敛的。如上例的方程也可改写 得迭代格式 仍取x0 = 1算得： 显然，该迭代序列不趋向于某个定值，这种不收敛的迭代过程称为发散。那么迭代格式要满足哪些条件才能保证迭代收敛呢？下面我们来讨论这个问题。 * |en+1| ~ g’|en |p. p = 1 收敛;否则发散 * 迭代过程的收敛性 y=x y y = ) ( x j 1 ) ( 0 * ￠ x j P 0 P 2 P * Q 1 Q 2 x 1 x 0 x 2 x * x y P 1 * 迭代过程的收敛性 y y=x 0 ) ( 1 * ￠ - x j x 0 x x 1 x 2 x 3 x * y = ) ( x j ) ( x j P * * 迭代法的几何意义 1 ) ( * ￠ x j y=x y = ) ( x j P * x 1 x 0 x y x * x 2 * 迭代法的几何意义 y x 0 x x 1 x 2 x 3 x * y = ) ( x j ) ( x j P * y=x -1 ) ( * ￠ x j *