ANALISIS

NUMERICO
II
TAREA NO. 3
CAMILA SOTO ROSALES

INGENIER
IA CIVIL MATEMATICA

Date: 22 de junio de 2011.

1

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

Ademas la matriz Q resultante es:

0,6667 0,3333 0,6667

Q = 0,6768 0,1041 0,7288
0,3123 0,9370 0,1562


R
sabemos que: QA =
0
Luego la matriz R resultante y el vector c1 es:




0,6667
3,0000 2,3333
,
c1 =
R=
0,6768
0
2,1344
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

Como Rx = c1 tiene solucion exacta

||Rx c1 ||22 + ||c2 ||22 ||c2 ||22 + ||Rx c1 ||22 ,
donde x es la solucion exacta.
As ||Ax b||2 es mnima para x = R1 c, donde

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

Ademas la matriz Q resultante es:

0,1387 0,0277 0,9707 0,1941

0,0277 0,1387 0,1941 0,9707

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

Como Rx = c1 tiene solucion exacta

||Rx c1 ||22 + ||c2 ||22 ||c2 ||22 + ||Rx c1 ||22 ,
donde x es la solucion exacta.
As ||Ax b||2 es mnima para x = R1 c, donde

x =

0,5590
0,1180

Luego calculamos:
||Ax b||2 = 8,9095

2 1 1
A = 3 2 2 ,
1
1 5

iii)

3,7417 1,8708 0,2673

0
1,5811 5,3759 ,
QA =
0
0
1,0142

1
b = 4
2

4,2762
Qb = 0,6325
1,5213

Ademas la matriz Q resultante es:

0,5345 0,8018 0,2673

0
0,3162 0,9487
Q=
0,8452
0,5071 0,1690


R
sabemos que: QA =
0
Luego la matriz R resultante y el vector c1 es:

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

Como Rx = c1 tiene solucion exacta

||Rx c1 ||22 + ||c2 ||22 ||c2 ||22 + ||Rx c1 ||22 ,
donde x es la solucion exacta.
As ||Ax b||2 es mnima para x = R1 c, donde

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

Ademas la matriz Q resultante

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

Como Rx = c1 tiene solucion exacta

||Rx c1 ||22 + ||c2 ||22 ||c2 ||22 + ||Rx c1 ||22 ,

donde x es la solucion exacta.

As ||Ax b||2 es mnima para x = R1 c, donde

x =

1,4832
0,9441

Luego calculamos:
||Ax b||2 = 0,5034

Problema 18. Se desea resolver el problema de aproximacion

m
X

(yi (0 0 (ti ) + 1 1 (ti ) + 2 2 ti )2 = mn

0 ,1 ,2

i=1

m
X
(yi (0 0 (ti ) + 1 1 (ti ) + 2 2 ti )2
i=1

para los datos

i 1
2
3 4 5 6 7
ti 3 2 1 0 1 2 3
yi 12 8
5 3 5 7 13
a) Resolver el problema para (i) i (t) = ti ,i = 0, 1, 2, (ii) i (t) = exp(it),i = 0, 1, 2,
(iii) 0 1, 1 (t) = sin(t), 2 (t) = sin(2t). En cada caso, transformar la matriz
A R73 a forma triangular superior mediante la transformacion de Householder.
b) Graficar el resultado.
c) En cada caso, calcular tambien la matriz A A. Comparar condkk2 (A A) y condkk2 (R).

i) Con los datos de la tabla, formulamos la matriz.

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

Luego procedemos a transformar la matriz a forma triangular mediante la trasnformacion

de Householder (usando el programa anexado anteriormente).

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

Ademas la matriz Q resultante es:

Q=

0,3780 0,3780 0,3780 0,3780 0,3780 0,3780 0,3780

