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 

Simplification d'équation différentielle

 
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
bendesarts



Inscrit le: 27 Fév 2009
Messages: 234

MessagePosté le: 24 Fév 2010 16:52    Sujet du message: Simplification d'équation différentielle Répondre en citant

Bonjour,

Je dispose de 3 équations dont des équations différentielles. Je voudrais simplement imposer une variable comme étant constant: diff(psi(t),t)=Omega et du coup, j'aimerais que les dérivées secondes de psi se suppriment. Est-ce quelqu'un aurait une petite idée ?
En vous remerciant d'avance


Code:
   // d  / d        \\           
   ||--- |--- psi(t)|| Inertie[p]
   \\ dt \ dt       //           

        // d  / d        \\     / d  / d        \\  \    \               
      + ||--- |--- psi(t)|| R + |--- |--- psi(t)|| e| R m| beta(t) = L_1_2
        \\ dt \ dt       //     \ dt \ dt       //  /    /               
/            2              /             2                 2  \    \       
|/ d        \               | / d        \      / d        \   |    |       
||--- psi(t)|  Inertie[p] - |-|--- psi(t)|  R - |--- psi(t)|  e| R m| beta(t)
\\ dt       /               \ \ dt       /      \ dt       /   /    /       

     / 2               \ / d  / d         \\     1                     2     
   + \R  m + Inertie[p]/ |--- |--- beta(t)|| = - - gamma theta(t) Omega  I[p]
                         \ dt \ dt        //     8                           

     1       / d         \           
   + - gamma |--- beta(t)| Omega I[p]
     8       \ dt        /           
       / d  / d        \\           
       |--- |--- psi(t)|| Inertie[p]
       \ dt \ dt       //           

            // d  / d        \\     / d  / d        \\  \           
          + ||--- |--- psi(t)|| R + |--- |--- psi(t)|| e| R m = M_1_2
            \\ dt \ dt       //     \ dt \ dt       //  /           
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
prof_simplet



Inscrit le: 12 Sep 2006
Messages: 86

MessagePosté le: 24 Fév 2010 17:00    Sujet du message: Répondre en citant

Peut-être en faisant une substitution dans ton équation que je note eq:

Code:

> eq:=subs(diff(psi(t),t)=Omega,eq);
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
bendesarts



Inscrit le: 27 Fév 2009
Messages: 234

MessagePosté le: 24 Fév 2010 17:25    Sujet du message: Répondre en citant

Ok merci pour votre conseil;
Par contre, par la suite, l'équation n'est pas simplifiée. J'ai essayé de faire un "simplify" mais celà n'a rien donné!!!

Code:
eq6li := (diff(Omega, t))*(Inertie[p]+R^2*m+m*R*e) = M_1_2
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
bendesarts



Inscrit le: 27 Fév 2009
Messages: 234

MessagePosté le: 25 Fév 2010 11:12    Sujet du message: Répondre en citant

Pour plus de précisions, voici le résultat que j'obtiens:


Code:
// d       \              // d       \     / d       \  \    \               
||--- Omega| Inertie[p] + ||--- Omega| R + |--- Omega| e| R m| beta(t) = L_1_2
\\ dt      /              \\ dt      /     \ dt      /  /    /               
/     2              /      2          2  \    \       
\Omega  Inertie[p] - \-Omega  R - Omega  e/ R m/ beta(t)

     / 2               \ / d  / d         \\     1                     2     
   + \R  m + Inertie[p]/ |--- |--- beta(t)|| = - - gamma theta(t) Omega  I[p]
                         \ dt \ dt        //     8                           

     1       / d         \           
   + - gamma |--- beta(t)| Omega I[p]
     8       \ dt        /           
               / d       \ /              2          \       
               |--- Omega| \Inertie[p] + R  m + m R e/ = M_1_2
               \ dt      /                                   


Comment faire pour que les dérivées de omega par rapport au temps se simplifient ?
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
prof_simplet



Inscrit le: 12 Sep 2006
Messages: 86

MessagePosté le: 25 Fév 2010 12:33    Sujet du message: Répondre en citant

Utilisez eval pour forcer Maple à évaluer la nouvelle expression.
Exemple:

Code:

> eq:=diff(diff(psi(t),t),t)*Inertie[p];
                          / 2        \
                          |d         |
                    eq := |--- psi(t)| Inertie[p]
                          |  2       |
                          \dt        /

> eq:=subs(diff(psi(t),t)=Omega,eq);

                           /d       \
                     eq := |-- Omega| Inertie[p]
                           \dt      /

> eval(eq);

                                  0
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
bendesarts



