非線性方程的迭代法.pptVIP

非线性方程的迭代法 向 华 武汉大学数学与统计学院 (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越大收敛越快. 因此迭代法的收敛阶是衡量迭代法优劣的重要标志. X f(x) 收敛性与初始值的选取 inflection point local minimum or maximum near-zero slope 全局收敛性定理(定理3.3.1)：设 f (x)?C2[a, b]，若 f (a) f (b) 0； 在整个[a, b]上 f ?(x) ? 0； f ??(x)在 [a, b]上不变号 选取初始值x0 ? [a, b] 使得 f ??(x0) f (x0) 0； 则由Newton法产生的序列{ xk } 单调地收敛到 f (x)=0 在 [a, b] 的唯一根x*,且收敛速度至少是二阶的 故可用 作为停机准则 需要一个适当的误差估计作为停机准则。因为 其中 在 与 之间，从而有估计 如何判断收敛?(停机准则) 直观上希望计算残量很小: X f(x) 注意事项一: 检查连续两次迭代的差, 防止假收敛 X f(x) 注意事项二: 检查函数值 f(x), 防止假收敛. (1) 选定初值x0 ,计算f (x0) , f ?(x0) 计算步骤 (2) 按公式 迭代 得新的近似值 xk+1 (3) 对于给定的允许精度?，如果 则终止迭代，取 ;否则k=k+1，再转 步骤(2)计算 #include stdio.h #include stdlib.h #include math.h void main() { float epsilon, x, xb, y, y_derivative, error; void func(); printf( Initial guess ? ); scanf( %f, x ); epsilon=1.0e-8; error=1.0e10; xb = x; while( errorepsilon ) { func( x, y, y_derivative ); x = x - y/y_derivative; /* finds new x. */ error=fabs(x-xb); xb = x; } printf( \n-----------------------\n ); printf( Final solution = %g \n, x ); printf( -------------------------\n ); } void func(x, y, y_derivative) /*Computes y and y_derivative */ float x, *y, *y_derivative; { *y = pow(x, 3) - 5.0*pow(x, 2) + 6.*x; *y_derivative = 3.0*pow(x, 2) - 10.0*x + 6.; return; } 1.已知重数m,用修正的Newton迭代公式: 2.未知重数, 定义一个新的函数: 把Newton法用于 , 得到对重根仍然二阶收敛的格式: 重根的处理 3.自适应 Newton 法 2. Newton下山法 Newton 法的变形