Daniel Quishpe Tarea3

ESCUELA SUPERIOR POLITÉCNICA DE CHIMBORAZO
FACULTAD DE MECÁNICA
ESCUELA DE INGENIERÍA AUTOMOTRIZ
ESCUELA SUPERIOR POLITÉCNICA

DE CHIMBORAZO
METODOS NUMERICOS
Nombre Daniel Quishpe

Semestre Quinto B
Código 2128
Fecha: domingo 17 de diciembre 2017
Riobamba-Ecuador
Escuela Superior Politécnica de Chimborazo

Ingeniería Automotriz
Deber No. 3
Métodos Numéricos
Ejercicio 1
Use la eliminación gaussiana con sustitución hacia atrás y aritmética de redondeo a
dos dígitos para resolver los sistemas lineales siguientes. (La solución exacta de cada
sistema es x1 = 1, x2 = -1, x3 = 3).
𝟒𝒙𝟏 − 𝒙𝟐 + 𝒙𝟑 = 𝟖
a) 𝟐𝒙𝟏 + 𝟓𝒙𝟐 + 𝟐𝒙𝟑 = 𝟑
𝒙𝟏 + 𝟐𝒙𝟐 + 𝟒𝒙𝟑 = 𝟏𝟏
𝟒𝒙𝟏 + 𝒙𝟐 + 𝟐𝒙𝟑 = 𝟗
b) 𝟐𝒙𝟏 + 𝟒𝒙𝟐 − 𝒙𝟑 =
−𝟓
𝒙𝟏 + 𝒙𝟐 − 𝟑𝒙𝟑 = −𝟗
Ejercicio 2
Use eliminación de Gauss para resolver los sistemas:
𝟖𝒙𝟏 + 𝟐𝒙𝟐 − 𝟐𝒙𝟑 = −𝟐
a) 𝟏𝟎𝒙𝟏 + 𝟐𝒙𝟐 + 𝟒𝒙𝟑 = 𝟒
𝟏𝟐𝒙𝟏 + 𝟐𝒙𝟐 + 𝟐𝒙𝟑 = 𝟔
𝒙𝟏 − 𝒙𝟐 + 𝟑𝒙𝟑 = 𝟐
b) 𝟑𝒙𝟏 − 𝟑𝒙𝟐 + 𝒙𝟑 = −𝟏
𝒙𝟏 + 𝒙𝟐 =𝟑
Ejercicio 3
Dado el sistema siguiente de ecuaciones
−𝟑𝒙𝟐 + 𝟕𝒙𝟑 = 𝟐
𝒙𝟏 + 𝟐𝒙𝟐 − 𝒙𝟑 = 𝟑
𝟓𝒙𝟏 − 𝟐𝒙𝟐 =𝟐
a) Emplee la eliminación de Gauss con pivoteo parcial para obtener cuales serían los
valores de x.
b) Sustituya sus resultados en las ecuaciones originales para efectos de comprobación.

Ejercicio 4
Dado el sistema siguiente de ecuaciones
𝟐𝒙𝟏 − 𝟔𝒙𝟐 − 𝒙𝟑 = −𝟑𝟖

−𝟑𝒙𝟏 − 𝒙𝟐 + 𝟕𝒙𝟑 = −𝟑𝟒
−𝟖𝒙𝟏 + 𝒙𝟐 − 𝟐𝒙𝟑 = −𝟐𝟎
a) Emplee la eliminación de Gauss con pivoteo parcial. Efectué todos los pasos del
cálculo.
Obtenemos el siguiente sistema:
b) Sustituya sus resultados en las ecuaciones originales para comprobar sus respuestas.
Ejercicio 5
Emplee la eliminación de Gauss - Jordán para resolver el sistema siguiente:
𝟐𝒙𝟏 + 𝒙𝟐 − 𝒙𝟑 = 𝟏
𝟓𝒙𝟏 + 𝟐𝒙𝟐 + 𝟐𝒙𝟑 = −𝟒
𝟑𝒙𝟏 + 𝒙𝟐 + 𝒙𝟑 = 𝟓
Ejercicio 6
Resuelva el sistema:
𝒙𝟏 + 𝒙𝟐 − 𝒙𝟑 = −𝟑
𝟔𝒙𝟏 + 𝟐𝒙𝟐 + 𝟐𝒙𝟑 = 𝟐
−𝟑𝒙𝟏 + 𝟒𝒙𝟐 + 𝒙𝟑 = 𝟏
a) Por eliminación de Gauss simple.
b) Por el método de Gauss – Jordán.
c) Por eliminación de Gauss con pivoteo parcial.

