ELIMINACION GAUSSIANA CON PIVOTEO PARCIAL (SOLUCION POR ETAPAS)
ANALISIS NUMERICO COMPUTACIONAL
fprintf(' ELIMINACION GAUSSIANA CON PIVOTEO PARCIAL (SOLUCION POR ETAPAS)\n\n\n'); %fprintf me permite ingre!r "#ment!ri# $e m!ner! tet&! &e p&e$en %#rient!r ! &&!ri# en e &# $e pr#gr!m! %inp&t e &n "#m!n$# $e #i"it&$ $e entr!$! $e $!t# $e &&!ri# A*inp&t('Ingree ! m!tri+ A * \n'); ,*inp&t('\nIngree e -e"t#r ,. "#rrep#n$ite ! # termin# in$epen$iente ,*\n'); % L! m!tri"e A / , $e,en er ingre!$! entre "#r"0ete ep!r!n$# ! %"#&mn! me$i!nte "#m! '.' / ! fi! me$i!nte p&nt# / "#m! ';' %permite #,tener e t!m!1# $e ! m!tri+ A 2n.m3*i+e(A); C*2A.,3; % ! m!tri+ C. repreent! ! f#rm! $e ! m!tri+ !&ment!$! 2A,3 fprintf('\nL! M!tri+ C. &e "#rrep#nte ! ! m!tri+ !&ment!$! 2A,3 e * \n'); $ip(C); % ! f&n"i#n $ip n# permite imprimir &n! -!ri!,e en e ep!"i# $e tr!,!4# if n**m %"#mp!r! e n&mer# $e "#&mn! / fi!. p!r! #,er-!r i #n ig&!e f#r 5*67(n86) fprintf('\n ETAPA %g*\n\n'.5) m!/#r*9; %!ign! "#m# "er# e n&mer# m!/#r $e ! fi! fi!m*5; %!ign! ! fi! 5 "#m# ! "#&mn! &e tiene e n&mer# m!/#r f#r p*57n if m!/#r:!,(C(p.5)) %e ,&"! e n&mer# m!/#r en ! fi! ; m!/#r*!,(C(p.5));%"!m,i# $e m!/#r fi!m*p; %"!m,i# $e fi! en$ en$ if m!/#r **9 fprintf('\nE item! tiene infinit! #&"i#ne\n') ,re!5 %e interr&mpe e pr#gr!m! "#n ! intr&""i#n ,re!5. /! &e %i m!/#r*#. m! !$e!nte e #,tiene &n! $i-ii#n p#r
%"er# ee if fi!m <* 5
ANALISIS NUMERICO COMPUTACIONAL
f#r 4*67(n=6) !&*C(5.4); %p!r! p#$er inter"!m,i!r ! fi!. &tii+!m# &n! %-!ri!,e !&ii!r
C(5.4)*C(fi!m.4);
C(fi!m.4)*!&; en$ en$ en$ fprintf('\nL! m!tri+ "#rrep#n$iente ! et! et!p! !nte $e pr#"e#7\n')
$ip(C) fprintf('\nL# M&tipi"!$#re "#rrep#n$iente ! et! et!p! #n7\n') f#r i*(5=6)7n m(i.5)*C(i.5)>C(5.5); %f#rm&! m&tipi"!$#re
fprintf('\nm(%g.%g)*'.i.5)
$ip(m(i.5)); f#r 4*57(n=6) C(i.4)* C(i.4) 8 m(i.5)?C(5.4);%f#rm&! n&e-! fi! en$ en$ fprintf('\nL! m!tri+ "#rrep#n$iente ! et! et!p! $ep&e $e pr#"e#7\n')
$ip(C) en$ f#r i*n78676
&m!*9; f#r p*(i=6)7n &m! * &m! = C(i.p)?@(p); en$
@(i)*(C(i.n=6)8&m!)>C(i.i); %f#rm&! $e ! &ti"i#n regrei-! / #&"i#n $e ! -!ri!,e en$
ee %f&n"i#n !ign!$! $e if. en "!# $e &e ete e! f!#
ANALISIS NUMERICO COMPUTACIONAL
fprintf('\nERROR7 L! m!tri+ NO e "&!$r!$!\n'); en$ fprintf('\n\n SOLUCION7\n'); fprintf('\n\nL! m!tri+ A, fin!7\n'); $ip(C) fprintf('\n\nL! #&"i#n $e @6 0!t! @n e7\n'); %! "#ntin&!"i#n $e &tii+! &n! intr&""i#n f#r. p!r! m#tr!r e &&!ri#. %# re&t!$# $e &n! m!ner! m! #r$en!$! f#r i*67n @i*@(6.i); fprintf('\n@%g*'.i) $ip(@i); en$
MATLAB PIVOTEO PARCIAL ESCALONADO
%-Modulo: GPivParEsc %-Auor: Aior Cald!ro" Mari"!# $aiorc"av!&alia'co() %-Pro*osio: Calculo d! !cuacio"!s li"!al!s *or !l (!odo d! Gauss co" %
*ivo!o *arcial !scalo"ado
%-+u"cio"!s: %
-GPivParEsc ,u"cio" *ri"ci*al.
-I"!rca(/io+ ,u"cio" d!clarada local(!"!.
%-Co(!"arios: El codi&o u!"! !sa co(!"ado0 lo u"ico a r!salar !s !l uso % d!l
d! la u"cio" 1c*ui(!10 cu2a u"cio" !" !s! (odulo !s la o/!"cio"
% i!(*o ra"scurrido dura"! la !3!cucio" d!l *ro&ra(a $i(!*o4i"icial %
i!(*o45"al')
u"cio" R 6 GPivParEsc,A0B. 6 c*ui(!7 %O/!"!(os !l los s!&u"dos acual!s or(a ra7 %P!di(os al *ro&ra(a 8u! us! raccio"!s !" v!# d! "u(!ros $co" !so !l !rror s!ra (!"or) ANALISIS NUMERICO COMPUTACIONAL
i ,,"ar&i" 9 . ; ,"ar&i" < .. %si !l "u(!ro d! ar&u(!"os "o !s !l *!dido''' *ri",1El "u(!ro d! ar&u(!"os di5!r! d!l "!c!sario'="1.7 A 6 i"*u,1Por avor !scri/a a co"i"uacio" la (ari# d! i"co&"ias A:="1.7 B 6 i"*u,1A co"i"uacio" !scri/a la (ari# d! !l!(!"os i"d!*!"di!"!s B:="1.7 !"d >+0C? 6 si#!,A.7 %O/!"!(os !l ord!" d! la (ari# -< + 6 "u(!ro d! 5las C 6 "u(!ro d! colu("as >+B0CB? 6 si#!,B.7 %O/!"!(os !l ord!" d! la (ari# -< +B 6 "u(!ro d! 5las CB 6 "u(!ro d! colu("as i ,,+ @6 C. ; ,CB < . ; ,+B @6 +.. *ri",1Error: El ord!" d! la (ari# A 2o la (ari# B0 "o !sso" !llos corr!cos' =" Por avor r!vis!lo1.7 r!ur"7 %Mosra(os u" (!"sa3! 2 ali(os !"d S 6 >7?7 %La (ari# S AB 6 >A0B?7 %La (ari# r!sulado d! co"ca!"ar A 2 B or i 6 :+ %Para cada 5la''' M 6 a/s,A,i0..7 %To(a(os !l *ri(!r !l!(!"o or 3 6 :C %Para cada !l!(!"o d! la 5la''' $!l *ri(!ro 2a lo !(os o(ado) i a/s,A,i03.. < M %Si !l !l!(!"o acual !s (a2or 8u! !l (ai(o r!&isrado''' M 6 a/s,A,i03..7 %Es! *asa a s!r !l (ai(o !"d !"d S,i. 6 M7 %El !l!(!"o d! la 5la i d! S !s !l (ai(o r!&isrado i S,i. 66 %Si al&u" !l!(!"o d! S !s ''' *ri",1No !is! solucio" u"ica *ara !l sis!(a' +i"ali#ado SIN !io1.7 %No !is! solucio" u"ica r!ur"7 %Para(os !l *roc!so d! cr!acio" d! la (ari# S 2 sali(os !"d
ANALISIS NUMERICO COMPUTACIONAL
!"d or i 6 :,C-.% Para cada colu("a''' or 8 6 i:,+.%Ord!"a(os la colu("a ' 8' a,0i.S,.9,0i.S,.9'''9a,80i.S,8. Ma+ 6 87 %La 5la co" !l (ai(o $asi&"a(os !l *ri(!ro al i"icio) Ma 6 a/s,AB,80i..S,8.7 %El (ai(o $asi&"a(os !l *ri(!ro al i"icio) or 3 6 ,8F.:+ %R!corr!(os la colu("a !" /usca d! oro *osi/l! (ai(o i AB,30i.S,3. < Ma %Si !l !l!(!"o acual !s (a2or 8u! !l (ai(o r!&isrado'' Ma 6 a/s,AB,30i..S,3.7 %Asi&"a(os !l valor d! !s! (ai(o Ma+ 6 37 %Asi&"a(os la "u!va 5la (ai(a a su varia/l! !"d !"d i ,Ma+ @6 8. %Si !l (ai(o "o s! !"cua"ra !" la 5la d! i"icio AB 6 I"!rca(/io+,AB0Ma+08.7 %I"!rca(/ia(os la 5la d! AB *ara 8 coi"cida S 6 I"!rca(/io+,S0Ma+08.7 % Ta(/i!" s! d!/! ca(/iar !" S0 *u!so 8 !l ord!" a ca(/iado !"d !"d %Colu("a 2a ord!"ada i ,AB,i0i. 66 . %El !l!(!"o d! la dia&o"al NO *u!d! s!r 0 si !s c!ro "o !is! solucio" u"ica *ri",1No !is! solucio" u"ica *ara !l sis!(a' +i"ali#ado SIN !io1.7 %No !is! solucio" u"ica r!ur"7% Sali(os !"d or 3 6 ,iF.:+ %Modi5ca(os las 5las *ara ac!r 1s ( 6 ,AB,30i.AB,i0i..7 %( s!ra !l (uli*lica"do0 !s "!c!sario d!clararlo %a"!s d! i"!rca(/iar las 5las *or8u! !l valor d! los !l!(!"os %d! las (is(as variara cua"do a&a(os 1s0 2 *or co"si&ui!"! !l %valor d! ( or 8 6 :,CF. %Para cada !l!(!"o d! la 5la''' AB,308. 6 AB,308. - (AB,i08.7 %R!ali#a(os !l ca(/io ANALISIS NUMERICO COMPUTACIONAL
!"d !"d !"d or i 6 :+ i ,AB,i0i. 66 . %Ni"&u" !l!(!"o d! la dia&o"al NO *u!d! s!r 0 si !s c!ro "o !is! solucio" u"ica *ri",1No !is! solucio" u"ica *ara !l sis!(a' +i"ali#ado SIN !io1.7 %No !is! solucio" u"ica r!ur"7 %Sali(os !"d !"d H 6 >7?7 %(ari# co" los r!sulados H,+. 6 AB,+0+F.AB,+0+.7 %Asi&"a(os !l uli(o !l!(!"o or i 6 ,+-.:-: % va(os calcula"do los d!(as Su( 6 7 %i"iciali#acio" d! la varia/l! Su( $*ara calcular !l su(aorio) or 3 6 ,iF.:+ %/ucl! *ara !l su(aorio0 *ara cada 5la''' Su( 6 Su( F AB,i03.H,3.7 %Calcula(os !l r!sulado acual 2 lo aJadi(os al co"ocido !"d H,i. 6 ,AB,i0+F.-Su(.AB,i0i.7 %O/!"!(os H,i. !"d *ri",1R!solucio" d!l sis!(a co"cluida !": %d s!&u"dos' El r!sulado: ="10 ,c*ui(! - ..7 R 6 H7 %El r!sulado !s H %A co"i"uacio" s! !"cu!"ra la u"cio" I"!rca(/io+0 cu2a u"cio" !s %i"!rca(/iar 5las d! u"a (ari#' E" *ri"ci*io !sa o*!racio" s! *u!d! %r!ali#ar d! dos or(as0 la !*u!sa 2 sa/i!"do 8 A,"0:. d!vu!lv! la 5la %co(*l!a " d! la (ari# A' o ! !sco&ido la *ri(!ra or(a *or8u! !s u" %(!odo (as &ra5co 2 (!"os /asraco'
u"cio" K 6 I"!rca(/io+,V00. %+u"cio" i"!r"a *ara i"!rca(/iar 5las
ANALISIS NUMERICO COMPUTACIONAL
>+V0CV? 6 si#!,V.7 %O/!"!(os !l "u(!ro d! 5las 2 colu("as $+ 2 C r!s*!ci/a(!"!) CV 6 7 %A8ui s! u" valor or " 6 :CV %Bucl! *ara r!corr!r las colu("as CV 6 V,0".7 %CV 6 5la V,0". 6 V,0".7 %Asi&"a(os a V,0". 6 CV7 %Asi&"a(os a la 5la " la 5la CV !"d K 6 V7 %!l r!sulado !s la "u!va (ari# V
ANALISIS NUMERICO COMPUTACIONAL