Processing math: 100%

User Tools

Site Tools


calculo1:estabilidadedecomputacao

Erro de Aproximação, Estabilidade de computação

Problema: Considere a função f(x)=5x+3. Qual é a precisão na entrada desta função para que a saída tenha erro menor do que 104?

Observem que o problema acima pode surgir num problema aplicado. por exemplo quando aumentamos as dimensões de uma figura na tela de computador (5 vezes maior) e gostaríamos de que a imagem resultante não seja de má qualidade. Somar com 3 pode representar uma translação que não implicaria nada na resolução!

Bom, se inserirmos x1,x2 na função. Queremos saber quão pequena precisamos escolher |x1x2| tal que |f(x1)f(x2)|=|(5x1+3)(5x2+3)|104 ou seja

|5(x1x2)|104. portanto é claro que se tivermos |x1x2|104/5 então teremos a precisão desejada na saída da função.

Exemplo:

Considere a função quadrática f(x)=x2. Queremos saber qual é a precisão necessária de x para que o erro de aproximação de seu quadrado seja menor do que 103.

Neste exemplo não podemos dar uma resposta independente do x. Vamos analisar melhor: sejam x1,x2 duas entradas da função e queremos determinar δ tal que se |x1x2|δ então

|x21x22|103 ou equivalentemente |x1+x2||x1x2|103.

Observe que não existe algum δ universal para satisfazer a desigualdade acima. Suponhamos que por absurdo exista tal δ. Então, escolhemos x1,x2 dois números maiores que 1δ e que |x1x2|δ/2. por exemplo, x1=1/δ,x2=1/δ+δ/2. Assim, |x1x2|=δ/2<δ e |x1+x2|=2/δ+δ/2 e portanto

|x21x22|=(2δ+δ2)(δ2)=1+δ24. ou seja nunca |x21x22| será menor do que 103!

O que podemos fazer????

O ponto é que qualquer erro de aproximação de x acarreta um erro “perto” de dobro deste erro multiplicado por x. Vamos ver isto usando álgerba! (assim que álgebra entra no mundo de aproximação!)

(x+h)2x2=h2+2hx ou seja o erro de calcular o quadrado de x, usando (x+h)2 é h2+2hx. Observe que se h é um número pequeno (<1) então h2 é menor ainda (os físicos as vezes até ignoram este termo!!!!! nós também, mas vamos fingir que não! :-X). Porém 2hx pode ser um número grande dependendo do x.

Considere um problema específico: seja A=15,a1a2a3 que os dígitos an são calculáveis com certo custo. Queremos saber quantos dígitos são necessários para que o erro de computação de A2 seja menor do que 103.

Ou seja queremos achar n tal que se An=15,a1a2a3an tenhamos certeza de que |A2A2n|103.

Observe que |AAn|10n e já que A e An são menor do que 16,

|A2A2n|=|AAn||A+An|32×10n. Portanto basta que

32×10n103 e assim 10n32×103 e basta escolher n5.

Exemplo:

Considere a função f:{xR:x>0}R com regra f(x)=1x. Novamente gostariamos de saber qual é a precisão necessária para determinar x tal que o erro no cálculo de f(x) seja menor do que 103. Neste caso, como exemplo anterior, a precisão necessária vai depender do x.

Se a diferença entre dois números grandes for 10n então a diferença entre seus inversos é pequeno. Porém se dois números pequenos tiverem diferença 10n suas inversas podem ser muito diferente. vamos ver o que queremos dizer com isto:

Seja x1=10,x2=10+103 então f(x1)f(x2) é aproximadamente 0,000009999. Enquanto isto se x1=103,x2=103+103 então f(x1)f(x2)=500!

Olhem o gráfico da função f(x) para visualizar geometricamente o que descrevemos no parágrafo anterior.

O gráfica mostra que escolhendo dois números com a mesma distância resulta distantes diferentes entre suas inversas. Assim que, para números pequenos precisamos de ter precisão mais alta para ter um erro computacional fixo.

Por exemplo suponhamos que A=0,02a3a4a5 queremos determinar δ de tal maneira que se |AA|δ então |1A1A|102. Observe que

|1A1A|=|AA|AA e essa última fração mostra o cuidado que devemos ter na hora de escolher δ! se escolhermos o δ=A podemos escolher A muito perto de zero e assim a fração |AA|AA vai ser enorme! Por isto primeiro escolhemos um número δ1 menor do que A, tipo 0,02. Então se |AA|<δ1 já que A0,02 concluímos que A>0,01. Sendo assim, se |AA|<102 obteremos

|1A1A|<|AA|0,01×0,02=1042|AA|

Agora queremos que o lado direito da desigualdade acima seja menor do que 102 e portanto basta que |AA|<2×106.

8-O Afinal das contas qual será a escolha da precisão δ?

Bom, basta escolhermos δ=2×106 pois assim δ<δ1 e automaticamente todos os cálculos acima estarão corretos.

calculo1/estabilidadedecomputacao.txt · Last modified: 2022/03/08 14:22 by 127.0.0.1