Matlab - Graficar las iteraciones de Newton Raphson

 
Vista:
sin imagen de perfil

Graficar las iteraciones de Newton Raphson

Publicado por Toni (2 intervenciones) el 09/12/2014 23:25:12
¡Hola! Me gustaría poder plotear (graficar) las iteraciones que se van dando, es decir los resultados parciales en cada iteracion en la grafica.

El código es el de Newton Raphson:
%Algoritmo de Newton Raphson
%fun y der son funciones que se usaran como handle. fun está en el archivo
%f.m y der en df.m
function [x,it]=newton(fun,der,x0,tol)
fx=fun(x0);
%No hay restricción de iteraciones:
maxiter= 99999;
%Un contador para las iteraciones para la variable de salida "it"
it=0;
for iter=1:maxiter
fx_der=der(x0);
x=x0-fx./fx_der;
fx=fun(x);
if abs(x-x0)<=abs(x0)*tol
fprintf('\n Newton Raphson CONVERGE en %d iteraciones a %f',it,x)
return
else
x0=x;
end
it=it+1;
fprintf('\n %i %f %f',it,x,fx);
end
fprintf('\n Newton Raphson NO CONVERGE en %d iteraciones',maxiter)

end
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder