30.01.2015 Views

4 Método dos quadrados mínimos

4 Método dos quadrados mínimos

4 Método dos quadrados mínimos

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

Da mesma forma que no caso discreto, temos funções ortogonais com relação ao<br />

produto escalar, como mostrará o exemplo abaixo.<br />

Exemplo 4.7.3:<br />

Os polinômios de Legendre, defini<strong>dos</strong> por<br />

(k)<br />

1 d<br />

P 0 (x) ≡ 2 k<br />

1, P k (x) =<br />

[(x − 1)] , k = 1, 2, ...<br />

k (k)<br />

2 k! dx<br />

b<br />

são ortogonais em [–1, 1], com relação ao produto escalar p , q = ∫ p(x)q(x)dx .<br />

a<br />

Fica como exercício a verificação de que os três primeiros polinômios de Legendre<br />

P0(x) ≡ 1, P1(x) = x e P2(x) = 1 (3x<br />

2 − 1)<br />

são ortogonais entre si.<br />

2<br />

Uma observação interessante é que, em geral, polinômios ortogonais satisfazem<br />

uma fórmula de recorrência de 3 termos, ou seja, da<strong>dos</strong> P 0 (x) e P 1 (x), conseguimos<br />

construir P k (x), k = 2, 3, ...<br />

No caso <strong>dos</strong> polinômios de Legendre, a fórmula de recorrência é<br />

⎛ 2j + 1⎞<br />

⎛ j ⎞<br />

P j + 1 (x) = ⎜ ⎟xP<br />

j(x)<br />

− ⎜ ⎟Pj−1<br />

(x)<br />

, j = 1, 2, ...<br />

⎝ j + 1 ⎠ ⎝ j + 1⎠<br />

Exemplo 4.7.4:<br />

Vamos aproximar f(x) = 4x 3<br />

intervalo [a, b] = [0, 1].<br />

por um polinômio do primeiro grau, uma reta, no<br />

ϕ(x) = α 1 g 1 (x) + α 2 g 2 (x) = α 1 + α 2 x, α 1 , α 2 ∈ R<br />

(g 1 (x) ≡ 1 g 2 (x) = x).<br />

Pelo que vimos, (α 1 , α 2 ) é a única solução de Aα = b onde<br />

A =<br />

⎡a11<br />

⎢<br />

⎣a21<br />

a12<br />

⎤<br />

a<br />

⎥<br />

22 ⎦<br />

α =<br />

⎡α1<br />

⎤<br />

⎢ ⎥<br />

⎣α2<br />

⎦<br />

b =<br />

⎡b1<br />

⎤<br />

⎢ ⎥ , sendo<br />

⎣b2<br />

⎦<br />

a 11 =<br />

a 12 =<br />

b<br />

1<br />

2<br />

1 = 1<br />

∫ g (x)dx ∫1dx<br />

=<br />

a<br />

b<br />

0<br />

1 2<br />

1<br />

x 1<br />

1 a 21<br />

2 2<br />

0<br />

0<br />

∫ g (x)g2<br />

(x)dx = ∫ xdx = = =<br />

a<br />

103

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

Saved successfully!

Ooh no, something went wrong!