cuaderno computación cientifica
DESCRIPTION
métodos numéricos en matlab para ingeniería eléctrica y electrónicaTRANSCRIPT
-
Computaaon aentiflcaen IEE
Tv Jen te escribe aca'
desecindotelomejor !svs svs svs Svs SVS SVS svs svs
Fe Amo Juan Pablo Arevalo 08.02.14
-
27/01/15Corte
Antonetta 4 digitostyredondo
P= 0.54617 , 9=0,54601- Redondo
fl( PHO .s46z ; flttfasybo Valor exact :P -9=0,00016--0.168
fl(flA ) - fl( ofD= 0.0002=0,200010-3 Calcio del error relative
|2Xv0t6=
ozscsxiu' '
IQIBXIOT La aproxmauon Here
P* aproxlmac , on de p un Solo diyitoSignificantno .
IP .p* { snot H es el enter no - negative mas grande )La aproxlmquon Elene t digitos Significations
- por carte
fl / p) = 0.5461
fl(q)= 0.5460fl (0.5461-0.5460)=0.110-3
pp .pif = 0.375
-
XP =-62*62.06=-62.10IXZXTT = 2.610 "Xz= -62.0838384276 121
reescrib1endo1uexpresion-7InEIxnT-ybtkaty0yuknyitinxnyxnTiff.itaranxensn-zyIl.A-ln-la-x4zparan-liNatYny-itxtIZ-nIy-itxza-atHnjtenerencuenta1aslimitacionesenpuntoY-ytaf1otanteenMaqunasFlnparaffy-x3-G.lx2t3.2xtl.sXX2Xs6.1X23.ZXxnn.xnosxExoarotoeeesi.egeeieeImYsFYgIffHxsmredonJe0-l3.4errorrelativ.S.3KO.OD6.0Stt0.04fLx1-Hx-6.txt3.yxti.sfH-4.7D-t4.2fFnrnEatiEFaYHfageqIfminiIeeomaEaruneienentoOperauores-transpvestaiRnx1pRnxm-MvHiP1icaaonCntremutrices.EAt-7Cii-Cn.RnxPxNxn_yRnxn-SumaiRnxnxnmxn-RmxnC-AB-DCii-fklihbki@pnotuuonproouctocartesianol.C
--A+DCii=AijtBi ;- Multiplication porvnescalar : RRnmRn "
( =XACi;=XAij
-
Rnmxrm 'nRnm( = A .*B ( i ; =9iibi ;
C =A/BCij=ai%i ; ( bi ; to )
EFFETEfxxeifxitir[ Di , Xli ) Para derotnrun element del vector
IElR2I= . ,xn]Opera clones
REEN y FEINE=aI Zi - axi
zE=I+FZi=Xity : ( proJuutolnntemo
)C=Ity=XiYi
oy=5tax( actual , zauo 's del vector 5)Alyoritno para ( alcular txty
C=0n - Sumas
for i = tin n - multiplicauonesc=c = ctxlilyli )
o( n )29/01/15Sea Aern
"
yauerenoscalcdar 5 '=y+AI done IEINY TERM
+ ^ -^
-
Representation de la matriz por una Particia en filas o en column as
AnafFI ftp.#Hf.Etdiritdis=H
for iain A- ki , . . . .ci ] G=# , CzjfYi=yi+KtxFor j=1 . n
Stops = float operations [ gweiytxociNotauia
infEItEfg.ee#.aeilt:oiiiFnntaa
:*
for j= In
[en
5=5 tx ; AC : , ;)
A -. AtxytAtrnxixcrnsyem
*DNIF KeatenQmn)for Elim
mn+mxi ( [for ;= , ;n
5 '#i% )
Aij =AijtXiY ;for itlim end
[ . Ali , :) - Ali,:) +i5 end for Jain[ eaACiD=AEIHYHE
-
Opciones para multiplier Zmutnces
1 ;;]ftp.H.FI?*hY.YtE:D=kDtHakHtsHDAe1Rntr,BcRnny(enm.n#56ttf
, ][ 78 ]
C =CtAD for it
:m%O( mnr )
For J=i : n
( |;" tiiiitcaiittaliihnth 'sjki
fur j=l:ni# for k=lir5=5+a~f|yYdEiiHH#
=5+AI
Kij for k=i:rfor i=lim
( [ [roan anti.es?fttax( list . Clint ACIFDBCK'Dend endend
-
3/02/15Si sterna de ecuacioneshnea#a
tdizxzt . . . +ainXn=b ,
Az ,X,tAnXz+ . . . .tAznXn=bz >AI= : - =[AiD,< i ,i^a!ix,+anxz+
.
.+annxkb-ftp.xtekf( Ordain para que exist una solvcion vmcui Jet ( Altometodos de sowaonf Directors
[ Ltevativos
Alyvmscnsos : A = diagonal ( an ,An , . . . , Ann )oettthanau
. . .an=IIaii for ] aiixitbixisbilaiito- Matriz triangular superior
= = = Xn=bXamo
"
%o" =
an .mx . .
.tw#xn=bn..A = ~ , =b~ , - Un - i ,nnAmXn=bn -An - , ,n . ,MEIR [ - eslna ecuacioni Aiixitai , it , Xotit . . . . tai .^Xn= bi
Sustituuoi haaa atria. Xi=a.[bi - .,ait[ it EITXE ; - EzEz+aE ,Esttztx ,E ,Eliminator Gaussian a
Efate xx xEz= -31+22 +63=4 * , g }x ##EJ
.
SX , - 12+53=6 X X ,
At"=[ai ; 'D A" :[ aim ]
se seleuionaelximas arardede los En n=Al " =[aie
iaE+psEz
-
EFEZTAZE ,
Eziazx , tAuXatAbX}=bzEl : A , ,X , taizxzt U , }X}=b ,
Ez( autxa , DX, tlazztxiaidtzt ( ant XA , )X}= bztxb ,( Como selecuonar x ?
Autxa , ,=OXz=jolt As =- 9 Yan
-7 0 Paso 1}g|4D ec . wrote IOX , -72+03=7EzEz+AE ,
A , = Y, or = - t % = 0.3
(-3+0.310) X , + (2+0.371) Xzt ( 6+0.383 = 4+0.37OX , - 0.1 Xz +63=6 . )
3 E ] to E ,Xi - Yz = -0.s
( St to.s)io) , tftf as )f7Dzt( Has)o) ] = ( Gtfo .s)7 )OX , + 2.52+53 = 2.5
:O:He .EE#EHEtPaso z
'
.
el. p Note 2.52+53=2 . s
Est Eztx Ez
a = - to .D/z, ,
= 0.04
t0.1+0.04 . 2. s ) Xz t(6+0.0451=4.1+0.042.5)02+6.23=6.2
-
Sistena r7 0 '
73=1.002.5 S Xz 2. s Xz = - 1.0 0 6.2 X ] 6.2 X|
= 0
for ITI :( A - ' ) ; 1 . Lato sobne las columnsfor I=jti in . % law sabre las Filas por debate de j/
maultitli.is/Ali.il ;Iii ) '. Ali , :j+mut*Afi , :) Ali ,iinFAli,iiDtmuH*Als;jin )blilibliltmulttbj ) ;
"
endend
10/02/19Factorization LV
>.
A = LU Axstbt.ir#r s . superior LU=ps
Tibtriangular inferior At LUuI=5
A :[ ai ;] Pqqawytarlgselezentoueliltzt. . .tn=tlei ,i< n " Uiltzt
.. .+n=nh
z
Suma :n(ntD=nmer0
de incognita
Paratenernewuciones para incognita : li ; - - 1,5=42 , . . . ,A-1=013
, '
- 1=lzN+Uz3=j't
,
Ii =,#n?, 4=l]*e3wYV]}=3.75= 3.334=4 , -1=42
-1=1 , ,u , ,Lu=- Yvi
,-0,4 4=lziUztVz2Uzz= 4 - luUh4 - 0.253.3
- 1=lsiU , ,l},=Yu , ,= -0.25 - 1=bsiU ,ztl]zUzzl,z= -1 - (-0.2/1)/3.75=-1/3
-
j ; J
ic ;
i oh's in , , , , = it|'mint:DAiileaU,;f%4lAii=i,
link ;
K "C is;lai;=, link ;
Lpavacn.PpZardj.esim0Pus0queuonocenos.2.a.YHi.QUtiit.iiitjLlinij-DLfsttimilAijsdiiUliAiiiEnlieVkiAzj-luU@tlzzUzjiaii-HE.liaueftUii.Ui-siligiUnt.i.tli.i.ilt.s
; U ;j=A;j - jlitUt;
4in .ph/,uli:;.n,
, :p , ,)" Vliit 'Di" Ipa , ;= ; tr%
"
Ariel ,U , ;tlzzUz ;:
D4iaie'teen.ua#tYj'Iahauxi
bjj - A ; ; - xlikVk;
3) Lthl :n , ;) jliti .. a )
-
Paramh"IhjfY'ftp.n
, ;y , , , ,
Tareniescnbircodigoparaeltoevehfort=:
resolver sistenu triangular para Vllijt , ;):eternnar Visit )% Determiner Ltithin , ;)
e
:2/02/15Inkialhar L Como una matrlz identical y V Como una mwtrlz de ceros
for jF ja- Af , ;)
. e :=At , ;)resolver Uitj - h , 1 :b- , ) )z= oil :( ; . 1 ) )U( it ; - i ) , ;)=z
NCJ :nj=a( I ;^) - Lljin , 1 :( I - I ) ) .zend Ai ;=lh"
"tirades
UH ,jj=v( j ) Ai ;= =fiHki =Llitiin , ; )=vVHti:n)/vl ;) I
end lixukitti ;Uj ;r= ,iitiBjesfjH
. Bj!s
{ #I; . ?
-
Aliinitlljinrjtjulitsittlljindluli
AelpixnA es una mutriz oefnwa positions , ETAI >o para too Vector Jena :O ' ' of DIlxn nxn nxl
A=#-
, =L to
'][to ] [ caobjfo ]=[i][Y]=a>o
- I,=[;][' ]faY][ , ]=[i][Y]=J>o-53=1il[ ii ][cY][ ''Ifattobff',]=a+b+c+o>o- Ee
, :f,][it][eY][!]=Df
, ]=a . b- cto >o
AelpixnA es una mutriz oefnwa Semipositions,ETAIZO para too Vector Jen- D 6 es triangular inferior
A = GGT si Aes definida positive
6 " 0 0 G , 62 , 63 , 61,2=4 Z,
"
,yh], 622f.+ 6226] 62,2+6,22=4C 31 6 }z 63 } 0 0 633 6312+6322+635=4G GT 44+3,55+9332=4
Gh6z , = - /-1=421931+932922
+1/4 + 6222=4 Gz , = - 1/2 43 }= 7
- 1=+44+12.932 gnittgz=FYI GNG } , = - IGzz =ltfz 63 , = - 1/2
-
".or I =L : n
If j > ,
Atsin , ;) = Ali :n , ;) - Ali : al :b - Dalin :( It )"
e :Hindtatiinikfend
Para :tie;y.ae #H:*' " 'lBD7=&b#Gi=l"
'
bixbixawnihnyie( =D
't At , ;).io#6lInH6liItnoIfhift. :c D
17/02/15tunica , lterativas matrices granites con ( aractensticas especial es
Metodosde iteration de Punto Fido y=sX e > Un Punto fyosi fK)=
Xkti #*hk=o , , , .in MathConviction inicial dada 'w:X AE=EAxtx=btx
E- MITE# IHAITI - b
I=g , =k+Ik' el problem es resolvermx'tb
-
19/02/15(
erosdeUnpo1n0mi0WnixFXntaiXnt.i.an.iXtanAjERi-h.i.1nwrlxFotienerauesfYoFhpYesosco.uyadoDPorefemp1o-pdinomiodeorden2ax2tbxtc-oXi-tsTs4acX2-b-b1S-b2-4acza-DiscrinnantePHesunpohnomiooeoroennlDiuuenddH2tPxt9tpo1inomodeoruen2-fautorcuadraticolDms0DcogxefteIyoyPkktx2tPxtacHtra1@reenn-y0rdeneDWnixFXntAiXnt.ian.iXtan-k2tptt4Hxt2tb.xt3t.i.tbn.3Xtbn.DtCRxtSTAi-b.tp-7bilP1.a
,- p
xn xn ' ' "
's
. .. .
. X ' 1
bhpbztpbibstpbzbnzt pbnz+9 t4b ,+qPbbI4Otbns 9 bn -2
1 A , gz U3 Oh -2 an -1 An
az= bztpbt 9- =bz=az - Pb , - of =bz( hada ]=b}tPbzt4bib}=az - pbz - of bi al vespaldo
aiibitppi - itqbi -2 bEa[ pbi - rofbi - 2anztbnttpbnztotbn -4 bn -2=9^2 - Pbn - 3-Fbn -4
-
An -1 = pbnztqbn . 3tR R= an - itpbn . z - qbnz =DAn= Fbn - zts Stan . 9- bn -2=0
ftp.g/=an--phn-zlp,q ) - of bn -31kg )Slp ,q ) - an - efbn . 21h41
q '=ba, , P '=an9bJ Mitooo de Lnbn -2( p,q1b , , . . . , bn -2 , KS ) ( p ' ,q ' )
R( P ,qj=ob
's ... ,bm,
R ,s )
Slp ,gt=o Metodo de Newton - RaphsonHxfftxolth - xdf 'Kd+yYf" lxojtlxxo )3f
"
'h + . . . .ftp.oo-SH~fkdtcx-xojf 'kde . xdf 'Lxd= - fko )
""xx*#=xot
,
Ike , )
RIP ,H= Rlpo ,%l+(rPDpf+( qaolsftth
SH,H= SIR,%1+(rPDDo+KItEIDfII%:*
' on:#
peep.tl?sYpogtaottEajtgfp*=Po+tfhEosFat#tEy*.o+ttH.3
. Ho