- Primer Sistema de Ecuaciones
resuelto solamente por el Método de Gauss-
Seidel - Programa que resuelve un
Sistema de Ecuaciones únicamente por el Método de
Jacobi - Programa que resuelve un
Sistema de Ecuaciones únicamente por el Método de
Newton-Rapson - Programa que resuelve un
sistema de ecuaciones ya sea, por el método de Jacobi o
si el usuario lo desea, resuelto por el método de
Gauss-Seidel - Programa que resuelve la
raíz de un sistema de ecuaciones - Programa que resuelve una
ecuación en especifico por el método de
Jacobi
Introducción
En este archivo de
Texto les dejo
una serie de códigos hechos en Matlab, donde resuelve una
serie de cuestiones; como problemas
resueltos por los métodos
Gauss-Seidel, Jacobi, Newton-Rapson,
así como de un código
que les ayudara a resolver problemas de matrices,
etc., espero que les sirva estos códigos en Matlab que por
cierto me costo mucho
hacerlos así que espero me lo agradezcan.
Primer Sistema de
Ecuaciones resuelto solamente por el Método de Gauss-
Seidel
A)
function gauss
clear; A=[5 6 1; -5 5 3; 3 1 7]; b=[2 1
3]
XO=zeros(1,4); X=XO; K=0;
Norma=1;
fprintf (' K X(1) X(2) X(3)
Norman')
while Norma>0.0001
K=K+1; fprintf ('%2d', K)
for i=1:3
suma=0;
for j=1:3
if i ~= j
suma=suma+A(i,j)*X(j);
end
end
X(i)=(b(i)-suma)/A(i,i); fprintf
('%10.4f',X(i))
end
Norma=norm(XO-X);
fprintf('%10.4fn',Norma)
X0=X;
if K>20
break
end
end
B)
function gauss1
clear; A=[4.81 10.00 7.45; 3.33 4.81 -2.22;
-2.22 7.45 15.00]; b=[0 5 2]
XO=zeros(1,4); X=XO; K=0;
Norma=1;
fprintf (' K X(1) X(2) X(3)
Norman')
while Norma>0.0001
K=K+1; fprintf ('%2d', K)
for i=1:3
suma=0;
for j=1:3
if i ~= j
suma=suma+A(i,j)*X(j);
end
end
X(i)=(b(i)-suma)/A(i,i); fprintf
('%10.4f',X(i))
end
Norma=norm(XO-X);
fprintf('%10.4fn',Norma)
X0=X;
if K>13
break
end
end
C)
function gauss2
clear; A=[72 0 0 9 0 0; 0 2.88 0 0 0 -4.5;
0 0 18 9 0 0; 0 0 9 12 0 0; 0 0 0 0 33 0; 0 -4.5 0 0 0 33]; b=[2
0.5 1 0 1.2 5]
XO=zeros(1,6); X=XO; K=0;
Norma=1;
fprintf (' K X(1) X(2) X(3) X(4) X(5) X(6)
Norman')
while Norma>0.0001
K=K+1; fprintf ('%2d', K)
for i=1:6
suma=0;
for j=1:6
if i ~= j
suma=suma+A(i,j)*X(j);
end
end
X(i)=(b(i)-suma)/A(i,i); fprintf
('%10.4f',X(i))
end
Norma=norm(XO-X);
fprintf('%10.4fn',Norma)
X0=X;
if K>9
break
end
end
D)
function gauss3
clear; A=[4 -2 0; -2 4 -1; 0 -1 4]; b=[0
0.5 1]
XO=zeros(1,4); X=XO; K=0;
Norma=1;
fprintf (' K X(1) X(2) X(3)
Norman')
while Norma>0.0001
Página siguiente  |