Matriz Inversa - o1_grupo 5

UNIVERSIDAD UNIVERSIDAD INDUSTRIAL DE SANTANDER FACULTAD DE INGENIERÍAS FÍSICO -QUÍMICAS ESCUELA DE INGENIERÍA DE PETRÓLEOS Métodos Numéricos en Ingeniería de Petróleos

MATRIZ INVERSA

ACOSTA CREUS MILEIDY LORENA ARISMENDI RUEDA RUBEN DARIO CARRILLO QUIJANO MARCELA ANDREA MARTÍNEZ JAIMES JULIETH PAOLA SILVA ALONSO KATHERINE LORENA

PH.D. EDUARDO CARRILLO Docente

UNIVERSIDAD INDUSTRIAL DE SANTANDER ESCUELA DE INGENIERÍA DE PETRÓLEOS MÉTODOS NÚMERICOS EN INGENIERÍA DE PETRÓLEOS BUCARAMANGA 2010


MATRIZ INVERSA

FUNDAMENTACIÓN BÁSICA 1. Definición de matriz Una matriz A de m x n es un arreglo rectangular de mnnúmeros dispuestos en m renglones y n columnas.

« a11 ¬ ¬ a21 ¬: A ! ¬ ¬ ai1 ¬ ¬: ¬ am1 -

a12

...

a1 j

a22

...

a2 j

:

ai 2 :

am 2

...

...

a1n »

¼

a2 n ¼

¼ ¼ aij ... ain ¼ ¼ : : ¼ amj ... amn ¼½ :

...

...

:

Una matriz con m renglones y n columnas se llama una matriz de m x n. El símbolo m x n se lee ³m por n´.

( a11

El vector renglón a12 ... a1 j ... a1n )se llama

El vector columna columna j

se llama

renglón o fila i La componente o element o ij de A, denotado por aij , es el número que aparece en el fila i y la columna j de A. En ocasiones se escribirá la matriz A como A= ( a ) ij . Por lo general las matrices se denotarán con letras mayúsculas. Por ejemplo, los coeficientes de las variables X 1, X2, X3 en el siguiente sistema: 2X1 + 4X2 + 6X3 =18 4X1 +5X2 + 6X3 =24 3X1 + X2 -2X3 =4

(1)

Se puede escribir como elementos de una matriz A, llamada matriz de del sistema:

co

efi c ientes c ientes


         2. Tipos de Matrices

2.1.

Matriz Cuadrada

Si A es una matriz m x n donde el número de filas fil as es igual al número número de columnas; es decir, m=n. Entonces la matriz Aes de orden n o llamadamatriz cuadrada. Una matriz cuadrada A de orden n es singular si su determinante es nulo. En tal caso se dice que dicha matriz no tiene inversa. Ejemplo: Sean las matrices

«1 2 3» ¬ ¼ A ! 4 0 5 ¬ ¼ ¬- 3 1 2 ¼½

B!

« 2 3» ¬ 1 5 ¼ ½

Entonces, A y B son matrices cuadradas de orden 3 y 2 respectivamente.

2.2.

Matriz Transpuesta

Sea A = (a ij) una matriz de m x n. Entonces la transpuesta de A, que se escribe A t, es una matriz de n x m obtenida al intercambiar los renglones por las columnas de A. De manera breve, se puede escribir A t = (a ji). En otras palabras,

Si

«a ¬a A ! ¬ ¬: ¬ - am

11

a12 ... 12

21

a22 ...

1

: am 2 ...

» « a11 ¼ ¬ a n ¼ , entonces At ! ¬ a12 ¼ ¬: : ¼ ¬ amn ½ - a1n a1n

a21

...

2

a22

...

:

a2 n ...

am1 »

¼ ¼ ¼ : ¼ amn ½ am 2

Simplemente se coloca el renglón o fila i de A como la columna i de At y la columna j de A como el renglón o fila j de At .

Ejempl o:


2.3.

                  

Matriz Simétrica

La matriz (cuadrada ) Bde n x n se llama simétrica si B t = B Es decir, las columnas de A son también los renglones o filas de B . Ejempl o:

« 1 ¬2 B ! ¬ ¬4 ¬ -6

2.4.

2

4

7

3

3

8

5

0