Ejercicio 7
Resuelva los sistemas de ecuaciones por medio de la descomposición LU (Crout y
Doolittle).
𝟖𝒙𝟏 + 𝟒𝒙𝟐 − 𝒙𝟑 = 𝟏𝟏
a) −𝟐𝒙𝟏 + 𝟓𝒙𝟐 + 𝒙𝟑 = 𝟒
𝟐𝒙𝟏 − 𝒙𝟐 + 𝟔𝒙𝟑 = 𝟕
 Método de Crout.
Ejercicio 8
Las ecuaciones algebraicas lineales pueden surgir al resolver ecuaciones diferenciales.

Por ejemplo, la ecuación diferencial siguiente proviene de un balance de calor para
una barra larga y delgada (véase la figura):
𝒅𝟐𝑻
′
𝒅𝒙𝟐 + 𝒉 (𝑻𝒂 − 𝑻) = 𝟎 (𝟏)
Donde T = temperatura (°C), x = distancia a lo largo de barra (m), 𝒉′ = coeficiente de

transferencia de calor entre la barra y el aire del ambiente (𝒎−𝟐), y 𝑻𝒂 = temperatura del
aire circundante (°C). Esta ecuación se transforma en un conjunto de ecuaciones
algebraicas
lineales por medio del uso de una aproximación en diferencias finitas divididas para la
segunda derivada,
𝒅𝟐𝑻
= 𝑻𝒊+𝟏 + 𝑻𝒊 + 𝑻𝒊−𝟏
𝒅𝒙𝟐
∆𝒙𝟐
Donde 𝑻𝒊 denota la temperatura en el nodo i. Esta aproximación se sustituye en la ecuación
(1) y se obtiene
−𝑻𝒊−𝟏 + (𝟐 + 𝒉′∆𝒙𝟐) 𝑻𝒊−𝑻𝒊+𝟏 = 𝒉′∆𝒙𝟐𝑻𝒂
Se puede plantear esta ecuación para cada uno de los nodos interiores de la barra, lo que
resulta en un sistema tri diagonal de ecuaciones. Los nodos primero y último en los
extremos de la barra están fijos por las condiciones de frontera.
a) Para una barra de 10 m con 𝑻𝒂 = 𝟐𝟎, 𝑻 (𝒙 = 𝟎) = 𝟒𝟎, 𝑻 (𝒙 = 𝟏𝟎) = 𝟐𝟎𝟎 𝒚 𝒉′ = 𝟎. 𝟎𝟐.

Desarrolle una solución numérica, con el uso de una solución en diferencias finitas con
cuatro nodos interiores según se muestra en la figura (∆𝒙 = 𝟐𝒎).
−𝑇𝑖−1 + (2 + ℎ′∆𝑥2) 𝑇𝑖−𝑇𝑖+1 = ℎ′∆𝑥2𝑇𝑎
Para 𝑖 = 1
−𝑇0 + (2 + 0.02(2)2)𝑇1−𝑇2 = 0.02(2)2(20)
−40 + 2.08𝑇1−𝑇2 = 1.6
2.08𝑇1−𝑇2 = 41.6
Para 𝑖 = 2
−𝑇1 + 2.08𝑇2−𝑇3 = 1.6
Para 𝑖 = 3
−𝑇2 + 2.08𝑇3−𝑇4 = 1.6
Para 𝑖 = 4
−𝑇3 + 2.08𝑇4−𝑇5 = 1.6

