| 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
![[tex]y(x) = y_1 (x) + \frac{{\beta - y_1 (b)}}{{y_2 (b)}}y_2 (x)[/tex]](images/latex/b7f6dce8a01989a0186a379dd464ed2d7cc0f5d4_0.png "[tex]y(x) = y_1 (x) + \frac{{\beta - y_1 (b)}}{{y_2 (b)}}y_2 (x)[/tex]") (*)
En general si ![[tex]u_{1,i}[/tex]](images/latex/2295b2c0bd163954e35c428b9c94b1b34cc2a95b_0.png "[tex]u_{1,i}[/tex]") y ![[tex]v_{1,i}[/tex]](images/latex/941d949a2726ccb2ebbd0766e050bf7db706076c_0.png "[tex]v_{1,i}[/tex]") son aproximaciones de ![[tex]O(h^n )[/tex]](images/latex/2f01643f39a21fa099ee067466bd980a2d1c6eb1_0.png "[tex]O(h^n )[/tex]") para ![[tex]y_1 (x_i )[/tex]](images/latex/10e1621e9a877e2eec1ba520b180d254661628eb_0.png "[tex]y_1 (x_i )[/tex]") y ![[tex]y_2 (x_i )[/tex]](images/latex/5dffa1db2737978ead9388618fc1af1c393860e9_0.png "[tex]y_2 (x_i )[/tex]") , con ![[tex]i = 0,1,...,N[/tex]](images/latex/a729546ac794b2d5d863f368ba41dcdf8f006797_0.png "[tex]i = 0,1,...,N[/tex]") , entonces ![[tex]w_{1,i}[/tex]](images/latex/20540e0d207ea699ecc60b7b97341d2fb64ac877_0.png "[tex]w_{1,i}[/tex]") será una aproximación de para ![[tex]y(x_i)[/tex]](images/latex/ed4e5758902439a1506cab2a7620f91c1c937773_0.png "[tex]y(x_i)[/tex]")
En particular,
![[tex]\left| {w_{1,i} - y(x_i )} \right| \le K \cdot h^n \cdot \left| {1 + \frac{{v_{1,i} }}{{v_{1,N} }}} \right|[/tex]](images/latex/aec1d15e7a4b8b9bca41da3bf2bb1b6c0e69c1a7_0.png "[tex]\left| {w_{1,i} - y(x_i )} \right| \le K \cdot h^n \cdot \left| {1 + \frac{{v_{1,i} }}{{v_{1,N} }}} \right|[/tex]")
para alguna constante ![[tex]K[/tex]](images/latex/3c31f3404a91311bb4be7cec9d49bb510d592d34_0.png "[tex]K[/tex]") .
Esto de alguna manera justifica porqué en el tp usamos ![[tex]\left| {w_{1,i} - y(x_i )} \right|[/tex]](images/latex/a54aaa8eed6e394cf11319863fb4e39bb92aa529_0.png "[tex]\left| {w_{1,i} - y(x_i )} \right|[/tex]") 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 ![[tex]\left| {w_{1,i} - y(x_i )} \right| \sim O(h^n )[/tex]](images/latex/ff485da7dea4a73380ce70da545fe690b072a2a6_0.png "[tex]\left| {w_{1,i} - y(x_i )} \right| \sim O(h^n )[/tex]") . 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
|
|
| |
|
|
|
|
