Método de Newton en matlab PDF

Preview

Summary

Código para resolver ejercicios con el método de Newton en matlab...

Description

%% Newton-Raphson %{ Este codigo resuelve ecuaciones no lineales mediante el metodo de Newton Raphson. Es un metodo abierto, por lo tanto no necesitamos encerrar a la raiz, pero si necesitamos una aproximacion inicial . %} %% Datos %{ Programadores : Ghilino Ramiro Eujanian Amir Villarreal Juan Rojas Juan %} %% Variables %{ p= Vector en donde guardaremos los puntos que vayamos obteniendo f= Lugar de asignación de la funcion df= derivada de la funcion M= variable para uso del menu interno N= numero de iteraciones %} %% Desarrollo del Codigo function [pfinal]=newtonR() clear clc disp('Newton Raphson'); syms x %creamos la variable x como symbolic %% Ingreso de Datos p(1)=input('Ingrese Punto de Partida: ');%Punto inicial prompt='Introduzca la funcion a continuacion >> ';%solicitamos funcion f=input(prompt,'s'); f=str2sym(f); df=diff(f);%calculo la derivada de f(x) disp('Seleccione las opciones') %Menu de seleccion disp('1: Iterar una N cantidad de veces') disp('2: Iterar hasta un margen de error') M=0; %Inicializo M a 1 para obligarlo a entrar while(M~=1 && M~=2) M=input('Ingrese la opcion deseada:');%Este while evita que ingresen una opcion incorrecta end if M==1 N=input('Ingrese el numero de iteraciones ');%Numero de iteraciones for i=1:N %calculo por numero de n iteracion fx=subs(f,x,p(i)); dfx=subs(df,x,p(i)); p(i+1)=p(i)-fx/dfx; end elseif M==2 error=input('Ingrese el error absoluto deseado: '); %calculo por margen de error i=1; fx=subs(f,x,p(i)); dfx=subs(df,x,p(i));

p(i+1)=p(i)-fx/dfx; while(abs(p(i)-p(i+1))> error) %Consultamos para ver si en la primera iteracion ya cumplimos con el error %en caso contrario iteramos hasta lograrlo i=i+1; %Comenzamos a repetir lo de antes fx=subs(f,x,p(i)); dfx=subs(df,x,p(i)); p(i+1)=p(i)-fx/dfx; end end %%Resultados pfinal = p(end); %graficacion de los valores obtenidos en distintas iteraciones disp('El valor de la raiz es: '); plot(p) title('CONVERGENCIA DEL METODO') ylabel('valor de p') xlabel('N? de iteraciones') end...