数学实验报告样本.docxVIP

数学实验报告 实验序号： 3 日期： 2013 年 12 月 14 日 班级 应数一班 姓名 陈菲 学号 1101114209 实验 求代数方程的近似根 名称 问题背景描述： 求代数方程 f ( x) 0 的根是最常见的数学问题之一， 当 f ( x) 是一次多项式时， 称 f (x) 0 为线性方程，否则称之为非线性方程． 当 f (x) 0是非线性方程时，由于 f (x) 的多样性，尚无一般的解析解法可使用，但如果 对任意的精度要求，能求出方程的近似根，则可以认为求根的计算问题已经解决，至少能满足实际要求． 本实验介绍一些求方程实根的近似值的有效方法， 要求在使用这些方法前先确定求根区间 [ a,b] ，或给出某根的近似值 x0 ． 实验目的： 了解代数方程求根求解的四种方法：对分法、迭代法、牛顿切线法 掌握对分法、迭代法、牛顿切线法求方程近似根的基本过程。 １ 实验原理与数学模型： 1．对分法 对分法思想：将区域不断对分，判断根在某个分段内，再对该段对分，依此类推，直到满足精度为止．对分法适用于求有根区间内的单实根或奇重实根． 设 f (x) 在 [ a, b] 上连续， f (a) f (b) 0 ，即 f (a) 0 ， f (b) 0 或 f (a) 0 ， f (b) 0 ．则根据 连续函数的介值定理，在 ( a, b) 内至少存在一点 ，使 f ( ) 0 ． 下面的方法可以求出该根： x0 a b 令 2 ，计算 f ( x0 ) ； 若 f ( x0 ) 0 ，则 x0 是 f ( x) 0 的根，停止计算，输出结果 x x0 ． 若 f (a) f ( x0 ) 0 ，则令 a1 a ，b1 x0 ，若 f (a) f (x0 ) 0 ，则令 a1 ak bk ??，有 ak 、 bk 以及相应的 xk 2 ． f ( xk ) xk (3) 若 ( 为预先给定的精度要求 )，退出计算，输出结果 反之，返回 (1)，重复 (1)，(2)，(3)．  a1 b1 x0 ，b1 b x1 2 ． ； ak bk ； 以上方法可得到每次缩小一半的区间序列 {[ ak , bk ]} ，在 ( ak ,bk ) 中含有方程的根． ak bk 当区间长 bk ak 很小时，取其中点 xk 2 为根的近似值，显然有 xk 1 1 1 ak 1) 1 (b a) (bk ak ) (bk 1 k 1 2 2 2 2 以上公式可用于估计对分次数 k ． 迭代法 迭代法的基本思想： 由方程 f ( x) 0 构造一个等价方程 x ( x) 从某个近似根 x0 出发，令 xk 1 (xk ) ， k 0,1,2, 可得序列 { xk} ，这种方法称为迭代法． 若 { xk} 收敛，即 lim xk x* k ， 只要 (x) 连续，有 lim xk 1 lim ( xk ) (lim xk ) k k k 即 ２ 可知， { xk } 的极限 x* x* (x* ) 是 x( x) 的根，也就是 f ( x) 0 的根． 当然，若 xk 发散，迭代法就失败． 迭代过程 xk 1 (xk ) 收敛的常用判别标准： 当根区间 [ a,b] 较小，且对某一 x0 [ a, b] ， ( x) 明显小于 1 时，则迭代收敛 迭代法的加速： a) 松弛法： 若 (x) 与 xk 同是 x* 的近似值，则 xk 1 (1 k ) xk k ( xk ) 是两个近似值的加权平均， 其中 k 称 为权重，现通过确定 k 看能否得到加速． 迭代方程是： x ( x) 其中 ( x) (1 ) x ( x) ，令 (x) 1 ( x) 0 ，试确定 ： 1 k 1 1 k ( xk ) 当 ( x) 1 时，有 ( x) ，即当 1(xk ) ， (xk ) 时， 1 1 k 1 可望获得较好的加速效果，于是有松弛法： xk 1 (1 k ) xkk ( xk ) ， 1 (xk ) b) Altken 方法： x* ( x* ) ， x* 是它的根， x0 是其近似根． 设 x1 (x0 ) ， x2 (x1) ，因为 x* x2 [ x* x2 ] x2 [ (x* ) (x1)] x2 ( )( x* x1 ) ， x2 x1 ( x1 ) ( x0 ) 用差商 x1 x0 x1 x0 近似代替 ( ) ，有 x* x2 x2 x1 ( x* x1 ) x1 x0 ， 解出 x* ，得 * ( x x ) 2 x x2 2 1 x2 2x1 x0 由此得出公式 xk(1) ( xk ) ； xk