实验3_求代数方程的近似根.pdfVIP

Matlab 实验 求代数方程的近似根 HSU_HJW 问题背景 • 求代数方程f (x ) 0 的根是最常见的数学问 题之一 f (x ) 0 f (x ) • 当 是非线性方程时，由于 的多样 性，尚无一般的解析解法可使用 • 如果对任意的精度要求，能求出方程的近 似根，则可以认为求根的计算问题已经解 决，至少能满足实际要求 实验目的 • 介绍一些求方程实根的近似值的有效方法 –要求在使用这些方法前先确定求根区间[a,b] ， 或给出某根的近似值x0 • 在实际问题抽象出的数学模型中 x – 可以根据物理背景确定 0 – 可根据y f (x ) 的草图等方法确定 – 可用对分法、迭代法以及牛顿切线法大致确定 根的分布情况 实验目的 • 通过本实验希望你能 – 了解对分法、迭代法、牛顿切线法求方程近似 根的基本过程 – 求代数方程（组）的解 相关函数 • abs( ) ：求绝对值函数 • diff( ) ：求微分，f 为符号表达式 • roots( ) ：求解多项式的所有根 • solve( ) ：求表达式的解 • linsolve(A, b) ：求线性方程组A*x=b 的解 • fzero(fun, x0) ：在x0 附近求fun 的解 • subs(f, ‘x ’,a) ：将a 的值赋给符号表达式f 中 的x ，并计算出值 实验内容 • 对分法 • 迭代法 • 牛顿切线法 对分法 • 对分法思想：将区域不断对分，判断根在 某个分段内，再对该段对分，依此类推， 直到满足精度为止 • 对分法适用于求有根区间内的单实根或奇 重实根 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 迭代法 • 迭代法的思想：由方程f (x ) 0 构造一个等 x (x ) x 价方程： . 从某个近似根 出发，令 0 x (x ), k 0,1,2, k 1 k 可得序列{xk } {x } x * x (x ) • k 的极限 是 的根 x • 当然，若 发散，迭代法就失败