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,