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

Create successful ePaper yourself

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

2.4. Polinoame ortogonale 19<br />

face uz <strong>de</strong> proprietăt¸ile algebrice ale polinoamelor s¸i este sensibil la erorile <strong>de</strong> rotunjire.<br />

Fie {Q0,Q1,...,Qn−1} o familie ortonormală <strong>de</strong> polinoame, astfel încât gradul<br />

luiQi să fieis¸i fieQ n ⊥ Qi, i = 0,n−1.<br />

Să consi<strong>de</strong>răm polinomul<br />

Q n(x)−αxQn−1(x)<br />

Pentru o alegere convenabilă a lui α = 0, acest polinom are gradul ≤ n−1,<br />

<strong>de</strong>ci<br />

n−1<br />

Qn −αxQn−1 =<br />

i=0<br />

αiQi<br />

Dacă 〈Q n,Qi〉 > 0 pentru oricei = 0,n−1 trebuie să avem<br />

0 = 〈Q n,Qn−1〉 = α〈xQn−1,Qn−1〉+αn−1<br />

(2.13)<br />

0 = 〈Q n,Qn−2〉 = α〈xQn−1,Qn−2〉+αn−2<br />

Putem alege α = 1, <strong>de</strong>oarece înmult¸irea cu o constantă nu afectează ortogonalitatea.<br />

Deci αn−1 s¸i αn−2 se pot obt¸ine din ecuat¸iile <strong>de</strong> mai sus. Aplicând<br />

rat¸ionamente similare lui Qi pentru i < n−2 obt¸inem αi = 0 pentru i < n−2.<br />

Aceasta sugerează următoarea formulă <strong>de</strong> recurent¸ă pentru calculul lui Q n:<br />

s¸i<br />

Q n(x) = (x+an)Qn−1(x)+bnQn−2(x), n ≥ 2 (2.14)<br />

Qn = Q n<br />

Qn<br />

an = −〈xQn−1,Qn−1〉 (2.15)<br />

bn = −〈xQn−1,Qn−2〉 (2.16)<br />

Se verifică că pentru an s¸i bn astfel <strong>de</strong>terminate avem 〈Q n,Qi〉 = 0, i =<br />

0,n−2 s¸i căQ n cu an s¸ibn <strong>de</strong>terminate <strong>de</strong> (2.15) s¸i (2.16) este unic <strong>de</strong>terminat.<br />

Deci (2.14) ne dă o formulă <strong>de</strong> recurent¸ă pentru calculul polinoamelor ortogonale<br />

(ortonormale) în L 2 w [a,b]. Vom începe punând Q0 = b0, un<strong>de</strong> b0 este o<br />

constantă astfe încât Q0 = 1 s¸i luăm Q 1 = (x+a1)Q0. Din<br />

se <strong>de</strong>termină<br />

s¸i se continuă.<br />

〈Q 1,Q0〉 = 〈xQ0,Q0〉+a1 = 0<br />

a1 = −〈xQ0,Q0〉<br />

Exemplul 2.4.1 Pentru polinoamele Cebîs¸ev I aplicând (2.14)-(2.16) se obt¸ine<br />

Tn(x) = 2xTn−1(x)−Tn−2(x).

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

Saved successfully!

Ooh no, something went wrong!