Algoritmos101

Cursos Python

Etiqueta: matriciales iterativos

Métodos Iterativos – Sistemas de Ecuaciones

  • 1Eva_2021PAOI_T2 Atención hospitalaria con medicamentos limitados

    1ra Evaluación 2021-2022 PAO I. 6/Julio/2021

    Tema 2 (35 puntos) Durante el año 2020, ante el aumento de atención hospitalaria estatal en la región y el limitado acceso a medicamentos, como una primera estrategia de manejo de recursos se derivan el exceso de pacientes hacia la atención en hospitales privados. hospital dibujo con ambulancia

    En la tabla  se muestra la cantidad de los tres medicamentos (mg, ml) que se administran al atender a cada paciente clasificado por grupo etario: niños, adolescentes, adultos y adultos mayores.

    También se dispone del total de medicamentos existente en bodegas en cada semana.

    Niños Adolescentes Adultos Adultos Mayores Medicamentos /semana
    Medicamento_A 0.3 0.4 1.1 4.7 3500
    Medicamento_B 1 3.9 0.15 0.25 3450
    Medicamento_C 0 2.1 5.6 1.0 6500

    Es de interés conocer la cantidad de pacientes de cada grupo que se podría atender con los recursos disponibles.

    a.  Realice el planteamiento de un sistema de ecuaciones que permita determinar la cantidad máxima de pacientes de cada grupo etario que podrían ser atendidos usando todos los medicamentos disponibles.

    Una vez planteadas las ecuaciones, se le indica que la capacidad K para pacientes niños sea una variable libre, por consumir menos recursos y se podrían derivar al sistema privado.

    b.  Escriba el conjunto de soluciones posibles en función de la variable libre, considerando la cantidad de niños a atender como máximo de K=100.

    c. Determine la capacidad de atención usando un método Iterativo con una tolerancia de 10-2. Realice tres iteraciones completas y revise la convergencia del método. Se estima atender al inicio de semana al menos 100 pacientes de cada grupo.

    d. Suponga que la cantidad de pacientes en cada grupo para una semana dada es: [350, 1400, 1500, 1040]. ¿Hay suficiente cantidad de medicamentos para atender el promedio actual de pacientes? Analice y describa los resultados encontrados.

    e. Si se decide vacunar primero a todos los niños, entonces ya no requieren atención hospitalaria (K=0) ¿Cuál es el número máximo de pacientes de cada grupo que podría incrementarse dadas las condiciones actuales? Resuelva usando un método directo.

    Rúbrica: literal a (3 puntos), literal b (2 puntos), pivoteo por filas(5 puntos), iteraciones (10 puntos), análisis de convergencia (4 puntos), literal d (5 puntos) literal e (6 puntos)

    Referencias:
    - BBC News Mundo. El país que está vacunando contra el covid-19 primero a los jóvenes y no a los ancianos. 16/enero/2021. https://youtu.be/oo2itoBBwyY

    - Manejo clínico de la COVID-19, orientaciones provisionales 27/mayo/2020. https://apps.who.int/iris/bitstream/handle/10665/332638/WHO-2019-nCoV-clinical-2020.5-spa.pdf

  • 1Eva_IIT2019_T3 Circuito eléctrico

    1ra Evaluación II Término 2019-2020. 26/Noviembre/2019. MATG1013

    Tema 3. (30 puntos) El sistema de ecuaciones que sigue se generó por medio de aplicar la ley de malla de corriente al circuito de la figura. circuito eléctrico con resistencias y fuentes

    \begin{cases} 55 I_1 - 25 I_4 = -200 \\ -37 I_3 - 4 I_4 = -250 \\ -25 I_1 - 4 I_3 + 29 I_4 = 100 \end{cases}

    a) Use el método de eliminación de Gauss para calcular I1, I3, I4, I1 observando que
    I2 = -10

    b) Encuentre la norma infinita de la matriz de transición T en el método de Jacobi y comente.

    c) Con el método de Gauss-Seidel realice tres iteraciones comenzando con el vector cero. Además en la tercera iteración, encuentre una cota para el error relativo.

    Rúbrica: literal a (12 puntos), literal b (6 puntos), literal c (12 puntos)

    A = [[ 55.0, 0,  0, -25],
     [  0  , 0,-37,  -4],
     [-25  , 0, -4,  29],
     [  0  ,  1, 0,   0]]

