非线性方程及5种数值解法及其.pptVIP

非线性方程的5种数值解法及其探究 专业:信息与计算科学 摘要：本文首先分别介绍了二分法、newton迭代法、反函数法、求交法、反插值法等5种算法;再结合具体的实例,并在相同的误差精度下,用matlab语言对它们逐一实现;最后对这5种方法进行归纳总结. 关键词：非线性方程;二分法;迭代法;反函数法;反插值法; matlab程序; 论 文 结 构 框 架 1 引言 在实际问题中,求解非线性方程根的精确值很困难, 大部分的情况下,我们只需要求解出近似值即可.而数值解法,就是用数值迭代的方法来求解近似值的一种方法. 其中最早提出来的是二分法. 到了17世纪,牛顿又提出了一种迭代方法,这是目前为止 最基本,最核心,最经典的一种迭代法. 在20世纪,有学者提出了反函数法,并推出一系列具有高精度的求根公式,这对于那些不便于进行迭代的方程的求解具有较大的应用价值. 此外,还有学者在不动点迭代法的基础上,利用切线的几何意义来构造迭代公式,称为求交法. 同时,一些专家通过反插值来求解出自变量的值,即反插值法. 2 相关理论知识 · 二分法 首先确定方程的有根区间,然后再进行中点函数值判断,这样不断地接近零点,把零点所在的区间缩小一半,最终使区间的两个端点逐渐收敛于函数的零点. 那么,该方程的解就是: 当 4 算例分析 谢谢观看！ * * * 指导老师:魏春艳(讲师) 姓名:张萨 引言 相关理论知识 及算法步骤 算例分析 综合分析比较 研究意义 相关领域研究回顾 介绍了这5种方法的基本 原理及算法步骤 以方程 为例, 用matlab程序分别实现 分析比较,归纳其应用 范围和优缺点 取上式中的线性部分,并考虑 ,则有： 设 ,那么该方程的解就是： 再把 在 附近展开成泰勒级数,并且也取其线性部分有: ·netwon迭代法 由泰勒定理知： 这样,就可以得到该方法的一个迭代序列为: ·反函数法(该方法仅适用于高次代数方程,且最高次系数设为n ) 若函数存在单值的反函数,设为 ,并且有： 时,用泰勒展开式把 在点 处展开,有： 令 ,当 充分小时,将上式中的拉格朗日余项略去， 当 时,同理可得:迭代公式是： 那么可得迭代公式是： 把前面的3个式子代入到上式中,就可以得到： 有： 时,迭代公式是： 当 时可以得到更高精度的迭代公式. ·求交法 x y 0 k 首先把方程化成它的同解方程 ,然后令 , 方程的根,就等价于求这两个 ,求解 函数的交点k. 设初始值是 , x0 x x1 则可得到 ,过该点作直线的切线,得到 ,该点的横坐标 记为 ,可以得到三角形 ,那么可得: 因为 ,则 如此这样进行下去,就会得到一组迭代公式是: ·反插值法 给定方程的n个离散数据点: ,且 是介 , 之间的某个值 ,由介值定理可知:它所对 应的自变量的取值必在 之间,考虑在该区间内进行 线性插值,则有: 再令 ,就可以得到可以得到一个粗略近似值: 于 令 ,则可得到第一个近似值: 按照同样的方法进行二次插值 ,可以得到第二个近似值: 然后把 作为初始值, 考虑进行n次反插值,则有: 如此进行下去,就可以得到该方法的一组迭代序列是： 3 算法步骤 给定有根区间,初始值,并考虑在相同的容许误差精度要求 下,分别用上述的5种方法编写具体的matlab程序来求解 方程根的近似值. 以方程 为例,先用matlab程序求出它的有根 区间是 ,然后再选取初始值 和精确度 最后用matlab语言对这5种方法逐一实现,求解出该方程 根的近似值,并要求能得到每一步迭代的结果.(具体程序 5 综合分析比较 分别从基本原理,迭代次数,迭代时间,适用范围,优缺点 等方面,对这5种方法逐以分析,并举例说明,最后再进行 见附录). 综合归纳比较(见下表). -5.80383649910152 0.036秒 5 反插值法 -5.80383649910152 0.006秒 4 求交法 -5.80383649910152 0.01秒 3 反函数法 -5.80383649910152 0.01秒 4 netwon迭代法 -5.80383649934083 0.015秒 33 二分法 数值解 迭代时间 迭代次数 初始值 表1:最终的迭代结果比较 表 2:5种数值解法的综合比较 函数在有根区间上连续,且在区间端点处的函数值异号； ①收敛速度比较慢； ②只能求解奇数重根,不能求解偶数重根； ①算法简单,容易操作；②近似根序列一定收敛于方程的根； 二分法 ①在整个有