realizar rutina
Publicado por gonzalo (1 intervención) el 27/01/2010 22:13:03
Hola que tal.
Tengo una duda que me aqueja de hace bastante tiempo, tengo que hacer un algoritmo que realice lo siguiente,
necesito calcular el error global de una ecuación diferencial, la fórmula es:
e(n+1)=exp(-5*(n+1)*Delta)-y(n+1)
donde y(n+1)=y(n)-5*y(n)*Delta, con Delta paso entre cada paso del tiempo.
Necesito generar una fórmula que me genere el error global, cuando Delta=1, 1/2, 1/4, 1/8,..., 1/2^(13). Me complica generar los y(n), ya que a medida que cambio el Delta, el y(1) no es el mismo que el y(1) para Deltas distintos, me explico, por ejemplo,
Cuando Delta = 1, entonces el n=0 (exponente de 1/2^n ), así queda
y(1)=y(0)-5*y(0)*Delta=1-5*1*1=-4
pero cuando Delta = 1/2 , entonces n=1, y así tenemos que:
y(2)=y(1)-5*y(1)*Delta =y(1)-5*y(1)*(1/2),
pero
y(1) =y(0)-5*y(0)^*Delta=1-5*1*(1/2)=-3/2
el cual es claramente distinto al y(1) cuando el Delta = 1.
Es fácil generarlos uno por uno, pero la idea es que se generen todos de una vez.
Espero que me puedan ayudar con este problemita.
Saludos!!!
Gonzalo
Tengo una duda que me aqueja de hace bastante tiempo, tengo que hacer un algoritmo que realice lo siguiente,
necesito calcular el error global de una ecuación diferencial, la fórmula es:
e(n+1)=exp(-5*(n+1)*Delta)-y(n+1)
donde y(n+1)=y(n)-5*y(n)*Delta, con Delta paso entre cada paso del tiempo.
Necesito generar una fórmula que me genere el error global, cuando Delta=1, 1/2, 1/4, 1/8,..., 1/2^(13). Me complica generar los y(n), ya que a medida que cambio el Delta, el y(1) no es el mismo que el y(1) para Deltas distintos, me explico, por ejemplo,
Cuando Delta = 1, entonces el n=0 (exponente de 1/2^n ), así queda
y(1)=y(0)-5*y(0)*Delta=1-5*1*1=-4
pero cuando Delta = 1/2 , entonces n=1, y así tenemos que:
y(2)=y(1)-5*y(1)*Delta =y(1)-5*y(1)*(1/2),
pero
y(1) =y(0)-5*y(0)^*Delta=1-5*1*(1/2)=-3/2
el cual es claramente distinto al y(1) cuando el Delta = 1.
Es fácil generarlos uno por uno, pero la idea es que se generen todos de una vez.
Espero que me puedan ayudar con este problemita.
Saludos!!!
Gonzalo
Valora esta pregunta
0