» ¼ 5 ¼ 0 ¼ ¼ 4 ½ 6



B

t

« 1 ¬ ¬ 2 !¬ ¬ 4 ¬ ¬ 6 -

2

4

7

3

3

8

5

0

» ¼ ¼ 5 ¼ ¼ 0 ¼  4¼½ 6

Determinante

Los determinantes para una matriz de orden 2 se resuelven de la siguiente manera:

              

 

Determinantes de orden 3 son calculados con la regla de Sarros como sigue:

Dada una matriz cuadrada A de tamaño (n ) se define su determinante como la suma del producto de los elementos de una línea cualquiera de la matriz (fila o columna ) elegida, por sus correspondientes adjuntos.

 



det(A) = a i1 Ai1 + ai2 Ai2 +«+ +«+ a in Ain det(A) = a 1j A1j + a2j A2j +«+ +«+ anj Anj




Para una matriz cuadrada de orden n, A = (a ) ij se llama adjunto del elemento aij, y lo representamos por Aij al producto (- 1 )i+j aij , es decir: Aij = (- 1)



i+j

aij

La matriz cuyos elementos son los adjuntos de los elementos de una matriz cuadrada A se llama matriz adjunta de A y se denota por Adj(A ).

                                                                                                                                   

Los adjuntos de la matriz:

La matriz adjunta de A es



2.5.

Matriz triangular superior o inferior

Una matriz triangular es un tipo especial de matriz cuadrada cuyos elementos por encima o por debajo de su diagonal principal son cero. Debido a que los sistemas de ecuaciones lineales con matrices triangulares son mucho más fáciles de resolver, las matrices triangulares son utilizadas en análisis numérico para resolver sistemas de ecuaciones lineales, calcular inversas y determinantes de matrices.


      

Triangular Superior

2.6.

Triangular Inferior

Matriz aumentada

La matriz ampliada o aumentada se obtiene al combinar dos matrices tal y como se muestra a continuación. Sean las matrices A y B, donde

                                 

Entonces la matriz aumentada ( A | B) se representa de la siguiente manera:



2.7.

Matriz Bandeada

Matriz cuyo ancho de diagonal depende del sistema de solución de ecuaciones con que se está trabajando; puede ser bidiagonal, tridiagonal, etc. Su forma se observa en la siguiente imagen:

En otras palabras, es una matriz que tiene todos sus elementos cero, e xcepto los de la diagonal principal.


Ejemplo:

2.8.

Multiplicación de Matrices

Sea A= (a ij) una matriz m x n, y sea B = (b ij) una matriz n x p. Entonces el producto c ), de A y B es una matriz m x p, C= ( c en d onde: ij =(renglón i de A) A)·(columna j de B)

c ij

Es decir, el elemento ij de AB es el producto punto del renglón i de A y la columna j de B. Si esto se e xtiende, se obtiene: = ai1 b j1 + ai2 b2j +«+ ainbnj

c ij

Si el númer o de col umnas de A es ig ual al númer o de rengl ones de B, ent onc es es se di c ión. ce que A y B s on compatibles baj o la multipli c c ac ión.

Ejemplo: Dadas las matrices Q y W calcular la multiplicación entre ellas.


MATRIZ INVERSA 1. Definición Sean A y B dos matrices de n x n. Suponga que: A·B=B·A=I -1 Entonces B se llama la inversa de A y se denota por A A . De esta forma se tiene:

-1

-1

A A = A A A I n 



=

Si A tiene inversa, entonces se dice que A es invertible.

2. Existencia de la Inversa Una matriz cuadrada que no es invertible se llama singular y una matriz invertible se llama no singular o regular. Es decir, si el determinante de una matriz A es de la siguiente forma, se tiene : |A|  0 M atriz atriz no sing ular o reg ular |A| = 0 M atriz atriz sing ular En definitiva, sólo e xiste matriz inversa de una matriz cuadrada si ésta es regular. Otros criterios para reafirman cuando una matriz es invertible son: i) A es invertible si y sólo si A es equivalente por renglones a la matriz identidad I n; esto es, si la forma escalonada reducida por renglones de A es In. ii) A es invertible si y sólo si el sistema A x= b tiene solución única para cada n-vector b. iii) Si A es invertible, entonces la solución única de A x= b está dada por x=A-1b.


