《理学院计算方法课程设计报告》由会员分享，可在线阅读，更多相关《理学院计算方法课程设计报告（18页珍藏版）》请在金锄头文库上搜索。
1、计算方法计算方法课程设计报告课程设计报告学学生生姓姓名名:杜珊珊学学 号：号：2012309010105李南飞2012309010116尹明仁2012309010117学学 院院 :理学院班班 级级:信计 121题题 目目:牛顿法、弦截法及抛物线法求解非线性方程指导教师：指导教师： 周 硕 职称：职称： 教授 朱振菊 实验师 2014 年 12 月 30 日-目 录摘要1一、选题背景21.1 背景介绍 2二、算法设计22.1 理论基础22.2 算法流程图4三、程序设计及功能说明63.1 算法功能设计63.2 算法实例6四、结果分析. 74.1 结果分析74.1.1 例题 1 结果分析 74.1
2、.2 例题 2 结果分析 84.1.3 例题 3 结果分析 94.2 误差分析115、总结及心得体会. 125.1 总结125.2 心得体会12参考文献.13源程序.14- 0 -摘 要众所周知，代数方程求根问题是一个古老的数学问题。早在 16 世纪就找到了三次、四次方程的求根公式。但直到 19 世纪才证明了 n5 次的一般代数方程是不能用代数公式求解的，或者求解非常复杂。因此需要研究用数值方法求得满足一定精度的代数方程的近似解。在工程和科学技术中许多问题常归结为求解非线性方程的问题。正因为非线性方程求根问题是如此重要的基础，因此它的求根问题很早就引起了人们的兴趣，并得到了许多成熟的求解方法。
3、本文主要介绍非线性方程的数值解法是直接从方程出发，逐步缩小根的存在区间，或逐步将根的近似值精确化，直到满足问题对精度的要求，主要的方法有牛顿法、弦截法和抛物线法；牛顿法是方程求根中常用的一种迭代方法，它除了具有简单迭代法的优点外，还具有二阶收敛速度（在单根邻近处）的特点，但牛顿法对初始值选取比较苛刻（必须充分靠近方程的根） ，否则牛顿法可能不收敛。弦截法是牛顿法的一种修改，虽然比牛顿法收敛慢，但因它不需计算函数的导数，故有时宁可用弦截法而不用牛顿法，弦截法也要求初始值必须选取得充分靠近方程的根，否则也可能不收敛。本文还写出这几种非线性方程的数值解法的算法步骤和例题，最后用三种方法进行计算，得出
4、结果并进行比较。关键词 牛顿法；弦截法；抛物线法；非线性方程; - 1 -一、选题背景1.1 背景介绍在利用数学工具研究社会现象和自然现象，或解决工程技术等问题时，很多问题都可以归结为非线性方程 f(x)=0 的求解问题，无论在理论研究方面还是在实际应用中，求解非线性方程都占了非常重要的地位。迭代法是求解非线性方程 f(x)=0 根的一种最重要的方法，而迭代法的优劣对于非线性问题求解速度的快慢和结果的好坏都有很大的影响，所以从实际出发，进行高计算效能迭代算法的研究具有重要的科学价值和实际意义。二、算法设计2.1 理论基础(1) 牛顿迭代法牛顿迭代法是一种特殊的不动点迭代法,其计算公式为1(),
5、0,1,2,.()k kk kf xxxkfx(2.1)其迭代函数为 ()( )()f xxxfx(2.2)牛顿迭代法的收敛速度，当( *)0,( *)0,( *)0f xfxfx时，容易证明,( *)0fx， ( *)( *)0( *)fxxfx ，牛顿迭代法是平方收敛的，且 1 2( *)lim2( *)kkkefx efx 。 (2.3) 牛顿法几何意义方程 f(x)=0 的根 x*是曲线 y=f(x)与 x 轴交点的横坐标,设是根 x*的某个近似值,过曲kx线 y=f(x)的横坐标为的点引切线交 x 轴于 , 并将其作为 x*新的近似kx)(,(kkkxfxp 1kx值,重复上述过程,
6、可见一次次用切线方程来求解方程 f(x)=0 的根,所以亦称为牛顿切线法。图 2-1 牛顿法几何意义图- 2 -（2）弦截法 将牛顿迭代法中的()kfx用( )f x在1kx，kx处的一阶差商来代替,即可得弦截法：11 1()()()()k kkkk kkf xxxxxf xf x (2.4)弦截法几何意义弦截法也称割线法,其几何意义是用过曲线上两点，的割)(,(000xfxP )(,(111xfxP 线来代替曲线,用割线与 x 轴交点的横座标作为方程的近似根。再过点和点2x1P作割线求出,再过点和点作割线求出,余此类推，当收)(,(222xfxP 3x1P)(,(333xfxP 4x敛时可求
7、出满足精度要求的。kx图 2-2 弦截法几何意义图（3）抛物线法弦截法可以理解为用过11(,(),()kkkkxf xxf x两点的直线方程的根近似替 ( )0f x 的根。若已知( )0f x 的三个近似根kx，1kx，2kx用过1122(,(),(,(),(,()kkkkkkxf xxf xxf x的抛物线方程的根近似代替( )0f x 的根，所得的迭代法称为抛物线法,也称密勒(Muller)法。（4）确定有根区间的方法为了确定根的初值，首先必须圈定根所在的范围，称为圈定根或根的隔离。在上述基础上，采取适当的数值方法确定具有一定精度要求的初值。对于代数方程，其根的个数（实或复的）与其次数相
8、同。至于超越方程，其根可能是一个、几个或无解，并没有什么固定的圈根方法。求方程根的问题，就几何上讲,是求曲线 y=f (x)与 x 轴交点的横坐标。由高等数学知识知, 设 f (x)为区间a,b上的单值连续, 如果 f (a)f (b)0x3=x2-(2*f2/(w+sqrt(w2-4*c*f2);enderr=abs(x3-1.87938524);x0=x1;x1=x2;x2=x3;y=feval(f,x2);fprintf(k=%.0f x%d=%.8f e%d=%.6e y%d=%.8en,k,k,x2,k,err,k,y)if (erre)|(y=0)|(k=max)break;endend
