Questions marquées «ode»

J'ai plusieurs problèmes d'optimisation globale non convexe difficiles à résoudre. Actuellement, j'utilise la boîte à outils Optimization de MATLAB (en particulier, fmincon()avec algorithm = 'sqp'), ce qui est assez efficace . Cependant, la majeure partie de mon code est en Python et j'aimerais...

Je suis actuellement à la recherche de méthodes parallèles d'intégration ODE. Il existe de nombreux ouvrages anciens et anciens décrivant un large éventail d'approches, mais je n'ai trouvé aucun sondage récent ni article récapitulatif décrivant le sujet en général. Il y a le livre de Burrage [1],...

Cette question a été migrée à partir de Mathematics Stack Exchange car il est possible d'y répondre sur Computational Science Stack Exchange. Migré il y a 6 ans . Je suis sur le point de commencer à travailler sur une bibliothèque de logiciels de solveurs numériques ODE , et je n'arrive pas à...

Comment le problème gravitationnel des n-corps peut-il être résolu numériquement en parallèle? Un compromis précision-complexité est-il possible? Comment la précision influence-t-elle la qualité du modèle?

Dans ce commentaire, j'ai écrit: ... intégrateur SciPy par défaut, qui, je suppose, utilise uniquement des méthodes symplectiques. dans lequel je me réfère à SciPy's odeint, qui utilise soit une "méthode non rigide (Adams)" ou une "méthode rigide (BDF)". Selon la source : def odeint(func, y0, t,...

Considérons un IVP pour le système ODE , y ( x 0 ) = y 0 . Le plus souvent, ce problème est considéré comme rigide lorsque la matrice de Jacobi ∂ fy′= f( x , y)y′=F(X,y)y'=f(x,y)y( x0) = y0y(X0)=y0y(x_0)=y_0a à lafois desvaleurs propres avec une très grande partie réelle négative et des valeurs...

J'étais simplement curieux de savoir pourquoi les méthodes Runge – Kutta d'ordre élevé (c'est-à-dire supérieures à 4) ne sont presque jamais discutées / utilisées (du moins à ma connaissance). Je comprends que cela nécessite un temps de calcul plus long par étape (par exemple RK14 avec une étape...

Je voudrais consolider la résolution de systèmes d'ODE sur des GPU, dans un cadre «trivialement parallélisable». Par exemple, faire une analyse de sensibilité avec 512 jeux de paramètres différents. Idéalement, je veux faire une résolution ODE avec un solveur d'horodateur adaptatif intelligent...

quelles sont les méthodes de pointe pour la solution numérique des ODE à côté droit discontinu? Je suis surtout intéressé par les fonctions du côté droit par morceaux, par exemple le signe. J'essaie de résoudre l'équation d'un type suivant: X˙v˙= v= { ( | Fexterne| - | Ffriction| )signe(...

Je souhaite simuler le comportement d'un système à double pendule. Le système est un robot manipulateur à 2 degrés de liberté qui n'est pas actionné et se comportera donc principalement comme un double pendule affecté par la gravité. La seule différence principale avec un pendule double est qu'il...

En pratique, le temps d'exécution de la résolution numérique d'un IVP est souvent dominé par la durée de l'évaluation du côté droit (RHS) . Supposons donc que toutes les autres opérations sont instantanées (c'est-à-dire sans coût de calcul). Si le temps d'exécution global pour résoudre l'IVP est...

Je travaille assez fortement en biologie mathématique / épidémiologie, où la plupart des travaux de modélisation / science informatique sont encore dominés par des ensembles d'OD, certes parfois assez élaborés. L'un des avantages de ces modèles est qu'ils sont assez faciles à décrire et à...

L'essence de ma question est la suivante: j'ai un système de deux ODE. L'un a une contrainte de valeur initiale et l'autre une contrainte de valeur finale. Cela peut être considéré comme un système unique avec une contrainte de valeur initiale sur certaines variables et une contrainte de valeur...

Je me demande: quel est le meilleur algorithme pour résoudre réurét= A uréurét=UNEu\begin{equation} \frac{du}{dt} = Au \end{equation} OùUNEUNEAest unematricen × nn×nn\times nréelle. A n'est pas explicitement dépendant du temps, généralement clairsemé mais pas nécessairement en bandes. Ses valeurs...

Salut et merci d'avoir donné le temps de regarder ma question. Ceci est une version mise à jour de ma question que j'ai publiée plus tôt dans physics.stackexchange.com J'étudie actuellement un condensateur à excitons 2D Bose-Einstein Condensate et je suis curieux de connaître l'état fondamental de...

J'ai utilisé GSL comme base de plusieurs de mes simulations, mais c'est un peu exagéré pour mes besoins et il définit son propre type complexe pour des raisons héritées. Plutôt que de coder mon propre solveur Runge-Kutta ODE, qui ne serait probablement pas très efficace, existe-t-il des solveurs...

J'essaie de résoudre une équation du type: (−∂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) Où a un pôle simple à 0 , pour les N valeurs propres et vecteurs propres les plus petits . Les conditions aux...

J'ai un système linéaire inhomogène A x = bUNEX=b Ax=b où UNEUNEA est une matrice n × nn×nn\times n réelle avec n ≤ 4n≤4n\leq 4 . L'espace nul de UNEUNEA est garanti d'être de dimension nulle, l'équation a donc un inverse unique x = A- 1bX=UNE-1bx=A^{-1} b . Puisque le résultat entre dans le côté...

J'essaie d'implémenter la méthode Runge-Kutta de quatrième ordre pour résoudre un ODE de premier ordre en Python ie . Je comprends comment fonctionne la méthode, mais j'essaie d'écrire un algorithme efficace qui minimise le nombre de fois oùf(x,y)est calculé car cela coûte assez cher. On m'a dit...

Les algorithmes existants pour résoudre les ODE gèrent les fonctions , oùy∈Rn. Mais dans de nombreux systèmes physiques, l'équation différentielle est autonome, doncdyréyrét= f( y, t )dydt=f(y,t)\frac{dy}{dt} = f(y, t)y∈ Rny∈Rny \in \mathbb R^n,y∈Rn, sanst. Avec cette hypothèse simplificatrice,...