iv)

A es invertible si y sólo si su forma escalonada reducida por renglones tiene n pivotes.

Por ejemplo, un caso de matriz de 3 x3 que no es invertible se presenta a continuación:

      × Sea A =

. Calcule A-1 si existe.

                                                              Hasta aquí se puede llegar. La matriz A no puede reducirse a la matriz identidad, por lo que se concluye que A no es invertible.

3. Propiedades La matriz inversa de una matriz cuadrada, si e xiste, es única. La inversa del producto de dos matrices es el producto de las inversas cambiando el orden: (B · A)-1 = A-1 · B-1

Si la matriz es invertible, también lo es su transpuesta, y el inverso de su transpuesta es la transpuesta de su inversa, es decir: (AT)-1 = (A-1)T


Una matriz es invertible si y sólo si el determinante de A es distinto de cero. Además la inversa satisface la igualdad:

     

Donde |A| es el determinante de A y adj(A) es la matriz de adjuntos de A.

4. Métodos Método s para hallar la matriz inversa Por definición o Método directo Por el método de Gauss Por determinantes o matriz adjunta Por factorización LU 4.1.

Por definición o Método Directo

Si tenemos un número real, por ejemplo el 2, podemos interesarnos en buscar el inverso del 2 para el producto, es decir un número real x tal que 2· x = 1. Evidentemente, en el caso de los números reales es bien fácil despejar x para obtener que x =1/2, es decir, el inverso de un número real es otro número que multiplicado por él da el elemento neutro, el 1. Todo número real, salvo el 0, tiene inverso. Trasladando esto a las matrices, nos podemos plantear si dada una matriz cuadrada A de orden n cualquiera, e xiste su inversa X para el producto de matrices, tal que: A X = l n 

Es decir, el producto de A por su inversa produce el elemento neutro matricial, la matriz identidad I n.

 

Sin embargo, hay algunas diferencias con respecto al caso de los números reales es notoria cuando no podemos ³despejar´ la matriz X del modo hemos definido la división de matrices. Ejemplo:



, porque no


                                                    

Determinar la inversa de la matriz

El método directo consiste en determinar A-1 planteando un sistema de ecuaciones, es decir, se busca otra matriz de igual tamaño (orden 2 ) tal que -1

-1

A· A =I 2 2 y A · A=I 2. 2. Dicha matriz es -1

A · A=I 2 2



como:

Al resolver este sistema de 4 ecuaciones con 4 incógnitas, se obtiene que:

4.2.

Método de Gauss-Jordan Gauss-Jordan

Consiste en hacer transformaciones elementales en las filas de la matriz para llegar a obtener la matriz identidad. Realizando estas mismas transformaciones con la matriz identidad llegamos a la matriz A -1. Se llama transformación elemental en una matriz a: T1) Multiplicar o dividir una fila por un número real distinto de cero. T2) Sumar o restar a una fila otra multiplicada por un número real no nulo. T3) Intercambiar el lugar de dos filas entre sí.

 

Veamos cómo se realiza el método de Gauss-Jordan, realizándolo a la vez con la matriz

i)



.

Consideramos la matriz formada por A y la matriz identidad correspondiente. En nuestro caso:


     



ii)

Se hace la matriz triangular superior (es decir, hacemos ceros por debajo de la diagonal principal ) usando transformaciones elementales en filas. La mejor forma de realizar esto es hacer cero los elementos de la diagonal en la primera columna usando la fila 1. Luego, hacer cero los elementos por debajo de la diagonal en la segunda columna usando la fila 2 y así sucesivamente. En nuestro caso, basta sumar la fila 2 con la fila 1, y se obtiene:

         



iii)



Una vez hecha la matriz triangular superior, se hace la matriz triangular inferior, haciendo ceros a los elementos por encima de la diagonal. El proceso es parecido al anterior: Hacer ceros los elementos por encima de la diagonal en la última columna usando la última fila. Luego, hacer cero los elementos por encima de la diagonal en la penúltima columna usando la penúltima fila, y así sucesivamente. En nuestro caso:

 iv)

        

Ya tenemos una matriz diagonal. Lo único que falta es dividir cada fila entre el número adecuado para obtener unos en la diagonal principal, es decir, para obtener la matriz identidad en la parte izquierda:

 

 v)

        

  

