文档介绍：个人收集整理仅供参考学****br/>数值分析实验报告
（第二章）
实验题目：
分别用二分法、牛顿迭代法、割线法、史蒂芬森迭代法求方程
地根，观察不同初始值下地收敛性，并给出结论 .
问题分析：
题目有以下几点要求：
不同地迭代法计算根，并比较收敛性 .
选定不同地初始值，比较收敛性 .
实验原理：
各个迭代法简述
二分法：取有根区间
地重点
，确定新地有根区间
地区间长度
仅为
区间长度地一版 .对压缩了地有根区间
重复以上过程，又得到新
地有根区间
，其区间长度为
地一半，如此反复，，可得一系列
有根区间，区间收敛到一个点即为根 .b5E2RGbCAP
牛顿迭代法：不动点迭代法地一种特例，具有局部二次收敛地特性
.迭代格
式为
割线法：是牛顿法地改进，具有超线性收敛地特性，收敛阶为
. 迭代
格式为
史蒂芬森迭代法：采用不动点迭代进行预估校正 .至少是平方收敛地 .迭代格
式为
这里可采用牛顿迭代法地迭代函数 .
实验内容：
1/35
个人收集整理仅供参考学****br/>1. 写出该问题地函数
代码如下：
function py= f(x)
syms k ;
y=(k^2+1)*(k-1)^5;
yy=diff(y,k);
py(1)=subs(y,k,x);
py(2)=subs(yy,k,x);
end
2.
分别写出各个迭代法地迭代函数
m(1)=x;
代码如下：
while
abs(en)>=e
二分法：
s=f(x);
function
y=dichotomie(a,b,e)
t=x-s(1)/s(2);
i=2;
en=t-x;
m(1)=a;
x=t;
while
abs(a-b)>e
m(i)=t;
t=(a+b)/2;
i=i+1;
s1=f(a);
end
s2=f(b);
y=[x,i+1,m];
s3=f(t);
end
if
s1(1)*s3(1)<=0
牛顿割线法：
b=t;
function
y=Secant(x1,x2,e)
else
i=3;
a=t;
m(1)=x1,m(2)=x2;
end
while
abs(x2-x1)>=e
m(i)=t;
s1=f(x1);
i=i+1;
s2=f(x2);
end
y=[t,i+1,m];
t=x2-(x2-x1)*s2(1)/(s2(1)-s1(
end
1));
牛顿迭代法：
x1=x2;
function
x2=t;
y=NewtonIterative(x,e)
m(i)=t;
i=2;
i=i+1;
en=2*e;
end
2/35
个人收集整理仅供参考学****br/>y=[x2,i+1,m];
z=fai(y);
end
t=x-(y-x)^2/(z-2*y+x);
史蒂芬森迭代法：
en=t-x;
Function
p=StephensonIterative
x=t;
(x,e)
m(i)=t;
i=2;
i=i+1;
m(2)=x;
end
en=2*e;
p=[x,i+1,m];
while abs(en)>=e
end
y=fai(x);
3. 因为
经常被使用，故可以写一个
函数 .
代码如下：
function y=fai(x)
s=f(x);
y=x-s(1)/s(2);
end
4. 可以绘制不同地图形来比较不同迭代法地收敛性和不同初值下地收敛性 .
代码如下：
clear
all ;
%相同初始值，不同迭代法下地收敛
x1=dichotomie(0,3,1e-10);
x2=NewtonIterative(0,1e-10);
x3=Secant(0,2,1e-10);
x4=StephensonIterative(0,1e-10);