MÉTODOS NUMÉRICOS

CADERNO DE EXERCÍCIOS

Nuno Alexandre Cruz e Maria do Rosário de Pinho

Faculdade de Engenharia da Universidade do Porto

L.EEC– Métodos Numéricos

2023-2024
Conteúdo

Prólogo i
1 Preliminares . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . i
2 Problema . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . i
2.1 Uso do Excel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii
2.2 Usar MATLAB ou OCTAVE . . . . . . . . . . . . . . . . . . . . . . . . . ii
3 Outros Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iv

1 Aproximações e erros 1
1 Exercícios (aula de papel e lápis) . . . . . . . . . . . . . . . . . . . . . . . . . . . 1

2 Equações Não Lineares 3

1 Introdução . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2 Métodos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

3 Normas e Sistemas de Equações Não Lineares 7

4 Sistemas de Equações Lineares 9

5 Método dos Mínimos Quadrados 11

6 Interpolação Polinomial 13

7 Integração Numérica 15

8 Equações Diferenciais Ordinárias 17

3
4 CONTEÚDO
Prólogo

Problema Numérico: Descrição concisa de uma relação funcional entre os dados do problema
(input data) e o resultado procurado (output data). Dados e resultados são quantidades em
número finito.
Algoritmo para resolver um problema numérico: descrição de operações com as quais é
possível transformar os dados do problema no resultado propurado.

1 Preliminares

1. Responda às seguintes questões (se necessário, faça uma pesquisa (por exemplo, na na
bibliografia) que o ajude a responder):

(a) Qual a necessidade de elaborar métodos numéricos?

(b) O que entende por métodos iterativos?
(c) O que entende por convergência e divergência de um método iterativo?

No que se segue, propomos alguns problemas matemáticos cujas soluções aproximadas poderão
ser obtidas usando o Excel ou o Matlab (ou versão similar, o Octave, que pode ser obtida
sem custos, da internet). Grande parte dos programas usados no Octave são idênticos aos do
Matlab. Pretende-se que os estudantes se familiarizem com estas ferramentas computacionais
experimentando alguns casos simples mas relevantes.

2 Problema

1. Considere a função F (x) = −0.1x2 + 2.4. Pretende-se encontrar x > 0 tal que x = F (x).

(a) Verifique que x = 2 resolve a equação.

(b) Considere a fórmula de recorrência xn+1 = F (xn ) e x0 = 1. Calcule x1 e x2 .

De seguida, ilustra-se como poderá usar o EXCEL ou MATLAB (Octave) para determinar valores
aproximados para determinar o x > 0 tal que x = F (x).

i
2.1 Uso do Excel

1. Considere o problema anterior. Numa folha EXCEL escreva em 3 colunas A, B e C a

partir da linha 1, o seguinte fazendo fill down para preencher a tabela até à linha 23:

linha / coluna A B C
1 Iteração xn xn+1
2 0 1 = −0.1 ∗ power(b2; 2) + 2.4
3 1 = c2 = −0.1 ∗ power(b3; 2) + 2.4
.. .. .. ..
. . . .
22 ... ... ...
2. Indique quais os valores de x3 , x10 e de x20 obtidos.

3. O que pode dizer sobre a convergência ou divergência deste método?

4. Utilize método e tabela idêntica (mas só com 8 linhas!!!) para resolver o mesmo problema,
agora com F (x) = −x2 + 2.4. Lembre-se que tem que modificar a definição de F na tabela.

5. Comente os resultados nos dois casos anteriores.

2.2 Usar MATLAB ou OCTAVE

1. Considere o problema anterior com F (x) = −0.1x2 + 2.4. Como anteriormente pretende-se
calcular x > 0 tal que x = F (x), usando a fórmula de recorrência e considerando x0 = 1.
Para tal, crie o ficheiro problema1.m (por exemplo!) como o que está em baixo:

%%CRIAR um Ficheiro para armazenar os dados

fileID = fopen(’semana1_exerc1.txt’,’w’);

%% Definir função F para a fórmula

de recorrência x_{n+1}=F(x_n)

f = @(x) 2.4-0.1*x^2;

%% Escreve a primeira linha da tabela a ser apresentada

disp(’ it. x_n x_{n+1} |x_{n+1}-x_n| ’)

