第六章-非线性方程求根.pptxVIP

会计学;问题驱动：全球定位系统（GPS）;图卫星定位示意图;假设接收的信息如表所示。请设法确定R点的位置。 ;;1．根的存在性。方程有没有根？如果有，有几个根？;1.根的存在性;m重根;2. 根的搜索;;? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ;例1：考察方程;; 执行步骤;4、当;;用二分法求解，要使 ;①简单; ② 对f (x) 要求不高(只要连续即可) .; ? 问题 ? 虽然二分法计算简单，能够保证收敛，但是它对于方程单根存在区域信息要求太高，一般情况下很难实现，并且不能求??重根、复根和虚根。在实际应用中，用来求解方程根的主要方法是迭代法。;1．简单迭代法;解：从方程可以看出1是方程的一个根，在[0,1.5]之间。;依此类推,得 x3 = 0.9940 x4 = 0.9990 x5 = 0.9998 x6 = 1.0000 x7 = 1.0000;;定理6.1：如果? (x)满足下列条件 （1）当x?[a, b]时，?(x)?[a, b] （2）当任意x?[a, b]，存在0 L 1，使 则方程x = ? (x)在[a, b]上有唯一的根x*,且对任意初值 x0?[a, b]，迭代序列xk+1=? (xk) (k = 0, 1, …)收敛于x*。;由递推公式： 可以得到：;3．迭代法的误差估计;②;求方程;由迭代格式 （1） ;由迭代格式 （2） ; 由迭代格式（3） ;迭代格式（1）的迭代函数为 ; 迭代格式（2）的迭代函数为 ;知当;迭代格式（3）的迭代函数为 ;由;结论:;迭代法的计算步骤(P145页)：;4．迭代过程的局部收敛性; ;;可取 ;(2). 埃特金加速公式;;例5 试用简单迭代法和Steffensen加速法求方程;Aitken加速公式;§3 牛顿法;将 (x* ? x0)2 看成高阶小量，则有：;x*;定理：设f (x)在[a, b]上存在二阶连续导数且满足下列条件： （1）f (a) ? f (b) 0； （2）f’ (x) ? 0； （3）f? (x) 在区间[a,b]上不变号； （4）取x0? [a, b]，使得f? (x)?f (x0) 0 则由（）确定的牛顿迭代序列{xk}二阶收敛于f (x)在[a, b]上的唯一单根x*。;第50页/共95页;注：Newton法的收敛性依赖于x0 的选取。;2.迭代过程的收敛速度;3. 迭代法收敛阶的判别;4. 牛顿迭代法的局部收敛性与收敛速度;例6 用牛顿法求方程的根：;;用牛顿迭代格式;注：牛顿迭代法的优缺点:;6. 牛顿法应用举例;例8 求;②求倒数算法：;四、求m重根的牛顿法-修正牛顿法;修正格式二：;由于Newton迭代法的收敛性依赖于初值 ;2.牛顿下山法的计算步骤：;【例9】求方程 f(x) = x3 – x – 1 = 0 的根。 x0 =0.6;牛顿迭代法每迭代一次都需计算函数值 ;（1）割线法的几何意义;（2）割线法的局部收敛性与收敛速度;2、 双点弦截法 ：;双点弦截法的局部收敛性与收敛速度;计算结果;双点弦截法的计算步骤：;3.抛物线法;;（3）抛物线法计算公式的推导：;;例11 用抛物线法求解方程 f(x)=xex-1=0. ;抛物线法的迭代误差有下列渐近关系式; 如果f(x)是多项式，那么称 f(x)=0为代数方程。前面介绍的方程求根方法理论上适用于代数方程。但由于代数方程的特殊性，还有一些更为有效的方法。;我们可令;与右式比较可得：; 对于多项式 其中系数 均为实数。;§6 非线性方程组的数值解法; ;;图;;;使用Matlab求解得迭代4次就可以得到相当精确的结果。 ;历史与注记 ; 关于特殊的非线性方程求根问题的迭代法最早出现在古希腊、巴比伦和印第安人的著作中。牛顿法的只有一部分属于牛顿本人，1669年牛顿第一次提出了与现在牛顿法基本等价的方法，但令人惊讶的是该方法并没有使用导数，而是基于二项展开式，因此只适用于多项式。1690年，拉弗森对牛顿法作了简化和改进，称为牛顿—拉弗森法。在牛顿法中使用导数是由辛普森1740年首次提出的，并将其从一维空间推广到多维空间，这才是现在所称的牛顿法。1798年，拉格朗日提出了牛顿法简单而精巧的表达式，并在1831年由傅立叶作了推广。非线性方程组的大多数方法都采用了牛顿法的设计原理，并以牛顿法为度量标准。“牛顿法”已成为将不同类型非