−𝑇3 + 2.08𝑇4 − 200 = 1.6
−𝑇3 + 2.08𝑇4 = 201.6
Obtenemos el siguiente sistema:
A es una matriz banda, y es simétrica y definida positiva, entonces existe 𝐿 𝑡𝑞 𝐴 = 𝐿𝐿 𝑇

Resolvemos aplicando la variante del método de Choleski para este tipo de matrices.
Ejercicio 9
Determinar las corrientes de malla en el circuito de la figura, utilice uno de los

métodos vistos.
Ejercicio 10
Use el método de Cholesky para resolver los sistemas lineales siguientes:
𝟒𝒙𝟏 – 𝟐𝒙𝟐 =𝟎
a) −𝟐𝒙𝟏 + 𝟒𝒙𝟐 − 𝒙𝟑 = 𝟎. 𝟓
−𝒙𝟐 + 𝟒𝒙𝟑 = 𝟏
𝟓𝒙𝟏 + 𝒙𝟐 + 𝟐𝒙𝟑 − 𝒙𝟒 = 𝟏
b) 𝒙𝟏 + 𝟕𝒙𝟐 + 𝟑𝒙𝟒 = 𝟐
𝟐𝒙𝟏 + 𝟓𝒙𝟑 + 𝒙𝟒 = 𝟑
−𝒙𝟏 + 𝟑𝒙𝟐 + 𝒙𝟑 + 𝟖𝒙𝟒 = 𝟒
Ejercicio 11
Use el método de Crout y Doolittle para resolver los sistemas lineales siguientes:
𝟒𝒙𝟏 + 𝒙𝟐 − 𝒙𝟑 = 𝟖
a) 𝟐𝒙𝟏 + 𝟓𝒙𝟐 =𝟓
𝟑𝒙𝟏 + 𝟖𝒙𝟐 + 𝟗𝒙𝟑 = 𝟎
1. Método de Crout.
𝟒𝒙𝟏 + 𝟓𝒙𝟐 + 𝟐𝒙𝟑 − 𝒙𝟒 = 𝟑

b) 𝟓𝒙𝟏 + 𝟖𝒙𝟐 + 𝟕𝒙𝟑 + 𝟔𝒙𝟒 = 𝟐
𝟑𝒙𝟏 + 𝟕𝒙𝟐 − 𝟒𝒙𝟑 − 𝟐𝒙𝟒 = 𝟎
−𝒙𝟏 + 𝟔𝒙𝟐 − 𝟐𝒙𝟑 + 𝟓𝒙𝟒 = 𝟏
2. Método de Doolittle
Ejercicio 15
Resuelva los siguientes sistemas de ecuaciones lineales mediante los métodos de Gauss –
Seidel y de Jacobi, verificar la convergencia de cada método.
𝟓𝒙𝟏 − 𝒙𝟐 − 𝒙𝟑 = 𝟑 5 −1 −1
a) 𝒙𝟏 − 𝒙𝟐 + 𝟐𝒙𝟑 = 𝟎 𝐴= 1 −1 2
𝟑𝒙𝟏 − 𝒙𝟐 + 𝟐𝒙𝟑 = 𝟒 3 −1 2
Es convergente si y solo si, el radio espectral de la matriz 𝐾 −1𝑀 𝑒𝑠 < 1. → 𝜌(𝐾−1𝑀) < 1
Calculamos los valores propios de la matriz (𝐾−1 ∗ 𝑀).
det(𝐾−1𝑀 − 𝜆𝐼) = 0
Hallamos el radio espectral:
𝜌(𝐾−1𝑀) = max|𝜆𝑖| = 1.1546 > 1 por lo tanto no converge por Jacobi.