B = [-200,-250,100,-10]

  • 1Eva_IT2018_T3 Temperatura en nodos de placa

    1ra Evaluación I Término 2018-2019. 26/junio/2018. MATG1013

    Tema 3. (25 puntos). La temperatura en los nodos de la malla de una placa se puede calcular con el promedio de las temperaturas de los 4 nodos vecinos de la izquierda, derecha, arriba y abajo. nodos en una placa rectangular, metodo iterativo

    Una placa cuadrada de 3 m de lado tiene la temperatura en los nodos de los bordes como se indica en la figura,

    a) Plantee el sistema de ecuaciones y resuelva con eliminación de Gauss, encontrar a, b, c, d.

    b) Encuentre la matriz T de Jacobi y comente sobre la convergencia

    c) Con X[0]=[a=60, b=40, c=70, d=50], realice 3 iteraciones, estime el error, comente.

    d) Con la tercera iteración calcule el residuo y encuentre una cota del error absoluto y error relativo

    Rúbrica: Literal a (5 puntos), literal b (5 puntos), literal c (10 puntos), literal d (5 puntos).

    Referencia: Ejercicio 12.39 p339 Steven C. Chapra. Numerical Methods 7th Edition.

  • 1Eva_IT2016_T4_MN Conceptos

    1ra Evaluación I Término 2016-2017. 28/junio/2016. ICM02188 Métodos Numéricos

    Tema 4. (25 puntos) Responda las siguientes preguntas y justifique la respuesta

    a) Si la ‖Tj‖∞ > 1 , entonces el método de Jacobi no converge

    b) Si f ∈ C2[a,b] y p ∈ [a,b], tal que f(p)=0 y f '(p)≠0,
    entonces existe δ > 0 tal que el método de Newton converge para cualquier
    p0 ∈ [p - δ, p + δ]

    Rúbrica: literal a) falso  (6 puntos), Justificación (6 puntos), literal b) verdadero (6 puntos), demostración (7 puntos)

  • 3Eva_IIT2009_T3 Sistema de ecuaciones

    3ra Evaluación II Término 2009-2010. 23/Febrero/2010. ICM00158

    Tema 4. (25 puntos) Enunciar el teorema de convergencia del método iterativo para resolver un sistema de ecuaciones lineales AX=B.

    Exponer el método iterativo de Gauss-Seidel para sistemas ecuaciones lineales.

    Construir un ejemplo de un sistema de 3x3, cuya diagonal principal sea estrictamente dominante y realizar cuatro iteraciones con el método de Gauss-Seidel, comenzando con el vector cero.

  • 1Eva_IIT2017_T3 Circuito eléctrico

    1ra Evaluación II Término 2017-2018. 28/Noviembre/2017. MATG1013

    Tema 3. (25 puntos) El sistema de ecuaciones que sigue se generó por medio de aplicar la ley de malla de corriente al circuito de la figura.

    \begin{cases} 55 I_1 - 25 I_4 = -200 \\ -37 I_3 - 4 I_4 = -250 \\ -25 I_1 - 4 I_3 + 29 I_4 = 100 \\ I_2 = -10 \end{cases}

    circuito Electrico 01

    a) Use el método de eliminación de Gauss para resolver el sistema

    b) Use el método de Jacobi y determine el número de iteraciones para ε=0.01

    c) Si el coeficiente 55 se cambia a 54.9, encuentre el error relativo de la aproximación en el literal a.

    Rúbrica: Aplicación del método de eliminación de Gauss hasta 10%, Uso del método de Jacobi hasta 5% y determinación del número de iteraciones hasta 5%, Calculo del residuo y cota del error relativo hasta 5%.

    A = [[ 55.0, 0,  0, -25],
     [  0  , 0,-37,  -4],
     [-25  , 0, -4,  29],
     [  0  ,  1, 0,   0]]

