非线性方程的迭代法.pptVIP

#include stdio.h #include stdlib.h #include math.h #define TRUE 1 void main() { int it_limit, it; float a, b, c, epsilon, Y_a, Y_b, Y_c; float fun_f(); printf( Lower bound: a ? ); scanf( %f, a ); printf( Upper bound: c ? ); scanf( %f, c ); epsilon=1e-8; it_limit=100; it = 0; Y_a = fun_f( a ); Y_c = fun_f( c ); if( Y_a*Y_c 0 ) { printf( f(a)f(c) 0, break. \n ); exit(0); } else { while( TRUE ) { it = it + 1; b = (a + c)/2; Y_b = fun_f( b ); if( it it_limit ) break; if( fabs( b - a ) epsilon ) break; if( Y_a*Y_b = 0 ) { c = b; Y_c = Y_b; } else { a = b; Y_a = Y_b; } } } if ( itit_limit ) printf( Iteration limit exceeded.\n ); printf( Final result: Approximate root = %g \n, b ); } } float fun_f(x) /* Definition of Function */ float x; { float f; f = x*x*x - 3*x*x - x + 3; return( f ); 作业: 用二分法求解 [实验目的]掌握二分法,理解其特点. 2. 用不同形式的Newton法求解 假设不知道Newton法的收敛阶,如何用数值实验估计? [实验目的]掌握Newton法,对重根的处理,收敛阶的估计. 3. 用Newton法解非线性方程组 [实验目的]掌握非线性方程组的数值解法,内迭代要求解线性代数方程组. 4. 用Newton法解复数方程 在平面上把收敛到不同根的初始值分别标上不同颜色. [实验目的]了解Newton法收敛域的结构和其局部收敛性. Hamming’s Motto “The purpose of computing is insight, not number.” * * 非线性方程的迭代法 向 华 武汉大学数学与统计学院 (Sumerians, third millennium BC, n10实际中很少,但在computer algebra中n100,且系数高精度) Kepler方程: 取x0作为初始近似值，将f(x)在x0做Taylor展开: 作为第一次近似值 牛顿法及其几何意义 基本思路：将非线性方程f(x)=0 线性化 重复上述过程 ? Newton法每一步用不同的切线方程去逼近非线性方程,因此也称为切线法,它是一种将非线性方程线性化的方法. 若定义 ,则可以写成定点迭代的格式 何时收敛? 局部收敛性 思考 如何判断收敛? 收敛的快慢? 收敛的阶 停机准则 设{xk}为收敛于f(x)=0的解x*的序列,记 ek=|xk-x*|为第k次迭代误差.若存在某个实数 p=1及正数C,使 则称序列{xk}为p 阶收敛. 牛顿法的收敛性与收敛速度 当p=1且0C1时,称为线性收敛; 当1p2时,称为超线性收敛; 当p=2时,称为平方收敛. 数p的大小反映了迭代法收敛速度的快慢, p越大收敛越快. 因此迭代法的收敛阶是衡量迭代法优劣的重要标志. （局部收敛性定理） 设 f (x)?C2[a, b]，若 x*