A velocidade de ascensão de um foguete em vôo vertical próximo a superfície terrestre pode ser aproximada pela seguinte expressão:
Na qual:
velocidade de exaustão relativa do foguete
massa do foguete ao ser lançado
taxa de consumo de combustível
aceleração gravitacional
tempo (medido a partir do lançamento)
Considerando , , e , determinar o instante em que . Empregar o método da bisseção para determinar uma raiz aproximada para a equação, com intervalo inicial e tolerância de . Explicite o erro relativo a cada iteração.
Passo 1
Eita cacete. Que negócio é esse de física? Achei que estávamos em cálculo numérico!
Haha, relaxa cara, essas letras todas são só pra te complicar! Vamos começar nos livrando delas substituindo todos os valores dados e igualando a velocidade a :
Melhorou um pouco né? Agora é só uma equação, nada de foguetes haha. Mas ainda não podemos resolver isso usando bisseção né? Temos que dar uma mexida:
Opa, agora sim! Vai dar trabalho, mas pelo menos sabemos fazer haha. Só faltou um detalhe! Que negócio é aquele de ? E de erro relativo a cada iteração?? Então aqui ele está usando outro critério de parada! Ele quer que o erro relativo de uma iteração pra outra, dado por:
Seja menor que . Bom agora que já esclarecemos isso vamos partir em frente.
Passo 2
Ta pronto? Bebe uma água e toma coragem que esse é chatinho:
Iteração :
Logo, nosso novo intervalo será formado por e . Temos então que:
Critério de Parada:
Como não temos um chute inicial continuamos as iterações.
Iteração :
Logo, nosso novo intervalo será formado por e . Temos então que:
Critério de Parada:
Como não alcançamos o critério de parada continuamos as iterações.
Iteração :
Logo, nosso novo intervalo será formado por e . Temos então que:
Critério de Parada:
Como não alcançamos o critério de parada continuamos as iterações.
Iteração :
Logo, nosso novo intervalo será formado por e . Temos então que:
Critério de Parada:
Como não alcançamos o critério de parada continuamos as iterações.
Iteração :
Logo, nosso novo intervalo será formado por e . Temos então que:
Critério de Parada:
Essa foi quase em haha. Mas como não alcançamos o critério de parada continuamos as iterações.
Iteração :
Logo, nosso novo intervalo será formado por e . Temos então que:
Critério de Parada:
Beleza! Nosso erro relativo está menor que a tolerância, então podemos parar por aqui e dar como resposta a aproximação para a raiz .
Resposta
Instante segundos.
Erro relativo por iteração: