数值分析大作业2014.pdfVIP

课程设计 课程名称： 高等数值计算 设计题目： 数值计算 B 课程设计 学号： 姓名： 完成时间： 2014 年 10 月 20 日 题目一：非线性方程求根 用 Newton 法计算下列方程 （1） x3 x 1 0 x 1 x 0.45 x 0.65 ，初值分别为 0 ， 0 ， 0 ； （2） x3 94x 2 389x 294 0 其三个根分别为 1,3, 98 。当选择初值 x0 2 时给 6 出结果并分析现象 ,当 5 10 ，迭代停止。 一、摘要 非线性方程的解析解通常很难给出，因此非线性方程的数值解就尤为重要。 本实验通过使用常用的求解方法二分法和 Newton 法及改进的 Newton 法处理几 个题目， 分析并总结不同方法处理问题的优缺点。 观察迭代次数， 收敛速度及初 值选取对迭代的影响。 二、数学原理 构造迭代函数的一条很重要的途径是， 用近似方程来代替原方程去求根。 因 此，如果能将非线性方程用线性方程来代替的话，求近似根问题就很容易解决， 而且十分方便。 Newton 法就是把非线性方程线性化的一种方法。 在求解非线性方程 f (x ) 0 时，它的困难在于 f ( x) 是非线性函数，为克服这 一困难，考虑它的线性展开。设当前点为 x k ，在 xk 处的 Taylor 展开式为 f (x ) f (x ) f (x )( x x ) k k k 令 f ( x ) 0 ，可以得到上式的近似方程 f ( x ) f ( x )(x x ) 0 k k k 设 f ( xk ) 0 ，解其方程得到 f (x ) k xk 1 xk ( k 0,1,…) f (x ) k 这就是牛顿迭代公式。用牛顿迭代公式求方程 f (x) 0 根的方法称为牛顿迭 代法。 牛顿迭代法的几何意义为， 不断用切线来近似曲线得到方程的根， 我们知道 方程 f (x ) 0 的实根 x * 是函数 y f (x) 的图形与横坐标的交点， xk 1 是函数 f (x ) 在 点 (x k , f ( xk )) 处的切线与 x 轴的交点，此时就是用切线的零点代替曲线的零点，因 此，牛顿迭代法又称为切线法。 三、程序设计 基于 MATLAB 软件编写程序，先定义一个用 Newton 法求解的功能函数，然后 调用函数用于计算不同的方程。各变量定义见程序。 1、选取初值。 f (