gaussparcial - end % desplazar todas las filas...

Info iconThis preview shows pages 1–4. Sign up to view the full content.

View Full Document Right Arrow Icon
Sheet1 Page 1 % PROGRAMA: ELIMINACION GAUSSIANA CON PIVOTEO PARCIAL % Descripcion:El programa encuentra las raices de un conjunto de ecuaciones lineales de la forma % A x X = B, utilizando el metodo de eliminacion gaussiana con la estrategia de % pivoteo parcial. % curso: Calculo Numerico % programacion: Rony Yupanqui G. % fecha: Setiembre 2006 %------------VARIABLES---------% % numero de filas nf=3 nc=nf+1 ma=zeros(nf,nc) mb=zeros(nf,nc) % matriz A con los segundos miembros B ma=[1 3 -5 1 %-----ALGORITMO DEL PROGRAMA----------% disp(' PROGRAMA: ELMINACION GAUSSIANA CON PIVOTEO PARCIAL') for col=1:(nf-1), valmax=0 filmax=col % encontrar el valor maximo for a=col:nf, if abs(ma(a,col))>valmax valmax=abs(ma(a,col)) filmax=a end end % guardar una copia de la matriz for a=1:nf, for b=1:nc, mb(a,b)=ma(a,b) end end ma=zeros(nf,nc) % copiar hasta la fila anterior for nfil=1:(col-1), for ncol=1:nc, ma(nfil,ncol)=mb(nfil,ncol) end end % copiar la fila elegida for ncol=col:nc, ma(col,ncol)=mb(filmax,ncol)
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 2
Background image of page 3

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 4
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: end % desplazar todas las filas filaux=filmax while (filaux>col), for b=1:nc, mb(filaux,b)=mb(filaux-1,b) Sheet1 Page 2 end filaux=filaux-1 end % copiar todas las filas for a=col+1:nf, for b=1:nc, ma(a,b)=mb(a,b) end end % convertir a cero los coeficientes for a=col+1:nf, valor=ma(a,col)+0 for b=col:nc, r= double(ma(col,b))+0 ma(a,b)=double(ma(a,b))-double(ma(col,b))*double(valor)/double(ma(col,col)) end end end % eliminacion gaussiana con la matriz triangular x=zeros(nf,1) fila=nf while fila>0, x(fila,1)=ma(fila,nc)/ma(fila,fila) pos=nc-1 while pos>fila, x(fila,1)=x(fila,1)-ma(fila,pos)/ma(fila,fila) pos=pos-1 end for a=1:(fila-1), ma(a,fila)=ma(a,fila)*x(fila,1) end fila=fila-1 end disp('ARREGLO TRIANGULAR SUPERIOR DE LA MATRIZ A Y EL SEGUNDO MIEMBRO B DE LAS ECUACIONES') disp(ma) disp('SOLUCIONES X DE LAS ECUACIONES') disp(x) Sheet1 Page 3-2 -4 -1 3 3 1 -2 -1] Sheet1 Page 4...
View Full Document

Page1 / 4

gaussparcial - end % desplazar todas las filas...

This preview shows document pages 1 - 4. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online