Calcolo numerico: algoritmi

Esame Calcolo numerico

Facoltà Ingegneria dei sistemi

Dal corso del Prof. Antonietti Paola Francesca

Università Politecnico di Milano

Appunto

4,0 / 5 (2)

Scarica

L'appunto contiene gli algoritmi implementati durante il laboratorio informatico del corso di Calcolo numerico presso il Politecnico di Milano. Gli algoritmi presenti, utili alla risoluzione di molteplici problemi matematici, sono: bisezione, Newton, punto fisso, algoritmi per l'approssimazione di integrali, algoritmi per trovare autovalori e autovettori di una matrice, algoritmi per risoluzione di sistemi lineari e di equazioni differenziali

…continua

Anteprima

Vedrai una selezione di 1 pagina su 4

Disdici quando
vuoi

Acquista con carta
o PayPal

Scarica i documenti
tutte le volte che vuoi

Estratto del documento

Bisezione

xvect[tr, x0]=bisez (a, b, nmax, tol, ffun)
if(funz (a) * funz (b) > 0)
error (‘...‘)
end
tr = 0;
xvect[tr] = [a; b]; x1[tr] = b-a;
fa = funz (a);
xr = abs (b-a) * nmax
while (ens stol rem lt nmax)
x = (a+b)/2
tr = tr + 1; fx = funz (x)
if ((xvect Stoll x0)
fx = x0 (med);
else
exit
rest = abs XR
end
xvect [tr, med] = Vmed,
vecr[tr] = vece[x, x1]
if (fxvect == 0)
a = xov;
f = saum
end
fxb = x(aud);
end

Newton

xvect[ti] = myste ton[del, nocomax, tol, ffun
ngton (blo, dtonax, tol, funddown
lt = 0;
xevec[Let];
exit (fends)
while (onte xtol rext nmax)
kv = xveet(eud)
if (divv(red) == 0)
dnp (restore)
iv = tf 1k, xi,
xtc (break)
end
li xiv = kv - fun (xv) diffun (ru);
exit = b3 : (cv - nxl)
xvect (for \vect);
lsktr = xi
end
end

Ptofis

xevo, lt ei = ptofis (nn fonenv, enon, lot
0;
tr=0;
xvect[zeo]
ena=t=tol+i;
while (entot loz nmasyti)
x_vec v vect (ewd)
kin = hit (x yo) ;
xr1 = abs (n s-xv);
eirw
end

Predcomp

I = pmeadcomp (f, a, b, N)
fx = (b-a)/ N
x_med= = cib; b;
Yf > (fx mode);
g_env sum (g);
end

Simcomp

I=encomp (f, a, b, n)
_n= [0: g: N]
[I, conol n];
— >2 chol (chesis ( x mde) (audi)
_x_med= db; h; b;
→2 dougl; (ilt ewa2:d1)
eis g, extra (f(x_mod = )fmdeo (cmode))
+ 4 * zum (Ymed)
end

Eigpower / Power / Potenze simetri

<iambdain, ytei - eigpower (A, tol, nmajr, xo)
x for= where;; xj;; =
yt rei 40 (E);
X = inAv
right trip; _lambda =
while (xt stoli x nmax)
zte;
y= x/norm(x);
x -(Lj)
_{lambda(rev)\;}

_{&rambda.

fi=<
avdiv}

Dettagli

Publisher

cescamhr1

A.A. 2014-2015

4 pagine

SSD Scienze matematiche e informatiche MAT/08 Analisi numerica

I contenuti di questa pagina costituiscono rielaborazioni personali del Publisher cescamhr1 di informazioni apprese con la frequenza delle lezioni di Calcolo numerico e studio autonomo di eventuali libri di riferimento in preparazione dell'esame finale o della tesi. Non devono intendersi come materiale ufficiale dell'università Politecnico di Milano o del prof Antonietti Paola Francesca.