数值计算课程设计报告.docxVIP

数值计算方法课程设计报告 课程设计名称： 数值计算方法 课程设计题目： 非线性方程求根 年级专业： 组员姓名学号： 指导教师： 完 成时间： 非线性方程求根 一、 问题提出 随着科学技术，生产力经济的发展，在科学与工程计算中存在着大量方程求根问题， 例如贷款购房问题，工厂的最佳订货问题等都需要求解一类非线性方程的根，首先根据 实际问题列出数学模型，确定变量，给出各个条件及相关函数；然后对建立的模型进行 具体分析和研究，选择合适的求解方法；编写函数的程序，用计算机求出方程的解，通 过所求解分析具体情况. 求解非线性方程的问题有以下几种基本方法。二分法简单易行，但收敛较慢，仅有 线性收敛速度。而且该方法不能用于求偶数重根或复根，但可以用来确定迭代法的初始 值。牛顿法是方程求根中常用的一种迭代方法，它除了具有简单迭代法的优点外，还具 有二阶收敛速度（在单根邻近处）的特点，但牛顿法对初始值选取比较苛刻（必须充分 靠近方程的根），否则牛顿法可能不收敛。弦截法是牛顿法的一种修改， 虽然比牛顿法收 敛慢，但因它不需计算函数的导数，故有时宁可用弦截法而不用牛顿法，弦截法也要求 初始值必须选取得充分靠近方程的根，否则也可能不收敛。 二、 背景分析 代数方程的求根问题是一个古老的数学问题。理论上， n次代数方程在复数域内一 定有n个根（考虑重数）。早在16世纪就找到了三次、四次方程的求根公式，但直到 19 世纪才证明大于等于 5次的一般代数方程式不能用代数公式求解，而对于超越方程就复 杂的多，如果有解，其解可能是一个或几个，也可能是无穷多个。一般也不存在根的解 析表达式。因此需要研究数值方法求得满足一定精度要求的根的近似解。牛顿迭代法是 牛顿在17世纪提出的一种求解方程 .多数方程不存在求根公式，从而求精确根非常困 难，甚至不可能，从而寻找方程的近似根就显得特别重要. 而在各种科学和工程计算中往往要用到非线性方程组的求解，而牛顿法又是最基础 的迭代法，在各种计算力学、控制工程等领域中发挥了不可代替的作用.而在数值计算 中，非线性方程组的求解同样具有重要意义.随着计算机技术的成熟和高速发展，对于 非线性方程求根问题出现了大量的数学软件(如 MATLAB,SAS,SPSS等)，计算机已经成 为工程师应用数学解决工程问题的主要运算工具?同时，工程专业的学生对数学教育的 需求重点正在从手工演绎和运算能力的培养转变到结合计算机软件进行建模、求解和论 证能力的培养?我们采用Matlab数学软件平台，通过实例比较了二分法、 牛顿迭代法、 弦截法三种基本方法的优缺点。 三、基本算法思想与实现 二分法 单变量函数方程： f (x) =0 其中，f(x)在闭区间［a , b］上连续、单调，且f(a)*f(b)0, 则有函数的介值定理可知， 方程f (x) =0在(a, b)区间内有且只有一个解 x*，二分法是通过函数在区间端点的 符号来确定x*所在区域，将有根区间缩小到充分小， 从而可以求出满足给定精度的根 X* 的近似值。 下面研究二分法的几何意义： 设ai=i, bi=b,区间，中点xi= ―bL及f人，若f人=o,则/ =人，若 2 f( a1 )*f( x1 )0 ,令比=內，t2=x1，则根x* ［比，b2］中，这样就得到长度缩小一半 的有根区间［去卫］，若f( b1)*f( x1)0，令去二乂1 , t2 = b1，则根x* ［比亠2］中, 这样就得到长度缩小一半的有根区间 ［%, b2］,即f( a2)f( 3)0,此时鸟-笑=直 也 2 对有根区间［去，b2］重复上述步骤，即分半求中点，判断中电处符号，则可得长度有缩 小一半的有根区间［%,?2］, 如图所示： 新心十‘ 重复上述过程，第n步就得到根x*的近似序列xn及包含X*的区间套，如下： TOC \o 1-5 \h \z \o Current Document [ai,bi] [a2, b2] ?…[an,bn] ■■- * f(an)f(bn) 0, X [an,bn] an- bn」(an 1 bn i) =... = b a \o Current Document v 7 2* i ⑷ Xn 且 |X*-Xn| 瞎 (n=1,2,3 …■■) 2 2 显然lim Xn，且Xn以等比数列的收敛速度收敛于 X* ,因此用二分法求f (X) =0的实根X* 可以达到任意指定精度。 I-(u+A )/2 (2)牛顿迭代法 设方程f(x)=O在其根x*的某个领域U(x* ,)内有一阶连续导数，且 f (X*)工0。 求f(x)=0的根x*，首先要将f(x)=0转化为等价形式 x (x)，并使 (x)满足不动 点迭代的一般理论。 于是我们令 (x)=x+h(x)f(x),