Trabajo de laboratorio: Segunda parte

1- Dado el siguiente sistema de ecuaciones lineales en forma matricial: A. X = B, siendo

5.11 -4 3.33 -1.11 3.77

-2.01 1.41 0.5 2.95 5.4
A= 1 0.333 1.5 -0.33 B= 3
4.321 -1.95 0 2.8 0.13

-0.146383964 -0.21347254 0.39612992 0.21356453 -0.488466152

A = -0.30153483 -0.14819079 0.71880425 0.12130878
-1
a) X= 0.235166312
0.168028999 0.22498389 0.21864766 -0.14465592 2.485520052
0.015904354 0.22622958 -0.11071752 0.11204993 0.964013054

5.1 -4.0 3.3 -1.1 3.77

b) A1= -2.0 1.4 0.5 3.0 B1= 5.4
1.0 0.3 1.5 -0.3 3
4.3 -2.0 0.0 2.8 0.13

-0.1 -0.2 0.4 0.2 -0.5

A-1 = -0.3 -0.1 0.7 0.1 b) X1= 0.2
0.2 0.2 0.2 -0.1 2.5
0.0 0.2 -0.1 0.1 1.0

c) Condicionamiento (µ)

µ= Normi(A)xNormi(A-1)

Para A: 13.55 Su máximo: 13.55 Normi (A) Para A-1

6.87
3.163
9.071

µ= 17.47731379

Problema 2.- Considere la matriz A =

-5 0 1.73205 1.73205 1
 0 -5 1.73205 1.73205 -1
A= 1.73205 1.73205 -5 0 B= 0
1.73205 1.73205 0 -5 1

b) Resuelva por Jacobi y por Gauss Seidel para a = -5 y b = 3 trabajando con cinco cifras decimales hasta diez iteraciones.

Por JACOBI:

PASO X1 X2 X3 X4 Distancia
0 0 0 0 0 0.2
1 -0.2000 0.2000 0.0000 -0.2000 0.069282
2 -0.2693 0.1307 0.0000 -0.2000 0.04799996
3 -0.2693 0.1307 -0.0480 -0.2480 0.03325533
4 -0.3025 0.0975 -0.0480 -0.2480 0.02303996
5 -0.3025 0.0975 -0.0710 -0.2710 0.01596254
6 -0.3185 0.0815 -0.0710 -0.2710 0.01105917
7 -0.3185 0.0815 -0.0821 -0.2821 0.00766201
8 -0.3262 0.0738 -0.0821 -0.2821 0.0053084
9 -0.3262 0.0738 -0.0874 -0.2874 0.00367776
10 -0.3298 0.0702 -0.0874 -0.2874 0.00254803
-0.3298 0.0702 -0.0900 -0.2900 0.00176532
-0.3316 0.0684 -0.0900 -0.2900

Por GAUSS SEIDEL:

PASO X1 X2 X3 X4 Distancia
0 0 0 0 0 0.2
1 -0.2000 0.2000 0.0000 -0.2000 0.069282
2 -0.2693 0.1307 0.0000 -0.2000 0.04799996
3 -0.2693 0.1307 -0.0480 -0.2480 0.03325533
4 -0.3025 0.0975 -0.0480 -0.2480 0.02303996
5 -0.3025 0.0975 -0.0710 -0.2710 0.01596254
6 -0.3185 0.0815 -0.0710 -0.2710 0.01105917
7 -0.3185 0.0815 -0.0821 -0.2821 0.00766201
8 -0.3262 0.0738 -0.0821 -0.2821 0.0053084
9 -0.3262 0.0738 -0.0874 -0.2874 0.00367776
10 -0.3298 0.0702 -0.0874 -0.2874 0.00254803
-0.3298 0.0702 -0.0900 -0.2900 0.00176532
-0.3316 0.0684 -0.0900 -0.2900
 Para el producto de matrices:
1) Seleccionar el rango de la matriz resultado
2) Función "MMULT"
3) Selecciono todo el rango de la 1° matriz, y lo mismo con la 2°
4) Pulsar Shift + Ctrl + Enter

0.96955095 Su máximo: 1.28983866 Normi A-1

1.28983866
0.75631647
0.46490138

Despeje de variables. Paso N° 1

X1= (1/-5)(1- 0*X2 -1,73205*X3 - 1,73205*X4)
 X2= (1/-5)*(-1 + 0*X1 - 1,7325*X3 - - 1,7325*X4)
X3= (1/-5)*(0 - 1,7325*X1 - 1,7325*X2 - 0*X3)
X4= (1/-5)*(1 - 1,7325*X1 - 1,7325*X2 - 0*X3)

hasta diez iteraciones.

Por GAUSS - SEIDEL:

0.2

PASO X1 X2 X3 X4 Distancia
0 0 0 0 0 0.2
1 -0.2000 0.2000 0.0693 -0.2000 0.069282
2 -0.2453 0.1547 0.0000 -0.2000 0.03137229
3 -0.2693 0.1307 -0.0314 -0.2314 0.02173535
4 -0.2910 0.1090 -0.0480 -0.2480 0.01505869
5 -0.3025 0.0975 -0.0631 -0.2631 0.01043296
6 -0.3130 0.0870 -0.0710 -0.2710 0.00722816
7 -0.3185 0.0815 -0.0783 -0.2783 0.00500782
8 -0.3235 0.0765 -0.0821 -0.2821 0.00346951
9 -0.3262 0.0738 -0.0856 -0.2856 0.00240375
10 -0.3286 0.0714 -0.0874 -0.2874 0.00166537
-0.3298 0.0702 -0.0891 -0.2891 0.0011538
-0.3310 0.0690 -0.0900 -0.2900
z, y lo mismo con la 2°