Deber3 Ushiña Alex 3257

METODOS
NUMERICOS
ALUMNO: Alex Ushi

na
FECHA: 16-Enero-2017
NRC: 3257
DEBER No 2
SISTEMAS DE ECUACIONES LINEALES Y AJUSTE DE CURVAS
1. EJERCICIO: En el sistema siguiente, pruebe que Ax = B es equivalente al sistema triangular
superior U x = y que se da y halle su solucion.
4x1 + 8x2 + 4x3 + 0x4 = 8
x1 + 5x2 + 4x3 3x4 = 4

x + 4x2 + 7x3 + 2x4 = 10
1
x1 + 3x2 + 0x3 2x4 = 4
4x1 + 8x2 + 4x3 + 0x4 = 8
3x2 + 3x3 3x4 = 6

4x3 + 4x4 = 12
x4 = 2
SOLUCION:
Se desarrollo un programa para resolver el ejercicio, cuyo codigo es:
function x=sistem_gauss(A,B)
[n1 n2]=size(A);
if ni~=n2
disp(El sistema no tiene solucion unica);
end
n=n1;
Ab=[A;B];
1
for k=1:n
[LL bb]=max(abs(Ab(k:n,k)));
if LL==0
disp(El sistema no tiene singular)
end
mm=k+bb-1;
Ab=intercambio_filas(Ab,k,j,-Ab(j,k)/Ab(k,k));
disp(Ab);
end
end
fprintf(la matriz triangular superior es:\n);
disp(Ab);
x=triangular_superior3(Ab(:,1:n),Ab(:,n+1));
end
La solucion que nos arroja el programa es:
Donde se cumple Ax = B es equivalente al sistema triangular superior U x = y.

2. EJERCICIO: Resolver con calculadora (a mano) el siguiente sistema de ecuaciones, aplicando
la eliminaci
on gaussiana:
x1 + 8x2 5x3 = 3
3x1 2x2 + 3x3 = 1
2x1 + 3x2 x3 = 4
SOLUCION:
Como podemos observar al momento de realizar la eliminacion gausiana tenemos una inconsistencia 0 = 2, por lo que el podemos concluir que el sistema no tiene solucion.
3. EJERCICIO:Halle la soluci
on del siguiente sistema lineal, con calculadora y a mano,aplicando
Gauss Seidel.
x1 + x2 = 5
2x1 x2 + 5x3 = 9
3x2 4x3 + 2x4 = 19
2x3 + 6x4 = 2
SOLUCION:
Una de las condiciones principales para poder realizar el metodo de Gauss Seidel es que la matriz sea diagonalmente dominante. En este ejercicio se puede observar que la matriz A no es
diagonalmente dominante, por lo que el metodo no se puede aplicar ya que como se puede ver
hasta la iteraci
on 6 no se logro establecer los valores.
POR GAUSS:
RESOLUCION
1 1
0 0 5
2 1 5 0 9
0 3 4 2 19
0 0
2 6 2
Realizamos el intercambio de filas: F2 F1
2 1 5 0 9
1 1
0 0 5
0 3 4 2 19
0 0
2 6 2
M21 = 21
Realizamos la operaci
on: F2 -M21 F1
2 1 5
0 3 5
2
2
0 3 4
0 0
2
0 9
19
0 2
2 19
6 2
2 1 5 0 9
0 3 4 2 19
5
19
3
0
2 0 2
0 0
2 6 2
M32 = 12
Realizamos las operaciones: F3 -M32 F2
2 1 5
0 9
0 3 4 2 19
0 0 1 1 0
0 0
2
6
2
2 1 5
0 9
0 3 4 2 19
0 0
2
6
2
0 0 21 1 0
M43 =- 41
Realizamos la operaci
on: F4 -M43 F3
2 1 5 0 9
0 3 4 2 19
0 0
2 6 2
0 0
0 12 21
Formando el sistema de ecuaciones:

