Source: http://slideplayer.es/slide/1714769/
Timestamp: 2017-12-11 17:37:41+00:00

Document:
3. Métodos de resolución Ecuaciones algebraicas lineales - ppt descargar
3. Métodos de resolución Ecuaciones algebraicas lineales
Publicada porChicky Peralta Modificado hace 3 años
Presentación del tema: "3. Métodos de resolución Ecuaciones algebraicas lineales"— Transcripción de la presentación:
1 3. Métodos de resolución Ecuaciones algebraicas lineales
Ecuaciones algebraicas no lineales Métodos para una variable Métodos para multivariable Modelado y simulación en Ingeniería Química. Manuel Rodríguez
2 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
Ecuaciones Algebraicas Lineales No lineales Metodos Numericos Interval Halving (o bisection) False Position (o regula falsi) Ridder Succesive Substitution (o fixed-point) Secant Muller Wegstein Newton Raphson Metodos Analiticos Brent Broyden Homotopy Dogleg step Hook step Para problemas multidimensionales Modelado y simulación en Ingeniería Química. Manuel Rodríguez
3 Ecuaciones algebraicas lineales
Modelado y simulación en Ingeniería Química. Manuel Rodríguez
4 Ecuaciones algebraicas lineales
5 Ecuaciones algebraicas lineales
6 Ecuaciones algebraicas lineales
7 Ecuaciones algebraicas lineales
8 Ecuaciones algebraicas lineales
9 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
10 [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] Ejemplo 3x1+2x2=4 2x1+x2=1
Inversa de la matriz Ax=b x=A-1 b x1=-2 x2=5 Descomposición de Gauss [ ] [ ] 3 2 2 1 4 1 f1’=f1-2f2 [ ] [ ] -1 0 2 1 2 1 f1’=-f1 [ ] [ ] 1 0 2 1 -2 1 f2’=-2f1+f2 [ ] [ ] 1 0 0 1 -2 5 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
11 Ecuaciones algebraicas lineales
12 Ecuaciones algebraicas lineales
13 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
14 Ecuaciones algebraicas no lineales
Objetivo Sea f(x) una función no lineal en x. Hallar el valor de x, x*, tal que se cumple f(x*)=0. x* se suele denominar el cero o raíz de f(x) x* se puede determinar por medios analíticos (solución exacta) o por medios numéricos (solución aproximada) La elección del método numérico depende del problema a resolver (estructura del problema, tipo de ecuaciones, precisión requerida, rápidez del cálculo,....). Por tanto no existe un mejor método universalmente aplicable. Tipos de métodos Métodos acotados (bracketing methods) Métodos abiertos (open methods) Modelado y simulación en Ingeniería Química. Manuel Rodríguez
15 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
Métodos acotados Base: Una función cambia de signo en la proximidad de una raíz Una raíz está acotada en el intervalo [a,b] si el signo de f(a) es diferente al signo de f(b) Método de la bisección (o intervalo medio) Algoritmo Selecciona un intervalo [a,b] donde halla un cero Calcula el punto medio como nuevo punto Comprueba si hay cambio de signo en [a,p] o en [p,b]. Comprobación: f(a)*f(p). Si el producto es cero, entonces p es una raíz. Si no es cero volver al punto 2. [nuevo punto] [a,b] Modelado y simulación en Ingeniería Química. Manuel Rodríguez
16 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
Ejemplo Método bisección (Intervalo medio) Modelado y simulación en Ingeniería Química. Manuel Rodríguez
17 Método de la posición falsa
Ecuaciones algebraicas no lineales Método de la posición falsa [nuevo punto] [a,b] Algoritmo Selecciona un intervalo [a,b] donde halla un cero Calcula un punto intersección como nuevo punto Comprueba si hay cambio de signo en [a,p] o en [p,b]. Comprobación: f(a)*f(p). Si el producto es cero, entonces p es una raíz. Si no es cero volver al punto 2. Modelado y simulación en Ingeniería Química. Manuel Rodríguez
18 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
Ejemplo método de la posición falsa (Regula Falsi) Modelado y simulación en Ingeniería Química. Manuel Rodríguez
19 Comparación entre ambos métodos. Similaridades:
Ecuaciones algebraicas no lineales Comparación entre ambos métodos. Similaridades: Ambos métodos necesitan DOS valores iniciales Requieren un procedimiento para determinar el cambio de signo. Acaban convergiendo a la raíz con cierta tolerancia Diferencias: El cálculo del nuevo punto estimado se hace con diferentes estrategias En general el método de la posición falsa converge más rápido que el de la bisección. Modelado y simulación en Ingeniería Química. Manuel Rodríguez
20 Ecuaciones algebraicas no lineales
Métodos abiertos Emplean una aproximación funcional para obtener el nuevo valor estimado de la raíz (línea recta, cuadrática, polinomio) Métodos: Punto-fijo (sustitución sucesiva o directa) Newton-Raphson (línea recta empleando información del gradiente) Secante (línea recta empleando dos puntos) Muller (aprox. cuadrática empleando tres puntos) Modelado y simulación en Ingeniería Química. Manuel Rodríguez
21 Metodos acotados vs. Métodos abiertos
Ecuaciones algebraicas no lineales Metodos acotados vs. Métodos abiertos Métodos acotados La raíz está situada en un intervalo (necesita dos puntos). Acaba convergiendo dentro de una tolerancia. Métodos abiertos Sólo emplean un punto inicial (o dos puntos que no tienen por qué contener a la raíz) y una fórmula para encontrar la raíz. No siempre convergen, pero cuando lo hacen son mucho más rápidos que los métodos acotados. Modelado y simulación en Ingeniería Química. Manuel Rodríguez
22 Sustitución sucesiva Problema f(x)=0 Transformar a x=g(x)
Ecuaciones algebraicas no lineales Sustitución sucesiva Problema f(x)=0 Transformar a x=g(x) Seleccionar un punto inicial x0 Calcular nuevo valor xi+1=g(xi) Repetir hasta llegar a la tolerancia requerida Si: |g’(x)|<1 El algoritmo converge linealmente |g’(x)|>=1 El algoritmo diverge Modelado y simulación en Ingeniería Química. Manuel Rodríguez
23 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
Ejemplo Método punto fijo (Sustitución directa/sucesiva) ¿Converge? ¿Converge? Modelado y simulación en Ingeniería Química. Manuel Rodríguez
24 Newton Raphson g(xi) xi+1=xi- g’(xi) Problema g(x)=0
Ecuaciones algebraicas no lineales Problema g(x)=0 Seleccionar un punto inicial x0 Calcular g(xi) y g’(xi) Aplicar la tangente en ese punto y en el corte con el eje de abcisas tenemos el nuevo punto estimado Repetir hasta llegar a la tolerancia requerida Newton Raphson xi+1=xi- g(xi) g’(xi) Necesita conocer la derivada de la función Convergencia cuadrática (rápida) Puede no converger (depende de la función y de la estimación inicial) Modelado y simulación en Ingeniería Química. Manuel Rodríguez
25 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
Ejemplo Newton Method Modelado y simulación en Ingeniería Química. Manuel Rodríguez
26 Secante xi+1=xi- xi+1-xi g (xi+1)-g (xi) g (xi+1) Problema g(x)=0
Ecuaciones algebraicas no lineales Secante Problema g(x)=0 Seleccionar dos puntos iniciales x0,x1 Calcular la recta que pasa por esos puntos El corte con el eje de abcisas da el nuevo punto estimado. Volver a calcular la recta. Repetir hasta llegar a la tolerancia requerida xi+1=xi- xi+1-xi g (xi+1)-g (xi) g (xi+1) No Necesita conocer la derivada de la función (la aproxima). Necesita dos puntos iniciales. Puede no converger. Modelado y simulación en Ingeniería Química. Manuel Rodríguez
27 Ecuaciones algebraicas no lineales
La primera iteración da el mismo resultado, luego cada uno obtiene un nuevo punto estimado diferente Modelado y simulación en Ingeniería Química. Manuel Rodríguez
28 Ecuaciones algebraicas no lineales
* * Se computan obligando a que g(x) pase por los 3 puntos seleccionados. Modelado y simulación en Ingeniería Química. Manuel Rodríguez
29 Ecuaciones algebraicas no lineales
30 Sistemas de ecuaciones algebraicas no lineales
31 Sustitución sucesiva acelerada
Ecuaciones algebraicas no lineales Sustitución sucesiva acelerada La sobrerelajación es una técnica para acelerar la convergencia de métodos iterativos en la solución de ecuaciones lineales. La idea es aplicarlo al método de sustitución sucesiva. Se dan pesos a los valores anteriores y a los previos con el fin de dar pasos “mayores” hacia la solución. Xk+1=qXk+(1-q)g(Xk) q=0 sustitución sucesiva q<0 aceleración de la convergencia 0<q<1 estabilización de la convergencia por amortiguamiento Modelado y simulación en Ingeniería Química. Manuel Rodríguez
32 Ecuaciones algebraicas no lineales
33 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
En los problemas de diagrama de flujo las variables de iteración no están todas muy acopladas ni todas desacopladas, nos podemos encontrar: Todas las especies débilmente acopladas a través de equilibrio L-V (no ideal) Especies muy acopladas si participan en una reacción Si hay más de una corriente de rasgado, hay un fuerte acoplamiento entre los flujos de cada corriente de rasgado. Cuando las variables están desacopladas o débilmente acopladas es necesario un parámetro de aceleración diferente para cada variable. Modelado y simulación en Ingeniería Química. Manuel Rodríguez
34 Ecuaciones algebraicas no lineales
El método de Wegstein es muy bueno para particiones en las que hay una única corriente de rasgado. O cuando hay reciclo sin estar los componentes muy acoplados (presencia de reacción) Modelado y simulación en Ingeniería Química. Manuel Rodríguez
35 Ecuaciones algebraicas no lineales
Newton Raphson En lugar de la derivada emplea el jacobiano (matriz de derivadas parciales) La estimación del nuevo conjunto de raíces se computa mediante la siguiente ecuación: ¿Cómo resolverías la ecuación anterior sin tener que calcular la inversa de la matriz jacobiana? Jacobiano Modelado y simulación en Ingeniería Química. Manuel Rodríguez
36 Newton Raphson- procedimiento de resolución
Ecuaciones algebraicas no lineales Newton Raphson- procedimiento de resolución ¡Resuelve un sistema de ecuaciones lineales! Actualiza el valor hasta que ∆X es tan pequeño como se haya requerido Ventajas: a) Buena convergencia (cuadrática) b) Bueno para diagramas de flujo con mucha interacción, ya que esta interacción se tiene en cuenta en el Jacobiano. Desventajas: a) Requiere unas estimaciones iniciales buenas b) Como las funciones no se conocen explícitamente, el Jacobiano se aproxima de forma numérica. Modelado y simulación en Ingeniería Química. Manuel Rodríguez
37 Ecuaciones algebraicas no lineales
Método de Broyden No calcula el jacobiano, lo aproxima empleando valores previos de x y f(x). W es la aproximación a la negativa de la inversa del jacobiano. Es una extensión del método de la secante (o método quasi-Newton) Modelado y simulación en Ingeniería Química. Manuel Rodríguez
38 Modelado y simulación en Ingeniería Química. Manuel Rodríguez
Ventajas: Sólo requiere una “pasada” por el diagrama de flujo en cada Tiene en cuenta de forma aproximada la interacción entre variables. Bueno para diagramas de flujo con alta interacción. Desventajas: Convergencia más lenta que Newton. Necesita más pasos que el método de Newton pero el costo computacional de cada método es Muy utilizado si el número de ecuaciones no es muy grande (<100) Utilizado para convergencia de reciclos en diagramas de flujo. Para los métodos de Newton o Broyden es deseable escoger el mínimo número de variables de corriente que rasgan todos los lazos (es decir buscar el menor número de ecuaciones). Si todos los bucles están rasgados la elección de las corrientes de rasgado no influye mucho en la convergencia de estos dos métodos. Modelado y simulación en Ingeniería Química. Manuel Rodríguez
39 Ecuaciones algebraicas no lineales
40 fsolve Solve a system of nonlinear equations
Ecuaciones algebraicas no lineales fsolve Solve a system of nonlinear equations for x, where x is a vector and F(x) is a function that returns a vector value. Syntax x = fsolve(fun,x0) x = fsolve(fun,x0,options) x = fsolve(fun,x0,options,P1,P2, ... ) [x,fval] = fsolve(...) [x,fval,exitflag] = fsolve(...) [x,fval,exitflag,output] = fsolve(...) [x,fval,exitflag,output,jacobian] = fsolve(...) Modelado y simulación en Ingeniería Química. Manuel Rodríguez
Descargar ppt "3. Métodos de resolución Ecuaciones algebraicas lineales"
Programación entera En muchos problemas reales las variables sólo pueden tomar valores enteros Ejemplos: decisiones sobre inversiones, compras, arranques,
4. ANÁLISIS FACTORIAL Introducción Modelo factorial ortogonal
Funciones lineales Las funciones de la forma y = ax + b, donde a, b R se llaman funciones lineales. Recorrido: R Recorrido: R (0, b): ordenada en el.
Otra forma es representando lo que hay de cada lado
Problemas resueltos /Aplicaciones de la derivada /Método de Newton
Métodos de Análisis Ingenieril
DESCRIPCION DE SISTEMAS
Tema III: Solución de ecuaciones no lineales

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución

 resolución 
 resolución