|
Apprendre Maple Site dédié au logiciel de calcul formel Maple
|
Voir le sujet précédent :: Voir le sujet suivant |
Auteur |
Message |
Syracuse
Inscrit le: 01 Juil 2007 Messages: 14
|
Posté le: 23 Déc 2009 12:52 Sujet du message: système d'équations différentielles aux dérivées partielles |
|
|
Bonjour
A l'aide de la fonction pdsolve, j'ai réussi à résoudre un système d'équations différentielles aux dérivées partielles. Mais lorsque je souhaite ajouter des conditions initiales, j'obtiens le message d'erreur :
Error, (in pdsolve/numeric/process_PDEs) PDEs can only contain dependent variables with direct dependence on the independent variables of the problem, got {s(x, 0), s(0, t)}
voilà ce que j'ai tapé :
> eq1 := diff(s(x, t), t) = -.35*i(x, t)*s(x, t);
> cdt1 := {s(x, 0) = 5, s(0, t) = 0};
> eq2 := diff(i(x, t), t) = .35*i(x, t)*s(x, t)-.2*i(x, t)+.65*(diff(i(x, t), x));
> cdt2 := {i(0, t) = 0, i(x, 0) = 5};
> pdsolve({cdt1, eq1}, {eq2, cdt2}, numeric);
Voilà, merci pour votre aide |
|
Revenir en haut de page |
|
|
Syracuse
Inscrit le: 01 Juil 2007 Messages: 14
|
Posté le: 23 Déc 2009 13:28 Sujet du message: |
|
|
J'ai cherché un peu sur internet et finalement, en bidouillant un peu les équations précédentes, j'arrive à cela :
with(PDEtools);
eq := [diff(s(x, t), t) = -.35*i(x, t)*s(x, t), diff(i(x, t), t) = .35*i(x, t)*s(x, t)-.2*i(x, t)+.65*(diff(i(x, t), x))];
cdt := {s(x, 0) = x, i(x, 0) = 0, s(0, t) = .1}
sol := pdsolve(eq, cdt, [s, i], numeric, time = t, range = 0 .. 10)
et j'ai : module() export plot,plot3d,animate,value,settings; ... end module
Par contre, lorsque je souhaite visualiser les solutions, et que j'écris :
Splot := sol:-plot(s(x, t), t = 0, numpoints = 50, colour = red); Iplot := sol:-plot(i(x, t), t = 0, numpoints = 50, colour = red); plots[display]([Splot, Iplot])
et bien j'ai l'erreur :
Error, (in pdsolve/numeric/plot) unable to compute solution for t<0.:
matrix is singular
Error, (in pdsolve/numeric/plot) unable to compute solution for t<0.:
matrix is singular
Error, (in plots:-display) expecting plot structures but received: [Splot, Iplot]
Je comprends vraiment pas ce que Maple essaie de me dire :/ |
|
Revenir en haut de page |
|
|
zozo
Inscrit le: 03 Jan 2013 Messages: 125
|
|
Revenir en haut de page |
|
|
Syracuse
Inscrit le: 01 Juil 2007 Messages: 14
|
Posté le: 23 Déc 2009 17:48 Sujet du message: |
|
|
ok, merci pour le lien, je vais essayer de comprendre ce dont ils parlent. |
|
Revenir en haut de page |
|
|
|
|
Vous ne pouvez pas poster de nouveaux sujets dans ce forum Vous ne pouvez pas répondre aux sujets dans ce forum Vous ne pouvez pas éditer vos messages dans ce forum Vous ne pouvez pas supprimer vos messages dans ce forum Vous ne pouvez pas voter dans les sondages de ce forum
|
Développé par phpBB © 2001, 2006 phpBB Group
Traduction par : phpBB-fr.com
Apprendre Maple - ©
- Alain Le Stang - Navigation optimisée pour une résolution 1024 x 768.
|