|
Apprendre Maple Site dédié au logiciel de calcul formel Maple
|
Voir le sujet précédent :: Voir le sujet suivant |
Auteur |
Message |
fenec74
Inscrit le: 11 Fév 2007 Messages: 2
|
Posté le: 12 Fév 2007 15:27 Sujet du message: Resolution systeme de 4 equa-diff |
|
|
Bonjour,
mon problème est que je n'ai pas réussi à résoudre(numeriquement) un système du second ordre non linéaire de 4 équation.Or , sur le livre où je les ais pris , ils affirment :"easy to solve on a computer".J'aimerais savoir si c quand meme possible.Voici le systeme, les variables sont a1,a2,a3,r, elles dépendent uniquement du temps.
eq1 := a3' = 3 + 4 cos(a3) + 5 r + 5 | r |
eq2 := a1' = 3 + 2 cos(a1) - 2 cos(a3) - 2 | cos(a3) |
eq3 := a2' = 2, 5 + 2 cos(a2) - 2 cos(a1) - 2 | cos(a1) |
eq4 := 50 r'' + r' + r = 5 tanh(
1/2 cos(a1) + 1/2 | cos(a1) | - 1/2 cos(a2) - 1/2 | cos(a2) |)
Merci d'avance.
Fenec |
|
Revenir en haut de page |
|
|
ALS
Inscrit le: 11 Sep 2006 Messages: 647
|
Posté le: 13 Fév 2007 11:45 Sujet du message: |
|
|
Bonjour, le code suivant:
Code: |
> eq1 := diff(a3(t),t) = 3 + 4* cos(a3(t)) + 5* r(t) + 5*abs( r(t));
> eq2 := diff(a1(t),t) = 3 + 2* cos(a1(t)) - 2* cos(a3(t)) - 2*abs(cos(a3(t)));
> eq3 := diff(a2(t),t) = 2.5 + 2* cos(a2(t)) - 2* cos(a1(t)) - 2*abs(cos(a1(t)));
> eq4 := 50*diff(r(t),t$2) + diff(r(t),t) + r(t) = 5*tanh(1/2* cos(a1(t)) + 1/2*abs(cos(a1(t))) - 1/2* cos(a2(t)) - 1/2*abs( cos(a2(t))));
> sys:={eq1,eq2,eq3,eq4}: vars:={a1(t),a2(t),a3(t),r(t)}:
> dsolve(sys,vars);
|
n'a pas l'air de donner de réponse.
Essayer peut-être les différentes options proposées dans l'aide pour la fonction dsolve.
A+ |
|
Revenir en haut de page |
|
|
fenec74
Inscrit le: 11 Fév 2007 Messages: 2
|
Posté le: 13 Fév 2007 20:16 Sujet du message: |
|
|
Bonsoir(et merci d'avoir répondu),dans dsolve j'ai essayé avec l'option "numeric", mais maple me retourne un message d'erreur.Mon erreur vient peut etre des conditions initiales,car étant donné le système je ne sais pas combien de valeurs initiales je doit fournir (à cause du second ordre sur r(t)).
Aussi , je me suis aperçu qu'on pouvait surement négliger le "tanh" ( en mettant la fonction identité à la place), ce qui arrange un peu les choses meme si je n'y arrive toujours pas.
@+ |
|
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.
|