Inscrit le: 27 Fév 2009
Messages: 234

MessagePosté le: 25 Fév 2010 17:55    Sujet du message: Répondre en citant

Super effectivement.
La simplification fonctionne très bien.

Par contre, je fais les 6 lignes de codes suivantes :

eq1li:=eval(subs(diff(psi(t),t)=Omega,eq1li));
eq2li:=eval(subs(diff(psi(t),t)=Omega,eq2li));
eq3li:=eval(subs(diff(psi(t),t)=Omega,eq3li));
eq4li:=eval(subs(diff(psi(t),t)=Omega,eq4li));
eq5li:=eval(subs(diff(psi(t),t)=Omega,eq5li));
eq6li:=eval(subs(diff(psi(t),t)=Omega,eq6li));

En résultant, le résultat est réalisé uniquement sur les 3 dernières équations eq4;eq5,eq6:

/ 2 2 \
| / d \ / d \ |
|-|--- psi(t)| R - |--- psi(t)| e| m = X_1_2
\ \ dt / \ dt / /
// d / d \\ / d / d \\ \
||--- |--- psi(t)|| R + |--- |--- psi(t)|| e| m = Y_1_2
\\ dt \ dt // \ dt \ dt // /
/ 2 2 \
/ d / d \\ | / d \ / d \ |
-|--- |--- beta(t)|| R m + |-|--- psi(t)| R - |--- psi(t)| e| beta(t) m =
\ dt \ dt // \ \ dt / \ dt / /

/ d \
2 gamma |--- beta(t)| Omega I[p]
gamma theta(t) Omega I[p] \ dt /
-------------------------- - ------------------------------ - g m + Z_1_2
6 R 6 R
// d / d \\ / d / d \\ \
||--- |--- psi(t)|| R + |--- |--- psi(t)|| e| beta(t) R m = L_1_2
\\ dt \ dt // \ dt \ dt // /
/ 2 2 \
/ d / d \\ 2 | / d \ / d \ |
|--- |--- beta(t)|| R m - |-|--- psi(t)| R - |--- psi(t)| e| beta(t) R m =
\ dt \ dt // \ \ dt / \ dt / /
1 2 1 / d \
- - gamma theta(t) Omega I[p] + - gamma |--- beta(t)| Omega I[p]
8 8 \ dt /
/ d / d \\ // d / d \\ / d / d \\ \
beta(t) |--- |--- psi(t)|| + ||--- |--- psi(t)|| R + |--- |--- psi(t)|| e| R m =
\ dt \ dt // \\ dt \ dt // \ dt \ dt // /

M_1_2
0 = L_1_2
/ d / d \\ 2 / 2 2 \
|--- |--- beta(t)|| R m - \-Omega R - Omega e/ beta(t) R m =
\ dt \ dt //
1 2 1 / d \
- - gamma theta(t) Omega I[p] + - gamma |--- beta(t)| Omega I[p]
8 8 \ dt /
0 = M_1_2

Je ne comprends pas pourquoi la commande qui est effectuée sur les 3 premières lignes n'est pas effectuée.

PS:

Si mon message n'est pas clair, faites moi le savoir j'essayerais de le reformuler
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
prof_simplet



Inscrit le: 12 Sep 2006
Messages: 86

MessagePosté le: 26 Fév 2010 9:48    Sujet du message: Répondre en citant

Bonjour,
C'est peut-être la façon dont vous avez défini dans eq1, eq2,eq3, la dérivée seconde de psi qui pose problème: il faudrait écrire diff(diff(psi(t),t),t) et non pas diff(psi(t),t,t) ou diff(psi(t),t$2) pour que la substitution s'opère correctement.
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
bendesarts



Inscrit le: 27 Fév 2009
Messages: 234

MessagePosté le: 26 Fév 2010 14:19    Sujet du message: Répondre en citant

Ok merci pour les conseils mais en tout cas, ce qui est bizarre c'est que la façon de programmer les équations au niveau des dérivées de psi est la même pour les équations eq1,eq2, eq3 et les équations eq4,eq5,eq6.
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
bendesarts



Inscrit le: 27 Fév 2009
Messages: 234

MessagePosté le: 02 Mar 2010 0:18    Sujet du message: Répondre en citant

C'est bon j'ai trouvé ; J'ai refait mon fichier maple en faisant une execution de notation ">" à chaque ligne de calcul et j'ai pu trouvé!!! J'ai tout simplement l'expression qu'il faut mieux ne pas faire de trop grands paragraphes où l'on execute une dizaine de lignes de commande en bloc!!!
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.