数值分析实验报告之方程求根演练.docVIP

PAGE 14 PAGE 1 数学与计算科学学院 实 验 报 告 实验项目名称 方程求根演练 所属课程名称 数值方法B 实 验 类 型 验证 实 验 日 期 2013.12.3 班 级 学 号 姓 名 成 绩 一、实验概述： 【实验目的】 1.掌握迭代法和牛顿迭代法求解非线性方程组； 2.掌握方程求解的基本方法迭代法以及牛顿迭代法； 3.研究选取不同迭代函数对收敛性的影响； 【实验原理】 1.迭代公式收敛定理 方程 x = g ( x ）在（ a , b ）内有根 x*，如果： ① 当x?[a,b]时，g( x）?[a,b]; ② g ( x）可导，且存在正数 q 1，使得对于任意x[a,b]都有||q 1. 则有以下结论： ① 方程x = g ( x）在（a , b）内有唯一的根x*。 ② 迭代公式对（a , b）内任意初始近似根均收敛于x*。 ③ 近似根的误差估计公式为： 2.牛顿法收敛定理 设在区间上有二阶连续导数,且满足,在上不变号, 在上不等于0,令 有.则对任意,牛顿迭代格式收敛于在中的唯一实根,并且: ① ② ③ ,牛顿迭代法为2阶收敛。 【实验环境】 1.硬件环境： HP Microsoft 76481-640-8834005-23929 HP Corporation Intel(R) Core(TM) I5-2400 CPU @ 3.10GHz 3.09GHz,3.16GB的内存 2.软件环境： Microsoft Windows XP Professional 版本 2002 Service Pack 3 二、实验内容： 【实验方案】 方案一： 1．验证迭代法求解教材p141例题6.3 求方程在附近的根； 2.牛顿法求解教材P150例题6.7 求方程在取初值的根。 方案二： 1.用迭代法求方程的在初值根；分别选取迭代函数为和求解。分析比较迭代函数选取的不同对收敛性的影响。 2.用牛顿法求在和附近的根，迭代10次。 分析比较初值的选取对迭代法的影响。 【实验过程】（实验步骤、记录、数据、分析） 1方程在附近的根（迭代法） 知初值取1.5，首先我们通过MATLAB的solve函数，求得在附近的根为1.324718。通过迭代，得到下表： 表1 迭代次数与迭代值 迭代次数 迭代值 迭代次数 迭代值 0 1.5 10 1.35721 1 1.35721 11 1.33086 2 1.33086 12 1.32588 3 1.32588 13 1.32494 4 1.32494 14 1.32476 5 1.32476 15 1.32473 6 1.32473 16 1.32472 7 1.32472 17 1.32472 8 1.32472 18 1.32472 9 1.32472 19 1.35721 图 1 牛顿迭代法的迭代次数与迭代值的变化趋势 从表1和图1可以看出，如果仅取六位数字，迭代第7次时已经满足方程的根，而迭代次数与迭代值的变化趋势图表明迭代法在选取合适的初值时有很好的收敛性，而且迭代次数较少的情况下就能得到比较精确地解，所以验证比较成功，迭代法是一种收敛性效果很好的方程求根的数值解法。 2方程在附近的根（牛顿法） 知初值取0.5，首先我们通过MATLAB的solve函数，求得方程在取初值的根为0.567143。通过MATLAB程序的计算，牛顿法的迭代法计算的相应的结果见下表： 表2 牛顿法求解方程的迭代结果 迭代次数 迭代值 迭代次数 迭代值 0 0.5 10 01 011 02 012 03 013 04 014 05 015 06 016 07 017 08 018 09 019 0图2 牛顿迭代法的迭代次数与迭代值的变化趋势 从表2和图2可以看出，迭代次数较少的情况下就能得到比较精确地解，所以验证比较成功，迭代法是一种收敛性效果很好