Scribd
Cargar
1K vistas4 páginas

Matlab Ejercicios Resueltos

Este documento presenta cuatro ejercicios resueltos en Matlab relacionados con sistemas de ecuaciones y matrices. El primer ejercicio resuelve un sistema de tres ecuaciones mediante el método de eliminación de Gauss. El segundo ejercicio resuelve otro sistema de siete ecuaciones usando el mismo método. El tercer ejercicio diseña un algoritmo para calcular la determinante de una matriz. El cuarto ejercicio diseña un algoritmo para calcular e imprimir una matriz triangular superior.

Cargado por

Yosimar Suca Cespedes
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
1K vistas4 páginas

Matlab Ejercicios Resueltos

Este documento presenta cuatro ejercicios resueltos en Matlab relacionados con sistemas de ecuaciones y matrices. El primer ejercicio resuelve un sistema de tres ecuaciones mediante el método de eliminación de Gauss. El segundo ejercicio resuelve otro sistema de siete ecuaciones usando el mismo método. El tercer ejercicio diseña un algoritmo para calcular la determinante de una matriz. El cuarto ejercicio diseña un algoritmo para calcular e imprimir una matriz triangular superior.

Cargado por

Yosimar Suca Cespedes
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

Matlab ejercicios resueltos

Matrices y sistemas de ecuaciones.


01 resuelva el sistema de ecuaciones

-1.41X +2Y = 1

X 1.41Y + Z= 1

22Y -1.41Z = 1

Haciendo uso de eliminacin de Gauss, con estrategia de pivoteo

Realizamos un algoritmo.

Comprobamos el algoritmo.

>> A=[-1.41 2 0;1 -1.41 1;0 22 -1.41]

A=

-1.4100 2.0000 0

1.0000 -1.4100 1.0000

0 22.0000 -1.4100

>> B=[1 1 1]'

B=

>> eliminacion(A,B)

A=

-1.4100 2.0000 0 1.0000

0 0.0084 1.0000 1.7092

0 22.0000 -1.4100 1.0000

A=

-1.4100 2.0000 0 1.0000

0 0.0084 1.0000 1.7092

0 22.0000 -1.4100 1.0000

A=

-1.4100 0 -236.9748 -404.0420


0 0.0084 1.0000 1.7092

0 22.0000 -1.4100 1.0000

A=

1.0e+003 *

-0.0014 0 -0.2370 -0.4040

0 0.0000 0.0010 0.0017

0 0 -2.6081 -4.4545

A=

1.0e+003 *

-0.0014 0 0 0.0007

0 0.0000 0.0010 0.0017

0 0 -2.6081 -4.4545

A=

1.0e+003 *

-0.0014 0 0 0.0007

0 0.0000 0 0.0000

0 0 -2.6081 -4.4545

x=

-0.4895

0.1549

1.7079

02 haciendo uso del mtodo de eliminacin de gauss con pivoteo resuelva


el sistema.

Solucin
>> A=[0 4 5 7 8 8 7;9 0 0.783 1 7 9 6;5 3 2 -1 -7 8 10;-1 2 3 0 4 5 3;11 10 9 4 2 7
6;15 11 8 7 6 9 3;12 1 3 6 5 7 10]

A=

0 4.0000 5.0000 7.0000 8.0000 8.0000 7.0000

9.0000 0 0.7830 1.0000 7.0000 9.0000 6.0000

5.0000 3.0000 2.0000 -1.0000 -7.0000 8.0000 10.0000

-1.0000 2.0000 3.0000 0 4.0000 5.0000 3.0000

11.0000 10.0000 9.0000 4.0000 2.0000 7.0000 6.0000

15.0000 11.0000 8.0000 7.0000 6.0000 9.0000 3.0000

12.0000 1.0000 3.0000 6.0000 5.0000 7.0000 10.0000

>> B=[-1 0 1 2 1 6 3]'

B=

-1
0

>> eliminacion(A,B)

A=

0 4 5 7 8 8 7 -1

NaN -Inf -Inf -Inf -Inf -Inf -Inf Inf

5 3 2 -1 -7 8 10 1

-1 2 3 0 4 5 3 2

11 10 9 4 2 7 6 1

15 11 8 7 6 9 3 6

12 1 3 6 5 7 10 3

El resultado nos bota.

x=

NaN

NaN

NaN

NaN

NaN

NaN

NaN

03 disee un algoritmo y programe en Matlab para calcular la


determinante de una matriz.
Solucin
Para el clculo de la determinante de una matriz en el programa Matlab
utilizamos el siguiente algoritmo.
function d = determinante(A)
n=length(A);
if n==1;
d=A(1,1);
else
d=0;
sgn=1;
for j=1:n
M1j=[A(2:n,1:j-1) A(2:n,j+1:n)];
d=d+sgn*A(1,j)*determinante(M1j);
sgn=-sgn;
end
end
end
ejemplo
>> A=[1 2 1 4;2 0 4 3;4 2 2 1;-3 1 3 2]

A=

1 2 1 4

2 0 4 3

4 2 2 1

-3 1 3 2

>> determinante(A)

ans =

-180

04 disee un algoritmo y programa en Matlab que permita


calcular e imprimir una matriz triangular superior.
Solucin
el siguiente algoritmo nos permite calcular e imprimir una matriz triangular superior
% Matriz del sistema
A=[1 3 2 2 1; 0 1 2 4 1; 0 0 2 1 -1; 0 0 0 3 -2; 0 0 0 0 2];
% Trmino independiente
b=[-2;-1;-2;-9;6];
% Hallamos el nmero de ecuaciones del sistema.
n=length(A);
% Definimos un vector-columna para almacenar la solucin.
x=zeros(n,1);
% El bucle va desde la ltima ecuacin hasta la primera
for i=n:-1:1
x(i)= (b(i)- A(i,i+1:n)*x(i+1:n))/A(i,i);
end
disp('La solucin del sistema es ');
disp(x);

También podría gustarte