20.07.2013 Views

Culegere de probleme de Analiz˘a numeric˘a

Culegere de probleme de Analiz˘a numeric˘a

Culegere de probleme de Analiz˘a numeric˘a

SHOW MORE
SHOW LESS

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

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

Saved successfully!

Ooh no, something went wrong!