Autor |
Mensaje |
avogadro
Nivel 5
Registrado: 11 Ago 2008
Mensajes: 127
Carrera: Química
|
|
Los autovalores me quedaron 1+- Raiz de (-40,23 K).
O sea, como vos decias numeros complejos 1+- i Raiz (40,23 K)
Si le calculo el modulo al primer autovalor me da 1 + 40,23 k y esto es mayor que uno siempre, por lo tanto es inestable siempre.
Te quedo algo asi?
|
|
|
|
|
|
|
|
|
Granjero
Nivel 4
Edad: 33
Registrado: 06 Sep 2007
Mensajes: 115
Ubicación: poca
Carrera: Informática
|
|
Pastore escribió:
|
Si a mi me quedaron iguales las ecuaciones para R-K2.
Tenes idea el k o h optimo para el metodo Euler? Planteando me quedan autovalores complejos. Nose como se resolvera si sera incondicionalmente estable o que..
|
La condición es que los 2 autovalores sean menores a 1 en módulo, así que tenés que sacar el módulo de los nros. complejos, y ver cómo tiene que ser k para que eso sea menor a 1.
Ahora, tengo una duda con eso: yo lo resolví, y me queda que el método es incondicionalmente inestable.
Entonces ¿cómo elijo un paso de cálculo "adecuado" y cómo calculo experimentalmente el valor de paso máximo para el cuál Euler se hace estable, como pide el enunciado? No se supone que es siempre inestable?
|
|
|
|
_________________ Lo envenenaron con cianuro... para colmo... cianuro en mal estado...
|
|
|
|
|
Pastore
Nivel 6
Registrado: 06 Ene 2009
Mensajes: 283
Carrera: Informática
|
|
como les quedo la matriz para plantear los autovalores?
a mi algo asi pero creo que esta mal:
1-lamda k
-Lw^2k 1-lamda
porque me quedan siempre k<0 la condicion de estabilidad
|
|
|
|
|
|
|
|
|
Pastore
Nivel 6
Registrado: 06 Ene 2009
Mensajes: 283
Carrera: Informática
|
|
Por si no se entiende, 1-lambda son el A11 y el A22 de la matriz
|
|
|
|
_________________
|
|
|
|
|
Granjero
Nivel 4
Edad: 33
Registrado: 06 Sep 2007
Mensajes: 115
Ubicación: poca
Carrera: Informática
|
|
En donde dice -Lw^2k a mi me quedó -LCw^2k.
|
|
|
|
_________________ Lo envenenaron con cianuro... para colmo... cianuro en mal estado...
|
|
|
|
|
Pastore
Nivel 6
Registrado: 06 Ene 2009
Mensajes: 283
Carrera: Informática
|
|
gracias granjero, si me olvide el C...
otra cosa, el k optimo o el k donde se hace estable para Euler alguno lo encontro aproximadamente? A mi con un numero muy chico de k, por ejemplo 0,01 el z= 15 tengo un valor aproximadamente de 190Volts! MIentras que en z=0 son 10 volts..No creo que este tan bien el rdo..
|
|
|
|
|
|
|
|
|
Guason
Nivel 3
Registrado: 06 Feb 2010
Mensajes: 42
Carrera: Mecánica
|
|
A mi me da para h=0,01 138V con el Met de Euler y 5,3V con el Met de RK2, alguno de nosotros lo tiene mal (o capaz los dos)...
|
|
|
|
|
|
|
|
|
Kabron
Nivel 3
Edad: 37
Registrado: 20 Ago 2009
Mensajes: 49
Carrera: Electrónica
|
|
avogadro escribió:
|
Para el que quiera aportar el enunciado es basicamente resolver las ecuaciones d2 V(z,t)/dz2 = -LCw2V(z,t) y d2 I(z,t) / dz2 = -LCw2 I(z,t)
Por Euler explicito y por Runge Kutta 2 y comparar resultados y encontrar experimentalmente o analiticamente el h o k max para el cual el sistema es estable.
El que quiera aportar es bienvenido.
|
buenas
hay que resolver cada ecuacion por separado o las dos juntas?
|
|
|
|
_________________ Saludos
|
|
|
|
|
Guason
Nivel 3
Registrado: 06 Feb 2010
Mensajes: 42
Carrera: Mecánica
|
|
Las tenes que resolver por separado, pero el algoritmo que las resuelve es el mismo solo cambian las condiciones iniciales.
|
|
|
|
|
|
|
|
|
Kabron
Nivel 3
Edad: 37
Registrado: 20 Ago 2009
Mensajes: 49
Carrera: Electrónica
|
|
avogadro escribió:
|
Para el que quiera aportar el enunciado es basicamente resolver las ecuaciones d2 V(z,t)/dz2 = -LCw2V(z,t) y d2 I(z,t) / dz2 = -LCw2 I(z,t)
Por Euler explicito y por Runge Kutta 2 y comparar resultados y encontrar experimentalmente o analiticamente el h o k max para el cual el sistema es estable.
El que quiera aportar es bienvenido.
|
osea seria resolver ese sistema de ecuaciones diferenciales?
|
|
|
|
_________________ Saludos
|
|
|
|
|
Kabron
Nivel 3
Edad: 37
Registrado: 20 Ago 2009
Mensajes: 49
Carrera: Electrónica
|
|
Guason escribió:
|
Las tenes que resolver por separado, pero el algoritmo que las resuelve es el mismo solo cambian las condiciones iniciales.
|
ahora si!
y esa matriz de la cual sacan los autovalores que llegaria a ser?
|
|
|
|
_________________ Saludos
|
|
|
|
|
Guason
Nivel 3
Registrado: 06 Feb 2010
Mensajes: 42
Carrera: Mecánica
|
|
En verdad para el TP no te sirve para nada, es una forma teorica para saber el valor de k o h max de estabilidad pero no te da ninguna condicion sobre k o h, tenes que buscar el k o h en forma experimental como dice el TP........... tenes que ir probando para distintos pasos.
|
|
|
|
|
|
|
|
|
Kabron
Nivel 3
Edad: 37
Registrado: 20 Ago 2009
Mensajes: 49
Carrera: Electrónica
|
|
Buenas
tengo una pequeña duda !
en el apendice del enunciado dice v=66%c usando c=3x10^8
tengo v=1.98x10^10 esto es correcto?
luego tengo v=1/raiz(LC) usando v=1.98x10^10 tengo
LC=2.55076x10^-21
|
|
|
|
_________________ Saludos
|
|
|
|
|
Granjero
Nivel 4
Edad: 33
Registrado: 06 Sep 2007
Mensajes: 115
Ubicación: poca
Carrera: Informática
|
|
Guason escribió:
|
A mi me da para h=0,01 138V con el Met de Euler y 5,3V con el Met de RK2, alguno de nosotros lo tiene mal (o capaz los dos)...
|
A mi me queda 5,8 V con los dos métodos, usando un paso muy chico (k = 0,000001). Si uso k = 0,01, me quedan los mismos resultados que a vos. Lo que pasa es que 0,01 es un paso demasiado grande para Euler.
|
|
|
|
_________________ Lo envenenaron con cianuro... para colmo... cianuro en mal estado...
|
|
|
|
|
Kabron
Nivel 3
Edad: 37
Registrado: 20 Ago 2009
Mensajes: 49
Carrera: Electrónica
|
|
el (-LCw2) de la ecuacion V(z,t)/dz2 = -LCw2V(z,t) es una constante? osea no cambia en las iteraciones ?
si es que es una constante alguien sabe su valor pues a mi me da un valor extremadamente chico!!
(-LCw2) =-4.027998x10^-15
desde ya muchas gracias!!
|
|
|
|
_________________ Saludos
|
|
|
|
|
|