数值分析作业讲解.doc

数 值 分 析 作 业 ——非线性方程的求解方法与分析 学院： 学号： 姓名： 摘 要 本文主要阐述了五种非线性方程的求解方法，分别为二分法、简易牛顿法、牛顿迭代法、牛顿下山法与弦截法。并分别对五种求解方法的计算结果进行了相应地分析。二分法运用函数有根区间中点与端点的函数值，缩小根区间，从而得到较快的收敛速度。牛顿迭代法，是一种常见的求解具有单重零点的非线性方程的数值方法，具有局部二阶收敛性。简易牛顿法便是简化的牛顿迭代法，将迭代点的导数值固定为初始值点的导数值，从而简化计算次数。牛顿下山法，为避免初值选取不当而使得迭代不收敛而在牛顿迭代法改进的方法。弦截法，克服了牛顿迭代法需求零点处函数导数的缺点，使用两次迭代点的差商替代了函数的导数值。本文非线性方程的求解方法均运用MATLAB编程及实现。 关键词：非线性方程；二分法；牛顿迭代法；牛顿下山法；弦截法 目 录 第一章 非线性方程 1 1.1 非线性方程简介 1 1.2 非线性方程求解方法简介 1 1.2.1 二分法 1 1.2.2 牛顿迭代法 2 1.2.3 牛顿下山法 4 1.2.4 简易牛顿法 4 1.2.5 弦截法 5 第二章 计算机配置 7 2.1 处理器 7 2.2 存储设备 7 2.3 显卡 8 2.4 显示屏 8 2.5 操作系统 8 第三章 算法的MATLAB实现及结果分析 9 3.1 二分法 9 3.2 牛顿迭代法 12 3.3 简易牛顿法 14 3.4 牛顿下山法 17 3.5 弦截法 20 结 论 24 第一章 非线性方程 1.1 非线性方程简介 非线性方程，就是因变量与自变量之间的关系不是线性关系。 在永恒变化发展的自然界与人类社会中，在研究其内部规律的各个科学领域中，更深刻、更精确地描述其内部规律的数学工具之一，就是非线性方程。非线性代数是研究大规模离散数据的运算处理与内在性状的数学科学。科学技术离不开数据处理与数据分析，因此非线性代数具有非常广泛的应用，在力学、化学、生命科学、控制理论等众多科学领域中，非线性方程早已屡见不鲜。因此，非线性方程的求解就显得愈加重要。然而求解非线性方程有很多种方法，每种方法都有自己的优缺点。 1.2 非线性方程求解方法简介 求函数零解作为数学研究领域的一个热点已经延续了几百余年，所以已经建立了许多种方法，拥有比较完备的求解体系。本文中，主要介绍非线性方程求解方法中最常用也是比较简单的几种方法。 在解决实际问题的中，大都会遇到非线性方程或非线性方程组的数学模型，这类方程的求解用一般的代数方法求解是不可能实现的。所以，在解决这类问题的时候，多是将求零解转化为求近似解。 1.2.1 二分法 若是区间上的连续函数，且，则在内必有一个零点。因为，所以函数在区间上改变符号，因此它在这个区间内至少存在一个零点。二分法就是利用这一中值定理来求解非线性方程零解。 二分法求解的具体方法：若，则计算区间中点，并且检验是否为真。若为真，则在内有零点。因而把中点设为作为区间新的右极点。若检验为假，则在区间内有零点，因而把中点设为作为区间新的左极点。这样新的区间的宽度就为原区间宽度的二分之一。并在此区间中重复上述操作。 当然，若，则从而求出一个零点。然而由于舍入误差的存在，在计算机计算的过程中，精确为0是完全不可能存在的。因此，主卧室算法循环的停止判断准则不应该是是否成立，而必须提供一个合理的允许误差。当计算结果的值在误差范围内，便可停止运算。 1.2.2 牛顿迭代法 牛顿法迭代法是一种能在许多不同情况下应用的通用过程。特别地，当用牛顿法来求实值变量函数零点时，常常被称为牛顿-拉弗森迭代。通常，牛顿迭代法比二分法与弦截法获取答案的速度要快，这是因为它的收敛是二次的而不是线性或者超线性的。一旦二次收敛变得有效时，即牛顿法序列的值充分地接近根时，其收敛是如此之快以致于仅仅再需要几个数值即可。但是，牛顿迭代法并无法保证总是收敛的。所以牛顿法经常与其他较慢的方法结合形成一种数值上整体收敛的混合方法。 若存在一个函数，其零点由数值方法计算得出。设是的零点，而是的一个近似，若的阶导数存在并且连续，则由泰勒定理将函数在零点处进行展开可得： 其中。若较小（即在附近），则可以略去项，并且在余下的方程中求。由此可得到结果是。若使的一个近似，则应该是的一个更好的近似。牛顿迭代法从的一个估计开始，则归纳出迭代的格式为 下面叙述一下牛顿迭代法的几何意义。 是的根，选取作为的初始近似值，经过上的点做的切线方程：，求出与横轴焦点的横坐标，则称为的一次近似值。将作为下一次迭代的初值，重复上述过程可得到的二次近似值。如此循环，