Preguntas etiquetadas con ode

Tengo que resolver varios problemas desafiantes de optimización global no convexo. Actualmente uso la Caja de herramientas de optimización de MATLAB (específicamente, fmincon()con algoritmo = 'sqp'), que es bastante eficaz . Sin embargo, la mayor parte de mi código está en Python, y me encantaría...

Actualmente estoy buscando métodos paralelos para la integración de ODE. Existe una gran cantidad de literatura nueva y antigua que describe una amplia gama de enfoques, pero no he encontrado encuestas recientes ni artículos de descripción general que describan el tema en general. Existe el libro...

Esta pregunta se migró de Mathematics Stack Exchange porque se puede responder en Computational Science Stack Exchange. Migrado hace 6 años . Estoy a punto de comenzar a trabajar en una biblioteca de software de solucionadores numéricos de EDO , y estoy luchando sobre...

¿Cómo se puede resolver numéricamente el problema del cuerpo n gravitatorio en paralelo? ¿Es posible el equilibrio precisión-complejidad? ¿Cómo influye la precisión en la calidad del modelo?

En este comentario escribí: ... integrador predeterminado de SciPy, que supongo que solo usa métodos simplécticos. en el que me refiero a SciPy's odeint, que utiliza un "método no rígido (Adams)" o un "método rígido (BDF)". Según la fuente : def odeint(func, y0, t, args=(), Dfun=None,...

Considere un IVP para el sistema ODE , y ( x 0 ) = y 0 . Más comúnmente, este problema se considera rígido cuando la matriz de Jacobi ∂ fy′=f(x,y)y′=f(x,y)y'=f(x,y)y(x0)=y0y(x0)=y0y(x_0)=y_0tienedosvalores propios con gran parte real negativa y valores propios con parte real negativa muy pequeña...

Solo tenía curiosidad por saber por qué los métodos Runge – Kutta de alto orden (es decir, mayores de 4) casi nunca se discuten / emplean (al menos que yo sepa). Entiendo que requiere un mayor tiempo de cómputo por paso (por ejemplo, RK14 con paso incorporado de 12 ° orden ), pero ¿hay alguna otra...

Deseo simular el comportamiento de un sistema de doble péndulo. El sistema es un robot manipulador de 2 grados de libertad que no se acciona y, por lo tanto, se comportará principalmente como un doble péndulo afectado por la gravedad. La única diferencia principal con un péndulo doble es que está...

Me gustaría desarrollar sistemas de resolución de EDO en GPU, en un entorno 'trivialmente paralelo'. Por ejemplo, haciendo un análisis de sensibilidad con 512 conjuntos de parámetros diferentes. Idealmente, quiero resolver ODE con un solucionador de paso de tiempo inteligente y adaptable como...

¿Cuáles son los métodos más modernos para la solución numérica de EDO con lado derecho discontinuo? Me interesan principalmente las funciones del lado derecho liso por partes, por ejemplo, signo. Estoy tratando de resolver la ecuación de un tipo siguiente: X˙v˙= v= { ( | FexternoEl | - |...

Trabajo bastante en biología matemática / epidemiología, donde la mayor parte del trabajo de modelado / ciencia computacional todavía está dominado por conjuntos de EDO, es cierto que a veces son conjuntos bastante elaborados de ellos. Una de las ventajas de estos modelos es que son bastante...

En la práctica, el tiempo de ejecución de resolver numéricamente un IVP menudo está dominado por la duración de la evaluación del lado derecho (RHS). Supongamos, por lo tanto, que todas las demás operaciones son instantáneas (es decir, sin costo computacional). Si el tiempo de ejecución general...

Hola y gracias por darme tiempo para mirar mi pregunta. Esta es una versión actualizada de mi pregunta que publiqué anteriormente en physics.stackexchange.com Actualmente estoy estudiando un condensador Bose-Einstein Condensador de excitón 2D y tengo curiosidad sobre el estado fundamental de este...

He estado usando GSL como base de muchas de mis simulaciones, pero es un poco exagerado para mis propósitos y define su propio tipo complejo por razones heredadas. En lugar de codificar mi propio solucionador ODE Runge-Kutta, que probablemente no sería muy eficiente, ¿hay algún solucionador ODE de...

Me pregunto: ¿cuál es el mejor algoritmo para resolver donde es una matriz real . A no es explícitamente dependiente del tiempo, generalmente escaso pero no necesariamente en bandas. Sus valores propios tienen partes reales no positivas. A también es diagonalizable, pero puede ser demasiado...

La esencia de mi pregunta es la siguiente: tengo un sistema de dos EDO. Uno tiene una restricción de valor inicial y el otro tiene una restricción de valor final. Esto puede considerarse como un sistema único con una restricción de valor inicial en algunas variables y una restricción de valor final...

Tengo un sistema lineal no homogéneo A x = bAx=b Ax=b donde es una matriz real n × n con n ≤ 4 . Se garantiza que el espacio nulo de A es de dimensión cero, por lo que la ecuación tiene un inverso único x = A - 1 b . Dado que el resultado entra en el lado derecho de una EDO, que tengo la...

Estoy intentando resolver una ecuación del tipo: (−∂2∂x2−f(x))ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Donde f(x)f(x)f(x) tiene un polo simple en 000 , para los valores propios y vectores propios más pequeños...

Estoy tratando de implementar el método Runge-Kutta de cuarto orden para resolver un ODE de primer orden en Python, es decir, . Entiendo cómo funciona el método, pero estoy tratando de escribir un algoritmo eficiente que minimice el número de veces quese calculaf(x,y)ya que esto es bastante...

Los algoritmos existentes para resolver ODE manejan funciones , dondey∈Rn. Pero en muchos sistemas físicos, la ecuación diferencial es autónoma, entoncesdyreyret= f( y, t )reyret=F(y,t)\frac{dy}{dt} = f(y, t)y∈ Rnortey∈Rnortey \in \mathbb R^n,y∈Rn, con latfuera. Con este supuesto simplificador,...