%%Escreve a primeira linha na file semana1_exerc1.txt

fprintf(fileID,’%3s %11s %11s %11s\r\n’, ’ x_n ’, ’ x_{n+1} ’ , ’ |x_{n+1}-x_n|

%% Definir o n. max de iterações

ii
 itmax=20;
%% Definir um vetor onde armanezar
%% os valores das iterações e um vetor com |x_n+1-x_n|

x=[];
delta=[];

%% Definir o valor inicial x_0 e guardá-lo como x(1)

x_0=1; %Valor inicial

x(1)=x_0; %armaze o valor inicial como sendo o primeiro elemento do vetor x

%%%% MÉTODO ITERATIVO:

for k = 1:1:itmax+1;
%%A iteracao k=1 e aquela que
%% chamamos de iteração 0, ou seja, corresponde a n=0;
%% temos x(1)=x_0,

x(2)=x_1,x(3)=x_2, ..., x(k+1)=x_k

x(k+1)=f(x(k));

%% Calcula a diferença x_n+1

delta(k+1)=abs(x(k+1)-x(k));

%mostra os resultados dessa iteração

sprintf(’%6i %11.8f %11.8f %11.8f\r\n’, k-1, x(k), x(k+1),delta(k+1))

%Para imprimir resultados na file semana1_exerc1.txt

fprintf(fileID,’%6i %11.8f %11.8f %11.8f\r\n’, k-1, x(k), x(k+1),delta(k+1));

end
%%Se quer ver o que foi gravado na file semana1_exerc1.txt
%% retire o % do comando final na linha final

%type semana1_exerc1.txt

2. Indique quais os valores de x3 , x10 e de x20 obtidos e compare-os com os obtidos com o
Excel.

3. Considere agora F (x) = −x2 + 2.4. Modifique o ficheiro de forma adequada considerando
agora apenas 8 iterações. Liste os resultados obtidos. Verifique se os resultados são análogos

iii
 aos obtidos com o Excel. Comente a diferença dos resultados no caso em que F (x) =
−0.1x2 + 2.4 e F (x) = −x2 + 2.4.

3 Outros Problemas
1. (Traçar gráficos:) Trace o gráfico de f (x) = x e F (x) = −0.1x2 +2.4. Pode, por exemplo,
usar o ficheiro .m em Matlab como o da figura em baixo:

%%% Programa para fazer o GRÁFICO de uma função

%%%%Primeira coisa a fazer: Limpar a memoria e fechar janelas
clear all;
close all;
%% Fazer os gráficos de DUAS funcões, f e F, usando Matlab

%% PRIMEIRO: DEFINIR AS FUNÇÔES: f(x)=x e F(x)=-0.1x^+2.4

f=@(x) x;
F=@(x) -0.1*x.^2+2.4;
%% definir valores das abcissas; AQUI O X ESTÁ ENTRE 0 E 4
%%% CONSIDERA-SE VALORES DE X 0, 0.01, 0.02, 0.03.... ATÉ 4
x=0:0.01:4;

%%%%FAZER OS GRÁFICOS DAS DUAS FUNÇÕES NUMA SÓ IMAGEM

figure;
grid on;
xlabel=’x’;
ylabel=’y’;
yf=f(x);
plot(x,yf, ’b’);
hold on;
yF=F(x);
plot(x,yF,’r’)
legend(’y=x’,’y=-0.1*x^2+2.4’,’Location’,’southeast’);
title( ’Graph of y=x and y=F(x)’,’interpreter’,’latex’);

2. Pretende-se calcular o valor de x tal que x2 = 2.

(a) Considere a sucessão definida por recorrência da forma

xn+1 = F (xn )
x0 = 1.5
1 2 √
onde F (x) = (x + ). Mostre que lim xn = 2.
2 x n→+∞
(b) Elabore programas
√ de EXCEL e de MATLAB (ou OCTOVE) para calcular apro-
ximações de 2 usando a fórmula de recorrência anterior: usar como x0 os valores
x0 = −1, x0 = 1 , x0 = 1.5 e x0 = 3 e comente os resultados.

iv
 (c) Repita as alíneas anteriores usando a fórmula de recorrência com F (x) = 2x.
(d) Comente os resultados.

3. Suponha que se quer calcular o valor da função p(x) = a1 x3 + a2 x2 + a3 x + a4 num dado

ponto x = x0 . Quais deverão ser os dados e qual deverá ser o resultado desejado?

4. Seja x0 = 2. Pretende-se calcular p(x0 ) onde p(x) = x3 + 2x2 + 3x − 1. Verifique que

p(x0 ) = ((x0 + 2)x0 + 3)x0 − 1.

Será que a partir desta observação consegue elaborar um algoritmo para calcular p(2)?

v
vi
Capítulo 1

Aproximações e erros

1 Exercícios (aula de papel e lápis)

1. Considere
9.004, 9.000, 0.049999, 617500, 618500, 1.0000, 4.352 × 10−4
(a) Arredonde para 3 dígitos significativos
(b) Supondo os números dados com todos os algarismos significativos, determine erro
máximo absoluto e relativo de cada.
x∗  0
9.004
9.000
0.049999
617500
618500
1.0000
4.352 × 10−4
(c) Arredonde os números para os algarismos significativos correspondentes aos erros in-
dicados.
x  0 arredondado as
9.004 – 5× 10−4
9.000 – 5 × 10−4
0.049999 — 5 × 10−4
617500 – 5 × 10−3
618500 1.2 × 102 –
1.0000 5 × 10−3 –
4.352 × 10−4 7.8 × 10−7 –
2. Considere que mediu o valor da tensão elétrica num circuito e que o valor medido foi
v̄ = 630 mV. Sabendo que o aparelho de medida tem um erro relativo máximo de 1%,
indique o intervalo de valores que contém o valor exato.

1
 3. Considere o valor aproximado de x̄ = 82.30.

(a) Supondo que todos os algarismos representados são signicativos, determine majorantes
do tipo 5 × 10−k para os erros absoluto e relativo.
(b) Arredonde x para os algarismos signicativos correspondentes a um erro máximo ab-
soluto de 0.5.
(c) Arredonde o número x para os dígitos signicativos correspondentes a um erro máximo
relativo de 3%.

4. Sabendo que x∗1 e x∗2 são valores aproximados de x1 e x2 com erros relativos r1 e r2 , verifique
x∗
que x∗1 x∗2 tem erro relativo r1 + r2 . Verifique ainda que o erro relativo de x∗1 é r1 + r2 .
2

5. Considere x1 = 0.2344 ± 0.5 × 10−4 e x2 = 0.2343 ± 0.5 × 10−4 . Calcule o valor aproximado
para x1 − x2 e o erro absoluto desta diferença.

6. Calcule um valor aproximado de y = sin(x) e o correspondente erro máximo absoluto

quando x∗ ≈ 0.57 (i.e., x ∈ [0.57 − 0.005, 0.57 + 0.005]).
R1 R2
7. Considere a resistência equivalente do paralelo de R1 e R2 , isto é, .
R1 + R2
(a) Sendo R1 = 1 kΩ e R2 = 2.2 kΩ tolerâncias de 5% e 10%, respectivamente, determine
um valor aproximado do paralelo e a sua tolerância (erro relativo máximo).
(b) Considere agora que R1 = 1 kΩ com tolerância de 10%. Qual deve ser a tolerância
máxima de R2 = 1 kΩ para que a tolerância do paralelo não seja superior a 7.5%?
π
8. Considere a função y = 2 sin(x1 ) cos(x2 ), com x1 , x2 = π4 ± 12 . Indique a expressão para o
valor aproximado de y e do seu erro máximo absoluto, no formato y = ȳ ± y .

9. Sabe-se que sempre que o erro relativo é multiplicado por 10, perde-se um algarismo sig-
nificativo. Verifique se há ou não perda de algarismos significativos no cálculo de y = ex
quando x ∈ [−1, 0].

10. Considere que a variável w é obtida a partir de x, y, z pela relação w = 5xy + 2z, onde os
valores x, y, z ∈ [0, 3]. Determine o erro máximo absoluto no cálculo de w, w , a partir dos
erros máximos absolutos de x, y, z, respetivamente x , y , z .

11. Considere que a variável w é obtida a partir de x, y, z pela relação w = 3xy + 4z, onde os
valores x, y, z ∈ [0, 2]. Determine o erro máximo absoluto no cálculo de w, w , a partir dos
erros máximos absolutos de x, y, z, respetivamente x , y , z .

12. Considere uma máquina de calcular que calcula os valores da função f (x) = cos(x) (lembrar
x2k
que a função f tem representação em série de Taylor, f (x) = nk=0 (−1)k (2k)!
P
) usando o
1
polinómio de Taylor de grau 4 centrado no ponto 0, P4 (x), quando x ∈ [0, 2 ]. Indique o
erro de truncatura da série.

2
Capítulo 2

Equações Não Lineares

1 Introdução

1. Determine os números de Rolle da função f (x) = ex − ex − e, no intervalo [0, 2]. E quais

os números são os números de Rolle da mesma função em [−1, 2]?

2. Determine os números de Rolle da função x(t) = −t2 − 2, em R.

3. Verifique que as funções dadas têm um único zero nos intervalos indicados.

(a) f (x) = 10e−x − 5 no intervalo [0, 4].

√
(b) x + 3 = x no intervalo [1, 3].

4. Obtenha intervalos de largura não superior a 0.5, cada um contendo uma só solução da
equação x3 + 2x2 − 3x − 1 = 0.

5. Obtenha intervalos de largura não superior a 0.5, contendo um ponto máximo ou de mínimo
da função dada:

(a) f (x) = x3 + 2x2 − x − 2.

(b) g(x) = x4 − 2x2 − 1.

2 Métodos

1. Considere a função f (x) = x3 + ex .

(a) Determine um intervalo de largura 1 contendo o único zero desta função.

(b) Efetue 2 iterações do método das bissecções sucessivas para obter uma aproxim-
ção do zero desta função.
(c) Apresente um majorante do erro máximo absoluto da aproximação obtida na alínea
anteior.

3
 (d) Determine o número de iterações necessárias para que o erro máximo da aproximação
fornecida pelo método das bisseções sucessivas seja 10−5 .
(e) Usando métodos computacionais (Excel ou Matlab/Octave, por exemplo), confirme o
resultado anterior.

2. (Extra.) Considere um objecto de massa m que é largado a uma certa altitude. Sabendo
que no instante t = 0, a sua velocidade é 0, sabe-se que a velocidade do objecto no instante
gm c
t é dada pela equação v(t) = (1 − e−( m )t ) onde g = 9.8, c é o coeficiente de atrito, e m
c
a massa do corpo.

(a) Pretende-se saber qual o valor de c para que um objecto com m = 68.1, atinja uma
velocidade de 40 para t = 10. Verifique que esse valor está entre 12 e 16.
(b) Calcule um valor aproximado para c calculado pelo método das bisseções sucessivas
com erro inferior a 5 × 10−2 (use métodos computacionais).

3. Considere a equação f (x) = x − cos πx

4 no intervalo [0, 1].

(a) Verifique que esta função tem um só zero nesse intervalo.

(b) Pretende-se aplicar o método iterativo simples, para determinar um valor apro-
ximado da solução da equação, com erro inferior a um valor dado. Verifique se as
condições que garantem a convergência do método iterativo simples são satisfeitas.
(c) Aplique o método iterativo simples para determinar um valor aproximado da
solução da equação, com erro inferior 10−3 .
(d) O mesmo exercício mas com um erro inferior a 5 × 10−5 .

4. Considere a equação (x + 2)2 − 2ex = 0.

(a) Obtenha intervalos de largura não superior a 0.5, cada um contendo uma só solução
da equação dada.
(b) Que condições devem ser satisfeitas para garantir a convergência do método iterativo
simples a este problema?
(c) Elabore um programa computacional em Matlab/Octave (ou use Excel) para imple-
mentar este método na resolução desse exercício.
(d) Aplicando o método iterativo simples, determine o valor aproximado da maior
solução da equação ao fim de 10 iterações quando x0 = 1. Qual é o majorante do erro
dessa aproximação?
(e) Determine o valor aproximado da maior solução da equação dada quando |xn −xn−1 | ≤
5 × 10−5 , onde xn é a estimativa mais recente calculada na alínea anterior.

5. A tensão num circuito RC tem como expressão v(t) = 10e−t − 5, para t ≥ 0. Verifique
que existe s ∈ [0, 1] tal que v(s) = 0. Aplicando o método de Newton, determine um
valor aproximado de com erro inferior a 10−5 . (use o cálculo automático para resolver este
problema). Sugestão: Verifique se as condições que garantem a convergência do método
de Newton são satisfeita em x0 = 0.

6. Considere a equação x3 + 2x2 − 3x − 1 = 0.

4
 (a) Pretende-se calcular a menor solução da equação usando o método de Newton. De-
termine um valor valor inicial para esse método de forma a garantir a convergência e
escreva expressões para o erro na iteração k.
(b) Utilizando o método de Newton, faça 10 iterações para determinar a menor solução
da equação.
(c) Determine o erro absoluto da solução obtida na alínea anterior.

7. (Extra.) Pretende-se calcular um zero da função f (que tem todas as derivadas) e sabe-
se que o único zero dessa função tem multiplicidade 2. Como poderia calcular esse zero
usando o Método de Newton? Será que pode generalizar esse método para o caso em que
s é zero de multiplicidade superior a 2?
Sugestão: Seja s o zero dessa função. Escreva a fórmula de Taylor centrada em s para f e
f 0 e verifique que g(x) = f (x)/f 0 (x) tem um único zero em s.

5
6
Capítulo 3

Normas e Sistemas de Equações Não

Lineares
 
0.11
1. Sabendo que x = −0.25 calcule ||x||1 .
0.80
 
−0.31
2. Sabendo que x =  0.05  calcule ||x||∞ .
0.90
 
0.36
3. Sabendo que x =  0.07  calcule ||x||1 .
−0.70
4. Determine a norma (induzida) infinita da seguinte matriz:
 
−0.72 3.13 0
 1.1 0.54 −1.3
1.53 1.23 0.26

5. Determine a norma (induzida) infinita da seguinte matriz:

 
−0.74 3.18 0
 1.2 0.51 −1.3
0.53 2.21 0.17

6. Determine a norma (induzida) infinita da seguinte matriz:

 
−1.74 1.18 0
 2.2 0.51 −1.3
0.23 1.61 0.57

7. O sistema de equações não lineares

 2
x1 − 10x1 + x22 + 9 = 0
x1 x22 + x1 − 10x2 + 8 = 0

7
 pode ser reescrito da seguinte forma, adequada para aplicação do método iterativo simples

x21 + x22 + 9


 1
 x = g (x ,
1 1 2 x ) =
10


2
 x2 = g2 (x1 , x2 ) = x1 x2 + x1 + 8



10
(a) Verifique que a função G = [g1 , g2 ]T de R2 em R2 tem um único ponto fixo no conjunto
D = {(x1 , x2 ) : 0 ≤ x1 , x2 ≤ 1.5}.
(b) Aplique o método iterativo simples para determinar a solução do sistema de equações
no conjunto D. Usando o Excel ou Matlab (ou Octave ou outro métos) termine o
método assim que a norma ∞ (i.e. a maior das componentes do vector) do erro de
aproximação seja inferior a 5 × 10−6 .

8. Utilize o método iterativo simples para determinar a solução do sistema de equações


10x1 + x2 − 5 cos(x1 ) = 0
3x2 + x1 x2 − 3 = 0

Parta do ponto (0, 0) e termine assim que a norma ∞ da diferença entre duas estimativas
consecutivas seja inferior a 5 × 10−5 .

9. Considere o sistema de equações


3x1 + sin(x1 + x2 ) − 5 = 0
x1 + 4x2 + cos(x1 ) + sin(x2 ) + 1 = 0

Aplique o método de Newton para determinar a solução deste problema. Parta do ponto
(0.5, 0.5) e termine assim que a norma ∞ da diferença entre duas estimativas consecutivas
seja inferior a 10−5 .

10. Considere a circunferência definida por x2 +y 2 −1 = 0 e a elipse definida por x2 +3y 2 −2 = 0.

Utilizando o método de Newton, determine o ponto de interseção das duas curvas situado
no 1o quadrante. Termine a aplicação do método assim que as variações nas estimativas
de x e y entre duas iterações sejam inferiores a 5 × 10−6 .

8
Capítulo 4

Sistemas de Equações Lineares

1. Considere o sistema de equações


 0.22x1 + 0.3x2 − 0.1x3 = 0
0.2x1 + 0.2x3 = 1.0
1.5x1 + 2.1x2 = 0.9


Resolva-o por eliminação gaussiana utilizando pivotação parcial, arredondando os resulta-

dos finais para 3 algarismos significativos.

2. Repita o exercício anterior, mas usando agora pivotação total.

3. Considere o seguinte sistema de equações.

    
−0.71 3.2 0 x1 1.5
 1.2 0.51 −2.3 x2  = 0.94
0.53 0.21 0.17 x3 1.06

(a) Reescreva-o numa forma adequada à aplicação do método de Jacobi, garantindo a

respectiva convergência.
(b) Partindo de x1 = x2 = x3 = 0, efectue 3 iterações do método iterativo de Jacobi para
determinar uma solução aproximada deste sistema de equações.
(c) Calcule o resíduo da solução aproximada determinada em (a).
(d) Estime o erro da solução aproximada determinada em (a).

4. Considere o sistema de equações


 −x1 + 10x2 = 9
−2x1 + 10x3 = 6
10x1 − x2 − 2x3 = 7


(a) Reescreva-o numa forma adequada à aplicação do método de Jacobi, garantindo a

respectiva convergência.

9
 (b) Faça 3 iterações pelo método de Jacobi a partir do ponto (0, 0, 0), para obter uma
solução aproximada do sistema.
(c) Calcule o resíduo da solução aproximada determinada em (b).
(d) Indique um majorante para a norma do erro da 3a iteração.

5. Considere o sistema de equações


 x1 + 8x2 − 3x3 = 12
x1 + 4x2 − 10x3 = 15
3x1 + x2 = 4


(a) Reescreva-o numa forma adequada à aplicação do método de Gauss-Seidel.

(b) Determine um solução aproximada deste sistema de por aplicação do método de Gauss-
Seidel. Inicie as iterações com x1 = x2 = x3 = 0 e termine a aplicação do método
assim que kx(k) − x(k−1) k∞ ≤ 10−3 .
(c) Calcule o resíduo da solução aproximada calculada em (b).

6. Utilizando eliminação de Gauss-Jordan com pivotação parcial, determine a inversa da se-

guinte matriz.  
−0.71 3.2 0
 1.2 0.51 −2.3
0.53 0.21 0.17

10
Capítulo 5
MIEEC: Mestrado Integrado em Engenharia Electrotécnica e de Computadores
Métodos Numéricos
Método dos Mínimos Quadrados 2013/2014

— Método dos Mı́nimos Quadrados —

Nota: Todos os problemas neste capítulo que envolvem cálculo numérico podem e devem ser
PROBLEMAS
resolvidos usando o Excel ou o Matlab ou o Octave ou .....

1 — Determine a aproximação dos mı́nimos quadrados aos pontos

1. Considere a tabela de pontos:
xx 11 33 44 66 88 99 11
11 14
14
yy 11 22 44 44 55 77 88 99

por:Determine a aproximação dos mínimos quadrados aos pontos por

(a)
i. uma
(a) Uma reta g(x) = a0 + a1 x;
recta.
ii. uma parábola g(x) = a0 + a1 x + a2 x2 ;
(b) Uma parábola.
iii. uma cúbica g(x) = a0 + a1 x + a2 x2 + a3 x3 ;
(c) Uma cúbica.
iv. uma reta x = g(y) que minimize o erro em x.
(d) Uma recta x = g(y) que minimize o erro em x.
(b) Para cada um dos três primeiros casos anteriores, determine o valor a soma dos desvios
Xn
quadrados |yi − g(xi )|2 e comente.
i=1
2 — (c)
ParaFaça
medirosa gráficos
capacidade
das(C) de um gcondensador,
funções calculados usou-se
nos trêsumprimeiros
circuito RC com
caso daum interruptor.
alinéa a) e dos
No instante
pontos da tabela. Comente as diferentes aproximações. Mais ainda, em cada) no
t = 0 abriu-se o interruptor, tendo-se registado os seguinte valores da tensão (V caso
indique qual o desvio máximo (i.e., o valor de maxi=1,...,n |yi − g(xi )|.
condensador:

2. Para medir a capacidade

t(s) 0 C de
0.1 um
0.2condensador, usou-se
0.3 0.4 0.5 0.6 um
0.7 circuito RC 1.0
0.8 0.9 com um interrup-
tor. V (volt) 10 7.8 6.2 4.9 3.5 2.7 2.1 1.6 1.3 1.0 0.7

t=0

R
U C V
(1M )

t
Considere que v(t) = v0 e RC , t 0 e que R = 1 ⇥ 106 (⌦). Aplique o critério dos mı́nimos
11
quadrados para estimar um valor para C. Considere que v(0) também está sujeito a erro e que
por isso v0 é desconhecido “a priori”.
 No instante t = 0 abriu-se o interruptor, tendo-se registado os seguinte valores da tensão
V no condensador:

t (s) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
V (volt) 10 7.8 6.2 4.9 3.5 2.7 2.1 1.6 1.3 1.0 0.7
t
Considere que v(t) = v0 e− RC , t ≥ 0 e que R = 1 × 106 Ω.

(a) Aplique o método dos mínimos quadrados para estimar um valor para C. Considere
que v(0) também está sujeito a erro e que por isso o v0 é desconhecido “a priori”.
(b) Usando o resultado anterior estime o valor de V no instante t = 0.25.

3. Considere a tabela de pontos

x 0 1.0 2.3 4.0 5.1 6.0 6.5 7.0 8.1 9.0

y 0.2 0.8 2.5 4.3 3.0 5.0 3.5 2.4 1.3 2.0
π
Determine a função da forma g(x) = a+b sin( 10 x) que melhor aproxima a tabela de pontos
no sentido dos mínimos quadrados.

4. Sabe-se que os pontos da tabela

x 0.5 1 2 3 4
y 10.4 5.5 3.3 2.4 2
 √ 2
a+ x
pode ser modelizado por uma função do tipo g(x) = √ . Use uma transformação
b x
que lhe permita obter de g uma função do tipo ψ(x) = c1 φ1 (x) + . . . + ck φk (x). Use essa
função para obter a e b que garante que o g correspondente é a melhor aproximação dos
pontos da tabela no sentido dos mínimos quadrados.

Outros exercícios:

5. Determine a melhor aproximação da forma y = a sin(x) + b cos(x), no sentido dos mínimos

quadrados, dos pontos da tabela.
π 3π
x 0 2 π 2
y −1.3 1.55 1.2 −1.45

1
6. Pretende-se determinar uma aproximação da forma y = ax+b aos pontos da tabela.

x 1.0 1.2 1.4 1.6 1.8 2

y 0.219 0.188 0.173 0.154 0.144 0.136

Utilizando uma mudança de variável em y determine a aproximação por aplicação do

método dos mínimos quadrados.

12
Capítulo 6

Interpolação Polinomial

1. Considere a seguinte tabela de pontos.

x 0.2 0.6 1.0 1.4

y −0.8 −0.2 −0.4 0.1

(a) Determine o polinómio interpolador destes pontos utilizando a forma de Lagrange.

(b) Determine o polinómio interpolador destes pontos, na forma de Newton.
(c) Repita a alínea anterior usando a tabela de diferenças divididas e a tabela de diferenças
finitas.

2. Considere a tabela de diferenças divididas.

i x y[·] y[·, ·] y[·, ·, ·]
(a) Determine os valores y1 , y2 e z. 0 3 2
(b) Obtenha o polinómio interpolador. 4
1 6 y1 1
(c) Considere o novo nó x = 4. Determine o valor nodal que z
lhe corresponde, sabendo que o polinómio interpolador 2 7 y2
dos 4 nós é de grau 2.

3. Considere a função f (x) = sin(x) + cos(x).

π π
(a) Determine o polinómio de menor grau que interpola esta função nos pontos 0, 6, 3 e
π
2.
(b) Determine, para o intervalo [0, π2 ], um majorante do erro de interpolação da função f
pelo polinómio obtido na alínea anterior.

4. Considere a função f (x) = ln(x) que se pretende aproximar por um polinómio de grau 8
que interpola esta função em 9 pontos igualmente espaçados no intervalo [1, 2]. Obtenha
um majorante do erro de aproximação da função por este polinómio interpolador.

13
5. A função f (x) = e−x − x tem um zero no intervalo [0, 1].

(a) Utilizando 4 valores de x igualmente espaçados neste intervalo e interpolação inversa

de f (x), obtenha um valor aproximado do zero desta função.
(b) Calcule o valor de f (x) no ponto obtido.
(c) Compare o valor do zero obtido por interpolação inversa com o valor calculado pelo
método das bisseções sucessivas, com erro inferior a 10−3 .

14
Capítulo 7

Integração Numérica

Z 1.2
2
1. Considere o integral definido ex dx.
0

(a) Utilizando a regra dos trapézios e um espaçamento h = 0.2 entre pontos consecutivos,
calcule um valor aproximado deste integral.
(b) Obtenha um majorante para o erro cometido na alínea anterior.

2. Considere os seguintes pontos de uma função y = f (x).

x −2 −1 0 1 2
y −0.32 −0.15 0.21 0.45 0.09
Z 2
(a) Determine o valor aproximado de f (x)dx pela regra de Simpson.
−2
(b) Calcule o erro de arredondamento do valor obtido em (a), considerando que os valores
y da tabela estão representados com todos os algarismos corretos.

Z 3
3. Considere o integral ln(x)dx.
1

(a) Determine o número de intervalos necessários ao seu cálculo com erro não superior a
5 × 10−3 usando:
i. a regra dos trapézios;
ii. a regra de Simpson.
(b) Considerando um espaçamento h = 0.2, obtenha um valor aproximado do integral,
usando o método que garante um erro menor.
(c) Compare o valor obtido na alínea anterior com o valor exato do integral.
R
(Nota: ln(x)dx = x ln(x) − x)

15
 2
1
Z
4. Considere o integral definido dx.
1 x

(a) Determine um valor aproximado do integral usando a regra de Simpson com h = 0.25.
(b) Estime uma majorante do erro cometido.
(c) Calcule o valor exato do integral e determine o erro efetivamente cometido na apro-
ximação da alínea (a).
1
dx
Z
5. Usando o método de Romberg obtenha um valor aproximado de , considerando
0 1 + x2
sucessivamente 1, 2 e 4 subintervalos para a regra dos trapézios.

16
Capítulo 8

Equações Diferenciais Ordinárias

1. Considere o problema de valor inicial

x0 (t) = t − x2 (t), t ∈ [0, 1]
x(0) = 1
(a) Obtenha uma solução aproximada do problema usando o método de Euler progressivo
com passo 0.05.
(b) Determine uma solução aproximada do problema pelo método de Taylor de ordem 2
com um passo 0.05.
2. Considere a equação diferencial x0 (t) = t2 − 1+x12 (t) , com condição inicial x(1) = 2. Deter-
mine uma solução aproximada no intervalo [1, 3], com passo 0.1, utilizando:
(a) O método de Euler modificado.
(b) O método de Runge-Kutta de 4a ordem.
3. Considere o problema de valor inicial
x0 (t) + 2x(t) = 6t + 1, t ∈ [0, 1]
x(0) = 0
(a) Obtenha uma solução aproximada do problema usando os métodos de Euler progres-
sivo e regressivo com passo 0.05.
(b) Obtenha uma solução aproximada do problema, usando os métodos de Euler modifi-
cado e Runge-Kutta de 4a ordem, com passo 0.1.
(c) Atendendo a que a solução exata é x(t) = e−2t + 3t − 1 (confirme!), determine o erro
ao fim de cada passo nos métodos anteriores.
4. Considere o sistema de equações diferencias
 0
x = (x − y)2
y 0 = y − 2x
com condições iniciais x(0) = 0, y(0) = 1. Determine uma solução aproximada da solução
do sistema no intervalo [0, 2] usando:

17
 (a) o método de Euler com passo 0.1.
(b) o método de Euler modificado com passo 0.2.

5. Determine uma solução aproximada de

q 00 (t) + 3q 0 (t) + 2q(t)|q(t)| = 4

q(0) = 0, q 0 (0) = 1

no intervalo [0, 3] usando o método de Euler com passo 0.1.

6. Determine uma solução aproximada de

y 00 (t) + y 3 (t) = 0
y(0) = 1, y 0 (0) = 0

no intervalo [0, 2] usando o método de Taylor de ordem 2 com passo 0.1.

18