Baixe Mini-Curso de MAPLE EDO & Modelos Matemáticos e outras Notas de estudo em PDF para Engenharia Civil, somente na Docsity!
MAPLE E.D.O E MODELOS MATEMÁTICOS
AULA Nº 1.
INTRODUÇÃO
O Maple V é um software capaz de resolver problemas em diversos ramos da matemática. Dentro dos recursos disponíveis, temos o Cálculo Diferencial e Integral, bem como Gráficos de 2 e 3 dimensões. Este material tem como finalidade apresentar ao aluno conceitos iniciais deste aplicativo para servir de subsídio para os alunos que já estudaram Equações Diferenciais. Apresentaremos os principais comandos para resolução de equações diferenciais dando soluções exatas e gráficas. Daremos também, a resolução através do Maple, para alguns Modelos Matemáticos sejam eles: Decaimento Radioativo, Resfriamento de Newton, Crescimento Populacional, Taxas de Variação, Circuito Elétricos, Modelos Geométricos entre outros.
CONCEITOS INICIAIS
Uma linha de comando do MAPLE é caracterizada por > (prompt)
Ao final de cada linha de comando escrita no MAPLE devemos sempre colocar : (dois pontos) e teclar ente, neste caso o Maple executa o comando, ou colocar ; (ponto e vírgula) e teclar enter e neste caso o Maple executa e exibe o resultado do comando.
Para as operações devem ser utilizados os símbolos:
- (soma), - (diferença),* (produto),/ (divisão), ^(potência), x^(m/n)(para raiz n-ésima de x elevado a m-ésima potência), sqrt(x) (raiz quadrada de x), = (igual), <(menor),> (maior), <=(menor ou igual), >= (maior ou igual), <>(diferente). O número e é representado por exp(1) e o número pi é representado por Pi. Usa-se := para se definir uma expressão.
PACOTES
Alguns comandos do Maple são de uso específico e por isto são agrupados em pacotes. Para disponibilizar estes comandos é preciso carregar o pacote que é feito com o comando with. Pode-se carregar mais de um pacote em uma mesma aplicação. Faremos uso nesta apostila dos seguintes pacotes:
**> with(DEtools):
with(plots): with(student):**
Ou usando; no final de cada pacote que exibe o resultado.
COMANDOS MAIS USADOS PARA RESOLUÇÃO DE EQUAÇÕES
DIFERENCIAIS
**> dsolve(ode);
dsolve(ode,y(x),extra_args); dsolve({ode,ICs},y(x),extra_args) dsolve({sysODE,ICs},{funcs},extra_args)**
Os parâmetros tem os seguintes significados:
ode- uma equação diferencial ordinária y(x)- qualquer função indeterminada de uma variável ICs- condições iniciais de uma equação diferencial {sysODE}- um conjunto com um sistema de equações diferenciais ordinárias {funcs}-um conjunto com funções indeterminadas extra_args- opcional, dependendo do tipo de problema a ser resolvido, como por exemplo: "implicit"-quando se deseja a solução na forma implícita "explicit"- quando se deseja a solução na forma explícita.
dsolve é capaz de resolver diferentes tipos de problemas de E.D.O que incluem desde simples equações como sistemas de equações diferenciais, envolvendo ou não condições iniciais, bem como soluções numéricas dentre as quais estudaremos soluções por séries, alem de soluções usando transformadas de Laplace e de Fourier.
> DEplot(deqns, vars, trange, inits):
DEplot - este comando traça as curvas solução de uma EDO de primeira ordem ou de um sistema de equações de primeira ordem. inits - condição inicial para o traçado de curvas solução num campo direcional.
OUTROS COMANDOS
rhs - Expressa o lado direito da expressão
Exemplo:
impliciplot - para traçado de gráfico na forma implícita.
Exemplo:
with(plots): implicitplot3d(z=x^2+y^2,x=-3..3,y=-3..3,z=0..6);
OBSERVAÇÃO :
Se ao final dos comandos plot ou implicitplot colocarmos: color=cor(altera a cor do gráfico para a cor especificada) title= 'nome'(escreve o nome do gráfico) scaling=constrained(coloca escalas iguais nos eixos horizontal e vertical)
Exemplo:
with(plots): implicitplot(x^2+y^2=1,x=-3..3,y=-3..3,color=blue,title='modelo');
simplify - simplifica a expressão. Exemplo:
simplify(4^(1/2)+3);
simplify(sin(x)^2+cos(x)^2);
unapply - expressa uma função de R em R.
Exemplo:
p:=x^2+sin(x)+1;
f:=unapply(p,x);
q:=x^2+y^3+1;
t:=unapply(q,x);
showtangent - gráfico da função e da reta tangente.
Exemplo:
with(student): showtangent(x^2+5,x=2);
Veremos a seguir alguns exemplos que mostram como usar o comando dsolve na solução de simples equações diferenciais.
Exemplos:
Resolver as seguintes equações diferenciais ordinárias:
Exemplo 1:
g:=2y(x)/(1-x^2);*
ode:=diff(y(x),x)=g;
ans1:=dsolve(ode);
ic:=y(0)=y0;
soln:=dsolve({ode,ic},y(x));
phi:=rhs(soln);
plot({subs(y0=-1,phi),subs(y0=-3,phi)},x=1..10);
Exemplo 2:
h:=y(x)ln(y(x))/sin(x);*
edo:=diff(y(x),x)=h;
ans2:=dsolve(edo);
ic:=y(Pi/2)=exp(1);
soln:=dsolve({edo,ic},y(x));
cal:=rhs(soln);
plot(cal,x=-1..1);
Exercício 1: Encontre a solução geral da equação 4 ( ) 6
x y x
x y x dx
dy
Vamos traçar as curvas soluções:
inits:={[-2,0.5],[1,3]}: p1:=DEplot(ode1,y(x),x=-2..2,inits): p2:=dfieldplot(ode1,y(x),x=-2..2,y=-2..2): display({p1,p2});
Exercício 1:
Trace o campo direcional e esboce as curvas soluções para y dx
dy = −.
Apresentaremos a seguir, as soluções, através do Maple, para alguns Modelos Matemáticos.
with(plots): with(student):
with(DEtools):
Exemplo 1.
Encontre a equação da curva sabendo-se que a inclinação da reta tangente num ponto
qualquer é y
− 2 x .
ode1.1:=diff(y(x),x)=-2x/y(x);*
eq1:=separablesol(ode1.1,y(x));
implicitplot({1/2y^2+x^2 =1,1/2y^2+x^2 =30,1/2y^2+x^2 =3,1/2y(x)^2+x^ =8,1/2y(x)^2+x^2 =13},x=-10..10,y=-10..10,color=blue,scaling=constrained);*
Exercício 1:
Determine a equação da curva sabendo-se que a 2ª derivada em cada ponto (x; y) é igual a 1/x 3. Sabe-se ainda que a curva passa pelo ponto (2 ; 4) e a equação da reta tangente neste ponto é 2y = 3x + 2.
Exemplo 2: A aceleração de uma partícula é diretamente proporcional ao tempo t. Para t = 0 seg. a velocidade da partícula é v = 9 m/s. Sabendo-se que ambas, velocidade e coordenada de posição são iguais a zero quando t = 3s, escreva a equação do movimento para a partícula.
eq3:=diff(v(t),t)=kt;*
eq3.1:=dsolve({eq3,v(0)=9},v(t));
14 é longa, aproximadamente 5745 anos, quantidades mensuráveis de carbono 14 estão presentes após milhares de anos. Libby mostrou que, se aproximadamente 0,002 ou mais da quantidade original ainda estão presentes, então pode-se determinar precisamente a proporção da quantidade original de carbono 14 que resta, por dosagem de laboratório adequada. Em outros termos: Se Q(t) é a quantidade de C14 no tempo t e Q0 é a quantidade original, então a razão
Q(t)/Q poderá ser determinada, desde que esta quantidade não seja tão pequena. a) Supondo que Q(t) satisfaça à equação dQ/dt = kQ, determine a constante k, de decaimento, para o carbono 14. b) Encontre a expressão Q(t) em qualquer tempo t, se Q(0)= Q0. c) Suponha que se descubram certos restos arqueológicos em que a quantidade residual de C14 seja 20% da quantidade original. Determine a idade desses restos.
Solução:
eq5:=diff(Q(t),t)=kQ(t);*
eq5.1:=dsolve(eq5,Q(t));
eval(solve(subs(Q(t)=Q[0],t=0,eq5.1),{_C1}));
eq5.2:=subs(%,eq5.1);
solve(subs(Q(t)=Q[0]/2,t=5745,eq5.2),{k});
eq5.3:=subs(%,eq5.2);
solve(subs(Q(t)=0.2Q[0],eq5.3),{t});*
soln:=rhs(eq5.3);
plot({subs(Q[0]=6,soln),subs(Q[0]=.9,soln),subs(Q[0]=2,soln),subs(Q[0]=.5,soln),subs( Q[0]=-.6,soln),subs(Q[0]=1.5,soln),subs(Q[0]=3,soln)},t=-10000..10000);
Exercício 2:
Num castelo inglês existe uma velha mesa redonda de madeira que muitos afirmam ser a famosa Távola Redonda do rei Arthur. Por meio de um contador Geiger (instrumento que mede a radioatividade) constatou-se que a massa M existente na mesa, no ano de 1999, é 0,894 vezes a massa M0 de C!4 que existe num pedaço de madeira viva com o mesmo peso da mesa. M0 é também a massa de C!4 que existia na mesa quando esta foi feita há t anos. A mesa pode ser a famosa Távola Redonda?
Exemplo 4: Uma certa cidade tinha uma população de 25.000 habitantes em1970 e, uma população de 30.000 em 1980. Assumindo que sua população continuará a crescer exponencialmente a uma taxa constante, que população podem seus planejadores urbanos esperar para o ano 2010?
eq7:=diff(P(t),t)=aP(t);*
eq7.1:=dsolve({eq7,P(0)=25000},P(t));
solve(subs(P(t)=30000,t=10,eq7.1),{a});
eq7.2:=subs(%,eq7.1);
subs(t=40,eq7.2);
evalf(%);
eq10.2:=subs(%,eq10.1);
solve(subs(T(t)=70,t=15,eq10.2),{k1});
eq10.2:=subs(%,eq10.2);
solve(subs(T(t)=40,eq10.2),{t});
plot(30+70.exp(-.3730771920e-1t),t=-10..10);**
Exercício 3:
O corpo de uma vítima de assassinato foi descoberto. O perito da polícia chegou à 1:00 h da madrugada e, imediatamente, tomou a temperatura do cadáver que era de 34,8º C. Uma hora mais tarde ele tomou novamente a temperatura e encontrou 34,1º C. a temperatura do quarto onde se encontrava a vítima era constante a 20º C. use a lei de resfriamento de Newton parar estimar a hora em que se deu a morte, admitindo que a temperatura normal de uma pessoa viva é de 36,5º C.
AULA Nº 2.
Apresentaremos a seguir, uma série de exercícios que devem ser resolvidos de maneira análoga aos exemplos vistos na aula anterior.
Determine a equação da curva tal que a inclinação da reta tangente em ponto qualquer (x;y) da curva é igual a metade da inclinação da reta que liga a origem ao ponto de tangência.
Um ator de cinema precisava fazer um regime para emagrecer, em virtude do seu papel num novo filme. O diretor exigiu que ele perdesse a terça parte do seu peso, que era de 120 kg, seguindo uma dieta que o emagrecesse proporcionalmente ao peso de cada instante. Nestas condições, sabendo-se que iniciada a dieta o artista emagreceu 20 kg em 40 dias, quanto tempo será necessário para que ele comece atuar no filme?
Um jarro de leite, inicialmente a 25ºC, é deixado para esfriar na varanda onde a temperatura é 0ºC. Suponha que a temperatura do leite tenha caído para 15ºC após 20 minutos. Quando a mesma será de 5ºC?
Circuito em série:
Em um circuito em série contendo apenas um resistor e um indutor, a segunda lei de Kirchhoff diz que a soma da queda de tensão no indutor (L(di/dt)) e da queda de tensão no resistor (iR) é igual à voltagem (E(t)) no circuito, ou seja Ldi/dt + Ri = E(t)
- Uma bateria de 12 volts é conectada a um circuito em série no qual a indutância é de 1/ henry e a resistência, 10 ohms. Determine a corrente i se a corrente inicial é zero.
Circuito em série L-R-C:
A segunda lei de Kirchoff diz que a soma da queda de voltagem através do indutor, resistor e capacitor é igual à voltagem no circuito, ou seja : L di/dt + Ri + q/C = E(t) Como a carga q(t) no capacitor está relacionada com a corrente i(t) por i(t)=dq/dt temos: L dq/dt + R dq/dt + q/C = E(t)
- Encontre a carga q(t) no capacitor em um circuito em série L-R-C quando L=0,25 henry, R=10 ohms, C=0,001 farad, E(t)=0, q(0)=q coulombs e i(0)=0.
2) Um ator de cinema precisava fazer um regime para emagrecer, em virtude do seu papel
num novo filme. O diretor exigiu que ele perdesse a terça parte do seu peso, que era de 120 kg, seguindo uma dieta que o emagrecesse proporcionalmente ao peso de cada instante. Nestas condições, sabendo-se que iniciada a dieta o artista emagreceu 20 kg em 40 dias, quanto tempo será necessário para que ele comece atuar no filme?
eq6:=diff(P(t),t)=rP(t);*
eq6.1:=dsolve({eq6,P(0)=120},P(t));
solve(subs(P(t)=100,t=40,eq6.1),{r});
eq6.2:=subs(r = 1/40ln(5/6),eq6.1);*
solve(subs(P(t)=80,eq6.2),{t});
evalf(%);
plot(120exp(1/40ln(5/6)t),t=-1000..1000);*
3) Um jarro de leite, inicialmente a 25ºC, é deixado para esfriar na varanda onde a
temperatura é 0ºC.Suponha que a temperatura do leite tenha caído para 15ºC após 20 minutos. Quando a mesma será de 5ºC?
eq10:=diff(T(t),t)=-h(T(t)-T[a]);*
eq10.1:=subs(T[a]=0,eq10);
eq10.2:=dsolve({eq10.1,T(0)=25},T(t));
solve(subs(T(t)=15,t=20,eq10.2),{h});
eq10.3:=subs(h = -1/20ln(3/5),eq10.2);*
solve(subs(T(t)=5,eq10.3),{t});
evalf(%);
T:=unapply(rhs(eq10.3),t);
plot(T(t),t=-100..100);