𝑖=1
Además, matriz A no es simétrica definida positiva, entonces no converge por Gauss – Seidel.
𝟐𝒙𝟏 + 𝒙𝟐 + 𝒙𝟑 = 𝟒 2 1 1
b) 𝒙𝟏 + 𝟐𝒙𝟐 + 𝒙𝟑 = 𝟒 𝐴= 1 2 1
𝒙𝟏 + 𝒙𝟐 + 𝟐𝒙𝟑 = 𝟒 1 1 2
VERIFICAMOS CONVERGENCIA
Calculamos los valores propios de la matriz (𝐾−1 ∗ 𝑀).
det(𝐾−1𝑀 − 𝜆𝐼) = 0
Hallamos el radio espectral:
𝜌(𝐾−1𝑀) = max|𝜆𝑖| = 1 por lo tanto no existe convergencia por Jacobi.

𝑖=1,2
La matriz A es simétrica definida positiva, entonces converge por el método de Gauss – Seidel.
APLICAMOS GAUSS – SEIDEL.
𝑥1𝑘 = 4 − 𝑥1𝑘 − 𝑥3𝑘−1 4 − 𝑥1𝑘 − 𝑥2𝑘

4 − 𝑥2𝑘−1 − 𝑥3𝑘−1
𝑥2𝑘 = 𝑥3𝑘 =
2 2 2
� � � � ‖ 𝑥⃗𝑘 − 𝑥⃗𝑘−1 ‖ <

� 0.01
� � � ∞
1 2 3
0 0 0 0 ---
1 2.0000 1.0000 0.5000 2.0000
2 1.2500 1.1250 0.8125 0.7500
3 1.0313 1.0781 0.9453 0.2187
4 0.9883 1.0332 0.9893 0.0449
5 0.9888 1.0110 1.0001 0.0222
6 0.9944 1.0027 1.0014 0.0083
𝑥⃗ = ( 0.9944 , 1.0027 , 1.0014 ) ← SOLUCIÓN

Ejercicio 13
Use el método de Cholesky para resolver los sistemas lineales siguientes:
𝟒𝒙𝟏 – 𝟐𝒙𝟐 =𝟎
a) −𝟐𝒙𝟏 + 𝟒𝒙𝟐 − 𝒙𝟑 = 𝟎. 𝟓
−𝒙𝟐 + 𝟒𝒙𝟑 = 𝟏
𝟓𝒙𝟏 + 𝒙𝟐 + 𝟐𝒙𝟑 − 𝒙𝟒 = 𝟏
b) 𝒙𝟏 + 𝟕𝒙𝟐 + 𝟑𝒙𝟒 = 𝟐
𝟐𝒙𝟏 + 𝟓𝒙𝟑 + 𝒙𝟒 = 𝟑
−𝒙𝟏 + 𝟑𝒙𝟐 + 𝒙𝟑 + 𝟖𝒙𝟒 = 𝟒
ES EL MISMO EJERCICIO N° 10
Ejercicio 14
Use el método de Crout y Doolittle para resolver los sistemas lineales siguientes:
𝟒𝒙𝟏 + 𝒙𝟐 − 𝒙𝟑 = 𝟖
a) 𝟐𝒙𝟏 + 𝟓𝒙𝟐 =𝟓
𝟑𝒙𝟏 + 𝟖𝒙𝟐 + 𝟗𝒙𝟑 = 𝟎
𝟒𝒙𝟏 + 𝟓𝒙𝟐 + 𝟐𝒙𝟑 − 𝒙𝟒 = 𝟑
b) 𝟓𝒙𝟏 + 𝟖𝒙𝟐 + 𝟕𝒙𝟑 + 𝟔𝒙𝟒 = 𝟐

