牛顿法与弦截法的上机报告.doc

Word文档，可编辑修改！ PAGE 理学院数值分析实验报告 实验名称 牛顿法和弦截法 实验时间 2013年3 月28 组长 11级数学与应用数学班 学号 11107100112 成绩 组员 学号 11107100107 学号 11107100108 学号 11107100109 学号 11107100111 学号 11107100119 一．实验目的，内容 二．相关背景知识介绍 三．代码 四．数值结果 五．计算结果分析 六．计算中出现的问题，解决方法及体会 一 实验目的，内容 1.了解牛顿法和弦截法的原理和主要思想； 2.掌握牛顿法和弦截法的计算步骤。 牛顿法的基本思想 把非线性函数f(x)在 处展开成泰勒级数 f(x)=f( )+(x- )f′( )+(x- ) + … 取其线性部分，作为非线性方程f(x)=0的近似方程，则有 f( )+(x- ) f′( )=0 设f′( )≠0?，则其解为x = - (1) 再把f(x)在x 处展开为泰勒级数，取其线性部分为f(x)=0的近似方程，若 f′(x ) ≠0，则得x = - 如此继续下去，得到牛顿法的迭代公式：x = - ...(n=0,1,2,…) (2) 几何意义 方程f(x)=0的根就是曲线y=f(x)与x轴交点的横坐标x*，当初始近似值x0选取后，过( x0,f(x0))作切线，其切线方程为：y- f(x0)=f′(x0)(x-x0) 它与x轴交点的横坐标为x 一般地，设 是x*的第n次近似值，过( x,f(x))作y=f(x)的切线，其切线与x轴交点的横坐标为：x = - 即用切线与x轴交点的横坐标近似代 曲线与x轴交点的横坐标 牛顿法正因为有此明显的几何意义，所以也叫切线法。 弦截法是一种求方程根的基本方法，在计算机编程中常用。 他的思路是：任取两个数x1、x2，求得对应的函数值f(x1）、f(x2）。如果两函数值同号，则重新取数，直到这两个函数值异号为止。连接（x1,f(x1））与（x2,f(x2））这两点形成的直线与x轴相交于一点x，求得对应的f(x），判断其与f(x1）、f(x2）中的哪个值同号。如f(x）与f(x1）同号，则f(x）为新的f(x1。将新的f(x1）与f(x2）连接，如此循环。体现的是极限的思想。 二 相关背景知识介绍 牛顿法最初由艾萨克·牛顿于1736年在 Method of Fluxions 中公开提出。而事实上方法此时已经由Joseph Raphson于1690年在Analysis Aequationum中提出，与牛顿法相关的章节《流数法》在更早的1671年已经完成了。 要求熟练掌握牛顿迭代公式，并能针对不同的非线性方程构造其牛顿迭代公式及判别其收敛性；能正确应用所学方法求出给定的非线性方程满足一定精确度要求的数值解。 三 代码 1、主要迭代函数 clear %迭代的本次要求 esp=10^-5; b=-2; a=？;（取1为牛顿 取2为截） %迭代的程序核心 b=b1-esp-0.01; i=0; while(abs(b-b1)esp) i=i+1; n=b; b=f(b,b1,a); b1=n; end f(b,b1,a) i 2、迭代的方法部分 %牛顿法与弦截法的迭代函数 function y=f(x,k,a) switch a case 1 y=x-(g(x)/gd(x)); %当a=1的时候应用函数f为牛顿法 case 2 y=x-(g(x)*(x-k))/(g(x)-g(k)); %当a=其他时候应用函数f为弦截法 end end 3、目标函数 function y=g(x) y=x^4-2*x^3-4*x^2+4*x+4; end 4、目标函数导数 function y=gd(x) y=4*x^3-6*x^2-8*x+4; end 四 数值结果 i = 6 ans = -1.4142 i = 5 ans = -0.7321 五 计算结果分析 计算结果，但因为老师在给题目的时候已经告诉我们所给的方程中所要求解的函数肯定是收敛的，所以我们组再一开始就没有再次去验证函数的收敛与否。计算的结果告诉我们，牛顿法和弦截法分别适用于不同的情况。相对于二分法，牛顿法的迭代次数较少，它的收敛速度较快。应用牛顿法要合理选取它的初值。 六 计算中出现的问题、解决方法和体会 在此次牛顿法和弦截法的实验中也遇到了一些问题。我们开始编写程序的时候没有将文件名正确保存，导致了无法运行，后来经过仔细查找原来是这方面出现了问题。 通过这次实验我也充分体会到合作的重要性，往往某