0,5669 0,3780 0,1890
0
0,1890
0,3780
0,5669
0,5455
0
0,3273 0,4364 0,3273
0
0,5455
0,0276
0,1154 0,5946 0,7453 0,2361 0,1278 0,0702
0,0026
0,3043 0,5101 0,2512 0,7322 0,2075 0,0703
0,1600 0,4726 0,2447 0,1979 0,2652 0,6993 0,3043
0,4601 0,6205
0,2017 0,0949 0,2283 0,4072 0,3683

Luego la matriz R resultante y el vector c1 es:

2,6458
0
10,5830
,
0
5,2915
0
R=
0
0
9,1652

20,0321
c1 = 0,1890
9,0560

As la solucion al problema es:

3,6190
x = 0,0357
0,9881

Luego procedemos a calcular la matriz A A.

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

ii) Con los datos de la tabla, formulamos la matriz.

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

Luego procedemos a transformar la matriz a forma triangular mediante la trasnformacion

de Householder (usando el programa anexado anteriormente).

QA =

2,6458 11,9988 176,3477

0
17,9611
356,9766

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

Ademas la matriz Q resultante es:

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

Luego la matriz R resultante y el vector c1 es:

2,6458 11,9988 176,3477

0
17,9611
356,9766 ,
R=
0
0
85,2089

20,0321
c1 = 5,1547
3,4600

As la solucion al problema es:

7,2234
x = 0,5200
0,0406

Luego procedemos a calcular la matriz A A.

0,0001 0,0003 0,0047

A A = 1,0e + 005 0,0003 0,0047 0,0853
0,0047 0,0853 1,6579

Tambien calculamos:
condkk2 (A A) = 5,1065e + 004
condkk2 (R) = 225,9760

iii) Con los datos de la tabla, formulamos la matriz.

A=

1,0000 0,1411 0,2794

1,0000 0,9093 0,7568

1,0000 0,8415 0,9093

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

Luego procedemos a transformar la matriz a forma triangular mediante la trasnformacion

de Householder (usando el programa anexado anteriormente).

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

Ademas la matriz Q resultante es:

Q=

0,3780 0,3780 0,3780 0,3780 0,3780 0,3780

0,0800 0,5156 0,4772
0
0,4772
0,5156
0,1646
0,4531 0,5173
0
0,5173 0,4531
0,3719 0,0600 0,1283 0,8965 0,1019 0,1253
0,6405 0,0121
0,4900 0,1643 0,5421 0,1165
0,3911 0,6102 0,2401 0,1260 0,2202 0,5613
0,3496 0,0992 0,2183 0,1030 0,0955 0,2040

0,3780
0,0800
0,1646
0,1091
0,1229
0,1941
0,8712

Luego la matriz R resultante y el vector c1 es:

2,6458
0
0
0
1,7634 0,0426 ,
R=
0
0
1,7186

20,0321
c1 = 0,4356
0,2886

As la solucion al problema es:

7,5714
x = 0,2511
0,1679

Luego procedemos a calcular la matriz A A.

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

Problema 19. Resolver el problema de aproximacion.

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

para los datos

i
1
2
3 4 5 6
7 8
9
ti 2 3 5 0 1 4 3 1
7
yi 182 78 630 0 2 252 84 0 2408
Con los datos de la tabla, formulamos la matriz.

A=

1 2 4
8
16
32
1 3 9 27 81
243

1 5 25 125 625 3125

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

Luego procedemos a transformar la matriz a forma triangular mediante la trasnformacion

de Householder (usando el programa anexado anteriormente).

QA = 1,0e + 004

0,0003 0,0003 0,0038 0,0137 0,1154

0
0,0010
0,0028
0,0296
0,1488
0
0
0,0035 0,0154 0,1519
0
0
0
0,0093 0,0565
0
0
0
0
0,0251
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0

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

Luego la matriz R resultante y el vector c1 es:

0,0003 0,0003 0,0038 0,0137 0,1154 0,6313

0
0,0010
0,0028
0,0296
0,1488
1,1618

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

As la solucion al problema es:

x =

0
1,0000
0
0
1,0000
0