Algebra Lineal
Algebra Lineal
MLGEBMffi
LffiMEffiL
m flo$
Segunda Ediciór'¡
MatrEces
Determinantes
l
Sistemas de Ecuaciones
Factorizacióm
trV{as'de 'B O0 Eiercici«¡s
FLesu¡eItos
§c¡§ucionanlo de Exánreiles
2$$6-20{ { itrcluyendo
ct¡rsos interlsivos
rl I
IlIl I
I
1,,) I i r.rll I
lrll I
-'.,- l i lli I I
["]'ll
\ /t
t=
!
I
[.iB]= I l lLr I
I
l^l- lr¡+2\ I
i
I
[ , ]l I llt I
,l I
*z¡l
¡2,,
[ ,,-t )l
/i
, CONTENIDO DEL TEXTO
Mf,AT H[GtrS
INTRODUGCION
por su facilidad de
extensos, mayormente se utiriza matrices
cuando ros sistemas de ecuaciones rineares son de sistemas de
datos y ." ur.l no solo en la resoluciÓn
manejo. Las matrices son ordenamientos de
nurnéri"o, en la resoluciÓn de sistemas de ecuaciones
ecuaciones (lineales), sino ademá= "1 "i cánuto en geometrfa'
tas matrices tamoi¿n aparecen de forma natural
diferenciares y de derivadas parciares.. Además
i.ii"", eóonom ía, informática, física, etc'
".üá que myghas
de ecuaciones rineares. sin embargo, puesto
Er árgebra matriciar puede ser aplicada a sistema y otras pueden ser convertidas a
mediante ecuaciones rineares
reraciones económicas pueden ser aproxrmadas
parte evitada'
réiác¡ones lineales, esta limitación puede ser en
DEFINICION I. MATRIZ
o retras u objetos (*amados erementos) dispuestos en
nx
una matriz es un anegro rectangurar de números
(columnas) de la forma:
lÍneas horizontales (nlai) y rl líneas verticales
(Notación de LEIBNITZ)
Í\ »txn -
a¡l ü¡2 a¡3 oü d i¡t
o,rrn
lml Qm2 am3 o,nj
tnxn
ElconjuntodeelementosoComponentes(objetos)deuna.m,atrizpuedenSerre.presentadosentrecorcheteso de la
lás filas de la matriz y los verticales son sus columnas
entre paréntesis. Los tárminos horizontales son
misma.
alt alZ aÉ aln <_
a2z 023 a2n <_
<- Filas de A
t:
<-
antl 0nt2 am3 ottttt +-
t t t tt
Columnas de A
Observaeiones.
; Los objetos a,7 pueden ser números reales, números complejos o cualquier objeto no
numérico como
por ejemplo las fichás de un ajedrez o los apellidos de personas cuando estas se codiflcan en orden
alfabético.
. Una matriz es un aneglo y como tal no un valor numérico
Una matríz de rz fifas y ll columnas se dice que tiene dimensión o que es de orden mxn, y al conjunto de
todas las matrices de orden mxn lo denotaremos por:
T-1
A**n = [ay _] ittotación de KRONECKER)
Con 1<i <m ,l l< i <tt
Esto es:
ORDEN DE
I\4ATRlz
f\lATRrz -___LA
,l f- ELEivIENTOS
vl
i
'l -V
tt
¡V tt
1\
r-\ ¡v7 v7 I LÍ, I
¡- IL^
¡l i
L I',.tl-
)
coLUt,tNA
--',
CCLUf\'1NA- i FIL¡,
El conjunto de matrices de orden mxn, con coeflcientes en IK, donde IK Puede o complejo se
Ejemplos matrices:
L,o oj ;
l;", -rl
2abc2 L4
a
rt9.. Jtll=lo UBERHUAGA C-
IYG, §C. INC.{AFAEL VALET{CUI GOYZUETA !
ALGEBRA LINEAL Y TEORIA ¡ÚATR¡GIAL
a)Escribimoslascomponentesdecadamatrizsegúnladefinicióndada:
? +2=3 lz 3l
4r
ázr§
= 1-!
+ u=l t .l
1l
Llo 1
3 - lndicar el
a).
¿f - En Ia siguiente matriz
{-
L
TTPOS DE MATRICES
grupo
Las matrices pueden clasifiearse en dos grandes
I'
Lanú
dm2 ,,rrr)
Aun:lotj-] = A€iR''^"
ft.r,
Se denomina, matriz flla a
FILA que consta de una única
fu | atz ái3 'rry
AeR"''l
=
ffi,?'l
Se denof
COLUMNA ntxl
lri
*r
EÉ; :l
{
ffi L
untl
l<i<tt
e IR'*' con
1l i<rt
ffi A, -t
Ir,
I
l:
azt
r ür2
azz
a1¡t
ü2tt
att00 0
K o o ol
Es una matriz diagonal (Y en a,, -K i= i
I 0 K 0 0l
bonsecuencia, una matriz =]A;: -0 i+j' o o rl
u
Ar, An=
ESCALAR cuadrada) con todos los [, I
lo¡j
por el símbolo In .
00 o 1_l
tNG. JULIO UBER}IUAGA C.
MG¡ SC. ING. RAFAEL VALET{CIA GO}T¿UETA
ALGEBRA LINEAL Y TEORTA MATRrcIAL
:
Es aquella que todos su§
ESENTAGIO
I-9
0
:ltl
elementosson0Ys'e hrrr,, =o =lo¡i = o] v i,7 A=l:
representa Por 0. [oo ool
ffiadaquetienenulostodosloselementosqueestánaunmismo
lado de la diagonal principal. Las matrices
lorr aU a-1
Si los elementos que es
v i<j 0 ctzz
"rrl
TRIANGULAR
SUPERIOR
por debajo de la dit
i> j ^'-l-l o 0
.l
: lo 0 aru )
0
[,,,
Si los e
I nzt azz
TRIANGUI-AR -t t.
INFERIOR
i.]
I
atn2
Lantl
: at" 4l Qlrt
Q3,,
n,1,,,,
0
Se llaman a
equivalente
ESCALONADA rectangulares
triangular
ü,;
a2n
Qi,,
o,t,,l
TRANSPUESTA colurnna
de A
¿e AT, la segunda fila
es la segunda columna de 'r,. ' *,,r 1 l*', 'Z'' '""")
AT y asÍ sucesivamente. De la
deflnición se deduce que si A es
INVOLUTIVA
*
.f: l,*
Y';$
{
tB",:!
ffi''
L cumoleAí.FAi=J
.éi
F,"f,si, q
T
- I&-
ORTOGONAL
1 t..riE
i.i
ir
ki¡
azl
o2z
az3
Qnl
ant2
a»t3
10000
01000
00100
00010
ol¡ a2r, Qnut .00001
**..,.oLftrelB-An=|ou)siestaseigualaalatranspuestadesu
HERMITANIA
conjugada, esto es: A:|tr)
t'"-1 s¡ esta se iguala al negativo de la
Una matriz cuadrada y compleja A, = [ ,; -t
HEMIHERMETICA
i- +-^ñ-^,,^cta ocln o" A= - AT
conlu|faua utr §u .lrqrrr¿ ''-""Hry
(o sea con el mismo numero
Dadas ra matriz o,,,*n --lr,r7 r la matriz 8,,,*p= tur]
AMPLIADA de filas), se denomina matriz ampliada a la matriz
de la forma [eln] =lo,¡lb,¡) o"
GOYZUETA ING, JULIO UBERHUAGA C'
MG. SC. ITG. RAFAEL VALEHGiA
ALGEBRA LINEAL Y TEORTA I'ATRICIAL
\,
tlAl) es distinto de
es regular, si su determinante2 es
:
adjunta3 cuando sus elementos son
los
ADJUNTA
ES ESPECIALES
PROPIEDADES Y
CUADRADA
onal princiPal.
,.e=Ir,'A'7n=A
2. I+I+I+
IDENTIDAD
3. 1'?=I, V ¡
4. Sea A una mal A'I:I'A:A
Sl k es el menor número !ue satisface la condiciÓn Ak:I se dice que
A es una mairiz de Periodo Ak+l:A , Ak+2:A2 . Ak*3:43 y así
sucesivamente.
IDEMPOTENTE A2:Aentoncet@Si
Si k es un numero entero Positivo qu
NILP.TENTE
i e*-t:g, Ak*2:0, Ak*3:0 yasísucesivamente
1. A"=A si l¡es
a su existencia'
, a""a*,. ," d.*-*nte lo desarrollaremosadelante, en el capitulo correspondiente, aqui solo.hacemos refere.ncia
3 adjunta O".arrotl"do mas aquí solo lo mencionamos para definir el tipo matriz'
El concepto de solo lo mencionamos para identificar el tipo
=ur¿
. Er concepto de matriz i.r"*. *rv amprio y será deiariolrado en detalle mas adllante, aqui
=
de matriz.
con mas profundidad, aquí simplemente la nombramos como un tipo
de
5
El concepto formal de submafiz se desanollara mas adelante
matriz que se presenta en el est'¡dio de matrices'
ING. JULIO UBERHUAGA C.
ALGEBRA LINEAL Y TEORIA MATRICIAL MG. SC. ING. RAFAEL VALENCIA GOT¿gETA
t.
{o')': o
2. (*)': k'Ar
TRANSPUESTA 3. (e+e)r: 6r+8r
4. (a-s)': BrAr
5. (ro)t:ln
1. Si A essimétrica, ent k.A también -s simétrica V k escalary k* 0
: k.A también esantisimetrica V k y k*0
É,4+AT esd*Efüa@$qlétrioa ¡
F&{dl;i¡; rxjr;: - t
riamente normal.
CONJUGABA
isxl:H esherméticaVkeiR
imbién es hemehermetica V k e IR
, e*eT es hermética
HERMETICA ces: A-AT es hemihermetica
HEMIHERMETICA elementos son números complejos en general
puros) se puede expresar como la adiciÓn de
hemihermetica, así:
", c:l(e-n')
1. Si A eso l, su determinante es distinto de mas menos uno, esto es
lal+ tt
2. Si A y Bson matrices ortogonales, entonces A'B estambién ortogonal
ORTOGONAL 3. Si A es una matriz ortogonal, entonces AT = A-1'
4. Si A es una matriz ortogonal, entonces AT también es ortogonal
s. Si A es una matriz ortogonal, entonces A-ltambién es ortogonal'
1. A'A-l =I
INVERSA z. (a-r)-t = e
3. (e e)-':s-1.4-t
-at-
ING. JULIO UBERHUAGA C.
MG. SC. IiIG. RAFAEL VALENCIA GOYZUETA
ALGEBRA LINEAL Y TEORTA MATRICIAL
4. (o')-':("-')'
5. AR=BA:I = B:A-i
6. matriz es única
La inversa de una
1. A:B
2- A:B
: sumade Ay B a otra
Adición de
1. clausuraVAB € l]("rxl¡ j
2. Conmutativa :A,B = ¡6//rxl?
j
3. AsociaiivaVAB,C aTKtttxtl => ,{
4. Elemento neulro aditivo VA e lKtnxtl 'r,*n é [*$:S*,{:A
opuesto (inverso) aditivo
5. Elemento
(e,e)-+A+B
-9-
ING. JULIO UBERHUAGA C.
ALGEBRA LINEAL Y TEORIA MATRICIAL MG. SC, II.¡G. RAFAEL VALENCIA GOYZUETA
VA,B e IK
A
Elemento
5.
b.
ffi
"':g'ffi
eo,. t.nto (Ñ'i tL* nj;I de los números reales
Multiplicación de matrices' Para poder atrices, Ia primera debe tener el mismo número de
columnas que filas la segunda. La matriz re§ producto quedará con el mismo número de filas de la
primera y con el mismo número de columnas d
-10-
ING. JULIO UBERHUAGA C.
MG. §C. IT{G. RAFAEL VALEIGTA GOT¿UETA
iLGEBRA LINEAL Y TEORIA MATRICIAL
4r
4t
B: 4r matriz columna de orden P x
ary)matnzfiladeorden 1xp y
1
sean A:[a11 an aB
ui^
y B (en ese orden)' es Ia matriz lvl Oe orden 1x 1 :
Propiedades de la
1. Asociativa VA,
+B.C
2. DistributivaVA"B,
4. No cancelativa. Si A'B:A'
para dos matrices de orden
5. No conmutativa- En general es
A,,,*p=[nr] y B p*r=[ar] tt
6. Si A.B:0 no necesariamente imPli
7. Si A:B = AC=BC v CA:CB
8. Si AB:BA entonces A Y B son matrices les'o conmutables
A'A:A2
A.A2:A2 .A:A3
9. Para una matnz cuadrada A : A'A3:A3 ' A:A1
:
A.Ar-l-An-t. \:§n
10. Para una matriz cuadrada A Ap 'Aq = A't 'Ap Yp,q e Z+
11. En general (ete)2 :A2 + B2 +2'(a'g)
A'B*B'A va que
- 11-
MG. SC. !i¡G. R.AFAEL VALENCIA GOYZUETA ING' JULIO UBERHUAGA C'
ALGEBRA LINEAL Y TEORTA MATRICIAL
tj,
1. rr (e+n) = r,'(A) + rr (e
z. *(k-A):&./,'(A)v¿
3. rr(a.B) = ¡r(e'
5. rr (e) = 0,,
6. rr(I)=,r '
De forma Pract
b=)'
d =v'
.,:,=,7
Escalar por v A e ¡7L2x2
matriz p escalar
7
Los siguientes ejemplos se desarrollaran como máximo para matrices
cuadradas de orden 3
- tL -
ING. JULIO UBERHUAGA C.
ALGEBRA LINEAL Y TEORIA MATRICIAL MG. SC. ING. RAFAEL VALENCIA GOT¿UETA
Unitaria
^'[itl] ^'[l :il
A3=[l -4 tl
=Ll
4
[eool
0 o.=lo 6 ol
[ooe]
l-r o ol
Triangular suPerior inferior =l-, r ol
[-o -1 ,]
[r o ol
a.=lo , ,l
[ooo]
:
Ir
o r-l
=10 24ll
[, 4 3]
[o t z1
ANTtstMETRIcA i A: =
.,=l-, o ,l
l-, -3 ol
[r ll -sl lz
a
-+l
+i
IDEMPOTENTE o,=17
- tr 'l s I e', =l-1 1
J
It
Lz z) -+l Ir -L
'.)
-:]
Ii -j -+] Il -j -tl
N¡LPOTENTE o.- = [' -rl o. = [-',
L4 -z)
o,=l-,, ol
t -3 -+_]
=l-'' ol
[_o I L1 -3 4)
1 rl [-: -6 z) l+ 3 :l
INVOLUTIV,A -lo ol
o,=[-' rt o,=[l -1
:] ^'=L : :;] =| -r 0
*4
-rl
-rl.l
lo 0 [-o
T
II i-3i slt-t l-¡ 1+ 3i s.l
CONJUGADA A¡ =l 2 i-1
I
_.1
I
I
[-o n ,, 0I l4-3¡ -3 0
ING. JUUO UBERHUAGA C.
ALGEBRA LINEAL Y TEORIA MATRICIAL MG. SC. ING. RAFAEL VALENCIA GOY¿UETA
[+ -i t+zif 4 -i
AI=i i -3 +-til
-f
3+2il
HERMATICA er=l i 4 +-til
ls-z¡ 4+7i 6 _l li-zi 4+7i 6 I
iz , -ol lz -3 4-',
[+ -i 2)
6 3l lqs ol
-3 6l Lo 4sl
ft -ll "j''
i11 -3 i
A,=l
' 13 1lI #.*l ie a_l
I
PERMUTABLE f¡
IT
Br =l
l-
I
,ul
10 5 I
12j
.l
4 r) 1-2 [+ -1 -t)-l
ó _L a 5-5 l=l -z -11 o
I L, -1 6l
I
,\
7-ü 0-1
rl
I
:
lc.t., lar A.B
l
primero el numero de filas x columnas- Sitenemos que una matriz
es 3
ar si
GOYZUETA ING' JUL¡O UBERHUAGA C'
ALGEBRA LTNEAL Y TIT RIA MATR]CIAL MG. SC. !I.IG' RAFAEL VALENCIA
-----{- -.
E[ lamaño de la
respuesta es
3x2
Si lm números cenlrales
9 son iguale§ entOnces se
puede multiplicar Y el
tamaño de la resPuesta
son los nÚmeros de los
Debe ser, extremos 3 x 2
Ahora la multiplicaclón:
Ir +0 0I I 4+31 l+ 71
A.B =lo 0+r0+61=1,, ,ul
Io + 2+o L, 2)I
1 Demostrarque A.(B.C)=(A.
Sea las matrices AeIK e IK""
Estas definidas por u=lo,f c = [tr,]
si B.C =1d,,]entonces, d
¡r = Z(arr X
A'= i
fit
Si A. B =fn,r] entonces, e¡r =L(r¡)(r'ir)
J=l
En consecuencia si A-(B.c)=[.lr],- (a.e)'c=[&,]
-15_
GOYZUETA ING' JULIO UBERHI"'AGA C' '
ALGEBRA LINEAL Y TEORIA MATRICIAL MG. SC. ING. RAFAEL IJALENCIA
('")
Por lo tanto
Luego
Entonces
Q
portanto(,t.¿) = ¿(e'n)
Demostrar que li'(A'
Sean las matrices con I U<t cron A Brr*,r=
['r]
Asítenemos: =C,,*,., = ,i =Z(o¡¡)(a¡o )= c¡¡ =Z@,r)(br,)
Á-=l k=l
= Il¡=1\Á.=t
U bnn¡t'|
)
t't
= I(a*o)
i=l
= lr (D) Por tanto tr ( A.B) = tr(n' A)
-16-
MFAEL VALENCIA GOYZUETA ¡NG. JULIO UBERHUAGA C.
ALGEBRA LINEAL Y TEORIA MATF }IAL MG. SC. II*G.
I'l
Multiplicando
k.n'(A)
n. A^,,)=1(A- AT)
Ahora sl escribimos a-
"2=l(a+ar)
Luego A. es simétrica y Ao es antisimetrica, en consecuencia, hemos expresado la matriz A como [a
suma de una matriz simétrica v otra antisimetrica.
A'B
=lo,r7 " p
Si hacemos = C es una matriz de orden t?lx
-17 -
j
ft
)(a'e) = t
Ie).4-t
r).¿-1
A.A-t
I
Para b) e-l .(e-1 e).e
n-i .(i).e
B.B-l
Por lo tanto de a) Y b) se concluYe
I (a .e)-1 : g-l ' 4-t
Se debe demostrar
Desarrollando
r.\,:'
Para b) Ay B-1 Multiplicamos Por B-l tanto porizquierda como por derecha a la matriz:
B-l .(A.n)-n-t = B-l .(e.a)'e-l
(e-' axe .e-t ) = (e-t 'e)(a e-t )
(e-' .e)(r) = (r)(a'e-t)
B-1 -A = A.B-l PortantoAv B-1 son Permutables
-19-
GOYZIETA ING' JUUO UBERHUAGA C'
ALGEBRA LINEAL Y TEOEIA MATRICIAL MG' SC. !i¡G. RAFAEL VALENCIA
14SeanAyBmatricessimétric,,,@tricasisolosiA.B=B.A
La demoslración üene dos Partes:
primero: ffi*tilXJi*:H,*rcs simétrica' rueso B es A= Ar n. B = Br
La segunda hipótesis, puesto que A'Bes (a.g)T = A.B,
simétril:a es:
Por otm lado (a.g)r 'Ar = B' A
=BT
-20 -
ING. JULIO UBERHUAGA C.
v ALGEBRA LINEAL YTrcRtA }'IATRIC¡AL MG. SC.II.¡G. RAFAEL VAL=NCIA GOYZUETA
Entonces -.'
Pero 42 = 0, eato i
Ahora, operando
BT 'A.B es idempotente
r2
Se debe demostrar que B) = B' .A'B
Como Aesidempotentey B ortogonal A B.BT:I
Desarrollanda¡.2/tr
(et.A.-eJ = (B' .A B)
)a.s ¿
=BT
=BT
=BT
Sean A matriz simétrica entonces dem métrica cualquiera que sea la
20
matriz B de orden lt t-¡
Se debe demostrar que
'}í
r?'
Como A es
Sean4' demostrar que
21
(n+»)(
Por la
#É
pro4Ec
distributiva de la suma
(A i+D) D
"i¡'J $ry
_D2 Portanto (a+P)(a-D) = A2 -Dz
cy p son escalares, demostrar que A
22 ySi A y B:aA + /l son matrices del , donde
B son
conmutables
Se debe demostrar que A'B = B'A
Entonces: A.B = A.(aA+ §I)
= A.(aA)+e.(pt)
= a(x.A)+É(A.I)
= (aA+ pI)-A Portanto A y B son conmutables
EQUIVALENCIAS
DEFtNtctóN 4. TRANSFoRMACIONES U OPERACIONES ELEMENTALES.
Transf
I de Este efecto se
Sustituir la fita la
lntercambiar la colurnna i por la columna fr de la matriz A ( ci (-- c1. ). Este efecto se Produce al
multiplicar, por la derecha, la matnz A por la matnz elemenial E7A- , siendo esta el resultado de intercambiar
efecto se produce al multiplicar, por la derecha, la matriz A por Ia matriz E.¡(a) d*A, siendo esta el
resultado de multiplicat pot d la columna .¡ de la matnz I,,
'23'
GOY¿UETA ING' JULIO UBERHUAGA C'
ALGEBRA LTNEAL y TEoRIA IIATRICIAL MG. SC. t!'lc. RAFAEL VALEI,IC¡A
3. Sustituir una columna de la matriz por lá ¿¡"ta cálumna con un múltiplo escalar
de otra columna en Ia matr¡z.
"um;?.
A su columna * multiplicada por a*O(ac¡r+c¡ c; )' Este
Sumar a la columna ,r de ta matriz
A por la matriz ng{a), -)
siendo esta la resultante de
efecto se produce al multiplicar, por la derecha, la matriz
Sumar a la columna j de la matriz I, su columna * mulüplicada por d *0
la b cl
Ejemplo. Sea la matriz A=l p q r i
operaciones elementales:
OPERACION
F.lfob,1
)q'l:l * , ,l
,,,1 L, s ,) -
u,'):¡:,:,;l
iT:jf;
o o riL, .! r) L'. .t' z)
7 [r + ol;" t cj la-zr t-z,t "-»1
.A=loiollrs,i-p q ,l
Lroil['.,,,]l_* )' z)
COLUMNA
cl7+ c: [rool
i:lo r oi o *,,=[í
lo o I-l 1 ;]ti : i]:[; ?"')
[r o ol '?f
l, , .l[t o ol la b r,l
o 2c)
ol
r:loiol +K3(2)=lo,ol A.K:{il=ln t.li o t nl:l , ,
lo o r_l [o o 2J l, .r, ,_lL, o 2) l, v 2z)
t-r o ol h o 3l fo t clft o ¡l [a b ¡"l.c)
1:l orol =Kr:(¡)=lorol A.Kr¡(:)=l r rilo , ol:l p q tu*,'l
Lorrl Lootl [, .v ,l[o o t-] [-. )' 3.r+zl
DEFINICIÓN 6. MATRIZ EQUIVALENTE.
partir de la otra por
Dos matrices A y B del mismo tamaño, se dicen equivalentes si una se puede obtener a
medio de un número finito de operaciones elementales de línea {es decir ya sean de fila
o columna).
Observaciones.
¡ Se utiliza la notación A = B para indicar que "la matriz A es equivalente a la matriz B ".
. Las matrices equivalentes han de tener el mismo orden'
rr z -, 4r rr z -,;l ll : I ],-l=ll i ; ¿l
-¡ l=s
o=l ; ; l tl=l o o s -rl*[o o s
-
-z 6 ],1=L:' -; -;l-L;; ; -'l lo
o 5 -r]
[-r
REDUCIDA O CANOITIICA'
DEFINICION 7. MATRIZESCALONADA
tre.s condic'j:::
o canéni3¡i cumple las si¡uientes
Una matriz se dice escalonada reducida que todas los
filas son distintos de cero' mientras
1. Uno o mas elementos de cada una
elementos de las ñlas restantes
no nula es siemPre Ia unidad'
2. El primer elemento no nulo de
\2,3,...,r) es e[elemento 1 de la fila i
3. El único elemento distinto d
en la
iguales a cero) se encuentra siempre
4. El conjunto de filas
parie inferiorde la
a la unidad crece aritméticamente
En cada una
de fila a fi
, de alguna fila, tienen ceros en
Todas
todas
Entonces una
FRlh'lER Er-Ef\rEfo
Observación.
2 e matrices escalonadas
. Las matrices que no verifican la condición
I
[r o osl tE 0l iF;T r
1 o.l
^=lotr : 0 nl Io o 0 :l
lo o I -3 1 L I
IL
Ejemptosdematricesescalonadasreducidasotambiénllamadasmatricescanónicas
sl
Ir o ol [o 1 o ol [t o 0 -l
:l o=lo I ol a=lo o ol o=lo I 0 -tl
A=[:
^-Lool L;oi] '
[ooor] foo 1
;t
-rl
3 -rl
[o i + o ,-l -r 3 orI
0
[o o ['l0 0
o
0 +l
I
I
o=l: : : I ;i o=lo o I 5 0l A=l
0I 1l
o o o o] Lo o o o ,l
-, ln 0 0
0
0 ol
L; L0
-25-
ING, JULIO UBERHUAGA C.
MG. SC. ING. RAFAEL VALENCIA GOYZUETA
ALGEBRA LINEAL Y TEORIA MATRICIAL
equ¡valente Por filas)'
PARA EL TEOREMA I
Las rnatrices elementales
2t34 -2f1+f2-+f2
00 1 ool matriz F = E¡ 'Ez'Et
4215 ,]; 10 -2
1023 01 ;?l F=
::li;?l
Lo 0
_z.tt
Ir
fz-+
00
li
'1
,i: t0 ',111
o -5 -31 "fz¡+h it¡
01.
_1
.2
0
:I =
i+ : ilr ,:l
,.§
213 4
F;
l;;:l
L-2 r ol
v o-| + r
Ir o ol[z r 3 +l
0 0 -5-3 ,F. A=l-l o rllr 2 t,l:u
l-zro-] Lt 023)
\ t: i ril:'
-!l:-.fz
operaciones
H1,H2,H3,"',H, todas :stas resultado de efectuar
2. Sean las matrices elementales
hallar la matriz equivalente, las matrices
elementalesdefilaenlamatrizidentidad.Entoncespara
elementalesasíobtenidaspremu¡tiplicanalamatrizA
Hn""'H:'Hz'H1'A=B
Siendo H1 es la primera operación de fllas' H2 la segunda' H3'"'
eSas resultado de efecluar oPeraciones
3. Sean las matrices elementales K1,K2,K3,"',K¡ todas J
para hallar la matriz equivalente, las
elementales de columna en la matriz idenüd1!'if-l
a la máriz A t,
matrices elementales así obtenidas pos
¿
A.K I=B
.-.¡
Siendo Ki la Primera oPeración de
trices elementales
4. Cualquier matriz se Puede
H¡ ',: "'F
De donde:
P=H¡r J
Portantolacond§i3nnecesariaysuficrentepaia::Í:*g,*ffffi.;"#]:H:T.":fil""iJ]:i
:
de manera que se u"t¡f¡q1H$:§ffE v ¡!*'\r
' '
dos matrices P *Egulares,
Observaciones.
deduce'.t' una matriz resular
puesto or","tfi¡¡rquivalente -fg.ffiffi
uu p'on'q
: ":
.
"
la matriz A '
. Las operac¡ones elementales no modif¡can ni d orden n¡ el rango de
matriils
"on
Propiedades det rango de una matriz
Las operaciones o transformaciones elementales aplicadas
a una matriz A no modifisan ni el
1.
orden ni el ranEo de la matriz
z. sea la **J=ü] ,; nu,. o" otden ntx" entonoe§ 0'p(A) I min {*'u\
g. Sea la matñz *rCoO" A:[or;] , no nula de orden n, entonces O < p(A) < rl
4. sea A y AT entonceu P(A)
r)
5. La condición nmaria Y que dos malrices A YB de orden
'?'x'?
sean
equivalentes es que
b. SiAyBsot respecto de la §uma A+B entonces
! p(a +n)
de orden p x /, es tal que A'B =0'
8. Si el
! el.
! 9_
ser superior al rango o
ue su determinante es no
10. T
nulo
11. El rango
ngo /" , siendo 0 la
\. 12. Si A tiene
matriz nula.
13. Sila matriz amPli §o menor o igual
ar á una matriz canónica
14. Toda matiz A dé orden
n'a matriz cualquiera A
también de orden Íltxtt .
también Podría hacerse con la A. Por lo tanto las oPeraciones
elementales de fila que se re ción de una matriz escalonada
canon¡ca
15. Si A
es una matriz regular, el rango la matriz A'B (Y también de B.A) es
l-z 1 zr -21
Ejemplos Hallar el rango de la matriz A=l-l -r ol
L, | 2 -rl
en la unidad y los elementos que estén por debajo
El objetivo es que los elementos de la diagonal se conviertan
sean ceros.
Primera iteraciÓn. Muftiplicar por -'t la -.fz t fz f-z. t 2 -2
fila2 Li-r o
[-rr:-1
Seg{.rnOa iteración. lnterc¿mbiar las Ir r -r o
filas 1 y 2 ¡<- fz
[-r t 2 -1
lo
l-'> 1 )
-2
L- -'l
GOYZUETA ING. JULIO UBERHUAGAC'
MG. SC. ING. RAFAEL VALENCIA
BRA LINEAL Y TEORIA MATRICIAL
¡2la Zfi+fi-+f3
primera fila Y sumarla a la tercera ¿
La matriz que resrha de
la
f,, -1 ol septima)'
J
- lz iitin',. itutrcion 6a
lo I o -rirl
po' /A Itz ,-)
Quinta iteración' Multiplicar ^ .. eouivatente a la matiz A '
[o 3 o -tJ eita'matriz tiene dqs ñlas
no
¿
fr o -l t,3l nriát, luego su rargo es 2 ¿
@Por.-l ala -fz+ft+ rt
fila Y sumarla a la fila lo r o -url p(A)=2 J
Lo-e o -r
1
r"ér.0"
l ¿
@iPlicar.Por-3a -tfz+h-+fi ¿
f, fila y sumarla a la tercera
lugunO" J
)a P Y Q, tales que
Ir 1
las matrices
EjemPlo. Reducir la matriz A = I
2'
t3
P.A'q = }[,, ..:,'i
ar/
Como la matriz A es de orden-3x4
¿
de filas Y columnas
r' o.
z -)
00
1000 0 l'l 0i.
10
0100 00 1.,1
00 0 i
0010 00 0
10 0 0100
r3
N- 0001 3
ó-o-57-zio
t 2 3-z 100 -7-
-7
0--6-57-i01
2-2 i 3 010 -30
-3 0 -{ I 00 i -zct
-fz+h+f¡
-Zfi+ fz -+ Íz -3c1+33
-3fi+ -h -+ ft 2ci+34 -+
r1-?')
| -L
t )'i !.; -l ''
-;'1 3
100
) L* -l A p(A)=2
1000 100 10 00
0-6-57 -2 l0 00 -2 10
01 -1-11
00 00 -1-i1 00 00
-%,, -) c2 c"
-%,, *cac3- 5
c3
I a
P_= -.
I
1 tl A Q= I'lo -i¿
o I
/u I
ol
A
00 0,1
I
L--1 -1 1l
[oo 0 r]
I'G. SC.ING. RAFAEL VALE}ICIA @YAJETA
ffitrUNEAMAW
INTRODUCCION
Los determinantes hícíeron su aparición áticas más de un siglo antes que las matrices. El término
mabiz, precisamente, fue creado Por Sylvester con el propósito de distin-guÍr entre matrices y
determinantes. De hecho, su ide ?nadre de los determinantes'.
::jl
Entre otras cosas dicho numero proporcibna criterios de invertibilídad de A y resulta ser también igual al
volumen de un paralefepípedo P en el espacio n-dimensional, donde las aristas de P vienen de las filas de A .
Eldeterminante se denota por alguna de las siguientes formas (nosotros adoptaremos la primera)
oz2 a2n
Esto es. Sea Ia matriz A=l'?' I
.lrl I
I on, orn)
qtt at2 aln
0¡.I t A¡'t a2n
lal: LL
, este número se obtiene operando de cierta forma con los elementos de la matriz
::
aln ün2 Qnn
-31 -
MG. SC. INg. RAFAEL VALENCIA GOTZUETA ING. JUUO UBERHUAGA C' '
ALGEBRA LINEAL Y TEORTA UATRICIAL
An lel
. Eldeterminante es una funciÓn de R'"'r
en IR
. La correspondencia entre eldeterminante y los elementos de tra matriz es única, luego eldeterminante
de una matriz cuadrada es único.
'#.l
Lr
N eñ N los elementos
Estas biyecciones de
sideran ahora
de N, que inicialmente se su
elernentos. El
formando Ia nueva sucesión (
conjunto de todas ellas se deno
22
En este ejemplo con los tres elementos se han obtenido
Transposición. Se llaman trasposiciones a aquellas perm las que, salvo dos elementos de N,
que vamos a llamar i y i. Todos los demás p decir, coinciden con su imagen); los
elementos que varían se transforman uno en el otro (es decir y 7 en i).
indice. Se llama Índice de una permutación y se denota por ¡(o), al numero de trasposiciones que presenta
dicha permutación
-
Y DEFINICION I5. TERNIM) DE UNA MATRE CUADR.ADA
Dada una matriz cuadrada O, =lo¡) de orden r se denomina término de esta matriz a:
v
(-rld = q
üfzjza3jt
--- ary,
Que contiene como fuchres solamente a un elemento de cada fila y uno de cada columna, además
jt, jZ, j2,...,7, representan a cada uno de los números (1,13,...,n) no necesariamente en este orden y d
es el número de decrementos de la permutación jy j2,ft,...,irrfurmada por los segundos subíndices.
att
: Eiemplo. Sea Ia matriz A = dos térmirps de esta malriz:
v
-
(-l)d : ot segundos subíndices forman la permutación
h' jz-j
: ( subíndices forman la permutación
Y
t
-33-
MG, SC. IHG. RAFAELVALENCIA GOT¿UETA ING. JUUO UBERHUAGA C.
ALGEBRA LINEALY TEOR]A HATR'CIAL
t a,l= orrorr-anazl
Si n=2 = P2 ={(iz), (Zt)) v por tanto se tiene lAl=lo,
azzl
lozt
Desarrollando el sudatorio: lAl=l',
t anl:(-,)r, at§zz +Gid' otZa¡l = ollaz2 - alZal1
l'zt a22l
Donde:
-
dt=0la permutaciÓn (tZ) no tiene decremento -
dz:1 la permutaciÓn (21) t¡ene un decremento -
-
";:l = .(-'iB
lel== lr, t o.t.t o-¡.
"
Es decir:
[]l C
v
"::, \zoz:ast
Desarrollando
Luego:
lAl = a, f22s33 - a12a2§33 * aga2p32 - a¡a23a32 * a1ta23o31 - aga22a31
l¿.a"1
: É" lAl '
orrl: .l,.l
tc2
azzl
-r
8- Una matriz es
:
9. Si la matriz A posee constituida por ceros, su determinante es
: nulo
l3 ;l]
10. Si se intercambian dos líneas (fllas o has) del determinante, eldeterminante cambia de signo.
bl=o¿-rn ol=ru-ad=-Gd-cb)
l" dl
lc = lr ld cl
es cero
12. S¡ la matriz A tiene dos líneas paralelas (filas o columnas) iguales, entonces su determinante
lal= o
f ol
la "l=ab-bo=0 v \a bl
l- 1l=rb-bo=0
ib bl la bl
13. Si la matriz A posee dos líneas paralelas proporcionales, Es decir tiene una colurnna o fila que es
múltiplo de otra columna o fila, su determinante es nulo lal:0
-i!:-;r+
4t nz 4n nn
dl*
dnl dn2
en dos
en en esa
fila o columna los Primeros Y s los mismos
e¡ehenrile er oetermrnar
D e+'
) ,l
Ít t
+
{-
I
lb d f l=o+o:o
I
ls.
t,
pd pA
Aplicando propiedad 13 Aplicando propiedad 13
lineas paralelas lineas paralelas
17. Si un determinante es nulo entonces, las filas consideradas como vectores son linealmente
dependientes. Si es distinto de cero serán independientes'"'
1g. El determinante de un matriz cuadrada no cambia si se le suma a una línea (fila o columna) cualquiera
una combinaciÓn linealde otras líneas paralelas a ella'
--
tn
Esta propiedad es bastante util para los problemas en que nos pidan estudiar si un conjunto de vectores son dep€ndientes o
independientes sin tener que hacer combinaciones lineales.
-36-
v ALGEBRAUI{EALYTEORTAHATRICIAL MG. SC.ING. RAFAELVALENCTA@YZUETA ING' JUUO UBERHUAGA C'
v
APlicando ProPiedad 13
Iíneas paralelas
lg. Todo deteminanb de una matríz cuadrada se puede convertir en otro del mismo valor que el dado, tal
que todos los elemenbs de una fila, previamente elegkla, sean ceros, excepto uno de ellos.
a. Es consecr¡er¡cía de la anterior, bmh aplicarla de forma reiterada y conveniente para
: conseguir etre cierta ñta de la matriz, sin que var[e su determinante'
Frecuentemente en ceros los elementos de las columnas sucesivas, a partir
v de Ia primera, diagonat principal. De.este modo se obtiene una matriz
triangular y fácilmente.
: dsl
v tJl
\? ser una matriz diagonal) es siempre
t,
21. Elvalorde I
sr lAl + o entonces
lo 'l = Él
Si A es inversible se cumPle
De donde O"rp"¡rfu-
CALCULO DE DETERMINANTES
hace mucho mas Pesada a
Regla de SARRUS. La definiciÓn de determi
meá¡da que aumenta el orden de la matriz'
el cálculo de dichos determinantes-
En el caso de las matrices cuadradas de orden 3,
"tqri."trP"ri"rl% -38-
U
ING, JULIO UBERHUAGA C'
ALGEBRAuNEALYTEoRAIATRICTALIG'sc.ING.RAFAELVALENGIAGoYa,ETA
Y entonces
d'14
\-
y{
üú
:
v
v
v
Asi pues:
2t'
una matriz de orden 2 tiene 4 menores complement"'"tf
rP¡trrrr'rrta,uDL ?lf ?')
r JIl J
-- |e-'' t)t tt(z-l)t
=2'2 = 4
anl
a=[''t azt A
lozt ) = ",=ilillil lil;ll lill=;l lT^',:l=:::;
l:\l¡) 3t 3!
=9
una matriz de orden 3 tiene e menores complementrt'"t
[r.J[, )= fb,."6-ú.=3'3
-
p al proceso de obtener un menor de onden
p+i v
DEFINICION 19. Se lfama orlar un menor de orden '
a
añadiándole una fila y una columna al menor original'
menores de orden trz+1, que Pueden
Si un menor de orden h de una matriz A es distinto de
s que nc frguran en el menor son
formarse orlando este con la frla p de la matriz y cad3
n en ei menor.
nulos, entonces la ñla p es combinaión lineal de
Así J
Los colactores de una matriz de
ü' ''r
- 4rl ,#
. lr¿i
^
A =I liI
IL o¡t o11'
jt L+\
-.¡,
, at2 o,r
O=lo¡[', a2z ort
l
= Ai= Y:"1
32 1;ll
u33 I
;tl l;l ::i[ I
Lo: r a¡2
l
azl lerrl ' lon ar:
)
*ln, I
lo,, '''l
orrl *l'" '''l I
DEFINICION 2I. FO
v
Eldetermi suma de los productos de los elementos de
una línea nto). Es decir supuesta la fila p el
de
a
O¡¡'Ap*
:
t7 Para la malriz elementos de una línea
(fita o columna) por los línea (fila o columna)
paralela siempre es igual
Si se consideran Ia p -esima
Finalmenle podemos anotar que un determi esarrolla con cualquiera de las formulas:
Ir , -'l ll 1,
L¿I3X}I
1-catcutarlel=l: + sl. (6) (5) (1)
3 - Calcular e[ determinante
,-l
Fri¡-¡ar.¡ i.l,f V l--.ril;i¡= :-r:ltii'llnn f.afA
-t
5l
I _t
--'lJI
l
=1l r il=(,)[.i]
\ :,1).(,)[_li j,l).(,)[.li ;l)
lr 2 -31
lal=a{(+)(-:)-(sz)}*(-5){{z)(-r)-{:)(r)}+(:){(z)(z)-{t)(+)}
= a(-zz)-5(-lt)+:(o)
. = -"'.tt = -"
=
' ING' JULIO UBERHUAGA C'
ALGEBR,A LINEAL Y TEORIA MATRICIAL MG. SC. ING. RAFAEL VALENCIA GOYZUETA
lr o o rl
l0 -1 0 4l
4 - Calcule el determinante de: ¡.+l= l. j, 0 I
V 0l
li s -z ol
ú
Desarrollamos el determina
la pnmera flla
¡¡=3 e
-) *an'(-r)1*3lrrt,rl)
= lal: f
j=l
qr.( trrl ,? iorr-irt,:,)
;¡+(:) lnt,,l)
:ii'§
Vliá:*tíera columna y segundo de eliminar la primera
Los menores resultan de, primero eliminar la primera
fila y tercera columna:
ffi+
s=l€ -1 4l =
}
l-
}irr =l-'
1 ItJ
= \[rr=lt')l'
-11
il
l_i
L§3ol
Remplazando Y oPerando:
Y
ING. JULIO UBERHUAGA C.
v ALGEBRA Uñ¡EAL Y TEORTA HATRICIAL
MG. §C.II.¡G. RAFAEL VALET{GI,A COYZUETA
v 2 5
-1 3 0
| 2. -l 3 -5
5 - Calcular eldeterminar,t*, lAl= 2 0-1 11
10 2 L -l
2 t -l 0 0
es el elemento osz=1.
para Ia primera itenación: La fila elegida es la 5, el elemento ¡rivote seleccionado
es:"*',,j:i.'§i
elementales y evaluamos
Realizamos las operaciones erementa
qofadot: el.
-, l-zr 'N\ z
.:r,.T:
1 2 -r, -tl -+ l-r2A-1 z 1\si :-s l-z 3 2 sl
lel= lz o -l I ll
11
-2cz+q = tAr=(1)(-)l; j, i -il
021-1
r 0 2l l, zt -,1
zE-r.$S. llD.oo
pivote seleccionado es el elemento a-11 = 1'
Para la segunda
Realizamos las
-
{rl-
,,, 3 (ln
lei=-i i,
! ,-) I t-
I
-l r,rr-l l, ?_fl
-l I
\ l-\
r- -r
,lJ
[-ln z
!,qt= 1'; : \
z-il
t_ I t,l
_,l--i .fl 3
Finalmente:
23 )\2
=
fr
l) t
t_)
r
6 - Hallar el rango de la matriz -l
l) 4 6
L-
Tomamos las submatrices mas grandes (de orden 2) en la matriz
[; 6l [¡ -elI [6, -91
\.{
.,, I
=1, +) I
,\ lvl .=i
" 12 -61 ¡\ NI.
' =l 14 -6-]
I
BEUMffi
DEFINICION 23. MATRIZ INVERSA
Sean tas matices O, =lo, _ ,, =[Ar], matrices :-=:-ajas je orden n, tales :-e A'B:B'A:I
o tambié:r s::3e que A es la
entonces se dice que B es la rne:-z irversa de A y se ie-::: ;:: B:A-1,
inversa de B y se denota Por A = B-i
11: qJ. una matrlz cuadra:= :':se= ir'versa es
que su
Teorema 23.1 Una condiclón ne:=sa'a y suficiente para
forma escalonada c¿nórica sea ia n:="-u 'ridad'
Demostración. Si su fci;::a escaio:=:a carÓnica es I, , entcnces existe la mati-iz B = IK"t''' tal que'
= B=A'
1
B.A = I,i
Siexiste A-l tal que:
A-1 'A=Ir'É..1 B=A-1
Tal que B'A: I,, y por tanto, I, es la forma et"abtdL canónica de A '
:
t
Teorema 23.2: \ es inversible si, y solo si, lA!* '
0
Pordefinición tenemos:
A.B:r:B.A A A'c:I:c'A
De estas dos igualdades se deduce que
A.B:A.C
1t elementales (filas o columnas' pero no
E.t" teorema nos permite c¿lcular la matriz inversa, de una matriz dada, mediante operaciones
ambas simultáneamente).
_{U_
ING. JULIO UBERHUAGAG.
TllG. src. I}¡G. RAFAEL VALEI{CI'A GOY'¿UETA
ALGEBRA LINEAL Y TEORIA MATR¡C¡AL
O sea:
A.B-A-c=o -? A'(B-c):e
Como existe A-l entonces:
A+0=B-C=OeB=C
única'
En consecuencia la inversa de una matriz es
si, trlene una línea combinación lineal
de las paralelas'
Teorema 23.4: tJnamatriz m singular (lAl = 0 ) si, solo
y
Demostración:
forma escalonada canónica tiene una fila
de ceros'
a) s¡ lAl = 0 algún pivote es nulo' por lo que su de las
fita era necesar¡amente combinaciÓn lineal
Deshaciendo las transformaciones efectuadas' esa
demás. ', :\
por propiedades
a'
4.
CALCULO DE LA
or definición se
lnversa de una matriz trian
tiene A. B = I,, , esto es:
i
i¿]tt
I
i,r:r
t:
I
ouu | | b,;
| ünl
I
-t7-
ALGEBRA LiNEAL y rEoRrA MATRTcTAL 'MG. sc.llrrc- RAFAET- vALENctA GovaErA ING- JULIo uBEU{ulAGAi
,lllr:
Reiterando este procedimiento hasta la r¡-sima filade A podemos concluirquesi una matriztriangularinferior
A es inversible, entonces se cumPle:
. Todos los elementos de la diagonal principal de A deben ser distintos de cero.
,tI I 00 0
lot, 0 0
:l:l
o
_l 0 10 0
I
Si A es una matriz triangular superior su #tamUen es una matriz triangular superior y para
:
calcutarla se debe resolver la ecuaciÓn
0
I dtz
I
fr
t:
az2
L; ;
rersa, de una matriz
lnversa de una matriz Por
dada, mediante transformacio¡ mente). El método
's
consiste en lo siguiente: jii'{+
s'
,In
1. Para la matriz t*l ] ou oro"n ttx2tt
an
2. o; de openaciones elementales
si la l,[ÁlI,] ru"so
Adi(A):Cr = [r7i] =
-48-
ING. JULIO UBERHUAGA C.
MG. SC.Iil¡G. RAFAEL VALET{GIA GOY¿I,ETA
ALGEBR,A LINEAI. Y TEOfiTA HATRICIAL
1. Adi(Ín):rt
2. dat(N',):(rd(A))'
3. daL (a.r)=(ea¡1e¡)r
4. ea¡ (e,.e):,a4(n)'Aü (L)
5. ea¡(tc' A,',):kn-rAdj(Ar)v fr escalarv Aoe orden n
6, laa¡(x,,)l:le,l"-'
7. ,lar(a'):fr
8. A.Adi (¿):lel.t,
9. aa¡(aa¡ (t',))= |
rc. l,ta¡¡eat*(A
rr. ray(¿ffi
oesarrouáSi
Partamos de una
Sea la matriz A-
Luego se cumPle
nos permiten
Ahora los Productos escalares de
establecer las siguientes ecuaciones:
A^-,-l [r iFl-
l=l
l^l
]r,_l q2L 0u
I
= A-1 -
II ct22 -r,rr)
[: I
I
I
lolL-"rr orr -l
L-FT hI ]
Por otro lado hallemos la matriz de cofactores de A
-{9-
ING. JULIO UBERHUAGAC.
ALGEBRA LINEAL Y TEORIA MATRrcHL MG, SC. If.¡G. RAFAEL VALENCIA GOY¿UETA
-or,1
=l*",
A,,
+at ) = ' =ll?,. .'1r",)= Adj(^)
"i L-ozt
' l-ntz
que:
Comparando tos dos últimos resultados concluimos
_r Adj(A,)
A=1.!
^
v
*
fue desanollado Por Le Verrier Y
lnversa de una matriz por el esguerna Fadevva. Este método originalmente v
se tunda en las formulas de Newton Para la
iimptiRcaoo por el matemático soüético Faddewa, este método
sumadelosgradosdelasraícesdeunaecuaciónalgebraica. -
propuesio a conünuaciÓn:
Sea la matriz A, =[a¡] craOrada. Elaporte de Faddewa es eta§oritmo
Bt = Al -qt'l
,l ,$'
.É"+ r' (41; ) Bz=Az-qz'I
También:
M_ I
Donde:
AelKp"p A Ben{P*Q A ceIK{*P A DeIK{'{
lnicialmente se supone que la matriz D es una matriz no singular
que esta particionada en la misma forma que la matriz original, esto
es:
La inversa de la matriz M , es M-l
ING. JUUO UBERHUAGA G.
ALGEBRA LINEAL Y TEOHA IIATRIC¡AL HGt SC. ING. RAFAEL VáLEIICIA GOT¿UETA
X
ht -1- -.
_:
---. -i------
l
Donde:
XeIKP,P A Y=IKP,{ A Z elK4*P ,r U
Luego se cumPle:
q l.x lrl
M.M-r = I = Dl Lz lul
Desanollando el Producto
Obtenemos elsistema:
Suponiendo que D es no
Y = -X'B'D-I (s)
De (5), (6), C0, (8) tenemos la matriz inversa buscada'
lnversa de una matriz Por Orlado. ¿
eIKr nosingular,particionamoslamatrizA como n-(n-1)+1 ysuponerlmslamatriz
SeanlasmatrizA J
A como una submatriz de la siguiente forma: J
9
r/
t- ./
a
t ¿
:{?
Donde:
g etr(("-r)x(n-r) A C .n<(,,-j ¿x(r-t) ". dr*, €E<1*1
J
\/
;;r; '- i¡.1- :i
::,,.. '.i:..:; :-S
'{i::. .
:i,l.i:::*F1'
¡1r:i: -':,t¡'.-
"'iii'
!:
Donde:
Y e IK an e [(l'I
Luego se
.jr 1-l II el
A
,+.fri' . .{* ¡:-.1:¡..¡¡ =l
-
?,¡j
::::r
tdr, ]
lrl
L- 1l
,.
'
.,-..j ..5É
jF]:::-l 3F
Desarrollando el Prod
l
i"rreiE,{'§
. r ', f-: :: l): I
i;ffi[¡!: .r-] I
t
'q,1 l=[; I
I
?l
De (3) desPejamos Y
B.Y+ du'C = 0
B.y = -C.a, /le-t (izquierda)
y = -F^-l.C.an
(s)
- 52'
GOTZUETA ING' JULIO UBERHUAGA C'
ALGEBRA LINEAL Y TEORIA IIATRICIAL ftIG. §C. ING' RAFAEL VALETTICI¡A
(o*' ll(ar*r-D'B-1'.f'
r-1
-D-B-1 .c)-'r.I
,
(o*,- D'B-r c) (A
r-l
. B-1 c) (o)
De la misma forma,
CombinamQ
Como:
También:
alu
= lal A
= Ad¡ (a ))' uL(x ): (ur,(a))r
Por Io tanto
-53-
MATRIcIAL GOTZUETA lNG. JULIO UBERHUAGA c-
ALGEBRA L¡NEAL y TEoRIA MG. SC. ltt¡G. RAFAEL VALENCIA
eat(rt) =
[e-rl(a-r)-1
[a¡-r
(o-r¡-t
ft^f' u-')-'.
¿
aa¡ (n-')= 1,147{a)}-t J
consecuencia (aa¡1e¡)-'
En h, -: J
3 Demostrar que ,la¡ {d.g) = Adj (U) ,Aa}.(+)
J
le!*o A lel*o
;d¡@ n) : la.elaaT(a-e)-r
= lal lel-(e-r.e-l)
¡e¡ e-'){1al o-')
./j(B) .u¡ {x) Ad¡{A.s) = Adi (e) ,a,:7(e¡
'
la"i(a")
!a[" {a-})"
ítenemos:
(lul o-')'
nai(;, )= (,+,1{a))"
{aa¡{x))'
l* al(*.e)-1
r, lAl.r-t A-1
kn_t.(lol o_,)
k'-r . ,+a¡ (e) u.¡(tr.A) = t'',ld¡(l)
Hallemos la adjunta de la
adjunta lel"-'' [l^-'l (u-')-')
i
l^1"-' t- ,'A
lAl
lal"-' A
-5J-
v
ING. JULIO UAERHUAGAC.
v I|G. SC. ¡ñIG. RAFAEL VALETTCIA GOTZUETA
ALGEBRA LINEAL Y TEOH'T TATRICIAL
v
Hallamos el determinanle
laa¡(ea¡r.e))l = lnft ol
= (t"r-f nl
Finalmente
Adi(Adi(1l)) = lof'-')'
Si AyB son matrices
!
v Si A üene inversa
También
Desanollando
= lo-^-'- B A-')(A r B)
= io-r,'A-r'(e+a)
del producto A por B es igual al
si A y B son matrices cuadradas no singulares, entonces la inversa
producto de las inversas de A y B en orden permutado
Si A es no singular, entonces A tiene inversa de la forma: A-1 -,' A.A-1 = A-1 'A = I
_ B-l .B = I
Si B es no singular, entonces B tiene inversa de la forma: B-t t,B. B-t
Del enunciado tenemos (e g)-' = g-1 ' 6-t
GOYZUETA ING' JULIO UBERHUAGA C'
ALGEBRA LINEAL Y TEoR¡A MATR¡CIAL ¡t4G' SC' lltG' RAFAEL VALENCIA
:t
Por lo tanto B-l'A-1 es la inversa deifr,
Si A.B=B.A=I {+ +A3 +"'+A'-l
=-:I+A+A2
tT_ * g$igl* o1-'
El problema pide demostr \,
r3
Se debe ( [:+A
'
"rrP"t'
.;:'.ir.i.l ::::
'¿\
rambién'""ffiH ( ! A-+
1
'.i¡
"4
! i:¡
;f"i .nr }Élt;.
F.é."
"ttt:*! ..,, '1¡ u)(, +A+ k+A
^.
T-A,?
I-e
I
10
Desarrollando la
segunda cuestión
Efectivamente la inversa oe
Demostrar
(I-A)
que si B. A = A'C y A
,il'.-r'
t es
'A"'rrl
I_A,,
i-e
I
l0
= f =t-,1 [?
l0ll=
2->c2
-+.1 ^l 6 üi -+l
-(:z-eo) =zg
.l-;' :l
1l l-z 4 zl
e=l-10
Ir , ,l - a 68 lo|
= A,¡=l
\-31 lsi
L,o
l7
Segundo la
,tj
Luego la inversa de la matriz A i 4l
7
1 s el
J
L
3 5 -1
Calculamos eldeterminante de A
l¡
.1,
rl -1,
lz 1l *1,
*lr rll
-,1 -rl ,l
s
4 sl sl sl li [-s 1 1l
[; l+ *1, 4l
:) to :l zo -18 -:
-+l 'l: -il -1, ,l I
I
s -1]
l+ sl _l: sl *1, ol
[,
I
-57 -
GOYZITETA ING' JULIO UBERHUAGA C'
ALGEBRA LTNEAL y rEoRtA r,tATRtCtAL MG. SC. ll¡G. RAFAEL VALENCIA
Segundolaadjuntaesalatranspuestadelamatrizdecofactores
zs -l il
Adi(^)=(^r)' = Adi(^)=[[;X -:, :-|l' -18 tl
[L-rr , ')] LJ 7 t) :J
Luego Ia inversa de la matriz A es: ¿
t_t Zg _l zs -t tl
A-r = o-t = , -rB 7
'l -18 7
¿
¿
=
fiza71e¡
lAl 1l
'[r 7.t [j 71 ¿
lr o 0
o:l-, 0
l,L'
2
3 - Detenninar la matriz inversa de
)
principal, luego
La matriz A ¡nvers¡ul
planteamos la
= 0+1+3\t-O = b:z=-713
Luegoramatrizinversaes:
Ir 0 0l [ 1 0 0
,l-ll
A-r =u=lrr, 1i2 o l=lliz t'2 o
an
[.,
lo a22
Para by
[: T t]L':'
se efectúa el produdo escalar de
t:'
[::?] íil
nda fila de A por la tercera columna de B , esto es:
o l-{-:¡lo r ol=lo , ol
|
-L
Ar=l 0 0t =
I
clz --
-2-2-2 ,,- =l O -L
a
[o i -rl
I
[o 1 -z) [o o t.] [o 1 rl
A; = A'Bz ,,, =Ir'(ar)
J = B;=A3-qt'l
lr o ol [t o ol [t ool
o.,=lo , ol
l+1+1
-i = e,=lo or ol-(r)l
rl
I 0l=0
0
o
lo o 1.] lo [o 1.1
-59-
J
GOTAETA ING' JULIO UBERHUAGA C'
ALGEBR.A LINEAL Y TEOBIA MATRETAL MG. SC. Ii¡G. R/AFAEL VALENCIA
[t I ol J
¡ -l Bn-t -1 _8, =lo ol J
A'- A
^ = 1
=
Qn 8i
[o l ,l "¿
¿
l-s z ¿
6 - Determinar la matríz inversa de 5
^=l-o -Fadewa ¿
[-, 3 ;] ¿
Desanollando el algoritmo tenemos:
¿
rr(ar)
Al =A = ql=_
|
F¡^
l-) / o-l 10 ol
a, =l-+
I
s -rl OI ol
,-l
[-s 3 o.l
Az = A'Bt
ftt t o -+l -+l
\)
-el
'
o, =lro i -6
-r] -).
l,j
-r -61
I
A: =A'Bz
J
?$.u ,2 0
' o=.2
!
l-: o -2 I
=l s o -3 I
[], j, ;; -",1
7_Detenninartamarrizinversaou e=l | ; ; lo"rparrición
L, 3 -4 Io.]
Renombramos la matriz A y la particionamos, de la siguiente manera:
[1 3i-5:l
--1 o=[' :l
-l]l o=[-: ''1: .=[l ,
,ur=|"'i -?,1?
,",,:1.-? jlo i; I=o=l_, -sl "-1,,-r7l :] ,=[-oo ;]
Ir 3 i-4 10]
-60-
ING. JU LIO UBERHI'IGA G.
ALGEBRA L¡NEAL Y TEOATAtrATRICIAL HG. §G. ING. RAFAEL YALE EIA GOTAJETA
Calculamos'Y
Y : -X'B'D
|
I
-ztz
Loo \
[-l
I
sz :::-lIs
lt
¡o
-t l
-2J\
ie -t i
I+g -i8 i
l3 -l
-el
l-o ai
i
4 2)
L'
mat riz invers a es:
Calculamos U
U = D-l-D-l C.Y vl
+l isro -tr+l[t sl[+e -r*-l N,i-r = [x
ls
ro _1 I
IZ U-]
Lr/3 0 7
I -63 I
[+ 2-1 1l
Isro -1t
Lr/j 0 l-l 1t3 -r .l
[:
I
-21
l-r 1l
I
-6t-
\¿
i
GOTZI'ETA ING' JULIO UAERHUAGA C'
ALGEBRA LINEAL Y TEORIA IiATR¡CIAL MG' SC' ING- BAFAEL VALENCIA -l
,J
lr. I -l
e=12
J
B-Determinarla matrizinversaoe lootnartl"ion ,J
Ll -2 3l . --)
de la siguienle manera:
Renombramos la matriz A y la particionamos' ,.J
[r o .-tl D=[3]
*=1, ,i-,1 .J
¿
Ll -2,3) ¿
Como, lol = bl+ 0 , entonces D üene inversa
dada por: ¿
"#t -t)
Calcufamos X J
1-l 0, W tiene inversa
x= le-e.o-l
\
c)
Y
J
lll ol rlI t
|,\L- ,] rlil
n
[, [-r
llr
\L- 1
L-g
.J
al
(t, 0
ii¡
rl.3
'}¡
-o-l .c.x
i;i.'1fip,
-tir,6 til', -41.:i, :i',)
4l:]L-,
_1[-trzf _[ tio.l -!¡rsra
J
-rl
:i trz I i-tro-l l-sro t/31
matriz inversa es:
Calculamos U
U = D-l-D-l -l [x vl
= [i] [i]" -r1l 't61
lz rtl
'i -1l 6l Itto l/3 ttol
a ,¡vá1 = l-r,u zt3 -,tol
-L r[_uoJ
l-sre 1/3 rro )
J
r
_--i-t
rIt-l
J- JL¿)
-l ^ I
111
= ---=-
36 6 r
- 62'
@OT¿UETA ING' JULIO USERHUAGA C'
ALGEBRA LINEAL Y TBCNIAXATR¡C¡AL rG. SC. ING. RAFAE. VALEIIGTA
Ir -1
r=lz
t........
lr
0
= "=[l i] ' .=[;] ; D=[z -3] ; .,,,=[r]
L- i
B-l - f [ '-
-3.l_[-us 3/s
I
inversa dada pon
como, ld =ll 6
-sl'2 I I lzrs -1ls_l
Calculamos Z Calculamos ar
/5 l)-'
us))
[-srs
e t sl=[-s i:
*
-ít+li-+fi 100 luz o U2
¿
01 l0 r 0 h-f3 a
10 l-1t2 i 1t1
ttz o rtzi -
I
-fz+{t-fi I2 s=i -tiz t u2l
t
ol
zJi+-,f">f3 r 00
0 I
-
0
¿
|o I ol
ll
10 _t r
I
01 t
t
p
-t0 ol
)f:+fr-,.fr
01 0
1-
20 Ilt-+li
11 - Determinar la matnz I
m atriz[i,,lB]
Formamos la
22 \+fz-ft
-11
-t
I
ilt
-1 -5 ,,,:+rqij,f;,
fi+h'+f¡
2,f1+ [a-+ fa
) -ó
-i
",., ¡,lr:. i.I
I
i
-Í¡+rt 0 0 4
.fi, - f+
0
h+fz 0 1 -l :l
[r2i
1
_1 0 ¡. -1
lo f,l 2
r-l-12
r\-lrr 2
1
1
J 0
lo
10 o12i I 1 1 -l
Ir 0 J
z¡2* j: 00 nlr z 4 LI U -2 4
-1
.fzt-+ ft
01 o l, r
1 .1
_t - fz- »fz. 1 -1
00 t l, o I
t: 0
ft+.h-+fi -ó
lo 0 -={
02 1000 23 4
t'l0-1
lt l J t
_J
-2.ft+f¡-+f1 100 \1 2
-5
I r.,
01
--1
ln13 5
3ft+f2'-+f2 010 r\-r
3fi+fa-+fa 001 0-2\
1-3 loo-{ -i
['^l']
ELEMENTALES
INVERSA DE UN PRODUCTO DE MATRICES
Ya se conoce que:
P=]Hr""'H3'Hz'Ht ¡. Q=Kt'Kz'K:""'K"
*
ilG. §G. I'IG. RAFAEL VA¡'E}STA
TEOf,IA TTATECIAL
ALGEBRA LINEAL Y igual al Producto de las
de un pod'
y K posee irwe¡sa y la inversa
ffi ,ta de las mat¡ir-s H
: ¡nrátt"t"rd,
en orden contrario'
P-r - sll -Htl .H;l ...n;1 r, Q-l = r;t "'rt' 'rt' "r, '
Tenemos:
I .Htl'Htl"'H;l'n'r;l"'K31 'rtl "rll
A = p-l.B.e-l =Hr
:
v
de
También si B es una matriz
l'tas Rtas de la matriz
Luego el rango de A'B es igual
ntos nulos. Por tanto, el
N'Q-l 'B r
son las primeras filas de
caracteríslica de A ' De forma
rango de la metriz A B- no P'"9"
-1"-lma
,naiog, el rango de A'B no Puede
ser
puedenSernUmeroscualesquiera.Portantoe[rangodelamatrizQ-l.ey,enConsecuencia,la
que P - r '
de B , no Puecje ser mayor
reciprocas sobre la
inversas, se reariza las operaciones erementales
parra hallar las matnces erementales
matriz identidad original'
ALGEBRA LINEAL Y TEORTA IIATRICIAL MG. SC. tT{G. RAFAEL VALENCIA GOTZI.ETA ING. JULIO UBERHUAGA C' J
::IJ
OPERACION PARA HALI-AR LA ffiPARAHALLARI-A
MATRIZ ELEMENTAL MATRtrZhffERSA
-a.fz+ fi-+ -fi a' f2+ fv --> fy
t, It r I srl
,J; ? :l = ,,=il ; :l
t;
-i
;l = t; I :ll í :l +
ffi I u,,=[l ; :l
/
a)
loo,.l foo,J
loo,loo,l l-oo,loo,l foori ¿
¿
fi<+ fz ¡. fz¿+ ft 9
IIEiI ¿
irool lorol
,=1, , ol = e, =1, o ol
t00
l0 t.,ol
Eit =11 0 0l
[o o r.l Lo o r.j
I
= L. o rl
a'fi-+ Á / t
[r o ol
r=lo , ol Er=
o = tr:
[o '] r.iE iü.ü
a
l'! F,'Ll :il -
I
-a 'q +41 ",
g
-¡. : r
'{i
t- %5 l,T ci<
I
Fo
,lo
1la
-c2
I
r
0
|=
I
[r o ol
e¡'=lo I ol
L" u tl
rfr I .llEr'
tr o ol ff [oool
,=lo [r o olo o o'l
'
lo o
ol
rj
'ji.r,
'i. i lo r oh o rl =Efl =l'o rl
¡r:
--:'{
t.,,1,'ol [o r o.]
@- )
1
ct.\ ) c3 -'fi c3
a
.E1 | r rlEil
-66-
ING. JULIO UBERHUAGA G.
re. SC. tNG. RAFAE. YALEI{CIA @YZUETA
ALGEBRA LIHEALY TEOR1A IATRICIAL
an\ + + olnxn = bt
azfit t + oztzxtz : b2
: a3 + o3nxn = b3
§
ul,*
:
:r
Un sistema de
. Forma
. Matriz
de mxn'.
En elsistema de
Anxl -
Y
"t1
forma:
ot2
['' ' az.z
A.x:B + l'?t
t':
lan t ottL
--- Al adjuntar el vector columna B a la matriz A se determina una matriz de orden rnx(n+l) que
-67'
GOT¿UETA ING' JULIO UBERHUAGA C'
ALGEBRA LINEAL Y TEORIA MATRIC¡AL MG' SC, ING' RAFAEL VALENCIA
(),' '
I- ú
l.\-.]i
|
4-'- v
l-. - .'
La rnatriz aum sist
-J
')
LINEALES
DEFtNtCtoN 2G. SoLUCI9N DE UN SISTE¡¡*'DE EcUAE'oNES
incósnitas
se denomina sorución der sistema
': ';"':'"::- "'o
;:lif :';;"
Que verifican en forma simuitanea.todas las ecuaciones
del sistema.
la solución del sistema
El siguiente esquema grafico muestra forma de interpretar
D=TERfvlIJADO
SlSTFftlA$
f,enjlrilo
Saiuciór finito
CONSISTENTES
CO¡.,1PATIELES ;NDETERI,.IlNADO
S S-iElvlA
Tienen Solucitilr
-1E Cotr¡urrlo
ECUACiC:.JE S LiNEAL F.S SolLr ciór: i r-t [i n il.r,;
SISTEN.]AS
II..ICONSISTENTES
INCO[,lPATIBLES
No Exjste Solr-,ciÓn
ING. JULIO UBERHUAGAC.
,1'.QE3BI\ UNEAL Y TB§6UA HATRIC¡AL ilG. SC. ING. RAFAEL YALE¡{C¡A GOT¿UETA
dientes no es nulo
Es no homogéneo, es decir
DEF¡NlCloN 2e.Él§l
f 't "
El sistema lineati;-¡
Es homogéneo, es
+ Q¡¡tZ
tI¡ 1-Y1
+ tl12x1
+ {tZliz
I 'r,'t,
t:
I r¡,,, r.rr +üut,tf,=b,,,
La misma que puede ser expresada en su forma matricial o en su forma de matriz aumentada
1il- '
ING. JULIO UBERHUAGAG.
ALGEBRA LINEAL Y TEORIA MATRICIAL MG. SC. ING. NAFAEL VALENCIA GOY:AETA
J
P(A)= ,(o.)= P(ale) = r J
A* y rt es el nrúmero de
J
Si se cumple que el rango de !a matriz A es igual al ru¡go de la maüiz annentada
¿
incógnitas del sistema, estoes:
J
p(A)=r(o.)=p(alr)=r ,.. rr=# incosnitas ¿
Se presentan los siguientes casos: ú
. El sistema tiene una solución (solución única) cr¡ardo et niine¡o de
incÓgnitas dd $srerra e igual al ú
rango de la matriz aumentda ú
P(A)= r(o.) =P({ry= ¡ =tt -
\,
so-iuciones¡ cuando e{ n¡¡nero de incógnitas del
¿
. El sistema tiene mas de una solución (exisJen infints
sistema es mayor que el rafigo de la ma
á r<tt
De forma
-&,..
TIBLES
Tienen SoluciÓn
lnfir'ritas
U¡rrca Soitrciortes
Observaciones.
Para el sistema homogéneo:
. siendo el sistema de ,, incógnitas es fácil deducir que la matriz de coeficienles A y la matriz
aumentada tienen el mismo rango y por lo tanto el sistema es siempre compatible
[elg]
(llamada solución trivial) de la forma:
. El sistema homogéneo siempre tiene por lo menos una solución
Xl=X2=X3="'=f¡z=0
. Si p (A) = t, el sistema tiene solución única y esta es la soluciÓn trivial
la solución trivial (también se
. Si p (A) < ¡r el sistema tiene infinitas soluciones, donde una de ellas es
dice que elsistema tiene soluciones NO TRIVIALES)
GOT¿UETA ING' JULIO UBERHUA€A C'
ALGEBRA LINEAL Y TEO'RIA I'IATRICIAL UG. SC. ING. RAFAEL YALET{CTA
T
cttt al7 aB aln c1
ü1 at2
I a22 o.23 d2n Ct
A1 I u
0 a:; 4n c3
A* l
a3 1,,
ann ,"1
0 ol
t
rmaciones
Las propiedades enunc
elementales flla.
. Si se aplica a la al fila
h\\'\
(a.f¡ --) .í¡ , a-f¡, +f, rnatri2 resultanle
representa un sistema lineal equiva.
Método de Ia INVERSA. Sea el sistema de ecuaciones lineales de ¡r ecuaciones con il incÓgnitas (sistema
cuadrado)
Iar r.r'r + t\2r2 + + übrX,, = bL
I
l' trt
I'
:"
f
o,,, + at¡tzx2 + T Lrn¡'¡^Jl /.
url
A.X:B
pana 5espejarla matriz de varialles
Y al sistema asi expresadc !:,jemos aplicar operaciones ce rnairices
j a
Multiplicamos por la li¡¿;-sa de la n:ai;z de coeflcientes 1A-1
A-1 ,q..x:A-l .g
.
la ecuación (multiplice:"Jn por izor-:lJa)
7nx n
-a1 - !:
s el §istema tiene n esta dada por:
Diremos que es de Cramer
t"l
-11
I
b.l
'l
;l
'i
-l
,h..,
I
l^i
erüo¡ces nterior se tiene que X = A-l 'B es la solución
Demostracién. Si lqÉrt
única del sistema r++@ -á,A,
^ A lt7 -"'+4,Ar,, I
1+h2A12
2tt t hA:1 -á2A7r +"'+ bnAzn I
,1 '
'.1 I
'r=ff
Donde el denominador es desarrollo del determinante de la malnz
A, obtenid 2 e peftir de A sustituyendo la
[r,i
-l -l
l'r I [¿'
1",
I i Ii
r
.7 -esima cotumna
veaor
I not "t
-l
1
L *, L',, l
le,i
En consecuencia si j =1,?,3,"',tt se tier-' '- -l
re' 'tl -
r I
lAl
\,
GOYZUETA ING' JULIO IIBERHUAGA C'
\2 ALGEBRA L¡NEAI. Y TEOHA MATRICIAL ¡IG. SC. ING. RAFAEL VALET{CIA
:
Ir
¡,=lz
I
l1
!)
L
\-
lr -z :l['l [zl
el sistema en su forma matncial A'X:B
Expresamos = l,l', .;:lL:l=L;l
L-
Calculamos el determinante. De la matriz de coeflcienies:
tsl
fU -2 3[ -)
+l
r 0 0l
lal= lz -3 ,l 2c1+c-t + c-t z t -rl = lel = +(r)
ll ;l = -7 +25:tB
|, -t al --rrr +c-r --, c-i 3s-71
La solución por Cramer es:
-'13 -
C.,
ALGEBRA LINEAL Y TEORIA TATRIGIAL
MG.. SC. {ilG. RIFAEL VALENCIA GOY:4ETA ING. JULIO UBERHUAGA r/
7\J
J
ru
stz atl J
a2z az3 -3fr+fr-+f¡
xl=
a3z a-33
-2f2+fi-+f3
--1 LX 18
-5-49 =J
18
l8
lr,, \ or:
lr z 3l ¿
l, , E-l l-s -ll {
I
xz= lat
4 onl_1, sz I
-Zf.r+ f3-+ f3 ,l
-I818 =--i!-
--,1-,
t
=z
r,
¿
lal I
lru arz 4l
[j J
?
E -2
lor., azz árl 2 -3 h -31
l"r a:z 4l : J -1 lt 3l
-l-" 3+r5
---l
!= iÁi --'ll"
18 18
al
Luego el conjunto solución
v
\r/
prccedemos a escelona¡"
-2.f1+f2'->t2
-fi*"Ít-lt.,ili
+ 'l*fz-Í2,-'?r' 10 013
',"4 --a;:. .. "i,{ |-r-r--t (^
27 \ - : a!.. ¿r:a '
01 olz l.t=J
00 r ir t^
26
ffifl;,-,
o 1 0 Ízz- l2l I¡ tt rl
I
= 1)',=',
[z=l
ool: r} z
r ol,
za
l-z -2f2+f1-+f1 0r l,
1
ll -.fz+fz-ft
1
113
§
ING. JULIO UBERHUAGAC.
I/IG. SC. ING. RAFAEL UAI§¡{CIA @YZUETA
ALGEBRA L¡NEAL Y TEOHA }JIATRICIAL
INTRODUCCION
En el momento de resolver un si§ema de ecuaciones lineates
de n ecuaciones con ll incógnitas' podemos
es er método de eriminación de Gauss el cual
recunir a diferentes métodos- uno de los métodos más utirizdos
del sistema de
consiste en convertir la matriz arrnentada ,a.. = , donde A es la matriz de coeficientes
[Alf ]
ecuaciones, el cual esta expresada en ta forma escalonada reduci'a'
son difáíóntes de cero, entonces el sistema
§i U es una matriz triangular superior cuyos eleme
lineal U-X=B puede ser resuelto sin la matriz aumentada [Ulg] a la forma
escalonada.
,
[,,,
io
i
o
LulB_l= |
I'
v
L0
es de cero,
De forma parecida, si L es una matriz triang
entonces el sistema lineal L' X = B puede ser re
Y iguiente
La matriz aumentada tiene la forma:
(su
i4r 4z 4: It, br
I| o t-,,
-' t.¡., lz, b2
[rlu] = lo o /¡¡
I
t.
l:n b3 'T2 = /\I
t22 .
I' o :
Lo o J
,,,,,
b,, ,.I
b¡ -lt¡rxr
A'=l
'rj=-L ,J]
matrizdecoeficientes.Aquí,primeramentecons¡oer,,aremoslaEc-toñ'a"ionLUsinintercambiobasadaen
que da la factorizacion PA = LU'
matrices elementales conocida como de Doolitfle y luego ,uo*o.'ul-atgoritmo
Suponga que la matriz A es una ntatriz mxnse puede escribir como el producto de dos mairices:
A:L'U
Donde L es una matriz triangular inferior ntxnt y U es una matriz escalonada o triangular superibr mxn ' J
Entonces para resolver el sistema: J
A.X:B J
Escribimos
¿
A'X: (L' u)x:L(t-t'x)
9
una oosible estrategia de soluciÓn consiste en tomar Y:U'X y resolver para Y:
L'Y:B
Como la matriz L es triangular superior este sistema vérse mediante sustitución hacia edelante'
inicial se resuetve despejando X de
Una vez con los valores encontrados de y, las incógnitas
\.'á,:¿.-,fl'
Propiedades de la fa i'r-u -
i i¡-
de A v
Si se posee una facto,
;A:L'U de una matri entonces la inversa
puede calcularse medi
'*#i
.?
.1
A-r.R \/ qF entonces puede aPlicarse
Si
eli enla reducción
[rle] -+
.B = A-' 'B
lol-+ [ti!
:,1
Además, normalmente no ieñte realizar el producto P'B sino mías bien realizar el
movimiento de filas Y mías que realizar el movimiento de filas, se hacen trucos de
programación Para evitar talep rientos teniendo un vector que refiere a las filas de diferentes
posiciones.
Ejemplo. Determine una factorizacion LU de la matriz
lz , -'l
a _t l-o
I
-6 _5 I
¡\-l
l4 l8
t?
I
-e '
-31
I
-76-
!
§ ING' JuLlo UBERHUAGA c'
ALGEBRA LtnEAL y rEoFtA trATFuc¡AL x{L sc. lNG. RAFAELVALEIIcIAGoYZUETA
v del elemento all=2 pata ello debemos
E
Así para el ejemplo. Hacemos cero todos los elementos por debajo
reaf izar las sÍguientes operaciones
ones Elemental
1000 [t o o oil , , -, 2 3 -ll
_l
l¡ i o oil-6 4 5l 0 3 2l
I
3100 u'o=l;; t_
3f1+ f2 -+ f2 E,=
0010 iril * rB 6l= 4 18 6l I
Jr
itH
o
.fi+ fa --> ft E -lo 0128 012 I
K: :?
J
-2 -9 -3 0 -6 --+
Aquí, comenzarr¡os ia
matriz triangular inferior
cc
ir o o oi
t-3 I 0 0l
ry
diagonal pnncipal. Para -q
N :: ; ?l
it,=
l00u
0100
[rr]10
lF
"=L:
§F ,--, '-u.
lF:o l-lo o o
4 +ltttl0 0 0l
I
0l0l
[r o o ol
Colocamos los opuestos de
trril Lt-ll-: I o t, I
i
los
multiplicadores debajo del segundo elemento I
diagonal de , -+ 1,, I
J
I
L, Y obtenemos
[--' -' ,,t '
I
1 0 00 23-1
Por consiguiente una faclonzacion de TI -
1100 03 2
-I
^^-4,d.
Strt
I \-L_lv-i
-
2 4 | 0 00 0
-1 101 00 0
ALGEBRA LINEAL Y TEORÍA MATRICIAL
MG. SC; lNG. RAFAEL VALENCIA GOÍAEÍA ING. JULIO UBERHUAG.AA
J
J
DEFINICION 33. FACTORIZACION DE PA = LU J
En eslos casc se
Frecuentemente, no es posibre escaronar una matriz sóro con operaciones de eliminación'
intercambio de rengloner.
p*ra este tipo Oe niátriás m ex'ste la fadorizacion LU' Lo qle ¿
requiere realizar p ,n" mát'Ut de permutación' Estas rrÚim de
o ¿
aplica es la factorizacion PA = LU. Donde la matriz fadorizacion PA = LU se oHierc de
permutación se obtienán de la matriz i¿entiOaO intercamgand" ú;
¡ñ: l¡
registro ¿* l"t nc q¡e s9 inte.rcambian
y se efrctiat hs J
forma análoga a la fadorizacion LU p"ro ,u llára un
intercambios en una matriz que registra los inversos
oe las operaciones de eliminación' J
Algoritmo de PA -- LU ,
Entrada:
ú
. Matriz Ade orden ¡¡xm ú
Salida: ú
. P matriz de permutación de orden ll x ¡l ,
J
. L matriz triangular superior unilaria de orden fixtl '
. IJ matriz escalonada de orden rlx r¡
Que cumPlen:
Se toma:
P=In
Mientras que U no sea escalonada hac f v
Aplicar una oPeración ;f de I
L; ¿
a Si ./ es de la forma .fl I
a Si'.f es de la forma É' =k;
, Pz =Pt
ITATRICIAL MG. SC.ING. R,AFAE!.VALEI{CTAGOTZUETA ING' JUUO UBEETIUAGA C'
ALGEBRA UNEAL Y TEORI'I
si aplicamos sobre U2 la operación para intercambiar la tercera y la cuarta filas se obtiene Ia nueva matriz
U3 . Se aplica la operación de intercambio en L2 \ P2:
lr 2 -2 ll
[ ; : ;l r;.-'r, . ;ll
loorol
=
"[i i lil I "+ 1lil Ll
Puesto que la matriz U: Ia es escalonada, el procedimbntotermina y finalizamos haciendo L=L3+I y se
tiene:
?
tt
g=u:=lo
_J
::":j:::"":
f, independientes, entonces A Puede
La forma de hace
Demostración
Sean d1,or,ct-i,....ün al ortonormalizarlas
según el proceso de Gram-
Definamos
Donde R es la malriz cuyo elemento (i,7) es r¡ . Las matrices buscadas son las matrices Q V R: Q tiene
sus columnas ortonormales yR es triangular superior. Asimismo R debe ser invertible pues en caso contrario
GOY¿IETA ING' JULIO UBERHUAGA C.
MATRICIAL MG' SC' ING' RAFAELVALENCIA
ALGEBRA LINEAL Y TEORIA
R.X=0tendríainfinitassolucionesypoffi=0contradiciendoelhechodeque
las columnas de A son línealmente independientes'
R = eT-A
r'
LasformulaSparaelprocesodeortonormalizaciondeGrarn-schmidtson: -
r',)ti
'i=ffi
''3=m tt1 - (tt¡,ri ) vr - (rl3,
de la matri
Ejemplo. Detennine una fadorizacion QR
v
'aS"''
v
.T
2 _t1\
« 43)-
Ahora se tiene
1 00 t_
tl \ 0 0 t 1 ll
o il 1 11
0
Jz Je I , ll o 1 ol 0
I
,!7 --=
--QT A= Ell o o rl=
0
q = lqrlqzl qt)= 2.
U
'Jz J2
0--
Jo 1 2 , l[ o i r] 3
0 --F
_11 ----: {o G Gl- ./6
Jz Jo ],
Ir -2 tl
Determine una factorizacion QRoelamatrize=l-r
:
'7l
It -i -4]
C'
GOYzrl TA ING. JULIO UBERñT¡AGA
TEOR¡A MATRIC¡AL
ilG. Sc. NG. RAFAEL VAII}E'A
ALGEBRA LINEAL Y
,r2={-2,3,-tf ,3 =(Le-+)r
,4 =1t"-Lt)r
a las columnas de A oütener¡c
Al aplicarle el proceso de Grarn-Schmidt
r\r
,,=[;,-+'+,) n'=lo'i'#) *=[á'+' _-
J6)
I
Por tanto
i I
1
0
zl {J
-Í a
+",1i -:.6
JEI J, -zJl
I
+[
J3
J
q =[erlazle:1=
{J
1
Í
I
#l
0
2
E
I
-l o+6
1
+I
1
J: Í «
,,§
Algoritmo QR- Para una malri
que
TomarAO=A'
Para i = 0,1, 2r. ,
a
de A.
Aplicar el algoritmo QR a la
:l: '