Atividade 01

UNIVERSIDADE FEDERAL DE SÃO PAULO
Instituto de Ciências Ambientais, Quı́micas e Farmacêuticas

Departamento de Engenharia Quı́mica
CÁLCULO NUMÉRICO
Atividade 01 - 14/dezembro/2020
INSTRUÇÕES:
1. O programa deverá ter o seguinte cabeçalho:
/*******************************************************************
********** Atividade 1 - Cálculo Numérico **********
********** **********
********** Nome Completo **********
********** **********
*******************************************************************/
2. A atividade pode ser feita em dupla ou trio.
3. O programa deverá ser obrigatoriamente implementado no Scilab
4. Atividades atrasadas não serão recebidas.
5. Todo aluno deverá entregar o código e um relatório, tudo devidamente identificado. A

ausência de qualquer dos itens anteriores (inclusive identificação) poderá implicar na
invalidação da atividade ou na redução da nota final.
6. Data de entrega: até às 18:00 horas do dia 13/janeiro/2021 (via Moodle)
1
Curvas de Isoconcentração de Poluentes
1. Objetivo
O objetivo da atividade é apresentar um traçado das curvas de isoconcentração de um deter-
minado poluente na atmosfera. Os dados são concentrações medidas em várias localizações
e altitudes. Esses dados - coordenadas (xj , yj , zj ) e concentração correspondente fj (ppm)
- serão interpolados e teremos uma função de três variáveis f = f (x, y, z) que deverá ter
curvas de isoconcentração traçadas para várias cotas de altitude z.
2. Interpolação dos dados

Os dados (provenientes de medições) estão dispostos em uma tabela com N linhas:
x y z concentração f (x, y, z) em ppm

x1 y1 z1 f1 = f (x1 , y1 , z1 )
.. .. .. ..
. . . .
xN yN zN fN = f (xN , yN , zN )
A concentração do poluente em um ponto (x, y, z) é estimada segundo a seguinte média

ponderada
N
X fj
j=1
(x − xj )2 + (y − yj )2 + (z − zj )2
f (x, y, z) = N
X 1
j=1
(x − xj )2 + (y − yj )2 + (z − zj )2
Cuidado especial deve ser tomado na implementação da função f por causa de zeros no
denominador.
3. Aplicação do Método de Newton

Dada uma função de três variáveis f = f (x, y, z) (representando a concentração de poluen-
tes), com (x, y, z) ∈ [a1 , a2 ] × [b1 , b2 ] × [c1 , c2 ] e fixada uma cota de altitude z0 ∈ [c1 , c2 ], a
curva de nı́vel c ∈ R de f (x, y, z0 ) é o conjunto
{(x, y) ∈ [a1 , a2 ] × [b1 , b2 ] | f (x, y, z0 ) = c}

Para desenhar as curvas de nı́vel de f (x, y, z0 ) procedemos da seguinte forma:
Dividimos os intervalos [a1 , a2 ] e [b1 , b2 ] em n e m intervalos iguais (respectivamente)

de forma a obter uma malha com pontos (xi , yj ), onde xi = a1 + i(a2 − a1 )/n e
yj = b1 + j(b2 − b1 )/m para i = 0, . . . , n e j = 0, . . . , m.
2
Para cada xi aplicamos várias vezes o Método de Newton para obter um zero de
g(y) = f (xi , y, z0 ) − c. utilizando como chutes iniciais os valores y0 , y1 , . . . , ym . Veja
que poderemos encontrar várias raı́zes ou então nenhuma raiz.
Mutatis mutandis para cada yj e zeros de h(x) = f (x, yj , z0 ) − c
O procedimento acima nos fornece uma coleção de pontos que devem ser plotados na tela
do computador. O(s) gráfico(s) deve(m) ser salvo(s) e inserido(s) no relatório da atividade.
Na implementação do Método de Newton você deverá estabelecer um número máximo de
iterações e também calcular e utilizar a variação relativa entre iteradas consecutivas. Os
parâmetros para o critério de parada deverão ser establecidos por vocês. Não é má ideia
também calcular o valor de f nos pontos produzidos para saber ser estão mesmo na curva
de nı́vel procurada.
As derivadas parciais ∂f /∂x ou ∂f /∂y deverão ser calculadas analiticamente e depois imple-
mentadas no Scilab. O mesmo cuidado com zeros no denominador se aplica aqui também.
4. Dados do problema e exemplo de saı́da

