Apprendre Maple Index du Forum Apprendre Maple
Site dédié au logiciel de calcul formel Maple
 
  Page d'accueilPage d'accueil   FAQFAQ    RechercherRechercher    Liste des MembresListe des Membres    Groupes d'utilisateursGroupes d'utilisateurs 
S'enregistrerS'enregistrer    ProfilProfil   Se connecter pour vérifier ses messages privésSe connecter pour vérifier ses messages privés   ConnexionConnexion 

système d'équations différentielles aux dérivées partielles

 
Poster un nouveau sujet   Répondre au sujet    Apprendre Maple Index du Forum -> Maple en analyse
Voir le sujet précédent :: Voir le sujet suivant  
Auteur Message
Syracuse



Inscrit le: 01 Juil 2007
Messages: 14

MessagePosté le: 23 Déc 2009 12:52    Sujet du message: système d'équations différentielles aux dérivées partielles Répondre en citant

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
Voir le profil de l'utilisateur Envoyer un message privé
Syracuse



Inscrit le: 01 Juil 2007
Messages: 14

MessagePosté le: 23 Déc 2009 13:28    Sujet du message: Répondre en citant

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
Voir le profil de l'utilisateur Envoyer un message privé
zozo



Inscrit le: 03 Jan 2013
Messages: 125

MessagePosté le: 23 Déc 2009 14:56    Sujet du message: Répondre en citant

Salut,
aucune idée de ce comportement de Maple, mais essaie peut-être de voir sur mapleprimes par exemple ici: http://www.mapleprimes.com/forum/solvingsystempdes
A+
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
Syracuse



Inscrit le: 01 Juil 2007
Messages: 14

MessagePosté le: 23 Déc 2009 17:48    Sujet du message: Répondre en citant

ok, merci pour le lien, je vais essayer de comprendre ce dont ils parlent.
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
Montrer les messages depuis:   
Poster un nouveau sujet   Répondre au sujet    Apprendre Maple Index du Forum -> Maple en analyse Toutes les heures sont au format GMT + 2 Heures
Page 1 sur 1

 
Sauter vers:  
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


phpBB

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.