牛顿改进弦截法matlab,Newton迭代法、弦截法上机实验报告 2.doc

您所在位置：网站首页 > 海量文档

&nbsp>&nbsp高等教育&nbsp>&nbsp工学

Newton迭代法、弦截法上机实验报告 2.doc4页

本文档一共被下载：次,您可全文免费在线阅读后下载本文档。

下载提示

1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。

2.该文档所得收入(下载+内容+预览三)归上传者、原创者。

3.登录后可充值，立即自动返金币，充值渠道很便利

中国矿业大学(北京)理学院

数值分析实验报告

实验名称已知数据点,求Lagrange插值多项式.

实验时间2012.4.3组长签名龙纯鹏班级信息与计算科学(1)班学号11107200110成绩组员签名1110720010111107200102111072001031110720011911107200120一、实验目的,内容 二、相关背景知识介绍 三、代码

四、数值结果 五、计算结果的分析 六、计算中出现的问题，解决方法及体会一、实验目的、内容

实验目的：

内容：

已知数据点,求三次Lagrange插值多项式, 并求

相关背景知识介绍

1.牛顿迭代法原理:

设已知方程的近似根,则在附近可用一阶泰勒多项式近似代替.因此, 方程可近似地表示为.用表示的根,它与的根差异不大.

设,由于满足解得

重复这一过程,得到迭代格式

.

2. 牛顿迭代法的几何解析

在处作曲线的切线，切线方程为。令，可得切线与轴的交点坐标，这就是牛顿法的迭代公式。因此，牛顿法又称“切线法”。

3．牛顿迭代法的收敛性

计算可得,设是的单根,有,则

,

故在附近,有.根据不动点原理知牛顿迭代法收敛.

4.程序设计思路：

先用m文件先定义一个名为f1.m的函数文件。

function y = f1(x)

y = ;

再定义一个名为df1.m的函数文件

function y=df1(x)

y=

建立一个主程序prog1.m

clear

newton('f1','df1',1, 10^(-5), 10)

然后在MATLAB命令窗口运行上述主程序，即：

>> prog1

代码(Matlab)function [x2,k]=kkoo(x0,x1,esp)

esp=0.00005;

x0=2;

x1=3;

x2=x1-(x1-2^(-x1))/(1+log(2)*2^(-x1));

k=0;

while abs(x1-x0)>esp

x0=x1;

x1=x2;

x2=x1-f(x1)*(x1-x0)/(f(x1)-f(x0));

k=k+1;

end

x2

k

df1.m:

function [x1,k]=wu(x0,esp)

esp=0.00005;

x0=0.4;

x1=x0-f(x0)/f1(x0);

k=0;

while abs(x1-x0)>esp

x0=x1;

x1=x0-f(x0)/f1(x0);

k=k+1;

end

x1

k

四、

在matlab运算结果如下表：

k

k

k

1

2

3

4

5

0.5000

0.6250

0.6563

0.6484

0.6445

6

7

8

9

10

0.6416

0.6411

0.6412

0.6414

0.6412

11

12

13

14

15

0.6412

0.6412

0.6412

0.6412

0.6412

0.6412

K=7

五、计算结果的分析

方程,的近似根为0.6412

方程图如下:当k增大时，x 的值趋近于一个常数，即方程的根。

计算中出现的问题，解决方法及体会

二分法尽管其收敛速度不太快，但二分法简单直观，特别适合于为收敛的迭代公式提供好的初值。牛顿迭代法具有较快的收敛速度，但对初值选取要求较高。弦截法比牛顿迭代法收敛速度稍慢，但它的计算量比牛顿迭代法小，特别当都函数的导数的计算比较复杂时，弦截法更显示了它的优越性。

教 师 评 语

指导教师：

2013 年 3 月 20 日

1

发表评论

请自觉遵守互联网相关的政策法规，严禁发布色情、暴力、反动的言论。

用户名:

验证码:

匿名?

发表评论