Discussion :

[Nath31]

Bonjour,
je suis en thèse et j'ai un problème sur Matlab, je suis une utilisatrice novice et je veux résoudre un système de 3 équa différentielles et de 9 équa algébrique. Pour cela j'ai utilisé la méthode ODE15s mais la fenêtre de commande me dit que mon système doit avoir un "index" supérieur à 1... Nulle en math et en Matlab je ne sais pas ce que ça veut dire... quelqu'un peut il m'éclairer? Merci d'avance

[tug83]

bonjour la thése!!!
ce lien pourrait téclaircir sur ton message d'erreur:
http://www.mathworks.com/support/sol...lution=1-15S5M

[Nath31]

merci je fonce sur le lien... quelle rapidité!

[tug83]

j'espère que ça puisse t'aider parce que j'ai rien capté aux histoires d'index!

[Nath31]

ben en fait pour ce que j'en ai lu pour l'instant (mais je crois que je reprendrai ça demain à tête reposée) ça ressemble beaucoup à ce qu'il y a dans l'aide de matlab.... mais bon je me réserve le doc en pdf pour demain!

[DenisLorrain]

Pourrais-tu poster le message d'erreur exact ?

[Nath31]

voilà le message d'erreur :

??? Error using ==> funfun\private\daeic12
This DAE appears to be of index greater than 1.

Error in ==> ode15s at 394
[y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...

Error in ==> partetglobal6 at 175
[t,y] = ode15s(@f,tspan,y0,options);

[tug83]

ton message d'erreur est discuté ici:
http://www.mathworks.com/support/sol...lution=1-19Z26

[Nath31]

oulala... ça viendrait donc de mes conditions limites??? Il faut que je regarde ça de plus prêt... apparemment mon système est aussi mal conditionné c'est à dire que j'ai des valeurs pouvant avec de 1e3 à 1e-10... donc système très instable à résoudre...

[tug83]

autant dire que 3 ans pour faire la thèse c'est ce qu'il faut

[Nath31]

oui tu as tout compris.... en plus j'avoue le calcul numérique n'est pas ma tasse de thé et j'ai du mal à comprendre quelle méthode utiliser et comment ... alors la subtilité des messages d'erreur...
Mais ce forum m'aide beaucoup!

[DenisLorrain]

Un conseil, fais une petite étude biblio de ton problème. Tu trouveras sans doute des articles qui traitent d'un sujet equivalent.

[salseropom]

Salut, pour le peu que je connaisse des DAE, je crois que tout DAE d'ordre >=1 peut se transformer en DAE d'ordre 1 quitte à augmenter la dimension du sysème.

J'ai moi-même des ODE et des DAE mais grâce à mes conditions algébriques issues de mon DAE j'ai pu simplifier des équations de mon ODE et donc du coup diminuer la dimension de mon système d'ODE

[Nath31]

OK... mon problème est que je suis très très nulle en math et que je ne sais même pas identifier l'ordre de mon système... en réalité je pensais qu'il était d'ordre 1 et c'est l'erreur de Matlab qui m'a fait douter...

[Nath31]

bonjour je sais que ça fait très longtemps que j'ai posté mais je voulais vous dire que j'avais enfin compris ce que voulais dire le message d'erreur pour lequel j'avais posté au départ! En fait mon système d'équation présentais des équations "redondantes" et donc le déterminant était "nul". Voilà donc ça venait de mon système d'équations. Je n'ai plus qu'à trouver une autre équation et recoder pour voir si ça marche! en tout cas encore une fois merdi pour votre aide!