Scribd
Cargar
207 vistas7 páginas

PDF Metodos Numericos en Matlab 1 - Compress

Este documento presenta un manual sobre métodos numéricos en Matlab. Explica brevemente cómo utilizar el manual y luego procede a describir métodos para resolver ecuaciones algebraicas y trascendentes como la bisección, regla falsa, punto fijo, Newton-Raphson y secante. También cubre métodos de integración numérica como la regla del trapecio, Simpson 1/3, Simpson 3/8 y sus versiones compuestas, proporcionando código de Matlab para cada uno.

Cargado por

Abad Alejandro
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
207 vistas7 páginas

PDF Metodos Numericos en Matlab 1 - Compress

Este documento presenta un manual sobre métodos numéricos en Matlab. Explica brevemente cómo utilizar el manual y luego procede a describir métodos para resolver ecuaciones algebraicas y trascendentes como la bisección, regla falsa, punto fijo, Newton-Raphson y secante. También cubre métodos de integración numérica como la regla del trapecio, Simpson 1/3, Simpson 3/8 y sus versiones compuestas, proporcionando código de Matlab para cada uno.

Cargado por

Abad Alejandro
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
Está en la página 1/ 7

Métodos numéricos en Matlab 1

Mike Flores Theory

Sigue mis redes sociales

Mike Flores Theory Software para Ciencias e Ingeniería

Ingeniería |Tecnología |Educación (Tutoriales de Matlab y más)

https://www.facebook.com/Mike- https://www.facebook.com/Softwa
Flores-Theory-170750223735529/ reparacienciaseingenierias/

Mike Flores Theory

Ingeniería |Tecnología |Educación

https://www.youtube.com/channel
/UCgCpzENtk3DR1KCKjccDg2g?vie
w_as=subscriber

Preg úntate s i lo que es tás haciendo hoy te acerca al lug ar en el que quieres
es tar mañana. Walt Disney

Esperando que este manual sea de ayuda a todo aquel lector que se adentra
en el estudio de los métodos numéricos y el cómputo científico.

Página 1
Métodos numéricos en Matlab 1
Mike Flores Theory

Contenido

¿Cómo utilizar este manual? ........................................................................................................... 3

Solución de ecuaciones algebraicas y trascendentes ..................................................................... 4

Métodos cerrados (acotados) ......................................................................................................... 4

Método de la Bisección (Bipartición del Intervalo) ..................................................................... 4

Método de la Regla Falsa o de la Falsa Posición ......................................................................... 5

Métodos abiertos ............................................................................................................................ 6

Método del Punto Fijo................................................................................................................. 6

Método de Newton-Raphson ...................................................................................................... 7

Método de Newton-Raphson Modificado para el cálculo de raíces múltiples ........................... 8

Método de la Secante ................................................................................................................. 9

Integración numérica .................................................................................................................... 10

Regla del Trapecio ..................................................................................................................... 10

Regla del Trapecio Compuesta .................................................................................................. 11

Regla de Simpson 1/3 ................................................................................................................ 12

Regla de Simpson 1/3 Compuesta ............................................................................................ 13

Regla de Simpson 3/8 ................................................................................................................ 14

Regla de Simpson 3/8 Compuesta ............................................................................................ 15

Página 2
Métodos numéricos en Matlab 1
Mike Flores Theory

¿Cómo utilizar este manual?

Este manual es un compendio de códigos de métodos numéricos en Matlab, por


lo tanto no desarrolla la teoría correspondiente a cada método.

Se recomienda abrir un nuevo script en Matlab, allí copiar y pegar el código con el
cual se desea trabajar. Posteriormente correr el script.

Cada código contiene comentarios que indican el procedimiento a realizar, los


cuales inician con el símbolo (%) y están indicados en color verde. Como se
muestra en el siguiente ejemplo:

%Limpiamos la pantalla y mostramos el nombre del método


clear
clc
disp('Método de la bisección')

Estos comentarios pueden ser borrados por el usuario y no alterarán la ejecución


