数值分析实验五 非线性方程的求根4.docxVIP

PAGE

PAGE1

数值分析实验五 非线性方程求根

组号 班级 学号 姓名 分数

一、实验目的

1、掌握用MATLAB软件求解非线性方程和方程组的基本用法，并对结果做出初步的分析。

2、学会用二分法、牛顿法求非线性方程的根。

二、实验内容及基本知识介绍

1、二分法

假设f(x)∈C([a,b])，由连续函数的性质知，如果区间[a,b]内存在相异的两点x,x对

1 2

应的函数值异号。则方程f(x)?0在(x,x)内至少有一个根。在此我们假设f(a)f(b)?0称

1 2

(a,b)为有根区间。

yy

y

y? f(x)

a

1

b1

O

a

x?

x

0

b

为了找到[a,b]内的根α，令a?x,b?x

，则可以按下面的二分法计算

第一步：取中间中点c?

1 2

a?b

并计算f(c)

2

第二步：如果f(c)?0或c?a?0则认为a?c，算法终止；如果|f(c)|比较大，则f(c)或大于零或小于零，因而必与f(a)或f(b)异号，就是说要么区间[a,c]中存在方程的根，要么区间[c,b]中存在方程的根，如果[a,c]中存在方程的根，即f(a)f(c)?0,令b?c。返回

第一步。如果[c,b]中存在方程的根，即f(b)f(c)?0，令a?c，也返回第一步。

按二分法求根，只会出现两种情况：一种情况是在某一步f(c)?0成立；一种情况是不断更新的区间[a,b]的长度h逐次减半，直到趋于零，这时，因为??[a,b]，所以我们就可以认为??c，因此，二分法必有限步终止。

2、牛顿法 计算公式为

x ?x ?f(x)/f?(x) k=0,1,2,…

k?1 k k k

定理 设所要求解的方程为f(x)?0，且此方程有根x*。如果f(x*)?0，而且f?(x*)在x*

得某个邻域上连续，则牛顿法是局部收敛的并且至少是二阶收敛的。

设f(x)其单根x?附近有连续的二阶导数，将f(x)在x?x

处泰勒展开得

f(x)?f(x

k?1

)?f?(x?x )?

k?1

f?(?)2

(x?x

k?1

k?1

)2

从而 f(x?)?f(x

k?1

)?f?(x

k?1

)(x??x )?

k?1

f?(?)2

(x??x

)2

k?1

当x 充分靠近x?时，得

k?1

这样f(x?)?0导出

f(x?)?f(x

k?1

)?f?(x

k?1

)(x??x )

k?1

f(x )

x??x ???k?1

k?1

f?(x )

k?1

由此结果我们预期x

k?1

f(x )

? ? k

? ?

f(x )

因此x

k?1

精度更高的近似值。这样牛顿迭代法

k?1

x ?x

k k?1

f(x )

? ? k

? ?

f(x )

（*）

k?1

yL

y

L

x?

x

x

x

2 1 0

从几何上理解，牛顿迭代法每次迭代是用切线代替曲线，所以又称切线法。

f(x)

牛顿迭代法的映射为g(x)?x? 显见x?g(x)等价于f(x)?0，故迭代式（*）定

f?(x)

义合理，由于f(x?)?0而x?是单根时，f?(x?)?0

?[f?(x?)]2?f(x)f?(x) f(x)f?(x) f(x?)f?(x?)

?

g(x)?1? ? ?g(x?)? ?0

[f?(x)]2

[f?(x)]2

[f?(x?)]2

牛顿迭代法拒不收敛，且具有二阶收敛速度，计算结果条件一般用|x

k

x

k?1

|??，这是由于

g?(x?)?0，当x

k

与x?充分靠近时，L很小，一般能保证|x??x

k

|??。

三、实验问题及方法、步骤

1、二分法

例题：求方程f(x)?x3?2x2?4x?7?0在区间[3,4]内的一个实根，要求精确到小数点后的第二位。编程如下：

functionx=nabisect(fname,a,b,e)

%用途:二分法解非线性方程f(x)?0

%格式:x=nabisect(fname,a,b,