AN Alisis Num Erico Ii Tarea No. 3: Camila Soto Rosales Ingenier Ia Civil Matem Atica
AN Alisis Num Erico Ii Tarea No. 3: Camila Soto Rosales Ingenier Ia Civil Matem Atica
AN Alisis Num Erico Ii Tarea No. 3: Camila Soto Rosales Ingenier Ia Civil Matem Atica
NUMERICO
II
TAREA NO. 3
CAMILA SOTO ROSALES
INGENIER
IA CIVIL MATEMATICA
Problema 17. Para cada uno de los pares (A,b), usar el metodo de Householder para determinar una descomposicion QR de A, determinar un vector x = x que minimiza
||Ax b||2 , y calcular ||Ax b||2 :
Para la matriz:
i)
2 3
A = 1 1 ,
2
9
1
b= 0
0
Por medio del programa, realizado en Matlab, que se adjunta calculamos la descomposicion
QA y Qb.
Listing 1. Descomposicion QR mediante Householder
function [ A, b]= h o u s e h o l d e r (A, b )
[m, n]= s i z e (A ) ;
[ s1 , s 2 ]= s i z e ( b ) ;
bb=b ;
Q=eye (m) ;
for i =1:n
be =1/((norm(A( i :m, i ) , 2 ) ( abs (A( i , i ))+norm(A( i :m, i ) , 2 ) ) ) ) ;
w=[ sgn (A( i , i ) ) ( abs (A( i , i ))+norm(A( i :m, i ) ) ) ; A( i +1:m, i ) ] ;
P=eye (m+1 i )be ww ;
Pe=eye (m) ;
Pe ( i :m, i :m)=P ;
A( i :m, i : n)=PA( i :m, i : n ) ;
b ( i : s1 )=Pb ( i : s1 ) ;
Q=PeQ;
end
Q
R=A( 1 : n , 1 : n )
c1=b ( 1 : n , 1 )
x=inv (R) c1
min=norm( (Axb ) , 2 )
end
function a=sgn ( x )
i f x<0
a=1;
else
a =1;
end
3 2,3333
QA = 0 2,1344 ,
0
0
0,6667
Qb = 0,6768
0,3123
b||22
2
R
c
1
=
= ||Q(Ax
x
0
c2 2
Rx c1 2
= ||Rx c1 ||22 + ||c2 ||22
=
c2
2
b)||22
x =
0,0244
0,3171
Luego calculamos:
||Ax b||2 = 0,3123
ii)
10 5
2 11
A=
70 35 ,
14 77
5
9
b=
35
0
72,1110 49,9230
0
69,3375
,
QA =
0
0
0
0
34,4191
8,1818
Qb =
1,3199
8,8112
Q=
0,9790 0,1467 0,1399
0,0210
0,1467 0,9790 0,0210 0,1399
R
sabemos que: QA =
0
Luego la matriz R resultante y el vector c1 es:
72,1110 49,9230
R=
,
0
69,3375
c1 =
34,4191
8,1818
Observamos que:
||Ax
b||22
2
R
c
1
= ||Q(Ax
=
x
0
c2 2
Rx c1 2
= ||Rx c1 ||22 + ||c2 ||22
=
c2
2
b)||22
x =
0,5590
0,1180
Luego calculamos:
||Ax b||2 = 8,9095
2 1 1
A = 3 2 2 ,
1
1 5
iii)
1
b = 4
2
4,2762
Qb = 0,6325
1,5213
4,2762
3,7417 1,8708 0,2673
0
1,5811 5,3759 ,
c1 = 0,6325
R=
1,5213
0
0
1,0142
Observamos que:
||Ax
b||22
2
R
c1
=
= ||Q(Ax
x
0
c2 2
Rx c1 2
= ||Rx c1 ||22 + ||c2 ||22
=
c2
2
b)||22
4,0000
x = 5,5000
1,5000
Luego calculamos:
||Ax b||2 = 2,2204e 016
iv)
QA =
A=
1 1
2 4
3 9
4 16
5 25
7,4162 30,3390
0
7,6515
,
0
0
0
0
0
0
0,1348
0,4040
Q=
0,4308
0,5260
0,5970
R
sabemos que: QA =
0
Luego la matriz R resultante y
7,4162
R=
0
b=
2
7
13
21
31
Qb =
39,6430
7,2238
0,1356
0,2003
0,4415
||Ax
es:
0,2697
0,5465
0,3861
0,0195
0,6922
0,4045
0,4277
0,7909
0,1592
0,0497
el vector c1 es:
30,3390
,
7,6515
0,5394 0,6742
0,0475 0,5941
0,1765 0,0928
0,7356 0,3956
0,3668 0,1658
c1 =
39,6430
7,2238
Observamos que:
b||22
2
R
c1
= ||Q(Ax
=
x
0
c2 2
Rx c1 2
= ||Rx c1 ||22 + ||c2 ||22
=
c2
2
b)||22
x =
1,4832
0,9441
Luego calculamos:
||Ax b||2 = 0,5034
0 ,1 ,2
i=1
m
X
(yi (0 0 (ti ) + 1 1 (ti ) + 2 2 ti )2
i=1
A=
1 3 9
1 2 4
1 1 1
1 0 0 ,
1 1 1
1 2 4
1 3 9
b=
12
8
5
3
5
7
13
QA =
2,6458
0
10,5830
0
5,2915
0
0
0
9,1652
0
0
0
,
0
0
0
0
0
0
0
0
0
Qb =
20,0321
0,1890
9,0560
0,6451
0,4561
0,3415
0,9622
Q=
2,6458
0
10,5830
,
0
5,2915
0
R=
0
0
9,1652
20,0321
c1 = 0,1890
9,0560
3,6190
x = 0,0357
0,9881
7,0000
0
28,0000
0
28,0000
0
A A =
28,0000
0
196,0000
Tambien calculamos:
condkk2 (A A) = 68,0686
condkk2 (R) = 8,2504
A=
1,0000 0,0498
0,0025
1,0000 0,1353
0,0183
1,0000 0,3679
0,1353
1,0000 1,0000
1,0000
1,0000 2,7183
7,3891
1,0000 7,3891 54,5982
1,0000 20,0855 403,4288
b=
12
8
5
3
5
7
13
QA =
0
0
85,2089
0
0
0
,
0
0
0
0
0
0
0
0
0
Qb =
20,0321
5,1547
3,4600
5,1212
3,4621
2,5483
0,6816
Q=
0,3780
0,2497
0,2640
0,3672
0,4050
0,4771
0,4456
0,3780
0,2450
0,2442
0,2288
0,1411
0,1030
0,8091
0,3780
0,2320
0,1914
0,0814
0,2337
0,7480
0,3820
0,3780
0,1968
0,0541
0,8840
0,1262
0,1340
0,0075
0,3780
0,1012
0,2717
0,1161
0,7769
0,3940
0,0276
0,3780
0,1589
0,8072
0,1046
0,3750
0,1698
0,0014
0,3780
0,8658
0,3252
0,0141
0,0368
0,0157
0,0002
20,0321
c1 = 5,1547
3,4600
7,2234
x = 0,5200
0,0406
Tambien calculamos:
condkk2 (A A) = 5,1065e + 004
condkk2 (R) = 225,9760
A=
1,0000
0
0
,
1,0000 0,8415
0,9093
1,0000 0,9093 0,7568
1,0000 0,1411 0,2794
b=
12
8
5
3
5
7
13
QA =
2,6458
0
0
0
1,7634 0,0426
0
0
1,7186
0
0
0
,
0
0
0
0
0
0
0
0
0
Qb =
20,0321
0,4356
0,2886
5,6992
5,3346
1,0856
4,6177
Q=
0,3780
0,0800
0,1646
0,1091
0,1229
0,1941
0,8712
2,6458
0
0
0
1,7634 0,0426 ,
R=
0
0
1,7186
20,0321
c1 = 0,4356
0,2886
7,5714
x = 0,2511
0,1679
7,0000
0
0
0
3,1096 0,0751
A A =
0
0,0751 2,9553
Tambien calculamos:
condkk2 (A A) = 2,3934
condkk2 (R) = 1,5470
m
X
i=1
yi
5
X
l=0
!2
l tli
= mn
0 ,...,5
m
X
i=1
yi
5
X
!2
l tli
l=0
A=
1 2 4
8
16
32
1 3 9 27 81
243
1 0 0
0
0
0
1 1 1
1
1
1
,
1 4 16 64 256 1024
1 3 9 27
81
243
1 1 1 1
1
1
1 7 49 343 2401 16807
b=
18
78
630
0
2
252
84
0
2408
QA = 1,0e + 004
Qb = 1,0e + 003
1,1573
1,4981
1,5185
0,5655
0,2507
0
0
0
0
0,6313
1,1618
0,9099
0,5662
0,1807
0,0673
0
0
0
0
0
0,0035 0,0154 0,1519 0,9099
R = 1,0e + 004
0
0
0
0,0093 0,0565 0,5662
0
0
0
0
0,0251 0,1807
0
0
0
0
0
0,0673
c1 = 1,0e + 003
1,1573
1,4981
1,5185
0,5655
0,2507
0
x =
0
1,0000
0
0
1,0000
0