八 非线性方程及方程组的解法培训教程文件.ppt

第八章　非线性方程及非线性方程组解决 由此导出迭代公式 x y x* x0 解： 代入初值得： Newton法迭代公式为 由定理8.4，Newton法至少二阶收敛。 以上两种改进都是至少二阶收敛的。 在Newton公式中，用差商代替导数，即　　　　 即得迭代公式　　　　　　 3.2弦截法 * * a b x0 x1 a1 b2 x* §1　区间对分法（二分法） 1. 确定有根区间: 2. 逐次对分区间： 3. 取根的近似值: b1 a2 其误差为: 根的近似值: a b x0 x1 a1 b2 x* b1 a2 用对分区间法求根步骤： f=inline(x^3+10*x-20,x); [x,err]=bisection(f,1,2,0.00005,15) n x err 1.00000000000000 1.50000000000000 0.50000000000000 2.00000000000000 1.75000000000000 0.25000000000000 3.00000000000000 1.62500000000000 0.12500000000000 4.00000000000000 1.56250000000000 0.06250000000000 5.00000000000000 1.59375000000000 0.03125000000000 6.00000000000000 1.60937500000000 0.01562500000000 7.00000000000000 1.60156250000000 0.00781250000000 8.00000000000000 1.59765625000000 0.00390625000000 9.00000000000000 1.59570312500000 0.00195312500000 10.00000000000000 1.59472656250000 0.00097656250000 11.00000000000000 1.59423828125000 0.00048828125000 12.00000000000000 1.59448242187500 0.00024414062500 13.00000000000000 1.59460449218750 0.00012207031250 14.00000000000000 1.59454345703125 0.00006103515625 15.00000000000000 1.59457397460938 0.00003051757813 x = 1.59457397460938 err =3.051757812500000e-005 §2　简单迭代法 2.1 简单迭代法的一般形式及其几何意义　　　　 x y y = x x y y = x x y y = x x y y = x x* x* x* x* x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 　　 　　　　 显然，此迭代序列发散。 f=inline(x^3+11*x-20,x); iteration(f) enter initial guess x1=1.5 allowable tolerance tol=0.01 maximum number of iterations max=3 zero not found to desired tolerance step x 1.0000 1.5000 2.0000 -0.1250 3.0000 -21.3770 ans = 1.5000 -0.1250 -21.3770 f=inline(20/(x^2+10),x); iteration(f) enter initial guess x1=1.5 allowable tolerance tol=0.0000005 maximum number of iterations max=20 iteration method has converged step x 1.00000000000000 1.50000000000000 2.00000000000000 1.63265306122