数值计算设计(非线性方程求根)案例分析.docxVIP

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