correcta del Script

Figura 1.-Script en Matlab

Página 3
Métodos numéricos en Matlab 1
Mike Flores Theory

Solución de ecuaciones algebraicas y trascendentes

Métodos cerrados (acotados)

Método de la Bisección (Bipartición del Intervalo)

%Limpiamos la pantalla y mostramos el nombre del método


clear
clc
disp('Método de la bisección')
%Damos de alta la variable simbólica X
syms x
%Introducimos la función,los límites superior e inferior,así como
%porcentaje de error
f =input('Introduzca la función f(x):');
xa =input('Introduzca el valor de Xa :');
xb =input('Introduza el valor de Xb :');
err =input('Porcentaje de error :');
%Graficamos la función
ezplot(f), grid on
f =inline(f);
erro=100;
xr =0;
i =0;
%Comprobamos que la función cambie de signo en el intervalo
%Si f(xa).f(xr)<0 resuélvase xb=xr
if f(xa)*f(xb)<0
while erro>err
i=i+1;
ea=xr;
%Aproximamos la raiz con la fórmula correspondiente
xr=(xa + xb)/2;
%Si f(xa).f(xr) >0 resuélvase xa=xr
if f(xa)*f(xr)>0
xa=xr;
else
xb=xr;
end
%Calculamos el porcentaje de error
erro=abs(((ea-xr)/xr)*100);
end
else
fprintf('No existe raíz en el intervalo dado')
end
%Mostramos los resultados en pantalla (con 3 decimales)
fprintf('\nRaíz=%10.3f en %4d iteraciones\n',xr,i);

Página 4
Métodos numéricos en Matlab 1
Mike Flores Theory

Regla de Simpson 3/8

%Limpiamos la pantalla y mostramos el nombre del método


clear
clc
disp('Regla de Simpson 3/8')
%Damos de alta la variable simbólica X
syms x
%Introducimos la función,los puntos a y b,así como
%porcentaje de error
f =input('Introduzca la función f(x):');
a =input('Introduzca el punto a: ');
b =input('Introduzca el punto b: ');
%Graficamos la función en el intervalo [a,b]
ezplot(f,[a,b]);
grid,title('\bf Regla de Simpson 3/8')
f=inline(f);
h=((b-a)/3);
x=a;
sum=0;
for i=2:3
x= x + h;
sum=sum + 3*f(x);
end
%Aproximamos la integral con la fórmula correspondiente
int=((b-a)/8)*(f(a)+sum +f(b));
%Mostramos los resultados en pantalla (con 3 decimales)
fprintf('\nAproximación a la integral=%12.3f \n',int);

Página 14
Métodos numéricos en Matlab 1
Mike Flores Theory

Regla de Simpson 3/8 Compuesta

%Limpiamos la pantalla y mostramos el nombre del método


clear
clc
disp('Regla de Simpson 3/8 compuesta')
%Damos de alta la variable simbólica X
syms x
%Introducimos la función,los puntos a y b, número de subintérvalos así
como
%porcentaje de error
f =input('Introduzca la función f(x):');
a =input('Introduzca el punto a: ');
b =input('Introduzca el punto b: ');
n =input('Número de subintervalos (n): ');
%Graficamos la función en el intervalo [a,b]
ezplot(f,[a,b]);
grid,title('\bf Regla de Simpson 3/8 compuesta')
f=inline(f);
%Calculamos la longitud del subintérvalo
h =(b-a)/(2*n);
f0=0;
for i=1:n-1
x=a+h*(2*i);
f0=f0+f(x);
end
f1=0;
for i=1:n
x=a+h*(2*i-1);
f1=f1+f(x);
end
f0=2*f0+4*f1;
f0=f0+f(a)+f(b);
%Aproximamos la integral
aprox=(h/3)*f0;
%Mostramos los resultados en pantalla (con 3 decimales)
fprintf('\nAproximación a la integral=%12.3f \n',aprox);

Página 15

También podría gustarte