ESCUELA PROFESIONAL DE INGENIERIA MECÁNICA,
MECÁNICA ELECTRICA Y MECATRÓNICA
MECÀNICA COMPUTACIONAL II
PRACTICA N°: 05
GRUPO N° : 13
TEMA : Métodos Numéricos – Integración Numérica
PROFESOR : Juan Diego Cáceres.
ALUMNO : Ccopa Villalba Fredy
INDICE
ANALISIS OBJETIVO DEL PROGRAMA ................................................................................. 3
CUESTIONARIO FINAL ............................................................................................................ 6
ANÀLISIS FINAL ..........................................................................................................................8
OBSERVACIONES........................................................................................................................8
CONCLUSIONES ........................................................................................................................ 9
BIBLIOGRAFIA............................................................................................................................ 9
ANALISIS OBJETIVO DEL PROGRAMA
1. Localice gráficamente las raíces de f(x) = 0, siendo f(x) = |x| - ex.
%Metodo de Biseccion - Metodos Numericos Matlab
clc;
%Introducimos los datos propios del ejercicio
Fx=input('Ingrese la funcion: ','s');
a=input('Ingrese Limite Inferior a : ');
c=input('Ingrese Limite Supeior c : ');
e=input('Ingrese el error : ');
%Evaluamos la funcion Fx y lo ordenamos en una matriz de n x n
x=a;
Fa=eval(Fx);%Evalua la funcion Fx para su Limite Inferior
x=c;
Fc=eval(Fx);%Evalua la funcion Fx para su Limite Superior
fprintf('\n %6s %7s %8s %10s %8s %8s %8s \n
','A','B','C','F(a)','F(b)','F(c)','|c-a|');
while abs(c-a)>e
b=(a+c)/2;%Encontrar el punto medio entre los limites
x=b;
Fb=eval(Fx);
fprintf('\n %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f
\n',a,b,c,Fa,Fb,Fc,abs(c-a));
if Fa*Fb<=0
c=b;
Fc=Fb;
else
a=b;
Fa=Fb;
end
end
fprintf('\nEl resultado sera %.4f\n',b);
ezplot(Fx);%graficamos la funcion
grid on;
Imagen 01. Comprobación del método de la Bisección.
Este programa nos a encontrar anormalidades en las curvas de las funciones, tanto
si es cóncavo o convexo, ya que de acuerdo a los datos que nos muestra podemos
predecir su comportamiento.
2. Confeccionar en MATLAB un programa que:
Calcule la raíz de una ecuación algebraica no lineal o ecuación trascendente
disp('programa para resolver Ecuaciones Cuadráticas');
disp('=============================================');
% Módulo de introducción de los coeficientes a, b, c, desde teclado.
disp('Introduzca coeficientes de la ecuación a.x^2 + b.x + c = 0');
A = input('Valor del coeficiente a: ');
B = input('Valor del coeficiente b: ');
C = input('Valor del coeficiente c: ');
% Caso de ecuación (a = 0 y b = 0)
RootsFound = 0;
if A == 0 & B == 0,
fprintf('No es posible resolver la ecuación . Ecuación extremadamente
degenerada');
fprintf('%14.8g = 0.0\n', C);
RootsFound = 1;
end;
% Caso de ecuación (a = 0)
if A == 0 & B ~= 0,
Raiz1 = -C/B;
fprintf('Ecuación degenerada: Raiz = %5g\n', Raiz1);
RootsFound = 1;
end;
% Caso de ecuación (a distinto de 0)
if RootsFound == 0,
Discriminante = B^2 - 4*A*C;
if Discriminante >= 0
Raiz1 = (-B - sqrt(Discriminante))/(2*A);
Raiz2 = (-B + sqrt(Discriminante))/(2*A);
fprintf('Dos raíces reales: Raiz1 = %5g\n', Raiz1);
fprintf(' Raiz2 = %5g\n', Raiz2);
else % Caso contrario, raíces complejas
fprintf('Dos raíces complejas: Raiz1 = %5g\n + %5g i\n',-B/(2*A), sqrt(-
Discriminante)/(2*A));
fprintf(' Raiz2 = %5g\n + %5g i\n',-B/(2*A), -sqrt(-Discriminante)/(2*A));
end;
end
Imagen 02. Comprobación de las raíces de la función.
Este programa nos ayuda bastante para hallar las raíces de una función algebraica
que no necesariamente tenga raíces reales, sino que también pueden ser imaginarias
como es el este caso.
CUESTIONARIO FINAL
1. Use cuatro iteraciones del método del la Bisección para determinar las raíces de e2x
– 6x = 0, en el intervalo [0, 0.5].
¿Cuántas iteraciones son necesarias para obtener la aproximación a la raiz
redondeada a 5 cifras decimales?
clear all
format short;
a=input('Introduzca el valor de a: ');
b=input('Introduzca el valor de b: ');
cont=input('Introduzca el número de iteraciones cont: ');
fun=input('Introduzcal a funcion f(x)=','s');
f=inline(fun);
for k=1:cont
c=(a+b)/2;
e=abs((b-a)/2);
A(k,:)=[k a b c f(c) e];
if f(a)*f(c)<0
b=c;
else
a=c;
end
end
fprintf('\n \tk \ta \tb \tc \tf(c) \terror \n');
disp(A);
fprintf('Solución:\n c=%8.5f\n',c);
fprintf('f(c)=%8.5f\n',f(c));
fprintf('error=%8.5f\n',e);
Imagen 03. Comprobación del método de la Bisección.
Este programa nos ayuda a encontrar el número de iteraciones máximos hasta que
el error sea de 0.0001, dando una muy buena aproximación al grafico original, así
como ver el comportamiento de la función en todos los puntos (iteraciones).
ANÀLISIS FINAL
El desarrollo de este método de la Bisección es muy útil, ya que es un método fácil a
comparación de otros métodos que son un poco más precisos, pero por ende son más
complicados al momento de plasmarlo en un programa.
Además, este programa os ayuda a comprender, los fenómenos que ocurren en la
función dada, en todas sus iteraciones. Pudiendo hallar puntos de inflexión, como
partes cóncavas o convexas.
OBSERVACIONES
Se observo que para toda función hay un limite máximo de numero de
iteraciones las cuales dependen del error que nosotros decidamos ponerle en
la función.
Se observo la utilidad del programa para hallar el comportamiento de la
función, pudiendo ver todos los puntos mediante una matriz y analizar su
comportamiento.
CONCLUSIONES
Se aplico los métodos iterativos de intervalo, específicamente el método de la
bisección en la solución de ecuaciones no lineales algebraicas y trascendentes.
Se implemento programas en MATLAB que solucionen las raíces de
ecuaciones algebraicas y trascendentes por medio del método de la bisección.
BIBLIOGRAFIA
MANUAL BÀSICO DE MATLAB, M. Cristina Casado Fernandez. Primera
Edicion. Rescatada el 28/08/2018 de
[Link]
MANUAL DE INTRODUCCIÒN DE MATLAB, Manuel Lopez Martinez y Jose
Angel Acosta. 2004. Rescatada el 28/08/2018 de
[Link]
CREACION DE INTEFACES GRAFICAS, Gonzalo Fernandez Cordova.(2007).
Rescatada el 10/09/2018 de
[Link]
