Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- exercicios:calc1 [2012/05/11 17:11] – [Soluções Numéricas] adalardo
+++ exercicios:calc1 [2024/01/09 18:18] (atual) – edição externa 127.0.0.1
@@ Linha 104: / Linha 104: @@
 <code>
-# Lembre que dy/dt = f(y, t) no caso geral
+f <- function (N, t)
-# Vamos criar essa f:
+{
-f <- function (N, t) {
+        return (2*N)
-        return (2*N);
 }
 # No tempo inicial, N vale 20:
-N0 <- 20;
+N0 <- 20
 # O passo de tempo eh dt. Vamos rodar ateh tmax
 dt <- 0.1
-tmax <- 2;
+tmax <- 2
-euler <- function (f, N0, dt, tmax) {
+euler <- function (f, N0, dt, tmax)
+{
         # res vai retornar o vetor com todos os Ns
-        res <- NULL;
+        res <- NULL
-        N <- N0;
+        N <- N0
-        for (tempo in seq(0, tmax, dt)) {
+        for (tempo in seq(0, tmax, dt))
+        {
                 N <- N + f(N,tempo)*dt
                 res <- rbind(res, N)
@@ Linha 127: / Linha 128: @@
 # Examine a solucao numerica:
-numerica <- euler(f, N0, dt, tmax);
+numerica <- euler(f, N0, dt, tmax)
 x<- seq(0, tmax, dt)
 # A solucao correta da EDO:
-plot(x, 20*exp(2*x), typ='l', col='green');
+plot(x, 20*exp(2*x), typ='l', col='green')
 # Vamos comparar com a solucao numerica
-points(x, numerica, col='red', pch=4, ce=0.2);
+points(x, numerica, col='red', pch=4, ce=0.2)
 </code>
 A aproximação foi boa? Tente repetir o mesmo código com dt = 0.01 e 0.001 e compare.
 ====== Integração Numérica no R ======
-Vamos integrar numericamente algumas equações usando o pacote deSolve e a função //ode//.
+Não precisamos fazer todo o procedimento anterior para fazer a integração numérica no {{:exercicios:rlogo.jpg?20  |}}, existem soluções implementadas previamente que são muito mais eficientes e robusta que a nossa.
+Vamos integrar numericamente algumas equações usando o pacote //deSolve// e a função //ode//.
 Antes de tudo precisa instalar e carregar o pacote. Para instalar você pode usar o menu do RGui ou pela linha de comando, digite:
@@ Linha 216: / Linha 218: @@
   *1. $ \frac{dy}{dt} = y-y^2*f(t)$
-    * sendo: $f(t)= 2 + 0.2 sin(2\pi Mt)$;
+    * sendo: $f(t)= 0.01 + 0.01 sin(2\pi Mt)$;
     * M=15
   *2. $\frac{dn}{dt} = r(t)*n $
-    * sendo: $r(t)= 0 + 0.1 sin(2 \pi t)$
+    * sendo: $r(t)= 0.1 - 100 sin(2 \pi t)$
-Uma solução:
-<code>
-### primeiro caso
-fy3 <- function(time,y, parms)
-  {
-   n=y[1]
-   M=parms[1]
-   r=0.01+ 0.01*sin(2*pi* M *time)
-   dy.dt=n-n^2 * r
-   return(list(c(dy.dt)))
-  }
-y0 = 10
-prmt=c(M=1)
-st=seq(0.1,20,by=0.01)
-res.fy3= ode(y=y0,times=st, func=fy3,parms=prmt)
-plot(res.fy3[,1], res.fy3[,2], type="l", col="red",lwd=2, xlab="tempo", ylab="y")
-### segundo caso
-fy4 <- function(time,y, parms)
-  {
-   n=y[1]
-   rt= 0.1-100*(sin(2*pi*time))
-   dy.dt=rt*n
-   return(list(c(dy.dt)))
-  }
-y0 = c(.1)
-#prmt=c(M=1)
-st=seq(0.1,20,by=0.0001)
-res.fy4= ode(y=y0,times=st, func=fy4)
-plot(res.fy4[,1], res.fy4[,2], type="l", col="red",lwd=2, xlab="tempo", ylab="y")
-</code>