牛顿法与二分法的比较—matlab实现

刚学完牛顿迭代法，为了验证收敛的速率，用Matlab做了比较

首先是牛顿迭代法

%比较牛顿迭代法、
function [x,i]=newtonmethod(x0,f,ep,Nmax)%x0—初值，f—测试函数，ep—精度，Nmax—迭代的最大次数
i=1;
x(1)=x0;
while(i<=Nmax)
[g1,g2]=f(x(i));
if abs(g2)<=ep
%error('error');
disp('derivative is too smal')
return
end
x(i+1)=x(i)-g1/g2;
b=x(i+1)-x(i);
if(abs(b)<ep)
return;
end
i=i+1;
end
x

i

二分法

%二分法
function [x,k]=bisection(a,b,f,ep)%[a,b]区间，f—函数句柄，ep—最大二分次数
k=1;
eep=b-a;%eep区间长度
while(abs(eep)>ep)
x(k)=a+eep./2;
if(f(x(k)).*f(a)<0)
b=x(k);
else
a=x(k);
end
eep=b-a;
k=k+1;
end

函数f

function [fx,f1x]=f221(x)
% 课本p.42, 实验2.2, I中的函数及其导数, 含根区间为[1,2]
fx=x*x*x+4*x*x-10;
f1x=3*x*x+8*x;
ps:自己直接写的求导后的函数，为了方便

主脚本文件，即运行

%牛顿法选初始值1.5，二分法的含根区间为[1,2],
[x,i]=newtonmethod(1.5,@f221,0.00001,100);%@f221——获取函数句柄
[xx,k]=bisection(1,2,@f221,0.00001);
plot(x,'-rd');
hold on
plot(xx,'-gs');
xlabel('迭代序列');
ylabel('迭代结果');
title('shiyan');
legend('牛顿迭代法','二分法' );
hold off