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 

Resolution systeme de 4 equa-diff

 
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
fenec74



Inscrit le: 11 Fév 2007
Messages: 2

MessagePosté le: 12 Fév 2007 15:27    Sujet du message: Resolution systeme de 4 equa-diff Répondre en citant

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



Inscrit le: 11 Sep 2006
Messages: 647

MessagePosté le: 13 Fév 2007 11:45    Sujet du message: Répondre en citant

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



Inscrit le: 11 Fév 2007
Messages: 2

MessagePosté le: 13 Fév 2007 20:16    Sujet du message: Répondre en citant

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
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.