Guia matlab

GUIA MATLAB
SOLUCION DE ECUACIONES NO LINEALES Y SISTEMAS LINEALES
En este taller usaremos el programa MATLAB con el fin de resolver ecuaciones no lineales y sistemas de ecuaciones lineales, de manera rápida y fácil. Se usarán tanto las herramientas propias de MATLAB, como rutinas creadas por el usuario que nos llevarán paso a paso a la solución de problemas.

PRIMERA PARTE: SOLUCIÓN DEECUACIONES NO LINEALES DE LA FORMA F (X) = 0
Comencemos revizando unos fáciles ejemplos que nos mostrarán cuales son los pasos para resolver una ecuación no lineal empleando MATLAB, por medio de rutinas creadas por el usuario y que con anterioridad deben ser creadas por usted en la forma como se explicó en la inducción (y en la guía que aparece en la página del curso).

• Ejemplo: Determine valoresaproximados de las soluciones positivas de la ecuación ³x´ 1 exp − sin (x) = 0 2 3 Solución: Lo primero que debemos hacer es graficar la función f (x) =
1 2

exp

para identificar las raices de la ecuación, para ello empleamos la instruccion fplot que permite graficar funciones simbólicas. La sintaxis es: fplot( ‘funcion’,[xmin xmax]) = grafica la función para los valores de x en el intervalo[xmin,xmax]. fplot( ‘funcion’,[xmin xmax ymin ymax]) = grafica la función para los valores de x en el intervalo [xmin,xmax] y las imágenes comprendidas en el intervalo [ymin,ymax]. 1

¡x¢
3

−sin (x)

Ambas instrucciones permiten agregar color adicionando ’color’ antes de cerrar el paréntesis. Para distinguir las raíces adicionamos la instrucción gird on que activa una

cuadrícula a la gráficade la función. Ejecutamos las siguientes instrucciones fplot(‘0.5*exp(x/3)-sin(x)’,[-10 10 -1 1]), grid on) fplot(‘0.5*exp(x/3)-sin(x)’,[0 2 -0.3 0.3]), grid on) con la primera identificamos el intervalo donde se encuentran las raices positivas y con el segundo obtenemos la figura:

Para aproximar las raices podemos emplear el método de Bisección en los intervalos [0.5, 1] y [1.5, 2], clamamentela función es continua por ser la suma de funciones continuas (exponencial y trigonometrica). También podemos emplear el método de Newton (puesto que la función, la primera y la segunda derivada son continuas en todos los reales), pero debemos tomar el valor inicial muy cerca de la raíz para no caer en una zona de divergencia de la sucesión (es decir, que una aproximación sea el valor donde lafunción alcanza el mínimo). Para empezar, apliquemos el método de bisección con una tolerancia de 10−6 (es decir, 2

queremos que cuando la distancia entre dos aproximaciones sea menor a 10−6 paremos), para ello utilizamos la función "bisección" en MATLAB, que debe ser creada por el usuario de la misma manera como se indicó en la inducción (Una muestra sencilla de un programa típico se incluye alfinal del la guía). Digitando biseccion(’0.5*exp(x/3)-sin(x)’,0.5,1,0.000001) en la ventana Commad Window y oprimiento la tecla ENTER, obtenemos los siguientes resultados

3

Y así mismo para la otra aproximación

Para emplear el método de Newton tomaremos la misma tolerancia para la distancia entre las aproximaciones y emplearemos una función del método de Newton en MATLAB, que debe sercreada por el usuario (una muestra está al final la guía). ventana Command Window: newton(’0.5*exp(x/3)-sin(x)’,0,0.000001) Si digitamos en la

4

obtenemos los siguientes resultados para la primer raíz

y para la segunda raíz

• Ejemplo: El factor de fricción f para los fluidos turbulentos en una tubería está dado por 1 √ = 1.14 − 2 log10 f µ 9.35 e √ + D Re f ¶

llamada correlación deColebrook, donde Re es el número de Reynolds, e es la aspereza de la superficie de la tubería y D es el diámetro de la tubería. Resolver la ecuación para f utilizando el método de punto fijo para los siguientes casos: 1. (a) D = 0.1m, e = 0.0025m, Re = 3 × 104 (b) D = 0.1m, e = 0.0005m, Re = 5 × 106 Solución: Si queremos resolver el problema empleando el método de punto fijo debemos

5

llevar la...