Analisis Numerico Trabajo

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 11

1) Resolver el siguiente sistema de ecuaciones lineales por el método gráfico y

determinantes compare resultados

{12 x +2 y= 23 23 x− 13 y=2}
Solución:
Notemos que despejamos la variable y en cada ecuación

{
1 2 2 1 1 1
x +2 y = ⟶ 2 y= − x ⟶ y = − x
2 3 3 2 3 4
2 1 1 2
x− y=2 ⟶− y =2− x ⟶ y =−6+2 x
3 3 3 3

1 1
y= − x Procedemos a tabular
3 4

x y
-2 0,833
-1 0,58
-0 0,33
1 0,083
76/27 -10/27
3 -0,4166
4 -0,666

y=−6+2 x Procedemos a tabular

x y
-2 -10
-1 -8
-0 -6
1 -4
76/27 -10/27
3 0
4 2
Se solucionará el sistema anterior por el método de determinantes.
El sistema está dado por:

{12 x +2 y= 23 23 x− 13 y=2}
Luego,
| | ( )( )
2
2
3
−1 2 −1 −2 −38
2 −( 2 )( 2 ) −4
3 3 3 9 9 76

| |
x= = = = =
1
2 ( )( ) ( )
2
1
2
1
3
− (2 )
2
3
−1 4 −3 27
6

3 2
2 −1
3 3

| | ( ) ( )( )
1 2
2 3
2 1 2 2 4 5
2 ( 2 )− 1−
3 2 3 3 9 9 −10

| |
y= = = = =
1
2 ( )( ) ( )
2
1
2
1
3
−( 2 )
2 −1 − 4 −3 27
3 6 3 2
2 −1
3 3

76 10
La cual es la misma solución que se encontró con el método gráfico ( ,− ) .
27 27

2) Aplicar la regla de cramer para demostrar que x=1 , y=2 , z =3.

