Métodos numéricos/Exercícios computacionais

Alguns problemas computacionais.

1. Usando o Octave, calcule e explique os resultados das seguintes operações:

${\displaystyle 2^{1023}+2^{1023}-2^{971}=}$

${\displaystyle 2^{1023}-2^{971}+2^{1023}=}$

${\displaystyle 2^{1023}-2^{1023}=}$

${\displaystyle 2^{1023}-2^{1023}+0.1=}$

${\displaystyle 2^{1023}+0.1-2^{1023}=}$

2. O limite ${\displaystyle C=\lim _{n\to +\infty }c_{n}=0.57721566490\ldots }$ É chamada constante de Euler.

2.1 Escreva um programa que calcula ${\displaystyle C}$ com uma precisão de ${\displaystyle 10^{-6}}$ (Será que o consegue fazer?).

2.2 Verifique numericamente (para ${\displaystyle n\leq 10^{7}}$) que o erro em cada iteração satisfaz a relação

${\displaystyle \epsilon _{n}=c_{n}-C={\frac {1}{2n}}+O\left({\frac {1}{n^{2}}}\right)}$

Quando ${\displaystyle n\to +\infty }$.

Para encontrar as raízes de um polinómio ${\displaystyle p_{n}(x)=a_{0}+a_{1}x+a_{2}x^{2}+\ldots +a_{n}x^{n}}$, onde ${\displaystyle a_{1},a_{2},\dots ,a_{n}\in \mathbb {R} }$, pode-se desenvolver a factorização, onde ${\displaystyle z_{1},z_{2},\ldots ,z_{n}}$ são as raízes do polinómio,

${\displaystyle p_{n}(x)=a_{n}(x-z_{1})(x-z_{2})\times \ldots \times (x-z_{n})}$

Estabelecendo um sistema de equações não lineares com a forma

${\displaystyle a_{0}=a_{n}(-1)^{n}z_{1}\times \ldots \times z_{n}}$

${\displaystyle a_{1}=a_{n}(-1)^{n-1}\left(z_{1}z_{2}+z_{1}z_{3}+\ldots +z_{1}z_{n}+z_{2}z_{3}+z_{2}z_{4}+\ldots +z_{2}z_{n}+\ldots +z_{n}z_{1}+z_{n}z_{2}+\ldots +z_{n}z_{n-1}\right)}$

${\displaystyle \ldots }$

${\displaystyle a_{n-1}=-a_{n}\left(z_{1}+z_{2}+\ldots z_{n}\right)}$

Que tem uma única solução. Este processo leva a um método rápido e eficaz para se calcular todas as raízes de ${\displaystyle p_{n}(x)}$ se se aplicar o método de Newton à resolução deste sistema não linear.

1. Suponha que existem zeros complexos para um polinómio com coeficientes reais. Haverá possibilidade de convergência do método de Newton para a solução do sistema se considerar todas as aproximações iniciais reais? Por quê?

2. Para o caso de polinómios de grau três, com a forma ${\displaystyle p_{3}(x)=a_{3}x^{3}+a_{2}x^{2}+a_{1}x+a_{0}}$, escreva explicitamente o sistema não linear que deve resolver.

3. Aplique esse método para determinar aproximadamente as soluções de ${\displaystyle x^{3}+3x+1=0}$, após ter escolhido uma aproximação inicial para a solução do sistema anterior. Use como critério de paragem ${\displaystyle \vert \vert z^{(k+1)}-z^{(k)}\vert \vert \leq 10^{-3}}$.

Exercício sobre os polinómios de Berstein.

Esta página é um esboço de matemática. Ampliando-a você ajudará a melhorar o Wikilivros.