2x1 -x2 +5x3 =-9
3x2 -4x3 +2x4 =19
2x3 +6x4 =2
1
1
2 x4 = 2
Soluciones:
x1 = 2, x2 = 3, x3 = 2, x4 = 1
4. EJERCICIO:Demuestre que la inversa de una matriz triangular superior es una matriz triangular superior.
SOLUCION:
Sea la matriz triangular superior:
1 2 4
0 3 5
0 0 4
7
La matriz inversa adjuntando la matriz identidad es:
Por lo que se demuestra que la inversa de cualquier matriz triangular superior es otra matriz
triangular superior.
5. EJERCICIO:Resolver el sistema lineal de ecuaciones Ax = B, calculando a mano solo las
matrices correspondientes a la factorizacion P A = LU .
2 3 8
1
4
0
1 10
A=
16 4 2
1
0
7 1
5
1

1

b=

1

1
SOLUCION:
2 3 8
0
1
0
0
1 10
(1)
A =
P1 =
16 4 2
1
1
0
7 1
5
0
0 1 0
1 0 0
0 0 0
0 0 1
16 4 2
1
4
0
1 10
A(1) =P1 A(1) =
2 3 8
1
0
7 1
5
4
M21 = 16
2
M31 = 16
0
M41 = 16
16
(2)
(1)
A =L1 A =
2
16
0
16
(2)
(2)
A =P2 A =
0 0 0
1
4
2
1
16 4 2
16
1 0 0 4
0
1 10 0
1
3/2 41/4
0 1 0
1
7/8
2 3 8
0 7/2 33/4
0 0 1
0
7 1
5
0
7
1
5
0 0 0
4
2
1
4
2
1
16
16
0 0 1
1
3/2 41/4
7
1
5
0
0
1 0 0
7/8
1
3/2 41/4
0 7/2 33/4
0
0 1 0
0
7
1
5
0 7/2 33/4
7/8
M32 = 1
7
M42 = 7/2
7
1 0 0 016
4
2
1
16

0 1 0 0 0
0
7
1
5

(3)
(2)
=
A =L2 A =

0 1 1 0 0
1
3/2 41/4
0
7

0 12 0 1
0 7/2 33/4
7/8
0
(3)
(3)
A =P3 A =
7
1
5
0 19/14 267/28
0 31/4
27/8
0 0 0
16
1 0 0
0
0 0 1
0
0
0 1 0
16

0
7
1
5

=

0 19/14 267/28
0

0
0 31/4
27/8
4
38
M43 = 217
10
7
1
5
0 31/4
27/8
0 19/14 267/28
(4)
(3)
A =L3 A =
38
0 217
16
U=
7 1
5
0 31/4 27/8
0
0
10,12
P=P1 P2 P3
16

0
7
1
5

=

0 31/4
27/8
0

0 19/14 267/28
0
1
0
0
1/4
1
0
L=
1/8 1/7
1
0 1/2 38/217
1 1
L=L1
1 L2 L3
P=
0
16
0
0
0
0
1
0
1 0 0
0 1 0
0 0 0
0 0 1
2 3 8
1
4
0
1 10
A=
16 4 2
1
0
7 1
5
11
7 1
5
0 31/4 27/8
0
0
10,12
PA=LU
1 0 0
1
0
0
2 3 8
1
0 1 0 4
0
1 10
1
0
1/4
0 0 0
1
1
16 4 2
1/8 1/7

0 0 1
0
7 1
5
0 1/2 38/217
0
16
0
0
0
0
1
0
7 1
5
0 31/4 27/8
0
0
10,12
2 3 8
1
1
16 4 2
0
1 10
7 1
5
16 4 2
1
2
3
8
1
0
7 1
5
4
0
1 10
PROGRAMA:
>>SistemaLUPA(A,B)
La solucion del sistema es:
La solucion del sistema es:
ans=
0.0377 0.2182 0.2055 -0.0644
6. EJERCICIO: Considere la funci

on f (x) = x2 .ex . Se pide calcular un valor aproximado para
la integral de f (x) en el intervalo [2, 2] usando el polinomio de Lagrange, calculado a mano,
que interpola f (x) en los puntos x0 = 2, x1 = 1, x2 = 0, x3 = 1, y, x4 = 2.
SOLUCION:
12
7. EJERCICIO: Con el siguiente conjunto de nodos:
Obtener el valor de la funci

on para x = 90, con un polinomio de 2do grado, utilizando los
siguientes metodos:
a) Por interpolaci
on de Lagrange. (Aplicando el programa)
b) Construya solo la matriz de diferencias divididas para aproximar todos los puntos de la
tabla.
13
c) Evaluar el polinomio interpolador de Newton, de tercer grado, para x = 1,75.
SOLUCION:
a) Por interpolaci
on de Lagrange. (Aplicando el programa)
Como requerimos un polinomio de segundo grado tomaremos los tres puntos mas cercanos
al punto del cual queremos conocer el valor, en este caso son x = 40, x = 60 y x = 80
La ecuacion entonces es:

0,0001x2
+ 0,0252x 0,56
Y para x = 90, y = 2,62
b) Construya solo la matriz de diferencias divididas para aproximar todos los puntos de la
tabla.
c) Evaluar el polinomio interpolador de Newton, de tercer grado, para x = 1,75.

Para el polinomio interpolador de grado 3, tomaremos cuatro puntos de la tabla, entonces:
Podemos observar que evaluando el polinomio para x = 1,75 se obtiene y = 0,18857

8. EJERCICIO: La ecuaci
on x 9x = 0 tiene una solucion en el intervalo [0; 1]. Utilice la teora
de interpolaci
on polinomial en los nodos x0 = 0, x1 = 0,5, x2 = 1, para encontrar una soluci
on
aproximada x
de la ecuaci
on. (Aplicar los programas)
14

SOLUCION:
Para la solucion de este problema utilizaremos el siguiente codigo
function C=newtonInterpolacion2(x,f)
m=length(x);
A=zeros(m,n);
A(:,1)=f;
for j=2:m
for k=j:m
A(k,j)=(A(k,j-1)-A(k-1,j-1))/(x(k)-x(k-j+1));
end
end
C=A(m,m);
for k=m-1:-1:1;
C=conv(C,poly(x(k)));
mm=length(C);
C(mm)=C(mm)+A(k,k);
end
end
Corriendo el programa obtenemos:
Polinomio:
P (x) = 0,8889x2 + 2,778x 1
x=
2,778+
2,77782 4(0,8889)(1)
2(0,8889)
x1 =0.4151
x2 =2.7101
15

Deber3 Ushiña Alex 3257

Cargado por

Copyright:

Formatos disponibles

Deber3 Ushiña Alex 3257

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Deber3 Ushiña Alex 3257

Cargado por

Copyright:

Formatos disponibles

METODOS

ALUMNO: Alex Ushi

4x1 + 8x2 + 4x3 + 0x4 = 8

x1 + 5x2 + 4x3 3x4 = 4

4x1 + 8x2 + 4x3 + 0x4 = 8

3x2 + 3x3 3x4 = 6

Donde se cumple Ax = B es equivalente al sistema triangular superior U x = y.

Realizamos el intercambio de filas: F3 F2

Formando el sistema de ecuaciones:

La matriz inversa adjuntando la matriz identidad es:

A(1) =P1 A(1) =

6. EJERCICIO: Considere la funci

7. EJERCICIO: Con el siguiente conjunto de nodos:

Obtener el valor de la funci

c) Evaluar el polinomio interpolador de Newton, de tercer grado, para x = 1,75.

La ecuacion entonces es:

c) Evaluar el polinomio interpolador de Newton, de tercer grado, para x = 1,75.

Podemos observar que evaluando el polinomio para x = 1,75 se obtiene y = 0,18857

También podría gustarte