{ 2 y−z=1 ;−x+ y =1; x− y+ z =3

{
2 y−z=1
−x + y =1
x− y + z=3

El determinante de la matriz es entonces:

| |
0 2 −1
∆= −1 1
1 −1 1
0 =( 0 )
−1 1
−2 |
1 0 ( ) −1 0
1 1 | |
+(−1)
−1 1
1 −1 | | |
¿ ( 0 ) (1−0)−( 2 ) (−1+0)+(−1)(1−1)
¿ 0−( 2 )(−1 ) + (−1 ) ( 0 )
¿ 2 ≠0

| |
1 2 −1
∆ x= 1 1 0 =( 1 )
3 −1 1
−1 1 |
1 0 ( )1 0
−2
3 1| | |
+(−1)
1 1
3 −1 | |
¿ ( 1 ) (1−0)−( 2 ) (−1+0)+(−1)(−1−3)

¿ 1−( 2 ) ( 1 ) + (−1 )(−4 )


¿ 1−2+ 4=3

| |
0 1 −1
∆ y = −1 1 0 =( 0 )
1 3 1
3 1 | | |
1 0 ( ) −1 0
−1
1 1
+(−1) |
−1 1
1 3 | |
¿ ( 0 ) (1−0)−( 1 ) (−1+0)+(−1)(−3−1)

¿( 0)(1)− (1 ) (−1 ) + (−1 ) (−4 )


¿ 1+4=5

| |
0 2 1
∆ z = −1 1 1 =( 0 )
1 −1 3
|
−1 3
−2 | |
1 1 ( ) −1 1
1 3
+(−1) |
−1 1
1 −1 | |
¿ ( 0 ) (3+ 1)− ( 2 )(−3−1)+(−1)(1−1)

¿ 0−( 2 )(−4 )+ (−1 )( 0 )


¿8

| |
1 2 −1
1 1 0
∆x 3 −1 1 3
x= = =

| |
∆ 0 2 −1 2
−1 1 0
1 −1 1

| |
0 1 −1
−1 1 0
∆y 1 3 1 5
y= = =

| |
∆ 0 2 −1 2
−1 1 0
1 −1 1

| |
0 2 1
−1 1 1
∆z 1 −1 3
8
z= = = =4

| |
∆ 0 2 −1 2
−1 1 0
1 −1 1

3) Mediante la eliminación de incógnitas resolver

{ 4 x+ 4 y=−10; 8 x−5 y =5
Notemos que { a 11 x+ a12 y =b1 ; a21 x +a 22 y=b 2

(−5 ) (−10 )−(4)(5) 50−30 30 −15


Para la solución x 1= = = =
( 4 ) (−5 )−(4 )(8) −20−32 −52 26

( 4 )( 5 ) −(8)(−10) 20+80 100 −25


Para la solución x 2= ( 4 ) (−5 )−(4 )(8) = −20−32 = −52 = 13

4) Considere el sistema propuesto

A=[ 1−3 5−3−27 8 4 9 ]

Encuentre la solución del sistema asociado aplicando factorización


LU

Notemos que lo podemos expresar de la siguiente manera

( )
1 −3 5
A= −3 −2 7
8 4 9

Luego

( )( )( )
1 −3 5 1 0 0 1 −3 5
A= −3 −2 7 = a 1 0 0 d e
8 4 9 b c 1 0 0 f

Luego
( )( )( )
1 0 0 1 −3 5 1 −3 5
a 1 0 0 d e = a d−3 a 5 a+ e
b c 1 0 0 f b cd −3 b 5 b +ec +f

( )( )
1 −3 5 1 −3 5
a d−3 a 5 a+ e = −3 −2 7
b cd −3 b 5 b+ ec +f 8 4 9

Notemos que igualando miembro a miembro tenemos que:

a=−3

d−3 a=−2 ⟶ d−3 (−3 )=−2 ⟶ d+ 9=−2 ⟶ d=−11

5 a+ e=7 ⟶ 5(−3)+e=7 ⟶ −15+e=7 ⟶ e=22

Ahora bien, dado que tenemos que

b=8
−28
cd −3 b=4 ⟶ −11 c−24=4 ⟶ c=
11

−28
5 b+ ec+ f =9 ⟶ 5 ( 8 ) +22( )+ f =9 ⟶ f =25
11

Entonces

( )
1 0 0

( )
1 −3 5
−3 1 0
L= U = 0 −11 22
−28
8 1 0 0 25
11

Por lo tanto
)( )(
1 0 0

( )
1 −3 5 1 −3 5
−3 1 0
A= −3 −2 7 = 0 −11 22
−28
8 4 9 8 1 0 0 25
11
La factorización LU.

5) Utilice el método de jacobi para encontrar, a partir de cuatro


iteraciones, la solución a la matriz asociada a un sistema de
ecuaciones lineales dada por:

A=[ 0,20,1−0,5 0,7 ] , b=[ 0,75 0,35 ] x 0=[ 1 1 ]

Luego

{−0,5
0,2 x+ 0,1 y =0,75
x+ 0,7 y=0,35

Veamos que la matriz es diagonal estrictamente dominante

(−0,5 0,7 )
0,2 0,1
⟶ |0,2|>|0.1|verdad , ademas|0,7|>|−0,5|

Luego tenemos que:

(−0,5 0,7 ) ( 0
0,2 0,1
=
0,2 0
)(
+
0 0 0 0,1
+
0,7 −0,5 0 0 0 )( )
Donde se tiene que:

B= (0,20 0,70 ) , L=(−0,5


0
0)
0
, U =(0
0 0,1
0 )
Procedemos hallar la inversa de B
( )
5 0
−1
B =
0
10 Además, note C=L+U
7

( ) [(
5 0
T =−B C=
−1
0
10
7
0
+)(
0 0 0,1
−0,5 0 0 0
=)] (0
0,71428
−0,5
0 )
T= ( 0,71428
0 −0,5
0 )
Ahora se determina C de la siguiente forma:

()
15

( )(
5 0
C=B b=
−1
0
10
7
0,75
0,35
=
4
1 )
2

Entonces, para la primera iteración se tiene que:

( )( )
15 13
1
x= (
0
0,71428
−0,5 1
0 1
+
4
1
=
4
30357 )( )
2 25000

( )( ) (
13 15
2
x= (
0
0,71428
−0,5
0
4
30357
+
4
1
= )
3,14286
2.82141 )
25000 2

)( )(
15
3
x= (
0
0,71428
−0,5
0
3,14286
2.82141 1 )(
+ 4 = 2,339295
2.74488 )
2
( )(
15
x = 4 0
(
0,71428
−0,5
0 )(
2,339295
2.74488 )
+ 4 = 2,37756
1 2.17091 )
2

6) Aplicando el método de Gauss-seidel, resolver :

{2 x 1 +2 x 2−3 x3 =2−4 x 2−4 x 3=6 x 1 +3 x 2−10 x3 =6

Notemos que:

{
2 x 1 +2 x 2−3 x 3=2
−4 x 2−4 x 3=6
x 1 +3 x2 −10 x 3=6

Donde,

3
2 x1 +2 x 2−3 x3 =2 ⟶ 2 x1 =2−2 x2 +3 x 3 ⟶ x 1=1−x 2 + x
2 3

−3
−4 x 2−4 x 3=6⟶ −4 x 2=6+ 4 x 3 ⟶ x 2= −x 3
2

−3 +1 3
x 1+ 3 x 2−10 x 3=6 ⟶ −10 x 3=6−x 1−3 x 2 ⟶ x 3= x + x2
5 10 1 10

Lego tenemos que si x 3=0


−3 −3 −3
x 2= −x 3 ⟶ x 2= −0 ⟶ x 2= ≈−1,5
2 2 2

Lego tenemos que si x 1=0


x 3=
−3 +1 3
x + x 2 ⟶ x 3=
5 10 1 10
−3 1
+ ( 0) +
5 10 ( )
3 −3 −21
10 2
=
10
≈−1,05

Lego tenemos que si x 2=0

3
x 1=1−x 2 + x ⟶ x 1=1−( 0 )+
2 3 ( )
3 −21 −43
2 10
=
20
≈−2,15

Tenemos en vector

[ ]
x1 =−2,15
x2 =−1,5
x3 =−1,05

También podría gustarte