Lista de exercícios resolvidos do livro "Matrix Computations" do golub

Prévia do material em texto

eolucae LlTA DASJEMANAS 
AL6cbRA LwEAE COm PUTA UOnJAL 
fotorad 8 Cbv 
u..1 C6oUb) Moko qua 
e ma vmavi 
Vettvek A¬ C 
_nxn nica 
e thangular ade 
ponna A- LD Ua 2 D2 Uz 
Tria ngulo- fpeier 
in fero- om D ina diagonal, 
U: e Tria ngulo 
w AS ha dioprox E dagnol 
onten o p.vos). 
Ainvert vel e Li e 
tomben a gR 
oo 
veiueis. e De 
de+ (Li) -1 det Lu)) Segue 
triangulo 
L L =(D) ( bu) 
4haSiaana L e tionguo 
to ferio wr 
no dogoraX Gnguonto co (Da e (D»U2) 
ogo do igual dade 
= 
teio- 
triagul 
re peiore. enos m Matia acina 
tn ongulo- i9 Peo e sporo COm noa 
hdode 
diaganaR, 
(D2U2)(DU) E 
DU = D2U2.e 
DU =(D.) (U) e (DUe)=): Cone 
de U = DaVz que D= Dz 
-/ / 
->u,1 lbotub) 
Alaoritmo totoraa PA=LDU 
nxn 
Entrata Ae C Cinvevhrval) 
faa 0=A, rernutaao (42 i h) 
fora 
Cooha we taj+1,, nt Tal 9ue liol max |Uejl 
Peruta ao lio) lerrutaao CG) = S, ia 
Paro =L ) 
FIM 
L -,n 
oro j+i, -- n: 
PlM 
FM 
Fw 
Fao 
se ij 
D 
Ce i tJ 
Foa ()iy . 
O 
Ternutacav LE) Faya 
Ca Lontraio 
nxn nxn 4.2. 1 (bdub) H- A+iB EC hermitiona ASe2 
eHpoh va -d ttni do 
Siome tica e poS va de Pioi ta Aatrar C A 
C A 
= A- Logo H H () e = -B. A -B SAT C 
A 
Lom 
Te en taw 
Aye 
ime tioa. Agoa * ()e"Aa ulo 
tiy C na -hwlo 
oL aiy" aryy =-i )\AAB ) (2+i)= 
= A -aby +B 
o(Az -2by+Bz tyA) 
-B 
( 
ao lao A 
'Az2By + s4' Bz +fAu 
C>oeQ ausna Segue2 Logo pela 
b Alaoritmo pora 
roelver (A+iB) ariq) =(btic), 
4bc e 
Atid) (a+iy) 
= (Ae -b) t i\ B tAy) 
equa lente Veng 
aeolve (A+iB) 
atig) 
= b+ic 
SA Az bu b 
Be+Ay=t 
A 
C 
tse hmo 
istemaneo n x Zn 
de 
enyolve yto Ee 
2n = 8 n perawe, aritmeHas 
Computa onol de 
8n Q iebio de C) POro 
de no mati n x de e Orrmazena mentto 
Simetnua , 40 nn +i) n lZn ) ena tas
reoi 
O A Fortoroui de nma Matniz 
xe Smetna usta n perawes, 
agro matro que LAr, 2> 2T zj bæcr" ) oma 
aro +ontoi Preua co euinte Lema vomo 
onis më trica (5' = -5). a+ ma e sE 
Dep do Lena) 
frinmeiro 
invivA 
Lom bre e se e = -s 
ex evwuo pOra 2.5 Golub, exoia Como Lexeuio 
Ontreor na 2L9ta). Alem deco kalem 
y e 
= L94 = -
Sy> = - Zy S4S 
Lop 
-)(L+)g 4> = 
Agora, ze , Taue 4= LTz u la deiguol Agoro, e 
cace Qima 
-) +s) z> = <(5'Jz-)2i 2 May 
= [F] -
2 (I-) (L+5)z, (L+) z eanalogamete
Lema enerra Odemngtra eu co 
Agorai Ouo TE SDP, tone T PE SDP tol qu 
P rem co deoo p eparal de T P= Ja qe T SPP 
TSDP). Daz, e , 
F(E+ + PoP) F"*z ) P P(z+ PSP) Pz, z 
(zPs P) Px, P> Pe Pe, Pz> 
Psimétrtca PsP athetime trica 
Cpois PT=P eSzs) 
+Lemo 
sea 
T 2 
Conuo 
ie n QMoe> 
5ja A-LDu fatotaa LDU de A om d:= Dei 
Lono Le 0a tnongulves LOm s na daganol LeU 
in rerhivas Dai AL = DUC) Lonuo Le E inguo 
nteno wn no Capiok, tom benm e4ong (nPen m L 
iong na dogno e torto () (eu) a 
Dacona Sagpnal 
DUC)) Do (C) ex er> 
= L'A Ae A C)ex ()ek 
LA > LTY Antr) > 
See ecrenC peoR 
Mas 4= =()'er, -eMO nha de L 
quuR te 'ho diagn 
e min ullz nin T) ndn T)= 
Pes T SDP e 
L2.4 Loalub). onoicere ( ma 
A ) Veomo 
On too G A) 
9e re ioy Se 
(z 2 = zgt 
a+ 
e Poa bvo 
("C)U) Aogra Vea gi reok. 
1-1 -i 
Ax >o, Pal de mat ð Conno 
A nA 
t la deinida A NAD 
Ae T 
4.2. S 
lbolub)Spoha A¬ R 
(po Leohas e po loluRs) 
eA 
esmitole agnalle 
dom 
T Gne ioa 
dom nonte 
reaiS. atova ore 
Ti= A +A" 
stritole diagonale 
doni nonte 
2 |Tl 
T Tii 
} lajlt 
2lail <4 (lactl 
tlatel)
J4 
Qt 
T Ta-. 
Se AEZ e 
autovalo- 
de 
A-7ic ) 2e -
J# 
Vi Ai 
= T 
= T 
Dau 
Agpra J= 
Esolha io 
il= lall 0 Dad 
lill IA- Te lTo 
aa algun 
A-7yie2 erio
- tior
auloy a lore T AoaM one tricu 
e 
oi bva- de 
Anda 
COm peetro, 
tognoi uni ta fa, 
Vaa Agora, 
Ax ( 
A > 
ZTR( 2, A+A , A 
A 
tua - de hnita 
a nu Mootra 4.2. 6 l6olub) Ae 
A> o e Uma nomaa )=vCAz,AS> 
de tAoi da, to tq ( Se A ha oese Potn 
Az > L0. Da ( Pu) hao eta detni do u f) zo 
Poma Na sena novma em , (n z+O). e quolu 
A Simetca 
E) Vomo f inica mente 
aPonawCa4) e" CAz m oto latevno Nte qR 
hn AaAT) Y LAx AA, ¥S )glz+ A )= 
LAr t<AYS ag+ Agl214) 
uiAR = 4AreyS 
LA mi tia 
y() éAga = A'a 
20 eeo 
Pory A o) 
YLAZAS n na 
on a e og0 
horma 
n Laso Plz) Depa 
T A+A 
AUAD ime tia, 
agora 
ZAz z. ogp Pelo Loongtico 
T e 
norma 
ie 
a)= VTI 
42.8 6otb) 
AEC A)= 
Ara lzll: = 1f 
Porto rauas LO. A tem 
Ae: A(4: k, l: e) 
,VL kin, 
Atenn 
Moso gue 
OhlA) 
Vemo motra Jua 
emtouA 
nveiye 
C4 
moptro que A tem 
LO) 
most 
LU) 
Moatra 
nrertivel 
emto A ato va lore de 
Ar 
VOniO 
toro tonto CA ) autopr 
Ac, u 
ua E C 
zE C°|Hot 
e Ae Ar. 
A 
= AS = A 
a hwlos: Ae =A2. e 
AE C 
de Ak, 
A- Az, 
wA) , Onde 
PotS 
1z 
= ( eC". 
Ovn-e 
Ae tem auto ra loe 
Ac verti vel
.2 (bolo) 
Ae& m men. Grtau estem 
mug Pe C 
herm.4i ana sem -be 
Kinica ptivon)e e 
om nhas atur 
PA Logo P VAA 
Teo Cx 
A-PO
tnnas 
mwm e tats q 
tal mat Psathotay 
Mem dpo 
min 
vnito 
ono re 
a SVD redz.da de 
A:A= U 2, V" cC 
removemo THi mas Lem wlnos de h m nxm 
VieC e ortenorvmus 
21 6 VtenColunas 
2,eV) 
Daz. Uni toia 
e 
v IeC C*) e trna 
mxn 
C Cureva A= U: z, 
(UU,) V," 
= \u 2, U,") 
(Cv,). A 
Careva 
Dad Pe 
vote que 
Pebermitiana 
damoR reoR) e 
vote q 
P= U , U" 
e O= U,v,". 
seu-do Pni -
e reoR) e dagonoR
a poo tvo Pz 2 = 
2 U1E,U." z =<2, ue Ua 
Gnquonto co u"= UV," 
V, U= U.U. =I 
ten Enkos 
po VV=i 
e e un.tafa, ue maNa ue Utenm 
enhu 
Prora a eaiLtenta. OrtonoYmats 
Vote q U 2, u" 
U. 2, ." U z, U, 
A Or2, V" 2, U= U, 2,"=> p=AAH 
AA hermhona 
e senmi de Hnida Poh Voo 
Lano 
matiz nerm hana
em de mdo PovO ue ao qvodrado AA" 
ee de tecrea peal). c.o. 
toimo po VAA 
AA 
cop toR Pe nico erpreaen- 
( 
O Va teorema auma, eMtemn m matig UeC 
hermi tona emi-dok ido Pto) 
n ee C onnas ortonoimac) 
A 0f Vesee caso no entono, 
(A demorstrayao anloga) 
.. lo (Golub 
n) 
A= A t ey 
A=+ 
An t 
(ni ca) un 
enontror 
fator de dho leky6 
de se 
A= 6 e 6 
A se1 
Voma
tnonlo i0Pelo 
AT A A O 
maro not 
g+UU4rS 
deoo 
Ay oray hole e A 
nn nz- 
tennos 
aloraio Oholesu 
+ 
z 
V 
terno: pz 
Vno (rp) 
=2,--+ 2p 
2)(2341) -2 
22 Voz2- 
Eruptalar, 99zz 
a+ + 
( )aa) zazz 
A A aTxno e,2,P 
vAlha 
) 
VOmay 
atror 2 
Cnpletava a 
2 
ee p, pti
inua nu P 
hipote de md ( note mesrumeté 
eno 
V -+Icc t 
4re)u 322 Ge 
ete 
V24 
pelo 
= V 2+ 
L + nd A) teruo 
(A) 
anaente 
ipte 
grlo ra2 
Re + 
Daw 
= Gee 
Chd o elo alyOvtmo e 
p+1- 
2 
p 2 pi + p 
89,du 
+ +) -z 
+p* 
(AK) 2+2) p+)-72p 
exezQ p+ 
(A) 4 G229+u p+ 
. 
2 
992p )92 + 25)eR.). 
zz 
9uS72 23 
+ p+ 
(ree) -2 ** 
le 
Uma indua pde rwr 
sio fai /mente 
O 
(4AR) 
ovnula poa p, pL. 
S0 P+2 
Hgora, 0 gn,- 
2 
-
9p 
(z)z: ) 
(AaA) Rp+ 
Cha inuua 
pse pracur 
isso ruclreme 
U p+ 
p+ 
77-pp p,p+l 
enertO 
ula 
44.1 (bolub) Ae 7 Halus abm atn m�tni Ca 
Ako A = O 
uy I1 
d bmatnae upo 
nera meite noto 
Sa a,. A +o 
P wntra tuu) 
Arma 
Lom aj to. 
Cono A = A' 
orda Ce eneralade 
. 
A onad 
pela 
Enha e 
colnus 
submtiz Prn upa 
D tomo 
A 
mat J 
controria 
a h potese 
kouo A=s 
4.4.5 Ccolb) A11, Azz > 
DF)n+P) Smeio 
A axn 2 nxp 
A 
Az n Azz p- Snetfia9 
Comuo 
A'= A, seg A 
e Azz he tfiay 
Az A 
AT= 
T 
AT Az Az -Azz 
orkra LDO 
SDP An LDL PxP Onuo A 
(D.) o B- Azz 
+ Aa (LT) 8" L"Az f Dow 
Ontow ennoy 
D.L A 
A nn nxn 
B AC)D pe 
Pan 
Eo -A + 
hzi{L) LA 
Agora, vga 
71Otra 
4 Az24, +LD L Az LAzy qe 
eD. ndo que pohva de tin.da, ppiy f2t ois 
tomo e mpre o 
futrasa) LDU e 
Da onoidere 
A erer e 
wn ) >o. Daw 
D LA2 
T 
D2L A 
LL 
LD, LA 
L 
D A L 
-D AL L 
atoraop LDL' de A CouoD e a 
(V 
e stav 
OnuorrO