MODELOS MATEMÁTICO - SOLUCIONES
(Raíces de Ecuaciones).
Profesor: Jorge Ayala Mina
Modelos Matemáticos generalmente son no lineales
Velocidad de una burbuja
=
de galena.
=
Distribución partículas
Gaudin - Meloy
Modelo Gaussiano de la columna
de humo
Ecuación de transferencia de
calor:
SOLUCIONES MODELOS MATEMÁTICOS
(No lineales)
Resolver para
determinada Variable
Modelo
Matemático
Optimizar (Encontrar
valores máximos y
mínimos)
SOLUCIONES MODELOS MATEMÁTICOS
(No lineales)
% Valores iniciales
clc
clear
format long
x=0; % Valor de valor inicial.
error=eps;
k=1;
% cuerpo del programa.
CÓDIGO fx=exp(-x)-x;
MÉTODO DE dfx = - exp(-x)-1;
NEWTON xi = x -(fx/dfx);
RAPHSON errorcalculado= abs(xi - x);
while errorcalculado > error
fxi=exp(-xi)-xi;
dfxi=(-exp(-xi))-1;
xi1= xi - (fxi/dfxi);
errorcalculado = abs(xi1-xi);
xi=xi1;
k=k+1;
end
k,xi
SOLUCIONES MODELOS MATEMÁTICOS
(No lineales)
SOLUCIONES MODELOS MATEMÁTICOS
(No lineales)
Ecuaciones que no se pueden calcular por métodos analíticos como los
polinomios
O determinar la frecuencia de las oscilaciones transversales de una viga con
extremos empotrados y sometidos a un golpe
“ Cuando se trabaja con modelos matemáticos, generalmente al reemplazar
algunas constantes se convierte en un polinomio”
SOLUCIONES MODELOS MATEMÁTICOS NO LINEALES
UNIDIMENSIONALES
MÉTODO NUMÉRICO
- Método de bisección
- Método regla falsa
- Newton Raphson
- Método Secante
Modelo
Matemático MATLAB
FZERO
ROOTS
Raíces de ecuaciones no lineales.
(MÉTODOS NUMÉRICOS)
Métodos Cerrados
1. Método de Bisección
2. Método de la Regla Falsa
Métodos Abiertos
3. Método de Newton Raphson
4. Método de la Secante
1. Método Bisección
Aproximación a la raiz:
% Datos
clear
clc
format long
% Seleccionando el intervalo
xl=0;
xu=1;
error = 0.0001
error1= 1;
%Grafica de la formula o del modelo
x=xl:.1:xu;
% Modelo Matematico
Código de fx= exp (-x)-x;
plot (x,fx), grid
% la solución esta dentro del rango seleccionado?
método de
fxl=exp(-xl)-xl;
fxu=exp(-xu)-xu;
if fxl * fxu > 0
disp ('intervalo incorrecto, elija otros valores' )
bisección. end
xr= (xl + xu)/2;
k=1;
while error1 > error
% inicia las iteraciones
fxl= exp(-xl)-xl;
fxr= exp(-xr)-xr;
if fxl*fxr < 0
xu=xr;
end
if fxl*fxr > 0
xl= xr;
end
xr1 =(xl + xu)/2;
error1 = abs(xr1-xr);
xr = xr1;
k = k + 1;
end
k,xr
2. Método Newton Raphson
Determinar Pendiente:
Pasar al Jamboard
y Matlab
% Valores iniciales
clc
clear
format long
x=0; % Valor de valor inicial.
error=eps;
k=1;
% cuerpo del programa.
CÓDIGO fx=exp(-x)-x;
MÉTODO DE dfx = - exp(-x)-1;
NEWTON xi = x -(fx/dfx);
RAPHSON errorcalculado= abs(xi - x);
while errorcalculado > error
fxi=exp(-xi)-xi;
dfxi=(-exp(-xi))-1;
xi1= xi - (fxi/dfxi);
errorcalculado = abs(xi1-xi);
xi=xi1;
k=k+1;
end
k,xi
EJEMPLO DE TRANSFERENCIA DE CALOR POR CONDUCCIÓN Y RADIACIÓN.
T1(exterior) P T(interior)
= ?°K a = 625°K
r
e
d
Polinomio
Raíces de Ecuaciones (Matlab) - fzero
La función Valor inicial
El algoritmo, creado por T. Dekker, utiliza una combinación de bisección,
secante, y métodos de interpolación cuadrática inversa.
Ejemplo transferencia calor.
clear
clc
p = @(x) 4.53E-08*x.^4 +44*x-2.0602E+04;
x = 300:0.1:600;
plot(x,p(x)), grid
hold on
Z = fzero(p, 400)
plot(Z,p(Z),'r*')
Raíces de Polinomios (Matlab) - root.
r = roots (p)
Coeficientes
polinomio
Método de Muller ( Utiliza método Secante)
Método Bairstow ( Utiliza combinación método de Muller Newton Raphson)
CÓDIGO MATLAB
clc
clear
x = linspace(-6.1,5.1,100);
p = @(x) x.^6+3*x.^5-41*x.^4-
87*x.^3+400*x.^2+444*x-720;
plot (x,p(x)), grid
roots([1 3 -41 -87 400 444 -720])
Trabajo
1. Programar en Matlab el problema de transferencia de calor con el método de
Newton Raphson y comparar usando el comando Fzero del Matab.
2. La concentración de bacterias contaminantes c en un lago disminuye de
acuerdo a la ecuación
Determine el tiempo que se requiere para que la concentración de las bacterias
se reduzca a 15. Utilice la programación del Matlab y compare con los resultados
con el comando Fzero.
Libro de Chapra
Raices de
Ecuaciones
Utilizando
FZERO y
ROOTS
