【源码】牛顿法求解实值函数的根

示例：[ x, ex ] = newton( ‘exp(x)+x’, ‘exp(x)+1’, 0, 0.5*10^-5, 10 )

f：输入函数

df：输入函数的导数

x0：函数根的初值估计

tol：误差容忍度

namx：求解最大迭代次数

x：求解输出的近似根

ex：误差估计

MATLAB完整源码：

function [ x, ex ] = newton( f, df, x0, tol, nmax )

%

% NEWTON Newton’s Method

% Newton’s method for finding successively better approximations to the

% zeroes of a real-valued function.

%

% Input:

% f - input funtion

% df - derived input function

% x0 - inicial aproximation

% tol - tolerance

% nmax - maximum number of iterations

%

% Output:

% x - aproximation to root

% ex - error estimate

%

% Example:

% [ x, ex ] = newton( ‘exp(x)+x’, ‘exp(x)+1’, 0, 0.5*10^-5, 10 )

%

% Author: Tashi Ravach

% Version: 1.0

% Date: 16/04/2007

%

if nargin == 3

    tol = 1e-4;

    nmax = 1e1;

elseif nargin == 4

    nmax = 1e1;

elseif nargin ~= 5

    error('newton: invalid input parameters');

end



f = inline(f);

df = inline(df);

x(1) = x0 - (f(x0)/df(x0));

ex(1) = abs(x(1)-x0);

k = 2;

while (ex(k-1) >= tol) && (k <= nmax)

    x(k) = x(k-1) - (f(x(k-1))/df(x(k-1)));

    ex(k) = abs(x(k)-x(k-1));

    k = k+1;

end

end

关于牛顿法的MATLAB源码及PPT下载地址：

http://page5.dfpan.com/fs/blc4j2621129a165927/

更多精彩文章请关注微信号：