Autor |
Mensaje |
Juanse!
Nivel 5
Registrado: 27 Feb 2009
Mensajes: 180
Carrera: Mecánica
|
|
che en el tp q hay q programar RK 4.. el error con el paso 0,1 alcanza para la tolerancia q nos piden no?' cuanto les dio..?
|
|
|
|
|
|
|
|
|
Amintoros
Nivel 8
Registrado: 20 Mar 2008
Mensajes: 533
Carrera: Química
|
|
Sí, con el paso de 1 mm ya alcanza. No lo tengo a mano, en un rato te lo paso.
|
|
|
|
_________________
Elmo Lesto escribió:
|
Bistek escribió:
|
por qué pasa que a veces entro al foro y esta todo en aleman?
|
Ahí aplicaron la transformada de Führer
|
cuando la yerba mate
|
|
|
|
|
Juanse!
Nivel 5
Registrado: 27 Feb 2009
Mensajes: 180
Carrera: Mecánica
|
|
muchachos no me anda el runge kutta 4,, va para atras jaja es rarisimo xq el de orden 2 anda joya y el de orden 4 lo hice basicamente basandome en el anterior q anda y tiene orden 2 perfecto, la verdad nose q pasa,, subo el codigo a ver si alguien puede ayudarme,,
hice el cambio de variables
dT/dr=Z
dz/dr=-1/r . Z
cambiandoZ por u y T por V me queda el siguiente algoritmo:
r:=1;
for i:=1 to j do
begin
q1v:=paso*u;
q1u:=-paso*u/r;
q2v:=paso*(u+q1u/2);
q2u:=-paso*(u+q1u/2)/(r+(paso/2));
q3v:=paso*(u+q2u/2);
q3u:=-paso*(u+q2u/2)/(r+(paso/2));
q4v:=paso*(u+q3u);
q4u:=-paso*(u+q3u)/(r+paso);
u:=u+(q1u+2*q2u+2*q3u+q4u)/6;
v:=v+(q1v+2*q2v+2*q3v+q4u)/6;
vectorsalida[i]:=v;
r:=1+paso*i;
end;
nose cual es el problema, si alguien me ayuda le agradezco,..
GRACIIAAS!!!
|
|
|
|
|
|
|
|
|
Juanse!
Nivel 5
Registrado: 27 Feb 2009
Mensajes: 180
Carrera: Mecánica
|
|
sin palabras,, se puede hasta borrar el topic x la boludez q habia q corregir jaja
es mas.. borrenlo jaj
|
|
|
|
|
|
|
|
|
SaaS
Nivel 7
Edad: 34
Registrado: 17 Dic 2008
Mensajes: 310
Ubicación: San Martín
Carrera: Informática
|
|
era v:=v+(q1v+2*q2v+2*q3v+q4u)/6; ??
|
|
|
|
|
|
|
|
|
Juanse!
Nivel 5
Registrado: 27 Feb 2009
Mensajes: 180
Carrera: Mecánica
|
|
desde afueras es muy facil... yo no lo veia!!!!!!
|
|
|
|
|
|
|
|
|
MarianAAAJ
Nivel 7
Edad: 35
Registrado: 14 Ene 2009
Mensajes: 437
Carrera: Informática
|
|
con el euler con que paso da bien? con el de 1mm?
|
|
|
|
|
|
|
|
|
Amintoros
Nivel 8
Registrado: 20 Mar 2008
Mensajes: 533
Carrera: Química
|
|
No; no sé si lo habré programado bien, pero con Euler tengo que reducir casi 100 veces el paso inicial para que cumpla la condición.
Juanse, el error máximo que genero con RK4 es 1.9445e-004 ºC, a vos te da algo parecido?
|
|
|
|
_________________
Elmo Lesto escribió:
|
Bistek escribió:
|
por qué pasa que a veces entro al foro y esta todo en aleman?
|
Ahí aplicaron la transformada de Führer
|
cuando la yerba mate
|
|
|
|
|
rama_qac
Nivel 5
Edad: 35
Registrado: 26 Jul 2008
Mensajes: 142
|
|
gentee estoy atrasadisimo con este tp, me dieron un ultimatum para entregarlo hasta mañana, necesito ayuda de alguien q lo haya hecho.. la 1er duda q me surgio es si el error es respecto a la solucion exacta o a la discretizacion anterior??
gracias!
|
|
|
|
_________________ Rama!
|
|
|
|
|
Amintoros
Nivel 8
Registrado: 20 Mar 2008
Mensajes: 533
Carrera: Química
|
|
Tengo entendido que el error es respecto a la solución exacta.
EDIT: agrego esto que encontré en el Burden, a propósito del método del disparo lineal:
La solución que obtenés es de la forma
(*)
En general si y son aproximaciones de para y , con , entonces será una aproximación de para
En particular,
para alguna constante .
Esto de alguna manera justifica porqué en el tp usamos como expresión del error. El Burden dice "en general" y no "siempre" porque pueden haber casos en donde hayan cancelaciones de términos en (*), que hagan perder precisión y no se cumpla . Pero al parecer eso no pasa en el tp
Espero que te sirva, saludos
|
|
|
|
_________________
Elmo Lesto escribió:
|
Bistek escribió:
|
por qué pasa que a veces entro al foro y esta todo en aleman?
|
Ahí aplicaron la transformada de Führer
|
cuando la yerba mate
|
|
|
|
|
rama_qac
Nivel 5
Edad: 35
Registrado: 26 Jul 2008
Mensajes: 142
|
|
muchas graciasss!
me surgio una 2da duda, cuando disminuis el paso, hay q volver a calcular la pendiente de salida interpolando? xq con la q use para el 1er paso cuando lo disminuyo el ultimo valor no me queda en 20ºC
|
|
|
|
_________________ Rama!
|
|
|
|
|
SaaS
Nivel 7
Edad: 34
Registrado: 17 Dic 2008
Mensajes: 310
Ubicación: San Martín
Carrera: Informática
|
|
rama_qac escribió:
|
muchas graciasss!
me surgio una 2da duda, cuando disminuis el paso, hay q volver a calcular la pendiente de salida interpolando? xq con la q use para el 1er paso cuando lo disminuyo el ultimo valor no me queda en 20ºC
|
tal cual... la onda es... hago dos tiros (o uno), interpolo y resulevo con la pendiente obtenida... eso para cada vez que hagas algo...
|
|
|
|
|
|
|
|
|
rama_qac
Nivel 5
Edad: 35
Registrado: 26 Jul 2008
Mensajes: 142
|
|
|
|
|
rama_qac
Nivel 5
Edad: 35
Registrado: 26 Jul 2008
Mensajes: 142
|
|
ultima duda! puede ser q con el RK de orden 2 tenga q disminuir el paso al mismo valor q para Euler?? xq me dan igualeso bastante parecido! gracias!
|
|
|
|
_________________ Rama!
|
|
|
|
|
Amintoros
Nivel 8
Registrado: 20 Mar 2008
Mensajes: 533
Carrera: Química
|
|
Me parece que no, RK2 cumple con la condición mucho antes que Euler. Con dividir el paso inicial a la mitad alcanza. De todas maneras, entregalo como te haya salido, mientras que el tp esté completo es más importante entregarlo a tiempo que entregarlo sin errores...
|
|
|
|
_________________
Elmo Lesto escribió:
|
Bistek escribió:
|
por qué pasa que a veces entro al foro y esta todo en aleman?
|
Ahí aplicaron la transformada de Führer
|
cuando la yerba mate
|
|
|
|
|
|