A região para a análise da poluição está localizada no paralelepı́pedo [−3, 3] × [−3, 3] × [0, 3].
Queremos as curvas de isoconcentração para os nı́veis c = 15 ppm, c = 20 ppm, c = 25 ppm
e c = 30 ppm para as cotas z = 0.0, z = 0.5, z = 1.0, z = 1.5, z = 2.0, z = 2.5 e z = 3.0.
As concentrações estão na tabela abaixo:
x y z concentração (ppm)
0.0 0.0 0.0 5.0
0.0 0.0 1.0 35.0
0.0 0.0 2.0 25.0
0.0 1.0 0.0 10.0
0.0 1.0 0.5 5.0
0.0 1.0 1.5 45.0
0.5 0.0 0.8 0.0
0.5 1.5 0.8 20.0
1.0 2.0 1.0 50.0
1.0 1.5 1.5 40.0
Observações:
O gráfico abaixo apresenta o traçado das curvas de isoconcentração para as
concentrações de 15 ppm (todas as cotas), 20 ppm (cotas z = 0.0, z = 0.5), 25
ppm (cotas z = 0.0, z = 0.5) e 30 ppm (cotas z = 0.0, z = 0.5, z = 1.0).
Dependendo da forma como o problema é resolvido, para algumas concentrações

e cotas o programa pode demorar um certo tempo (acima de 1 h) para gerar o
gráfico.
3
• Você deverá obter coisas do tipo:

Atividade 01

Enviado por

Direitos autorais:

Formatos disponíveis

Atividade 01

Enviado por

Dados do documento

Descrição original:

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Atividade 01

Enviado por

Direitos autorais:

Formatos disponíveis

UNIVERSIDADE FEDERAL DE SÃO PAULO

Instituto de Ciências Ambientais, Quı́micas e Farmacêuticas

1. O programa deverá ter o seguinte cabeçalho:

2. A atividade pode ser feita em dupla ou trio.

3. O programa deverá ser obrigatoriamente implementado no Scilab

4. Atividades atrasadas não serão recebidas.

5. Todo aluno deverá entregar o código e um relatório, tudo devidamente identificado. A

Curvas de Isoconcentração de Poluentes

2. Interpolação dos dados

x y z concentração f (x, y, z) em ppm

A concentração do poluente em um ponto (x, y, z) é estimada segundo a seguinte média

3. Aplicação do Método de Newton

{(x, y) ∈ [a1 , a2 ] × [b1 , b2 ] | f (x, y, z0 ) = c}

Dividimos os intervalos [a1 , a2 ] e [b1 , b2 ] em n e m intervalos iguais (respectivamente)

Mutatis mutandis para cada yj e zeros de h(x) = f (x, yj , z0 ) − c

4. Dados do problema e exemplo de saı́da

Dependendo da forma como o problema é resolvido, para algumas concentrações

Você também pode gostar

Atividade 01

Enviado por

Direitos autorais:

Formatos disponíveis

Atividade 01

Enviado por

Dados do documento

Descrição original:

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Atividade 01

Enviado por

Direitos autorais:

Formatos disponíveis

UNIVERSIDADE FEDERAL DE SÃO PAULO

Instituto de Ciências Ambientais, Quı́micas e Farmacêuticas

1. O programa deverá ter o seguinte cabeçalho:

2. A atividade pode ser feita em dupla ou trio.

3. O programa deverá ser obrigatoriamente implementado no Scilab

4. Atividades atrasadas não serão recebidas.

5. Todo aluno deverá entregar o código e um relatório, tudo devidamente identificado. A

Curvas de Isoconcentração de Poluentes

2. Interpolação dos dados

x y z concentração f (x, y, z) em ppm

A concentração do poluente em um ponto (x, y, z) é estimada segundo a seguinte média

3. Aplicação do Método de Newton

{(x, y) ∈ [a1 , a2 ] × [b1 , b2 ] | f (x, y, z0 ) = c}

 Dividimos os intervalos [a1 , a2 ] e [b1 , b2 ] em n e m intervalos iguais (respectivamente)

 Mutatis mutandis para cada yj e zeros de h(x) = f (x, yj , z0 ) − c

4. Dados do problema e exemplo de saı́da

 Dependendo da forma como o problema é resolvido, para algumas concentrações

Você também pode gostar

Dividimos os intervalos [a1 , a2 ] e [b1 , b2 ] em n e m intervalos iguais (respectivamente)

Mutatis mutandis para cada yj e zeros de h(x) = f (x, yj , z0 ) − c

Dependendo da forma como o problema é resolvido, para algumas concentrações