Um teorema de aproximação. Interpolação linear

Seja f ( x ) uma função que possui pelo menos n derivadas f', f'', \dots, f^{(n)} em um intervalo. Seja x 0 um ponto deste intervalo. Dizemos que f ( x ) tem um zero de multiplicidade n de f ( x ) se f ( x 0 ) = 0 , f'(x_0) = 0, , f ( n 1 ) ( x 0 ) = 0 , mas f ( n ) ( x 0 ) 0 . Assim, se

f(x_0) = 0 = f'(x_0) = \dots = f^{(n-1)}(x_0),

podemos dizer que x 0 é um zero de f ( x ) de multiplicidade pelo menos n .

Se f ( x ) tem zeros (de multiplicidade 1 ) em x 1 e x 2 , então, pelo Teorema de Rolle, f'(x) tem um zero (de multiplicidade pelo menos 1 ) entre x 1 e x 2 . Se f ( x ) tem um zero de multiplicidade pelo menos 2 em x 1 e um zero de multiplicidade pelo menos 1 em x 2 e, digamos, x 1 < x 2 , então f'(x) tem um zero de multiplicidade pelo menos 1 em x 1 e um zero t de multiplicidade pelo menos 1 tal que x 1 < t < x 2 . Então f''(x) tem um zero de multiplicidade pelo menos 1 entre x 1 e t . Da mesma forma, se f ( x ) tem zeros em x 1 , x 2 , x 3 , e x 1 < x 2 < x 3 , então f'(x) tem um zero entre x 1 e x 2 e um entre x 2 e x 3 . Portanto, f''(x) tem um zero entre esses zeros de f'(x). Agora, façamos uma afirmação geral.

Seja f ( x ) com n zeros, contando multiplicidades, em um intervalo; mostramos que f'(x) tem pelo menos ( n 1 ) zeros no intervalo, contando multiplicidades. Se os zeros de f ( x ) são x 1 < x 2 < < x r , com respectivas multiplicidades m 1 , m 2 , , m r , m 1 + m 2 + + m r = n , então x 1 é um zero de f'(x) de multiplicidade m 1 1 , x 2 um zero de multiplicidade m 2 1 , e assim por diante. (Zeros de multiplicidade 0 são simplesmente não zeros.) Além disso, f'(x) tem um zero entre x 1 e x 2 , um entre x 2 e x 3 , etc., dando pelo menos r 1 zeros adicionais. O número total de zeros, contando multiplicidades, de f'(x) no intervalo é, portanto, pelo menos

\begin{aligned} (m_1 - 1) &+ (m_2 - 1) + \dots + (m_r - 1) + (r - 1)\\ &= (m_1 + m_2 + \dots + m_r) - r + (r - 1) \\ &=m_1+m_2+\cdots+m_r-1\\ &= n - 1 \cdot \end{aligned}

um teorema fundamental que nos levará a várias formas interessantes de aproximações.

Teorema 1

Suponha que f ( x ) tenha zeros em x 1 , x 2 , , x r , de respectivas multiplicidades pelo menos m 1 , m 2 , , m r . Seja x 0 tal que f ( x ) tem n = m 1 + m 2 + + m r derivadas no menor intervalo contendo x 0 , x 1 , x 2 , , x r . Seja

H ( x ) = ( x x 1 ) m 1 ( x x 2 ) m 2 ( x x r ) m r .

Então existe um ξ neste intervalo tal que

f ( x 0 ) = H ( x 0 ) n ! f ( n ) ( ξ )

PROVA. Suponha primeiro que x 0 é um de x 1 , x 2 , , x r . Então f ( x 0 ) = 0 = H ( x 0 ) , e podemos tomar ξ como qualquer ponto do intervalo. Resta considerar o caso em que x 0 é diferente de todos x 1 , x 2 , , x r . Então considere

g ( x ) = f ( x ) C H ( x ) ,

