équation différentielle raide
pom
-
SKZ -
SKZ -
Bonjour, je dois résoudre un système algébro-différentielle qui est extrêmement raide. Toutes mes méthodes numériques plantent ou ne prennent pas en compte la raideur.
Sur le net, j'ai vu des codes comme RADAU, DASSL etc...
En connaissez-vous les algos ou avez-vous des renseignements pour que je puisse les coder ?
Merci beaucoup.
Pom
Sur le net, j'ai vu des codes comme RADAU, DASSL etc...
En connaissez-vous les algos ou avez-vous des renseignements pour que je puisse les coder ?
Merci beaucoup.
Pom
A voir également:
- équation différentielle raide
- Editeur d'équation - Télécharger - Vie quotidienne
- Formule équation - Télécharger - Études & Formations
- Equation sinequanone - Télécharger - Calcul & Conversion
- Pourquoi je ne peux pas insérer une équation dans word - Forum Microsoft Office
- Télécharger microsoft equation 3.0 gratuitement ✓ - Forum Windows
1 réponse
Qu'entend tu par raideur ? Pente de la courbe, donc du df ? (dx étant le pas de temps).
Toutes les méthodes numériques plantent
A tu essayé TOUTES les méthodes
Essaye une méthode précise (RK4 par exemple) avec un pas de temps variable : quand la valeur d'un terme différenciel (un df) augmente, tu diminue d'autant (proportionellement) le dx (ou le dt, selon...)
En utilisant toujours le df max, tu garde ta précision, maintenant si tu dépasse les précisions admisent par les flottants, y'a toujours des bibliothèques à précisions arbitraire (la taille d'un flottant n'est limité que par la mémoire adressable), mais c'est ultra plus lent...
Et si ton graphe de dépendance entre équation (un noeud<->une équa ; un arc orienté<-> une dépendance) à des cycles à chier partout, t'es pas sorti de l'auberge.
Va voir eventuellement du coté de Maple ou soft de calcul formel, mais je promet rien, ces outils on leur limites aussi...
Après reste la bonne résolution à la main !!
Toutes les méthodes numériques plantent
A tu essayé TOUTES les méthodes
Essaye une méthode précise (RK4 par exemple) avec un pas de temps variable : quand la valeur d'un terme différenciel (un df) augmente, tu diminue d'autant (proportionellement) le dx (ou le dt, selon...)
En utilisant toujours le df max, tu garde ta précision, maintenant si tu dépasse les précisions admisent par les flottants, y'a toujours des bibliothèques à précisions arbitraire (la taille d'un flottant n'est limité que par la mémoire adressable), mais c'est ultra plus lent...
Et si ton graphe de dépendance entre équation (un noeud<->une équa ; un arc orienté<-> une dépendance) à des cycles à chier partout, t'es pas sorti de l'auberge.
Va voir eventuellement du coté de Maple ou soft de calcul formel, mais je promet rien, ces outils on leur limites aussi...
Après reste la bonne résolution à la main !!
