Exercício 1 A constante matemática de Euler-Mascheroni denotada por γ é dada por
1 1 1
γ = lim 1 + + +. . + − ln (𝑛)
→ 2 3 𝑛
Nota: γ =0,57721566 (primeiras 8 casas decimais).
a) Calcule γ para n=100 e n=1000. Nota: use format long para apresentar mais casas
decimais.
Exercício 2 Carregue o ficheiro [Link] que contem a velocidade do vento (m/s),
temperatura e pressão (médias de ½ em ½ hora) para cada um dos 30 dias do mês de
setembro de 2020.
Usando apenas os dados relativos à velocidade do vento calcule para cada instante
horário:
a) Média diária da velocidade do vento.
b) Média horária da velocidade do vento e aplique uma regressão polinomial com
um polinómio de ordem 6 a esta grandeza em função da hora do dia.
c) Represente graficamente, num gráfico de superfície e cor interpolada, a
velocidade do vento, e extraia conclusões sobre esta grandeza neste mês. Legende
todos os eixos de acordo com as variáveis correspondentes
d) Para clarificar as conclusões do ponto anterior faça um gráfico pseudo-color
(pcolor), sobreponha um de contorno e legende as curvas de contorno. Dê nomes
a todos os eixos.
e) Calcule o dia e a hora em que foi registada a maior velocidade do vento.
Exercício 3: Método de MonteCarlo (determinação aproximada do zero de uma função)
Considere a função y = ln(x) − 0.5x +2, a qual apresenta um zero no intervalo entre
x=1 e x=3.
a) Gere um vetor xtrial com 1000 valores aleatórios com distribuição uniforme no
intervalo entre 1 e 3.
b) Calcule ytrial para o vetor xtrial usando a função descrita e faça o respetivo
gráfico com pontos a vermelho.
c) Determine o menor valor absoluto de ytrial e o respetivo índice. Com o índice
determine a melhor aproximação, xnulo, do zero da fução. Coloque um círculo
preto nesse ponto.
Exercício 4: Método de Monte Carlo (integrais definidos)
Determine a área compreendida entre a função 𝑓(𝑥) = 𝑥 e o eixo dos xx entre x=0 e
x=1.
a) Para acompanhar melhor a solução faça o gráfico da função entre x=0 e x=1 e
ajuste os limites de ambos eixos para o intervalo [0 1].
1
b) Gere coordenadas aleatórias com distribuição uniforme para N=1E6 pontos
(xrand, yrand) dentro da área do gráfico e represente-os com pontos a vermelho.
c) Identifique os pontos que “caíram” entre a função f(x) e o eixo dos xx e
represente-com ‘.’ a verde. Calcule ainda a fração, Percent, do total que estes
pontos representam.
d) Calcule a área pretendida que será a área, Area, da figura multiplicada pela
fração de pontos determinada na alínea anterior. Nota: o valor exato até a 5ª casa
decimal é 0.78343.
Exercício 5 Resolução de sistemas lineares de equações
Num jogo de dardos a pontuação de cada seta depende da zona do alvo atingido:
circular A (central); área B (um anel concêntrico e confinante com A) e um anel C
(concêntrico e confinante com B).
a) Foram obtidos os seguintes resultados e pontuação de 3 equipas após o
lançamento de 18 dardos. Escreva as matrizes A e B, do sistema de equações
descrito por 𝑨𝒙𝑿 = 𝑩, que permitem determinar qual a pontuação atribuída num
acerto em cada uma das áreas [xA; xB; xC].
Área A Área B Área C Pontos
Equipa 1 3 5 10 40
Equipa 2 1 6 11 34
Equipa 3 2 3 13 32
b) Resolva a equação usando a matriz inversa de A não sem antes verificar que ela,
de facto, tem inversa.
c) Verifique a solução calculando 𝑨𝒙𝑿 − 𝑩 e comente.
d) Resolva a equação usando o método da divisão esquerda e verifique a solução.
Comente.
Exercício 6: Os planos Poupança-Reforma têm benefícios fiscais com dedução no IRS.
Estes beneficios estão detalhados na seguinte tabela.
Idade Dedução Limite Investimento para dedução máxima
Menos de 35 anos 20% € 400 € 2000
De 35 a 50 anos 20% € 350 € 1750
Mais de 50 anos 20% € 300 € 1500
a) Crie agora uma script CalcDedIRS.m que lhe permita calcular a dedução do IRS
para qualquer idade do investidor e da quantia entregue. Implemente os
seguintes passos:
1. Leia, com uma instrução input(), a Idade do investidor e o Montante
investido
2. Comece utilizando as instruções if elseif …. else end, por escrever na
Command Window as strings “Idade inferior a 35 anos”, “Idade entre 35
2
e 50 anos” e “Idade superior a 50 anos” conforme a Idade. Teste a sua
solução.
3. Acrescente as instruções necessárias para calcular as respetivas deduções
(Deducao).
4. Usando a função fprintf() no fim da scripts, escreva numa mesma linha o
valor do Montante e a Deducao em formato vírgula flutuante com duas
casas decimais.
b) Transforme a script numa function.
Exercício 7 A velocidade de um objecto de massa m, sujeito a uma força constante F e
uma força de atrito é descrita pela seguinte equação diferencial:
( )
𝑚 = 𝐹(1 − 0.025𝑣(𝑡)) com m=1 kg, F=10 N e 𝑡 em segundos. A velocidade inicial é
𝑣(0) = 0 (m/s).
a) Defina os símbolos/funções simbólicas necessárias ao problema.
b) Escreva a equação diferencial e obtenha a solução geral SolGeral (use a função
dsolve()).
c) Escreva a equação da condição inicial e obtenha a solução particular SolPart.
represente-a graficamente usando a função fplot() entre t=0 e t=60 s.
d) Calcule a distância percorrida nos primeiros 60 s. Nota integre a velocidade da
solução particular entre 𝑡 = 0 e 𝑡 = 60 s.
Exercício 8 A constante matemática de Euler-Mascheroni denotada por γ é dada por
1 1 1
γ = lim 1 + + +. . + − ln (𝑛)
→ 2 3 𝑛
Nota: γ =0,57721566 (primeiras 8 casas decimais).
a) Usando um ciclo while determine o valor de 𝑛 necessário para conseguir um γ
com um erro inferior a 1 milésima. Nota: Use o valor de γ acima como valor exato e
use format long para apresentar mais casas decimais.
Exercício 9 Considere a função 𝑦 = 𝑥 − 𝑐𝑜𝑠(𝑥) que tem um zero algures entre 0 e 𝜋/2.
a) Faça o gráfico da função entre 0 e pi para verificar a constatação acima.
b) O cálculo do valor do zero consiste em resolver a equação 𝑥 = 𝑐𝑜𝑠(𝑥) que não tem
uma solução algébrica, mas pode estimá-la usando sucessivas aproximações.
Usando um ciclo for repetido 20 vezes e um primeiro valor 𝑥 = 1 calcule uma
seguinte aproximação 𝑥 = 𝑐𝑜𝑠(1) e assim sucessivamente. 𝑥 = 1, 𝑥 = cos(1) =
0.5403; 𝑥 = cos(0.5403) = 0.8576; 𝑥 = cos(0.8576) = 0. 6543, 𝑒𝑡𝑐
c) Modifique o código para guardar agora os valores num vetor Nulos e, uma vez
concluído o ciclo represente graficamente os sucessivos resultados.