onde C é uma constante tal que g ( x 0 ) = 0 . Note que, para poder escolher tal C , é necessário que possamos resolver para C em f ( x 0 ) C H ( x 0 ) = 0 , ou simplesmente que H ( x 0 ) 0 , o que sabemos ser o caso. É claro que poderíamos escrever explicitamente o que é C , mas isso não é essencial.

Agora, g ( x ) tem x 1 como um zero de multiplicidade pelo menos m 1 , x 2 como um zero de multiplicidade pelo menos m 2 , , x r como um zero de multiplicidade pelo menos m r , e x 0 também é um zero de g ( x ) . Assim, g ( x ) tem pelo menos n + 1 zeros, contando multiplicidades, em nosso intervalo. Consequentemente, existe um ξ neste intervalo tal que g ( n ) ( ξ ) = 0 . Mas g ( n ) ( x ) = f ( n ) ( x ) C H ( n ) ( x ) . Observe que, quando expandimos H ( x ) , obtemos um polinômio com termo de maior grau x n e outros termos com expoentes menores. Mas, quando tomamos n derivadas, todos esses outros termos se tornam 0 e ( x n ) ( n ) = n ! Assim

g ( n ) ( x ) = f ( n ) ( x ) C n ! , 0 = g ( n ) ( ξ ) = f ( n ) ( ξ ) C n ! ,

e

C = f ( n ) ( ξ ) n !

Do fato de que g ( x 0 ) = 0 = f ( x 0 ) C H ( x 0 ) , vemos que

f ( x 0 ) = C H ( x 0 ) = f ( n ) ( ξ ) n ! H ( x 0 ) ,

o que era para ser demonstrado.

Como primeira aplicação, consideremos o problema de interpolar a partir de tabelas, por exemplo, a partir de tabelas de logaritmos. Se temos tabelas para a função f ( x ) e lemos f ( x 1 ) e f ( x 2 ) para x 1 < x 2 , queremos estimar o erro introduzido ao usar a reta que passa por ( x 1 , f ( x 1 ) ) e ( x 2 , f ( x 2 ) ) no lugar da curva para encontrar f ( x ) quando x 1 < x < x 2 . Seja y = a x + b a reta aproximadora. Então consideramos a função f ( x ) ( a x + b ) . Esta tem zeros de multiplicidade pelo menos 1 em x 1 e x 2 . Portanto, temos n = 2 , e

H ( x ) = ( x x 1 ) ( x x 2 )

Pelo nosso teorema,

f(x) - (ax + b) = \frac{(x-x_1)(x-x_2)}{2!} f''(\xi)

para algum ξ entre x 1 e x 2 . Agora, H ( x ) 0 , já que x x 1 e x x 2 têm sinais opostos, e a diferença entre f ( x ) e a x + b será, em valor absoluto, menor ou igual ao valor absoluto de H ( x ) em seu mínimo vezes o maior valor absoluto possível para f''(x) no intervalo, dividido por 2 . Mas H'(x) = (x - x_1) + (x - x_2), que é 0 no ponto médio x = 1 2 ( x 1 + x 2 ) . Este é o mínimo de H ( x ) , e o valor de H ( x ) lá é

( x 2 x 1 ) ( x 1 x 2 ) 4

Assim

|f(x) - (ax + b)| \le \frac{(x_2 - x_1)^2}{8} \cdot \text{Max } |f''(\xi)| \cdot \tag{*}

Para uma tabela de logaritmos, com x 1 = 1.01 , x 2 = 1.02 ,  o erro é no máximo

( .01 ) 2 8 Max  1.01 ξ 1.02 ( 1 ξ 2 ) = .0001 8 1 1.0201 < .00002

EXERCÍCIOS

Exercício 1.
  1. Enuncie uma regra da forma (*) para extrapolação (onde x está fora do intervalo [ x 1 , x 2 ] ).
Exercício 2.

log 70.0 = 4.24850 ; log 80.0 = 4.38203 . Extrapole para obter um valor para log 82.0 e estime seu erro.

Exercício 3.

log tan 1.00 = 0.1924 ; log tan 1.10 = 0.2933 . Estime o erro ao interpolar para qualquer ângulo entre 1.100 e 1.10.