B = [-200,-250,100,-10]

  • 1Eva_IT2017_T4 Componentes eléctricos

    1ra Evaluación I Término 2017-2018. 26/junio/2017. MATG1013

    Tema 4. (25puntos)

    https://es.dreamstime.com/tablero-electr%C3%B3nico-de-la-tv-image120402048
    Tablero electrónico de la TV. Sistema, tarjeta

    Un supervisor revisa la producción de tres tipos de componentes eléctricos.

    Para ellos se requieren tres clases de materiales como se indica en la tabla adjunta:

    Material 1 Material 2 Material 3
    Componente 1 5 9 3
    Componente 2 7 7 16
    Componente 3 9 3 4

    a) Si cada semana se dispone de un total de 945 gramos de material 1, 987 gramos de material 2 y 1049 gramos de material 3, ¿Cuántos componentes a lo sumo pueden producirse por semana? (Solo plantear)

    b) Si se resuelve con el método de eliminación de Gauss, ¿Cuántas multiplicaciones/divisiones como máximo se realizan?

    c) Si se resuelve con el método de Jacobi, encuentre la norma infinita de T y comente sobre la convergencia.

    d) Resuelva utilizando el método de Gauss-Seidel, realice tres iteraciones y estime el error de la tercera iteración.

    Rúbrica: Planteo hasta 5 puntos, Número de multiplicaciones hasta 5 puntos, ‖T‖∞ hasta 10 (con filas ordenadas), Iteraciones con Gauss Seidel con la estimación del error hasta 5 puntos.

  • 1Eva_IT2016_T2_MN Organismos patógenos en lago

    1ra Evaluación I Término 2016-2017. 28/junio/2016. ICM02188 Métodos Numéricos

    Tema 2. (25 puntos) Tres organismos patógenos decaen en forma exponencial en aguas de un lago de acuerdo con el siguiente modelo:

    p(t) = A e^{-1.5t} + B e^{-0.3t} + C e^{-0.05t}

    Estime la población inicial de cada organismo, dadas las mediciones siguientes:

    Tiempo
    (horas)    		 0.5 1 2 3 4
    Población
    (miles)    		 6.0 4.4 3.2 2.7 2.2

    a) Seleccione los tres primeros puntos y plantee un sistema de 3 ecuaciones.
    b) Con el método de Jacobi encuentre la matriz T y comente.
    c) Con el método de Gauss Seidel realice tres iteraciones y estime el error.

    Rúbrica: Ecuaciones (5 puntos), matriz (5 puntos), comentario (6 puntos), Iteraciones (5 puntos), estimación del error (4 puntos).

    Referencia: Cuales son los agentes patógenos del agua.
    https://www.ecomol.es/tratamientos/cuales-son-los-agentes-patogenos-del-agua/

     

  • 1Eva_IIT2014_T2 Componentes eléctricos

    1ra Evaluación II Término 2014-2015. 9/Diciembre/2013. ICM00158

    Tema 2. Un ingeniero eléctrico supervisa la producción de tres tipos de componentes eléctricos. componentes eléctricos

    Para cada componente se se requieren tres clases de materiales:
    metal 1, metal 2 y caucho.

    Gramos por componente Metal 1 Metal 2 Caucho
    Componente 1 15 0.25 1.0
    Componente 2 17 0.33 1.2
    Componente 3 19 0.42 1.6

    Se requieren disponer de componentes con el mismo desempeño, pero de menor tamaño y no se dispone de mas gramos de material que:

    materiales = [2.63, 0.0534, 0.202]

    a. Plantee el sistema de ecuaciones

    b. Utilice el método de eliminación de Gauss para resolver el sistema

    c. Encuentre la matriz de Jacobi y comente sobre la convergencia

    d. Realice tres iteraciones con Gauss Seidel y estime el error de la segunda iteración.

    e. Encuentre el número de condición y comente.

    A = np.array([[15, 0.25, 1.0],
              [17, 0.33, 1.2],
              [19, 0.42, 1.6]])
B = np.array([2.63, 0.0534, 0.202])

  • 1Eva_IT2012_T2_MN Modelo Leontief

    1ra Evaluación I Término 2012-2013. 3/Julio/2012. ICM02188 Métodos Numéricos

    TEMA 2. (35 puntos) La matriz insumo-producto propuesto por W. Leontief, es un modelo muy importante en Economía. Flores Ecuador 01

    En ésta matriz se describe la producción de los diferentes sectores económicos y la demanda interna para satisfacer a estos mismos sectores, expresada como una fracción de su producción.

    Ejemplo: Suponga que hay tres sectores
    A: agricultura,
    M: manufactura
    S: servicios
    y su demanda interna es:

    Matriz T Producción
    A M S
    Demanda A 0.40 0.03 0.02
    Interna M 0.06 0.37 0.10
    S 0.12 0.15 0.19

    Sea T el nombre de esta matriz.

    Para los datos propuestos, en la primera columna de la matriz T, el sector A requiere 0.4 de su propia producción, 0.06 del sector M, y 0.12 del sector S, etc.

    Sea D el vector de demanda externa de cada sector, y X el vector de la producción total de cada sector, requerida para satisfacer las demandas interna y externa:

    D = \begin{pmatrix} 80\\ 140\\200 \end{pmatrix} X = \begin{pmatrix} x_1 \\ x_2\\x3 \end{pmatrix}

    en donde x1, x2, x3 representan la producción total de cada sector.

    Entonces la ecuación X = TX + D proporciona la producción total X para satisfacer las demandas externa e interna.

    a) Formule un método iterativo en notación vectorial para usar la ecuación anterior. Indique cual es el nombre de la matriz T. Analice esta matriz y determine si el método iterativo es convergente.

    b) Comience con un vector inicial X = [200, 200, 200]T realice las iteraciones necesarias hasta que la norma de la diferencia entre dos vectores consecutivos sea menor a 1.

    Use la norma de fila.

    Referencia: Modelo Input-Output. https://es.wikipedia.org/wiki/Modelo_Input-Output, https://es.wikipedia.org/wiki/Wassily_Leontief

    T = np.array([[0.40, 0.03, 0.02],
              [0.06, 0.37, 0.10],
              [0.12, 0.15, 0.19]])

D = np.array([80.0, 140.0, 200.0],dtype=float)

Xa = np.array([200.0,200.0,200.0])