Culegere de probleme de Analiz˘a numeric˘a
Culegere de probleme de Analiz˘a numeric˘a
Culegere de probleme de Analiz˘a numeric˘a
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
104 Aproximări în medie pătratică<br />
Ecuat¸iile lui (7.2) se numesc ecuat¸ii normale. Determinantul lui (7.2) este <strong>de</strong>terminantul<br />
Gram al vectorilorg1,...,gn, G(g1,...,gn) = 0, căcig1,...,gn sunt<br />
liniar in<strong>de</strong>pen<strong>de</strong>nte.<br />
Deci g∗ există s¸i este unic.<br />
În cazul discret putem lucra analog cu<br />
m<br />
〈f,g〉 = w(xi)f(xi)g(xi).<br />
i=0<br />
Problema poate fi tratată s¸i astfel:<br />
Fie<br />
m<br />
<br />
G(a1,...,an) = w(xi) f(xi)−<br />
i=0<br />
n<br />
<br />
akgk(x)<br />
k=1<br />
Pentru a <strong>de</strong>termina minimul luiGvom rezolva sistemul<br />
∂G<br />
(a1,...,an) = 0, i = 1,n.<br />
∂aj<br />
Observat¸ia 7.0.5 Dacă funct¸iile gk, k = 1,n formează un sistem ortogonal<br />
coeficient¸ii λ∗ k sau a∗k se pot obt¸ine astfel<br />
a ∗ k<br />
Problema 7.0.6 Dându-se punctele<br />
= 〈f,gk〉<br />
〈gk,gk〉 .<br />
(0,−4),(1,0),(2,4),(3,−2),<br />
<strong>de</strong>terminat¸i polinomul <strong>de</strong> gradul I corespunzător acestor date prin metoda celor<br />
mai mici pătrate.<br />
G(a1,a2,...,an) =<br />
m<br />
i=0<br />
∂G<br />
∂ak<br />
= 2<br />
m<br />
<br />
i=0<br />
gj(xi) = g i j<br />
yi −<br />
m<br />
<br />
i=0<br />
yi −<br />
n<br />
ajgj(xi)<br />
j=1<br />
2<br />
n<br />
<br />
ajgj(xi) gk(xi) = 0<br />
j=1<br />
n<br />
ajgj(xi)gk(xi) =<br />
j=1<br />
m<br />
yigk(xi), k = 1,n<br />
i=0