Una vez se tiene la matriz identidad en la parte de la izquierda, la parte derecha es la matriz inversa, es decir, llegamos a:



                  


Si al realizar el método de Gauss-Jordan en algún momento alguna fila es de ceros, la matriz no tiene inversa. Cuanto mayor sea el orden de la matriz, mejor es este método frente al directo.

                                                                                                                                      

Veamos otro ejemplo:

Calcular la inversa de la matriz

por el método de Gauss-Jordan.

Siguiendo los pasos anteriores:





También se puede e xpresar sacando factor común:

Es la inversa de B.


            

Si calculamos por este método la inversa de resulta:









Como aparece una fila de ceros, la matriz A no tiene inversa.

4.3.

Determinantes Determinantes o matriz de adjuntos

Dada una matriz A, determinamos la matriz de adjuntos de su traspuesta. Si multiplicamos esa matriz por 1/|A| se obtiene la matriz inversa de A.

                               

Para calcular la inversa, primero calculamos el determinante:

Después calculamos cada uno de los adjuntos :

Y por tanto


        

    

Comprobación

4.5.

Descomposición Descompos ición LU para matrices inversas

La descomposición LU involucra solo operaciones sobre los coeficientes de la matriz [A], proporcionando un medio eficiente para calcular la matriz inversa o resolver sistemas de álgebra lineal.

                                      

L

Debido a que [A] = [L][U], al encontrar [L] y [U] a partir de [A] no se altera en nada la ecuación y se tiene lo siguiente:

PASOS

                               

1. Obtener la matriz triangular inferior L y la matriz triangular superior U. 2. Resolver Ly = b (para encontrar y ). 3. El resultado del paso anterior se guarda en una matriz nueva de nombre ³y´. 4. Realizar Ux = y (para encontrar x). 5. El resultado del paso anterior se almacena en una matriz nueva llamada ³ x´, la cual brinda los valores correspondientes a las incógnitas de la ecuación


EJEMPLO Obtener la inversa de la matriz [A]

Paso 1. Aplicar Gauss para obtener [U] y [L]


Paso 2. Obtener el vector de incógnitas {D} del siguiente sistema [L]{D}={B1}

De la figura 1 se resuelve el siguiente sistema [U]{ x1}={D}, aplicando sustitución hacia atrás:

Paso 3. Obtener el vector de incógnitas {D} del siguiente sistema [L]{D}={B2}, aplicando sustitución hacia adelante:



Paso 4. Obtener el vector de incógnitas {D} del siguiente sistema [L]{D}={B3}, aplicando sustitución hacia adelante:


5. Uso de la inversa de una matriz para resolver un sistema de ecuaciones Resuelva el sistema 2x + 4y + 3z = 6 y±z=-4 3x + 5y + 7z =3 olución. S olución.

Este sistema se puede escribir como A x = b, donde

                 


      

     

Así, la solución única dada por:



                    


BIBLIOGRAFÍA

GROSSMAN, Stanley l. Algebra lineal. 5ta Edición. Editorial Mc Graw Hill. 2008 http://www.ematematicas.net/matrices.php?a=6&tipo=6 http://www.terra.es/personal2/jpb00000/tmatrizinversa.htm http://es.wikipedia.org/wiki/Matriz_invertible http://personal.redestb.es/ztt/tem/t6_matrices.htm http://portales.educared.net/wikiEducared/inde x.php?title=Desarrollo_de_un _determinante http://sauce.pntic.mec.es/~jpeo0002/Archivos/PDF/T06.pdf http://www.investigacion-operaciones.com/Calculo%20matricial.htm http://books.google.com.co/books?id=zsQM4QfNVwoC&pg=PA67&lpg=PA6 7&dq=matriz+bandeada&source=bl&ots=tnAHbNsrMX&sig=2k62ChbPb1F0 M1vAiOg_jXP67PA&hl=es&ei=qWAgTL_CMoO8lQffuLSdAQ&sa=X&oi=boo k_result&ct=result&resnum=8&ved=0CDoQ6AEwBw#v=onepage&q=matriz %20bandeada&f=false

Matriz Inversa - o1_grupo 5

Recommend Documents