viernes, 4 de diciembre de 2009

Programa 3.1 (Sustitucion regresiva - pagina 135).

%Programa 3.1 (Sustitucion regresiva - pagina 135).
% Resolucion de un sistema triangular superior AX=B por el metodo de
% sustitucion regresiva. El metodo funciona solo si todos los elementos
% diagonales son distintos de cero. Primero se calcula X=b / a
% y luego se usa la regla
% X=( b - sumatoria(a·X) )/a ,para j=k+1,...,N y k=N-1,N-2,...,1.
%
%Sintaxis
% X=backsub(A,B)
%
%Entrada
% - A es una matriz triangular superior invertible de orden n x n (cuadrada)
% - B es una matriz de orden n x 1 (columna)
%Resultados
% - X es la solucin del sistema lineal

%Calculo de la dimencion de B e inicializacion de X
n=length(B);
X=zeros(n,1);
X(n)=B(n)/A(n,n);

for k=n-1:-1:1
X(k)=(B(k)-A(k,k+1:n)*X(k+1:n))/A(k,k);
end

No hay comentarios:

Publicar un comentario