Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

T1-Sistemas Lineales-P1 MAT 2106

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 5

MÉTODOS NUMÉRICOS II MAT – 2106 “A”

SISTEMAS LINEALES
PARTE I
1. CONCEPTO. – Un sistema lineal, es un modelo matemático que está representado por un
sistema con “n” ecuaciones lineales con “m” incógnitas. Como ser:

𝐚 𝐗 + 𝐛𝟏 𝐘 = 𝐂𝟏
{ 𝟏
𝐚𝟐 𝐗 + 𝐛𝟐 𝐘 = 𝐂𝟐

1.1. MATRIZ. – Es un arreglo rectangular de números reales dispuestos en filas y columnas


denotados como:
𝐀𝐧 ×𝐦 = [𝐚𝐢𝐣]
En forma extendida:

𝐚𝟏𝟏 𝐚𝟏𝟐 ⋯ 𝐚𝟏𝐦


𝐚𝟐𝟏 𝐚𝟐𝟐 ⋯ 𝐚𝟐𝐦
𝐀𝐧 ×𝐦 = ( )
⋮ ⋮ ⋮
𝐚𝐧𝟏 𝐚𝐧𝟐 ⋯ 𝐚𝐧𝐦

1.2. CLASES DE MATRICES


 Matriz Cuadrada: Nº filas = Nº columnas.
 Matriz Identidad: Matriz cuadrada con diagonal unitaria, el resto lleno de 0’s.
 Matriz Simétrica: Una fila es igual a una columna correspondiente.
 Matriz Triangular Inferior: La parte superior derecha contiene solo ceros.
 Matriz Triangular Superior: La parte inferior izquierda contiene solo ceros
1.3. OPERACIONES CON MATRICES
 Adición algebraica: suma de matrices.
−3 1 7 −1 4 0
( )+( )=( )
0 4 2 5 2 9
 Multiplicación de matrices: fila por columna.
(−3 ∗ 7) + (1 ∗ 2) (−3 ∗ −1) + (1 ∗ 5)
(−3 1) × (7 −1) = ( ) = (24 8 )
0 4 2 5 (0 ∗ 7) + (4 ∗ 2) (0 ∗ −1) + (4 ∗ 5) 8 20
 Determinante de una matriz: por pivote (existen otros métodos).
4 −1 3
|A| = (2 2 3 4 −1 |𝐀| = 𝟓𝟑
3 0) = 3 [ ] +4[ ] = −3 + 56 →
1 1 2 3
1 1 4
 Matriz Transpuesta: llevar las filas a columnas y viceversa
4 −1 3 4 2 1
A = (2 3 0) → AT = (−1 3 1)
1 1 4 3 0 4

SEM - II/2020
MÉTODOS NUMÉRICOS II MAT – 2106 “A”

 Convertir a Matriz Triangular Superior:


4 −1 3 4 −1 3 ->(+) 4 −1 3
A = (2 3 )
0 ->(+) → (2 3 )
0 *(-2) → (0 −7 3 ) ->(+)
1 1 4 *(-2) 0 1 −8 0 1 −8 *(7)
𝟒 −𝟏 𝟑
𝐀 = (𝟎 −𝟕 𝟑)
𝟎 𝟎 −𝟓𝟑
 Convertir a Matriz Triangular Inferior:
13 7 𝟓𝟑
4 −1 3 *(-4/3) − 0 *(-9/7) 𝟎 𝟎
A = (2 3 0) → ( 3 3 ) → ( 𝟕 )
2 3 0 ->(+) 𝟐 𝟑 𝟎
1 1 4 ->(+)
1 1 4 𝟏 𝟏 𝟒

2. FACTORIZACIÓN DE UNA MATRIZ. - Descomponer una matriz en dos factores.


2.1. FACTORIZACIÓN L.U.- Tiene la siguiente forma:

𝐀𝐟 ×𝐜 = 𝑳 × 𝑼
𝐷𝑜𝑛𝑑𝑒: 𝐿 ∶ 𝑚𝑎𝑡𝑟𝑖𝑧 𝑡𝑟𝑖𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟
𝑈 ∶ 𝑚𝑎𝑡𝑟𝑖𝑧 𝑡𝑟𝑖𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑠𝑢𝑝𝑒𝑟𝑖𝑜𝑟

2.1.1. L.U. DOOLITTLE: La matriz L tendrá la diagonal unitaria.

𝑨 = 𝑳 × 𝑼

4 6 −1 1 0 0 𝑈11 𝑈12 𝑈13


( 2 1 4 ) = (𝐿21 1 0) × ( 0 𝑈22 𝑈23 )
−3 3 0 𝐿31 𝐿32 1 0 0 𝑈33

Multiplicar las matrices L×U: fila por columna.

4 6 −1 𝑈11 𝑈12 𝑈13


( 2 1 4 ) = ( 𝐿21 ∙ 𝑈11 𝐿21 ∙ 𝑈12 + 𝑈22 𝐿21 ∙ 𝑈13 + 𝑈23 )
−3 3 0 𝐿31 ∙ 𝑈11 𝐿31 ∙ 𝑈12 + 𝐿32 ∙ 𝑈22 𝐿31 ∙ 𝑈13 + 𝐿32 ∙ 𝑈23 + 𝑈33

Igualar a la matriz A: Formando ecuaciones obtener los valores numéricos de los


componentes de L y U según su posición (color) en A. Los primeros valores que se tienen
directamente son:
𝑈11 = 4 ; 𝑈12 = 6 ; 𝑈13 = −1

𝐿21 ∙ 𝑈11 = 2 𝐿31 ∙ 𝑈11 = −3 𝐿21 ∙ 𝑈12 + 𝑈22 = 1


𝐿21 ∙ 4 = 2 𝐿31 ∙ 4 = −3 1
∙ 6 + 𝑈22 = 1
𝐿21 = 1/2 𝐿31 = −3/4 2
𝑈22 = −2

𝐿21 ∙ 𝑈13 + 𝑈23 = 4 𝐿31 ∙ 𝑈12 + 𝐿32 ∙ 𝑈22 = 3 𝐿31 ∙ 𝑈13 + 𝐿32 ∙ 𝑈23 + 𝑈33 = 0
1 3 3 15 9
∙ (−1) + 𝑈23 = 4 − ∙ 6 + 𝐿32 ∙ (−2) = 3 − ∙ (−1) + (− ) ∙ + 𝑈33 = 0
2 4 4 4 2

SEM - II/2020
MÉTODOS NUMÉRICOS II MAT – 2106 “A”

1 9 3 135
𝑈23 = 4 + ( ) −2 ∙ 𝐿32 = 3 + 𝑈33 = − +
2 2 4 8
𝑈23 = 9/2 15 1 129
𝐿32 = ∙ (− ) 𝑈33 =
2 2 8
𝐿32 = −15/4

Solución: Remplazar los valores obtenidos a la matriz L y U, entonces se tiene:

4 6 −1 1 0 0 4 6 −1
( 2 1 4 ) = ( 1/2 1 0) × (0 −2 9/2 )
−3 3 0 −3/4 −15/4 1 0 0 129/8
Sol.
Verificación: Multiplicar las matrices L×U y para comprobar debe ser igual a A.
4 6 −1 1 0 0 4 6 −1
( 2 1 4 ) = ( 1/2 1 0) × (0 −2 9/2 )
−3 3 0 −3/4 −15/4 1 0 0 129/8
4 6 −1 4 6 −1
( 2 1 4) = ( 2 1 4)
−3 3 0 −3 3 0 L.q.q.d
.
2.1.2. L.U. DE CROUT: La matriz U tendrá la diagonal unitaria.

𝑨 = 𝑳 × 𝑼

3 −2 0 𝐿11 0 0 1 𝑈12 𝑈13


(0 4 1) = (𝐿21 𝐿22 0 ) × (0 1 𝑈23 )
5 2 1 𝐿31 𝐿32 𝐿33 0 0 1
Multiplicar las matrices L×U: fila por columna.

3 −2 0 𝐿11 𝐿11 ∙ 𝑈12 𝐿11 ∙ 𝑈13


(0 4 1) = (𝐿21 𝐿21 ∙ 𝑈12 + 𝐿22 𝐿21 ∙ 𝑈13 + 𝐿22 ∙ 𝑈23 )
5 2 1 𝐿31 𝐿31 ∙ 𝑈12 + 𝐿32 𝐿31 ∙ 𝑈13 + 𝐿32 ∙ 𝑈23 + 𝐿33

Igualar a la matriz A: Formando ecuaciones obtener los valores numéricos de los


componentes de L y U según su posición (color) en A. Los primeros valores que se tienen
directamente son: 𝐿11 = 3 ; 𝐿21 = 0 ; 𝐿31 = 5 encontrar los demás:
𝐿11 ∙ 𝑈12 = −2 𝐿11 ∙ 𝑈13 = 0 𝐿21 ∙ 𝑈12 + 𝐿22 = 4
3 ∙ 𝑈12 = −2 3 ∙ 𝑈13 = 0 0 ∙ (−2/3) + 𝐿22 = 4
𝑈12 = −2/3 𝑈13 = 0 𝐿22 = 4

𝐿31 ∙ 𝑈12 + 𝐿32 = 2 𝐿21 ∙ 𝑈13 + 𝐿22 ∙ 𝑈23 = 1 𝐿31 ∙ 𝑈13 + 𝐿32 ∙ 𝑈23 + 𝐿33 = 1
2 0 ∙ 0 + 4 ∙ 𝑈23 = 1 16 1
5 ∙ (− ) + 𝐿32 = 2 5 ∙ 0 + ( ) + 𝐿33 = 1
3 1 3 4
10 𝑈23 = 16
𝐿32 = 2 + ( ) 4 𝐿33 = 1 −
3 12
𝐿32 = 16/3 𝐿33 = −1/3

SEM - II/2020
MÉTODOS NUMÉRICOS II MAT – 2106 “A”

Solución: Remplazar los valores obtenidos a la matriz L y U, entonces se tiene:


3 −2 0 3 0 0 1 −2/3 0
(0 4 1) = (0 4 0 ) × (0 1 1/4)
5 2 1 5 16/3 −1/3 0 0 1 Sol.
Verificación: Multiplicar las matrices L×U y para comprobar debe ser igual a A.
3 −2 0 3 0 0 1 −2/3 0
(0 4 1) = (0 4 0 ) × (0 1 1/4)
5 2 1 5 16/3 −1/3 0 0 1
3 −2 0 3 −2 0
(0 4 1) = (0 4 1)
5 2 1 5 2 1 L.q.q.d
.
2.2. DESCOMPOSICIÓN POR CHOLESKY: Matriz simétrica con diagonal dominante.

𝐀𝐟 ×𝐜 = 𝑹 × 𝑹𝑻
𝐷𝑜𝑛𝑑𝑒: 𝑅 ∶ 𝑚𝑎𝑡𝑟𝑖𝑧 𝑡𝑟𝑖𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟.
𝑅𝑇 ∶ 𝑚𝑎𝑡𝑟𝑖𝑧 𝑡𝑟𝑖𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑠𝑢𝑝𝑒𝑟𝑖𝑜𝑟
𝑦 𝑚𝑎𝑡𝑟𝑖𝑧 𝑡𝑟𝑎𝑛𝑠𝑝𝑢𝑒𝑠𝑡𝑎 𝑑𝑒 𝑅.

𝑨 = 𝑹 × 𝑹𝑻

4 3 −1 𝑟11 0 0 𝑟11 𝑟21 𝑟31


( 3 6 2 ) = (𝑟21 𝑟22 0 )×( 0 𝑟22 𝑟32 )
−1 2 4 𝑟31 𝑟32 𝑟33 0 0 𝑟33

Comprobar si es una matriz con diagonal dominante: El valor absoluto de los


valores de la diagonal principal deben ser mayores a la suma de los demás valores.
|4| ≥ |3| + |−1| → 4 ≥ 4 ¡ 𝐶𝑢𝑚𝑝𝑙𝑒!

|6| ≥ |3| + |2| → 6 ≥ 5 ¡ 𝐶𝑢𝑚𝑝𝑙𝑒!

|4| ≥ |−1| + |2| → 4 ≥ 3 ¡ 𝐶𝑢𝑚𝑝𝑙𝑒!

Multiplicar las matrices R×RT: fila por columna.

4 3 −1 𝑟11 2 𝑟11 ∙ 𝑟21 𝑟11 ∙ 𝑟31


( 3 6 2 ) = ( 𝑟11 ∙ 𝑟21 𝑟11 2 + 𝑟22 2 𝑟21 ∙ 𝑟31 + 𝑟22 ∙ 𝑟32 )
−1 2 4 𝑟11 ∙ 𝑟31 𝑟21 ∙ 𝑟31 + 𝑟22 ∙ 𝑟32 𝑟11 2 + 𝑟22 2 + 𝑟33 2

Igualar a la matriz A: Formando ecuaciones obtener los valores numéricos de los


componentes de R y RT según su posición (color) en A.

𝑟11 2 = 4 𝑟11 ∙ 𝑟21 = 3 𝑟11 ∙ 𝑟31 = −1


𝑟11 = √4 2 ∙ 𝑟21 = 3 2 ∙ 𝑟31 = −1
𝑟11 = 2 𝑟21 = 3/2 𝑟31 = −1/2

SEM - II/2020
MÉTODOS NUMÉRICOS II MAT – 2106 “A”

𝑟21 2 + 𝑟22 2 = 6 𝑟21 ∙ 𝑟31 + 𝑟22 ∙ 𝑟32 = 2 𝑟31 2 + 𝑟32 2 + 𝑟33 2 = 4


3 2 3 1 √15 1 2
( ) + 𝑟22 2 = 6 2
∙ (− ) +
2 2
∙ 𝑟32 = 2 (− ) + (1.420)2 + 𝑟33 2 = 4
2 2
9 √15 3 1
∙ 𝑟32 =2+
𝑟22 = √6 − 2 4 𝑟33 = √4 − − (1.420)2
4 4
11 2
𝑟32 = ∙
√15 4 √15 𝑟33 = 1.317
𝑟22 = 2
= 1.936
𝑟32 = 1.420

Solución: Remplazar los valores obtenidos a la matriz R y RT, entonces se tiene:

4 3 −1 2 0 0 2 3/2 −1/2
( 3 6 2 ) = ( 3/2 1.936 0 ) × (0 1.936 1.420)
−1 2 4 −1/2 1.420 1.317 0 0 1.317
Sol.
Verificación: Multiplicar las matrices R×RT y para comprobar debe ser igual a A.
4 3 −1 2 0 0 2 3/2 −1/2
( 3 6 2 ) = ( 3/2 1.936 0 ) × (0 1.936 1.420)
−1 2 4 −1/2 1.420 1.317 0 0 1.317
4 3 −1 4 3 −1
( 3 6 2) = ( 3 5.998 1.999)
−1 2 4 −1 1.999 4.001 L.q.q.d
.

SEM - II/2020

También podría gustarte