Autor |
Mensaje |
isma3l
Nivel 2
Registrado: 15 Ago 2007
Mensajes: 18
Carrera: Informática
|
|
Hola quisiera saber si alguien que esta desarrollando este tp, proceso de eutrofizacion, le pudo converger por cualquiera de los 3 metodos planteados.
Que no me converge Todas las cuentas estan bien hechas y nada.
Solo OSCILA....
Alguien
|
|
|
|
|
|
|
|
|
4WD
Administrador
Edad: 39
Registrado: 07 Sep 2006
Mensajes: 2430
Ubicación: Ingeniero
Carrera: Mecánica
|
|
Si estás tan seguro de que hiciste todo bien, ¿pensaste si la solución no puede ser oscilante? Capaz la solución es así... (igual es raro. En general los procesos térmicos y/o afines son monotónicamente convergentes)
Digo porque que converja a la solución implica oscilar si la solución es oscilante
|
|
|
|
_________________
|
|
|
|
|
Torbellino
Nivel 9
Edad: 37
Registrado: 29 May 2006
Mensajes: 1742
Ubicación: Congreso
Carrera: Electrónica y Informática
|
|
Mmmm los métodos de numérico 1, salvo que se especifique lo contrario, no son buenos para soluciones oscilantes. Al final del curso, recién, dan 1 o 2 que funcionan bien para esos casos. Así que no creo que la solución sea oscilante, a menos que la idea del TP sea comparar dos métodos que no sirven contra uno que sí.
Además, según wikipedia, eutrofización es un enriquecimiento de nutrientes, no me suena que la solución a eso sea oscilante.
Creo que la duda la resolvemos mejor si isma3l nos cuenta si oscila en el sentido en que oscila sen(x) o si el resultado en cada paso temporal da positivo y negativo alternadamente, o tiene un comportamiento similar. Supongo que se refiere a esto último, en ese caso el método no está convergiendo.
isma3l, si el método no converge, fijate el paso temporal que estás tomando en cada método. Si nos decís cuáles son los métodos que tenés que usar capaz te podemos ayudar mejor. También si posteas el código fuente que estás usando (pero sólo la parte que resuelve la ecuación!!! Por favor no nos postees el programa entero).
Saludos
|
|
|
|
_________________
No hay vuelta atrás...
Spike Spiegel escribió:
|
Por un lado se celebran las hazañas de San Martín, Bolivar y demases, la reforma de 1918, el cordobazo y otras tantas en Argentina, Latinoamérica y el mundo entero. No sé cuántos habrán llorado mirando Braveheart al grito de FREEDOM de Wallace y dicho "cuántos huevos, viejo", tenido ganas de cambiar el mundo cuando terminaron de ver V for Vendetta o celebrado toda la ficcionaria justicia que solía hacer El Zorro.
Y sin embargo...
"Ay, no, violencia no. Ay, no, corte de calle, no. Ay, no, piden democracia pero son antidemocráticos con sus métodos. Ay, no, a la facultad se viene a estudiar"
¡PERO QUÉ MANGA DE PUTOS!
|
|
|
|
|
|
isma3l
Nivel 2
Registrado: 15 Ago 2007
Mensajes: 18
Carrera: Informática
|
|
Si lo hice bien, cero errores, y ya no si deba converger, pero en el tp dice q debe alcanzar los 20ug/l y si parto de 4,21 dependiendo de los valores iniciales puede llegar a los 20 pero seguiria oscilando, y es que las constantes varian cada mes, cada año de de la misma forma.
Si alguien lo hizo y me responde le agradeceria muchisimo.
|
|
|
|
|
|
|
|
|
isma3l
Nivel 2
Registrado: 15 Ago 2007
Mensajes: 18
Carrera: Informática
|
|
Mira el paso tomado fue de 1/30 pues es por dia, luego le puse 1/300 y da lo mismo, luego 1/500 y es lo mismo.
Luego le puse 500 años con un paso de 1/30 y dio casi lo mismo.
Los metodos son euler implicto, euler explicito, rugen kutta 2.
Ahora los valores tomados en la ecuacion dependen del mes del año, para cada mes es diferente, quizas por eso oscila, no se.
Estan son las ecuaciones de cada variable, algunas son constantes q no las puse, tampoco puse los valores iniciales.
V(t)=Vo+(Qe-Qs)∆ t
Qe=Qm+∆Q.cos(wt) ,Qm y ∆Q son constantes
Qs=Qm+∆Q.cos(w(t-3)) w=(pi/6) [1/mes]
D: 0 si: P>=Pc; Do si: P<Pc Do es constante
Ld: Ld1 si t=9 , Ld2 si t=3 , 0 cualquier otro caso
Pe: Pe1 si: 3<t<=9 meses Pe2 si: 9<t<3 Pe1,Pe2 son constantes
H(t)= V(t)/areaHorizontal; //areahorizontal es area del lago
[B]Esta es la formula:[/B]
dP/dt= (1/V)*(Qe.Pe-Qs.P)- k.P/H + A.Ld/V + Ω.D.(Pc-P)^1/2
[B]Codigo fuente:[/B] hecho en java
[code]private void actualizar() {
this.actualizarQentrada();
this.actualizarQsalida();
this.actualizarVolumen();
this.actualizarAltura();
this.actualizarPe();
this.actualizarLd();
this.actualizarD();
}
private void actualizarVolumen(){
this.volumenLago = this.volumenLago + (this.Qe - this.Qs) * this.deltaT;
}
private void actualizarQentrada(){
this.Qe = this.Qm + this.deltaQ * Math.cos(this.W * this.tiempo);
}
private void actualizarQsalida(){
this.Qs = this.Qm + this.deltaQ *Math.cos (this.W * (this.tiempo -3));
}
private void actualizarLd(){
if ((this.tiempo >=8.0) && (this.tiempo <9.0)) //todoslos dias del mes 9
this.Ld = this.Ld1;
else if ((this.tiempo >= 2.0) && (this.tiempo <3.0))
this.Ld = this.Ld2;
else
this.Ld = 0;
}
private void actualizarPe(){
if ((this.tiempo > 3.0) && (this.tiempo <8.0))
this.Pe = this.Pe1;
else
this.Pe = this.Pe2;
}
private void actualizarAltura(){
this.H = this.volumenLago/ this.areaHorizontal;
}
private void actualizarD(){
if (Pc-P<=0)
this.D = 0;
else
this.D = this.D0;
}
private void eulerExplicito(){
P= P + this.deltaT *obtenerConcentracionDeFosforo(P);
this.tiempo += this.deltaT;
actualizar();
}
private double obtenerConcentracionDeFosforo(double P){
return (((1/volumenLago)*(Qe*Pe-Qs*P))-(k*P/H)+(areaHordeDrenaje*Ld/volumenLago)+
areaHorizontal*D*Math.sqrt(Math.abs(Pc - P)));
}
Ahi esta lo esencial del codigo, luego por otra parte le colocas el paso que es deltaT, y la cantidad de años.
|
|
|
|
|
|
|
|
|
Torbellino
Nivel 9
Edad: 37
Registrado: 29 May 2006
Mensajes: 1742
Ubicación: Congreso
Carrera: Electrónica y Informática
|
|
Lo que pasa es que yo no tengo el enunciado. Esa ecuación que te piden resolver, ¿te piden sacar la concentración en cada día? Porque si te piden sacar la concentración en cada día, tenés que tomar un paso mucho más chico que un día, en principio. Otra cosa es si ya hicieron un análisis de convergencia y resulta que alcanza con un paso de un día.
De todos modos, el Euler implícito, si no me acuerdo mal, converge le pongas lo que le pongas como paso. Tomalo con pinzas, capaz me falla la memoria.
Otra cosa, actualizarVolumen() también resuelve una ecuación diferencial y lo estás haciendo incondicionalmente con Euler Explícito. Cuando resuelvas con implícito o con RK-2, también deberías resolver eso con ese método.
No sé, la verdad que no veo más cosas, sin el enunciado no puedo ayudarte más. Me voy a morfar y a estudiar y si después me queda tiempo y ponés el enunciado, veo si te puedo ayudar en algo más.
Pero revisá eso último que te dije, porque puede ser el problema.
Saludos
|
|
|
|
_________________
No hay vuelta atrás...
Spike Spiegel escribió:
|
Por un lado se celebran las hazañas de San Martín, Bolivar y demases, la reforma de 1918, el cordobazo y otras tantas en Argentina, Latinoamérica y el mundo entero. No sé cuántos habrán llorado mirando Braveheart al grito de FREEDOM de Wallace y dicho "cuántos huevos, viejo", tenido ganas de cambiar el mundo cuando terminaron de ver V for Vendetta o celebrado toda la ficcionaria justicia que solía hacer El Zorro.
Y sin embargo...
"Ay, no, violencia no. Ay, no, corte de calle, no. Ay, no, piden democracia pero son antidemocráticos con sus métodos. Ay, no, a la facultad se viene a estudiar"
¡PERO QUÉ MANGA DE PUTOS!
|
|
|
|
|
|
Mr Nadie
Nivel 9
Registrado: 20 Dic 2007
Mensajes: 2885
Carrera: Civil
|
|
Mira, yo tuve unos problemitas con ese tp tambien, hasta q ayer, inspirado, creo haber alcanzado la solucion.
Si mal no recuerdo, no hay problema en q sea oscilante. Es decir, al menos entendi q dijeron eso en clase. A mi tmb me da oscilante. Sube en los primeros meses , baja en Marzo digamos, y vuelve a subir en Septiembre. Y si, tarda en alcanzar los 20 ug/l. Creo q lo tire para 10 años y no los alcanzo, lo cual, a mi entender, no esta mal por lo q dijeron en clase (no se si te acordas el grafico, q no era una sinusoide, sino una funcion que alcanzaba maximos y minimos, y los maximos, quizas no superaban los 20ug/l)
Igualmente, el lunes a Tarela lo voy a interrogar ni bien tenga la oportunidad, porque tengo un par de dudas tmb, entre ellas, esta, que para mi, segun entendi lo q dijeron en clase, no esta mal que oscile.
Saludos
|
|
|
|
_________________ Qué es registrar?
viedmense escribió:
|
PD: increible la capacidad de mantenerse en el mismo grado de pedo durante mas de 6 horas de mr nadie, ni mejoró ni empeoró
|
|
|
|
|
|
isma3l
Nivel 2
Registrado: 15 Ago 2007
Mensajes: 18
Carrera: Informática
|
|
Bueno ahi te va el link del enunciado:
Tp 2
|
|
|
|
|
|
|
|
|
isma3l
Nivel 2
Registrado: 15 Ago 2007
Mensajes: 18
Carrera: Informática
|
|
Ok Mr. Nadie. tambien oi a un compañero que le daba oscilante, pero si modifico el numero de padron hay veces q si alcanza los 20 pero sigue oscilando, luego si le doy otro padron no los alcanza.
Enfin, tendre presente lo que dices.
|
|
|
|
|
|
|
|
|
Torbellino
Nivel 9
Edad: 37
Registrado: 29 May 2006
Mensajes: 1742
Ubicación: Congreso
Carrera: Electrónica y Informática
|
|
Mr Nadie escribió:
|
Si mal no recuerdo, no hay problema en q sea oscilante.
|
De ahí venía mi pregunta original de si es oscilante en el buen sentido, como oscila un sen(x), o si es oscilante en el mal sentido, que la solución cambia de signo en cada paso y finalmente explota.
Si le da oscilante en el buen sentido, el método converge (otra cosae s si es la solución buscada, pero casi seguro que sí).
Saludos
|
|
|
|
_________________
No hay vuelta atrás...
Spike Spiegel escribió:
|
Por un lado se celebran las hazañas de San Martín, Bolivar y demases, la reforma de 1918, el cordobazo y otras tantas en Argentina, Latinoamérica y el mundo entero. No sé cuántos habrán llorado mirando Braveheart al grito de FREEDOM de Wallace y dicho "cuántos huevos, viejo", tenido ganas de cambiar el mundo cuando terminaron de ver V for Vendetta o celebrado toda la ficcionaria justicia que solía hacer El Zorro.
Y sin embargo...
"Ay, no, violencia no. Ay, no, corte de calle, no. Ay, no, piden democracia pero son antidemocráticos con sus métodos. Ay, no, a la facultad se viene a estudiar"
¡PERO QUÉ MANGA DE PUTOS!
|
|
|
|
|
|
isma3l
Nivel 2
Registrado: 15 Ago 2007
Mensajes: 18
Carrera: Informática
|
|
Si Torbellino oscila en el buen sentido, siempre se repite la misma onda cada año, y una pregunta que significa ese diodo de tenes de imagen.
|
|
|
|
|
|
|
|
|
Torbellino
Nivel 9
Edad: 37
Registrado: 29 May 2006
Mensajes: 1742
Ubicación: Congreso
Carrera: Electrónica y Informática
|
|
Ah, bueno, si te oscila en el buen sentido, es la solución. Yo pensé que no convergía. Cuando no converge, pasa lo que yo decía de que un resultado es máyor que el del paso temporal siguiente en módulo y de signo opuesto.
El diodo deja pasar al corriente en un sentido pero no en el otro, viene a que me pasé a electrónica estando en informática y al revés no pasa nunca.
Saludos
|
|
|
|
_________________
No hay vuelta atrás...
Spike Spiegel escribió:
|
Por un lado se celebran las hazañas de San Martín, Bolivar y demases, la reforma de 1918, el cordobazo y otras tantas en Argentina, Latinoamérica y el mundo entero. No sé cuántos habrán llorado mirando Braveheart al grito de FREEDOM de Wallace y dicho "cuántos huevos, viejo", tenido ganas de cambiar el mundo cuando terminaron de ver V for Vendetta o celebrado toda la ficcionaria justicia que solía hacer El Zorro.
Y sin embargo...
"Ay, no, violencia no. Ay, no, corte de calle, no. Ay, no, piden democracia pero son antidemocráticos con sus métodos. Ay, no, a la facultad se viene a estudiar"
¡PERO QUÉ MANGA DE PUTOS!
|
|
|
|
|
|
Mr Nadie
Nivel 9
Registrado: 20 Dic 2007
Mensajes: 2885
Carrera: Civil
|
|
Claro, a mi tmb, capaz q para Diciembre siempre me da algo cercano a 14, pero no pasa de ahi. Y si, si cambio el numero de padron,cambia todo. Es mas, creo q con 999 me dio como 40 en el primer mes, bajo dspues a 10, y volvio a llegar a 30, a 10 de nuevo, a 20 y x ahi se queda. Asi q no me sorprenderia q con algunos padrones ni llegue a 20.
Pero bueno, el lunes le preguntamos a Tarela bien.
|
|
|
|
_________________ Qué es registrar?
viedmense escribió:
|
PD: increible la capacidad de mantenerse en el mismo grado de pedo durante mas de 6 horas de mr nadie, ni mejoró ni empeoró
|
|
|
|
|
|
4WD
Administrador
Edad: 39
Registrado: 07 Sep 2006
Mensajes: 2430
Ubicación: Ingeniero
Carrera: Mecánica
|
|
Observen que el enunciado dice:
Enunciado escribió:
|
Resuelva el problema hasta obtener un estado estacionario (dinámico) como solución, o hasta un tiempo mínimo de 20 años.
|
En este caso podrías asumir que si achicando el paso de tiempo siempre llegás a la misma solución, tenés convergencia...
Si el método no convergiera deberías tener bastante variación al cambiar el paso de tiempo...
|
|
|
|
_________________
|
|
|
|
|
fede!
Nivel 3
Edad: 36
Registrado: 29 Feb 2008
Mensajes: 20
Ubicación: palermo
Carrera: Mecánica
|
|
che como hicieron para obtener el valor del paso?
decian de usar newton raphson, pero como lo aplican al problema?
|
|
|
|
_________________
|
|
|
|
|
|
Ir a página 1, 2 Siguiente
|
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.
|