Aplica

1
Aplicações de Métodos Matemáticos
Prof. Nestor Roqueiro Laboratório de Controle de
Processos EQA - UFSC e-mail nestor_at_enq.ufsc.br
2
O Problema

• Os modelos matemáticos que representam processos
  ou equipamentos são sistemas de equações
  algébricas, diferenciais e integrais que, em
  geral, não tem solução analítica.
• Para obter soluções, portanto, devem ser
  utilizados métodos numéricos.
• Para solução numérica de problemas de medio e
  grande porte é necessário o auxilio de
  computadores.
• E para que os computadores funcionem são
  necessários programas.

3
A Solução

• Comprar o software que resolva a classe de
  problemas que deseja tratar
• Desenvolver seus proprios programas
• Contratar alguem que desenvolva os programas
• Fazer combinações das três propostas acima

4
A Melhor Solução

• A melhor solução depende da verba disponível mas,
  em geral, o conhecimento de uma linguagem de
  programação de alto nível e programas que
  auxiliem na resolução de problemas numéricos
  costumam ser a solução mais eficiente, e mais
  barata.

5
Os Problemas Matemáticos

• Os problmas matemáticos que sugem na modelagem de
  processos e equipamentos são
• Solução de sistemas de equaçoes algébricas
• Solução de sistemas de equações diferencias
  ordinárias
• Solução de sistemas de equações diferenciais
  parciais
• Ajuste de curvas e interpolação
• Otimização

6
A Aula de Hoje

• Nesta aula abordademos
• Solução de sistema de equações algébricas
  (raízes)
• Solução numérica de EDO com condições iniciais
• Solução numérica de EDO com condições de contorno
• Solução numérica de EDP
• Interpolação

7
Software

• Usaremos o Matlab que combina uma linguagem de
  programação de alto nível com bibliotecas de
  programas de diversas áreas do conhecimento. Esta
  característica facilita a implementação de
  soluções numéricas.

8
Objetivo da Aula

• Resolver numericamente problemas matemáticos
  provenientes da modelagem de processos ou
  equipamentos.

9
Reator Continuo em Estado Estacionario

• O modelo é
• os programas usados para resolver são
• reator e cstr

10
Simulação de um tanque de nivel

• O tanque de nivel está representado pela equação
• os programas usados para resolver são edo,
  nivel e tanque.

11
Pellet de Catalisador

• O pellet está representado pela equação
• os programas para solução numérica são pell e
  pellet

12
Matemática simbólica

• Em uma calculadora a linha de comando
• ysin(x)
• dará erro se não for atribuído previamente
• um valor a
• x

13

• Em cálculo simbólico é possível operar expressões
  como sin(x) sem atribuir valores numéricos às
  variáveis.
• Ex Cálculo da derivada de sin(x)
• gtgtsyms x
• gtgtfsin(x)
• gtgtdiff(f)
• gtgt
• cos(x)

14
Outros exemplos

• Determinante de uma matriz simbólica
• syms a b c d
• M a,bc,d
• det(M)
• Integral
• fint(x3/sqrt(1-x),a,b)

15
Conversão de variáveis

• Para converter uma variável simbólica em numérica
  usa-se double
• Ex
• phisym((1sqrt(5))/2)
• double(phi)

16
Substituição de variáveis

• Para substituir uma variável em uma expressão
  simbólica usa-se subs
• Ex
• syms s
• fax2bxc
• subs(f,x,s)

17
Derivação

• Para derivação analítica utiliza-se a função diff
  
• Ex
• fax3bx2-xc
• diff(f,b,2)deriva em relação a b duas vezes

18
Integração

• Para integração analítica utiliza-se a função int
  
• Ex
• syms m n
• fsin(s2x)
• int(f,s,m,n)

19
Resolução de equações algébricas

• Um sistema de equações algébricas
• syms u c d v
• e2u-cdv-10
• e1d(cu)/2-v
• e3vd-uc/4
• e4vu-c8d-1
• resolve-se usando
• solve(e1,e2,e3,e4)

20
Resolução de equações diferenciais

• Uma equação diferencial pode ser resolvida usando
  a função dsolve
• Ex
• ydsolve('D2y-2Dy-3y0','y(0)0','y(1)1')
• tentar também
• ydsolve('D2y-2Dy-3y0')

21
Formatação e simplificação

• Para simplificar expressões usa-se a função
  simple
• f (1/x36/x212/x8)(1/3)
• ysimple(f)
• e para melhorar a apresentação usa-se
• pretty(y)

22
Representação gráfica

• Para representar graficamente uma expressão
  simbólica como
• syms t
• y-5t220t30
• ezplot(y)

23
Simulação de um tanque de nível

• O tanque de nível está representado pela equação
• A solução analítica pode ser calculada como
• hdsolve('Dhq/A-Cvh')
• e a solução pode ser observada usando
• pretty(h)

24

• Para representar graficamente usa-se a função
  ezplot, mas antes devem ser substituídos os
  valores de q,A,C1 e Cv
• hsubs(h, 'q',2)
• hsubs(h, 'C1',5)
• hsubs(h, 'Cv',1)
• hsubs(h, 'A',10)
• ezplot(h)

25

• Problema Proposto
• Resolver y-xy-y0
• Para resolver
• ydsolve('D2y-xDy-y0','y(0)1','Dy(0)2',
  'x')
• Para representar graficamente
• ezplot(y)

26
Pellet de Catalisador

• O pellet está representado pela equação
• E a linha de comando para resolver é
• Cadsolve('x2D2yxDy-x2yk0','Dy(0)0','y(R)
  Cas','x')
• Casubs(Ca,'Cas',1)
• Casubs(Ca,'k',1)
• Casubs(Ca,'R',1)
• plot(1/besseli(0,1)besseli(0,00.011))

27
Exemplos

28
Soluções

• caso1
• ydsolve('Dy1y2','x')
• caso2
• ydsolve('2xyDy-y2x20','y(1)1','x')
• ezplot('(-x22x)(1/2)')

29
Problemas sugeridos

30
(No Transcript)
31
(No Transcript)
32
(No Transcript)
33
(No Transcript)
34
(No Transcript)