gexianfashoulianjie.m

牛顿迭代法

function [x,k,t,z]=gexianfashoulianjie(f,x1,x2,dieci,xz)
% 求解非线性方程的割线法，其中
% f(x) --- 表示非线性方程函数
% x1,x2 --- 表示迭代初值
% dieci --- 表示指定迭代次数
% k --- 表示循环次数
% xz --- 表示不动点值,即x的真值
% t --- 表示lne(k)
% z --- 表示lne(k+1)
k=1;
y1=feval(f,x1);
y2=feval(f,x2);
x(k)=x2-(x2-x1)*y2/(y2-y1);
y(k)=feval(f,x(k));
k=k+1;
x(k)=x(k-1)-(x(k-1)-x2)*y(k-1)/(y(k-1)-y2);
while k<dieci
    y(k)=feval(f,x(k));
    x(k+1)=x(k)-(x(k-1)-x(k))*y(k)/(y(k-1)-y(k));
    e(k)=x(k)-xz;
    e(k-1)=x(k-1)-xz;
    t(k-1)=log(e(k-1));
    z(k-1)=log(e(k));
    k=k+1;
end
plot(t,z,'r*:');grid on;
xlabel('t轴，代表lne(k)');ylabel('z轴，代表lne(k+1)');
title('迭代之后(3)式图形');