Culegere de probleme de Analiz˘a numeric˘a
Culegere de probleme de Analiz˘a numeric˘a
Culegere de probleme de Analiz˘a numeric˘a
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).