Résolution d'une ODE

Bonjour à tous,

Je dois resoudre une ODE. J'ai rentrer mon chargement, ma matrice et mon fichier principal. Le probleme est que si je change les tolerances (abstol et reltol j'ai des résultats completement différents). Les composantes du vecteur X au final sont de l'ordre de la centaine pour les 7 premières composantes, et de l'ordre de 1e-3/1e-4 pour les 9 autres. Quelles sont les deux tolérances à prendre?

Le second soucis est dans la construction de ma matrice M. Je dois mémoriser la variable Max1 si lors d'un pas, la 7e variable du vecteur X dépasse l'ancienne valeur de Max1 et que le pas est positif . Quelqu'un peut m'aider sur ce point? (j'arrive à memoriser le max1 mais lorsque mon pas est négatif, max1 n'est pas à la bonne valeur).

Je vous donne mon code pour que vous voyiez à quoi cela ressemble.

Merci d'avance pour votre aide. Je galère depuis 1mois!

Bon courage à la lecture de mes codes.

Olivier

Salut, ta fonction chargement est indépendante de X. C'est normale ? Es-tu sûr que tu ne pourrais pas l'intégrer à la main et donc avoir une solution exacte ?

Pour tes pb de abstol et rtol, je n'ai pas exécuté ton prgm mais à ce que je vois, tu as des problèmes numériques. A cela, qu'une seule solution : changer de méthode d'intégration. As-tu un problème raide ? Essaye ode23tb, ode23s, ode15s (s pour stiff==raide)

Pour ma part, aucune des routines de matlab n'est capable de résoudre mon système d'équa diff (problème raide, matlab broutait à chaque fois). J'étais obligé de changer de solveur (j'ai chopé le solveur CVode, qui est extra. Il a même une interface avec Matlab si tu veux, sinon tu fais tout en C)

Je ne peux pas t'en dire plus...

Oui le chargement est indépendant de X.
J'ai déja fait un code avec la méthode d'euler. J'aimerai utiliser une fonction ode de matlab pour voir si les résultats sont proches. Avec certains paramètres Reltol et abstol cela se corrèle à peu près.

Pour le probleme de raideur, je ne pense pas que ce soit un probleme raide. Mes variables ne sont pas censer varier brusquement. Elle varie progressivement (mon chargement est linéaire). Par contre si j'utilise ode45, je n'arrive pas du tout au meme resultat qu'ode23. C'est assez bizarre non?

A +

Olivier