计算方法第二章new.pptVIP

浙江大学研究生学位课程 《实用数值计算方法》 第四章 第二章 非线性方程与方程组的数值解法 求解非线性方程的一般步骤 用扫描法或绘图对给定区间进行搜索，确定出各个仅存在单根的区间，此时区间内任意一点的值均可视为根的粗略位置或根的初值。 选用一种迭代方法使根精确化到所要求的精度。 方法：二分法、牛顿迭代法和弦截法、求方程全部根（实根和复根）的二次插值法。 导数变化对算法的影响图 作　业 * 2.1　引言 2.2　非线性方程的数值解法 　2.2.1　二分法 　2.2.2　简单迭代法 　2.2.3　牛顿法 　2.2.4　弦截法 　2.2.5　抛物线法 2.3　解非线性方程组的牛顿法 2.1　引言 非线性方程记为：f(x) = 0，包括： ⑴高次代数方程：含有变量的乘积、乘方，例如： x4?10x3+35x2?50x+24=0 ⑵超越方程：含有变量的三角函数、指数函数如sinx、cosx、ex、lnx等。例如： ex?sin(?x/2)=0 非线性方程的根（解或函数的零点）可能是实数或复数。 若f(p)=0，f’(p)?0，则p称为f(x)=0单根；若f(p)= f’(p)=…= f(m)(p)=0而f(m+1)(p) ?0， 则称p为f(x)的 m 重根。 求解方程的根一般需考虑三个问题： 在整个求根区间方程根的存在性。 　代数方程：n次方程的全部根共n个，当n为正奇数时，至少有一个根是实根； 　超越方程：根的个数可能只有一个，也可能有多个以至无穷多个，或没有根。 根的粗略位置。 知道根的足够好的近似值是求非线性方程解的重要条件。 如果要求确定出给定区间内方程的全部根，而且不知道根的分布情况，可采用以下两种方法对区间内存在的根进行搜索，确定出根的个数和其粗略位置，防止漏根。 (1)扫描法：从端点x=a出发，按预定步长h一步步地向右跨，每跨一步进行一次根的“扫描”，检查函数f(x0)和f(x0+h)的值是否同号。若f(x0)? f(x0+h)?0则说明所求根x*必在x0与x0+h之间，这时可取x0或x0+h作为根的初始近似值。 　　进行根区间扫描时遇到的困难是： (1)确定出合适的扫描区间（如果未给定）； (2)搜索步长h。 　　选择区间过宽，扫描计算工作量大；若区间不够宽或不合适，则可能漏掉一些根或找不到根。又如果方程有二个根靠得很近，选取步长过大也会使根被漏掉。为防止漏根，保守做法是使用较大的扫描区间和较小的步长。 根的精确化 已知一个根的近似值后，设法把根精确化，直到足够精确为止。 对于非线性问题，除少数情况，如二次三项式f(x)=ax2+bx+c=0，一般不能利用公式直接求解，而要采用某种迭代解法。 扫描法确定f(x)=0的根的粗略位置 　　为避免扫描法的盲目性，有效的做法是采用绘图法。 (2)绘图法：给定区间[a, b]和步长h，计算f(a), f(a+h),……, f(a+nh)，其中b=a+nh，并据此绘出f(x)的粗略曲线，观察曲线的走势和过零情况，或将f(x)＝0写成f(x)＝g1(x)?g2(x)=0的形式，观察g1(x)和g2(x)两条曲线相交的情况，交点处的x值即为方程的根。 迭代法：从给定的一个或几个初始近似值x0, x1,…, xr出发，按某种方法构造出一迭代序列（近似值序列）x0, x1,…, xk,…，使该序列收敛于f(x)=0的某个根p ，即 迭代法的收敛性 迭代序列的收敛性与初值的选取范围有关，有时也取决于方程的性态。 大范围收敛：从任何可取的初始值出发都能保证收敛。 局部收敛：初始值的选取必须充分接近于所要求的根。 迭代过程中，随着k的增加，|?k+1??k|连续增加，则说明此迭代过程发散。 2) 收敛速度一般由迭代方法所决定，方程的性态也会起一些作用。 收敛阶数?（反映收敛速度的快慢） 令 ek = xk ? p C1时， ?值越大，收敛速度越快。 最快的方法?=2，最慢的方法?=0.5。 当?=1，迭代法线性收敛 当?1，迭代法超线性收敛。 局部收敛方法比大范围收敛方法收敛得更快。 若存在实数?和非零常数C，使得 则称该迭代法的收敛阶数为?。 通常，非线性方程的根不止一个，而任何一种方法一次只能算出一个根。因此，在求解非线性方程时，要约定初始值或求解范围。 图 2.2 二分法方程求根 y 0 x a b x0 p 2.2 非线性方程的数值解法 2.2.1　　二分法 介值定理：对于连续函数f (x)，如果在x=a和x=b处异号：f (a)?f (b)0，则f (x)在[a, b] 内至少有一个根。 　　若|f (x0)| ? ?，且f (x0)?f (b)0，则用区间[x0, b]代替[a, b]；若f (x0)?f (b)0，则用区间[a, x0]代替