𝟑𝒙𝟏 + 𝟕𝒙𝟐 − 𝟒𝒙𝟑 − 𝟐𝒙𝟒 = 𝟎
−𝒙𝟏 + 𝟔𝒙𝟐 − 𝟐𝒙𝟑 + 𝟓𝒙𝟒 = 𝟏
Ejercicio 15
Resuelva los siguientes sistemas de ecuaciones lineales mediante los métodos de Gauss –
Seidel y de Jacobi, verificar la convergencia de cada método.
𝟓𝒙𝟏 − 𝒙𝟐 − 𝒙𝟑 = 𝟑
a) 𝒙𝟏 − 𝒙𝟐 + 𝟐𝒙𝟑 = 𝟎
𝟑𝒙𝟏 − 𝒙𝟐 + 𝟐𝒙𝟑 = 𝟒
𝟐𝒙𝟏 + 𝒙𝟐 + 𝒙𝟑 = 𝟒
b) 𝒙𝟏 + 𝟐𝒙𝟐 + 𝒙𝟑 = 𝟒
𝒙𝟏 + 𝒙𝟐 + 𝟐𝒙𝟑 = 𝟒

Daniel Quishpe Tarea3

Cargado por

Copyright:

Formatos disponibles

Daniel Quishpe Tarea3

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Daniel Quishpe Tarea3

Cargado por

Copyright:

Formatos disponibles

ESCUELA SUPERIOR POLITÉCNICA DE CHIMBORAZO

ESCUELA SUPERIOR POLITÉCNICA

Nombre Daniel Quishpe

Fecha: domingo 17 de diciembre 2017

Escuela Superior Politécnica de Chimborazo

b) Sustituya sus resultados en las ecuaciones originales para efectos de comprobación.

𝟐𝒙𝟏 − 𝟔𝒙𝟐 − 𝒙𝟑 = −𝟑𝟖

Obtenemos el siguiente sistema:

b) Por el método de Gauss – Jordán.

c) Por eliminación de Gauss con pivoteo parcial.

Las ecuaciones algebraicas lineales pueden surgir al resolver ecuaciones diferenciales.

Donde T = temperatura (°C), x = distancia a lo largo de barra (m), 𝒉′ = coeficiente de

a) Para una barra de 10 m con 𝑻𝒂 = 𝟐𝟎, 𝑻 (𝒙 = 𝟎) = 𝟒𝟎, 𝑻 (𝒙 = 𝟏𝟎) = 𝟐𝟎𝟎 𝒚 𝒉′ = 𝟎. 𝟎𝟐.

−𝑇𝑖−1 + (2 + ℎ′∆𝑥2) 𝑇𝑖−𝑇𝑖+1 = ℎ′∆𝑥2𝑇𝑎

−𝑇3 + 2.08𝑇4−𝑇5 = 1.6

A es una matriz banda, y es simétrica y definida positiva, entonces existe 𝐿 𝑡𝑞 𝐴 = 𝐿𝐿 𝑇

Determinar las corrientes de malla en el circuito de la figura, utilice uno de los

Use el método de Cholesky para resolver los sistemas lineales siguientes:

𝟒𝒙𝟏 + 𝟓𝒙𝟐 + 𝟐𝒙𝟑 − 𝒙𝟒 = 𝟑

Calculamos los valores propios de la matriz (𝐾−1 ∗ 𝑀).

Hallamos el radio espectral:

𝜌(𝐾−1𝑀) = max|𝜆𝑖| = 1.1546 > 1 por lo tanto no converge por Jacobi.

Calculamos los valores propios de la matriz (𝐾−1 ∗ 𝑀).

Hallamos el radio espectral:

𝜌(𝐾−1𝑀) = max|𝜆𝑖| = 1 por lo tanto no existe convergencia por Jacobi.

APLICAMOS GAUSS – SEIDEL.

𝑥1𝑘 = 4 − 𝑥1𝑘 − 𝑥3𝑘−1 4 − 𝑥1𝑘 − 𝑥2𝑘

� � � � ‖ 𝑥⃗𝑘 − 𝑥⃗𝑘−1 ‖ <

𝑥⃗ = ( 0.9944 , 1.0027 , 1.0014 ) ← SOLUCIÓN

Use el método de Cholesky para resolver los sistemas lineales siguientes:

b) 𝟓𝒙𝟏 + 𝟖𝒙𝟐 + 𝟕𝒙𝟑 + 𝟔𝒙𝟒 = 𝟐

También podría gustarte