20.07.2013 Views

Software matematic

Software matematic

Software matematic

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

3.3. Valori s¸i vectori proprii în MATLAB 103<br />

(b) Împărt¸it¸i fiecare linie cu maximul în modul din linia respectivă s¸i apoi utilizat¸i<br />

eliminarea gaussiană.<br />

(c) Rezolvat¸i sistemul folosind toolbox-ul Symbolic.<br />

Problema 3.12. (a) Dându-se o funct¸ie f ∈ C[a,b], să se determine ˆs1(f;·) ∈<br />

(∆) astfel încât<br />

S 0 1<br />

să fie minimă.<br />

b<br />

a<br />

[f(x)− ˆs1(f;x)] 2 dx<br />

(b) Scriet¸i o funct¸ie MATLAB care construies¸te s¸i rezolvă sistemul de ecuat¸ii normale<br />

de la punctul (a).<br />

(c) Testat¸i funct¸ia de mai sus pentru o funct¸ie s¸i o diviziune alese de dumneavoastră.<br />

Problema 3.13. Calculat¸i aproximat¸iile discrete în sensul celor mai mici pătrate ale<br />

funct¸iei f(t) = sin π<br />

2 t pe 0 ≤ t ≤ 1 de forma<br />

ϕn(t) = t+t(1−t)<br />

n<br />

cjt j−1 , n = 1(1)5,<br />

j=1<br />

utilizând N abscise tk = k/(N +1),k = 1,N. De notat căϕn(0) = 0 s¸iϕn(1) = 1<br />

sunt valorile exacte ale lui f înt = 0 s¸i respectiv t = 1.<br />

(Indicat¸ie: Aproximat¸if(t)−t printr-o combinat¸ie liniară a polinoamelorπj(t) =<br />

t(1−t)t j−1 ,j = 1,n.) Sistemul de ecuat¸ii normale are formaAc = b,A = [(πi,πj)],<br />

b = [(πi,f −t)], c = [cj].<br />

Ies¸ire (pentru n = 1,2,...,5) :<br />

• numărul de condit¸ionare al sistemului;<br />

• valorile coeficient¸ilor;<br />

• eroarea maximă s¸i minimă:<br />

emin = min<br />

1≤k≤N |ϕn(tk)−f(tk)|, emax = max<br />

1≤k≤N |ϕn(tk)−f(tk)|.<br />

Executat¸i de două ori pentru<br />

(a) N = 5,10,20,

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!