Autor |
Mensaje |
danie87
Nivel 5
Edad: 36
Registrado: 14 Feb 2010
Mensajes: 193
Ubicación: 34.5934°S 58.4445°W
Carrera: Informática
|
|
Hola tengo 2 dudas, en el tp2 de tarela.
1) Al usar Runge kutta 4, con h=1 el error en toda la solucion es menor que 0.05, si coloco h=0.5 el error se hace mas chico, pero al colocar h=0.25 el error crece, y con h=0.125 el error crece mas, entonces no sé que pasa.
Porq la discretizacion es perfecta y el algoritmo tambien, me inclino a pensar q es el error de redondeo, pero no recuerdo q mencionaran algo asi en clase.
2) Donde piden analizar la relacion de la magnitud del error y el tamaño de paso teniendo en cuenta los ordenes.
Note que con euler explicito el error se reduce a l mitad si reduzco el paso a la mitad, con RK2 el error se reduce 4 veces cuando reduzco el paso a la mitad, pero con Rk4 me dio mismo que Rk2, es correcto eso?
Entonces alguna sugerencia de como relacionar eso?
|
|
|
|
|
|
|
|
|
sabian_reloaded
Nivel 9
Edad: 34
Registrado: 18 Jun 2009
Mensajes: 2925
Ubicación: El bosque platense
Carrera: No especificada
|
|
Seguramente pase lo mismo que sucede con las formulas de derivacion, si planteas el error por definicion vas a ver que quedan dos terminos sumandose, uno que decrece con h y otro que crece (aunque mas lento que el primero). Derivando esa funcion vas a encontrar el h óptimo para cada algoritmo.
|
|
|
|
|
|
|
|
|
TuteQac
Nivel 2
Edad: 35
Registrado: 10 Jul 2008
Mensajes: 6
Ubicación: Ciudad Cervecera (Quilmes)
|
|
sabian_reloaded escribió:
|
Seguramente pase lo mismo que sucede con las formulas de derivacion, si planteas el error por definicion vas a ver que quedan dos terminos sumandose, uno que decrece con h y otro que crece (aunque mas lento que el primero). Derivando esa funcion vas a encontrar el h óptimo para cada algoritmo.
|
¿como quedaría el error por definición???
|
|
|
|
|
|
|
|
|
Amintoros
Nivel 8
Registrado: 20 Mar 2008
Mensajes: 533
Carrera: Química
|
|
Como te dijo Sabian, el error tiene un término debido al truncamiento (que es el que más estamos acostumbrados a analizar en esta segunda parte de la materia) y un término que corresponde al redondeo:
.
En el caso de RK, las ecuaciones en diferencias son complicadas y me parece poco práctico tratar de buscar una expresión para el error de redondeo (tendrías que conocer las cotas del error de redondeo en cada evaluación funcional). En el tp2, con RK el error cometido es muy pequeño, y si disminuís el paso aún más, la máquina representa cualquier cosa, por eso no tiene sentido estimar el orden de convergencia con todos los valores de paso que usaste con el método de Euler (que es lo que sugiere el enunciado). Con los primeros tres o cuatro pasos ya alcanza.
EDIT: me fui por las ramas así que borre unas cosas.
|
|
|
|
_________________
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
|
|
|
|
|
|
|
Ver tema siguiente
Ver tema anterior
Podés publicar nuevos temas en este foro No podés responder a temas en este foro No podés editar tus mensajes en este foro No podés borrar tus mensajes en este foro No podés votar en encuestas en este foro No Podéspostear archivos en este foro No Podés bajar archivos de este foro
|
Todas las horas son ART, ARST (GMT - 3, GMT - 2 Horas)
Protected by CBACK CrackerTracker365 Attacks blocked.
|
|
[ Tiempo: 0.2694s ][ Pedidos: 20 (0.2005s) ] |