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 

Procédure type "Dérivation formelle"

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



Inscrit le: 27 Oct 2006
Messages: 5

MessagePosté le: 27 Mai 2007 17:31    Sujet du message: Procédure type "Dérivation formelle" Répondre en citant

Bonjour !

j'ai bessoin de programer une procédure F un peu similaire a de la dérivation formelle.

elle prendait en argument des polynomes en D(y),... (D@@k)(y) et doit etre définit inductivement par les relation suivante :


si a et b sont des scalaire et u,v des polynome en D(y),... (D@@k)(y) alors:
f(au+bv)=aF(u)+bF(v)
F(u*v)=v*F(u)+u*F(v)


et l'initialisation ce fait par :
F(D(y))=0
et pour tous k>1 :
F((D@@k)(y)) = sum(binomial(n,k),(D@@k)(y)*(D@@(n-k))(y),k=1..n-1)


pour faire cela je pense que le plus simple serait de pouvoir manipuler un peu les arbres définissant les expressions... mais je ne sais pas si c'est possible avec maple ?


Merci d'avance !


NB : enfin... ceci dit si quelqu'un voit une expression de F directement comprehensible par maple je suis preneur moi, mon objectif etant d'étudié le noyaux de F.
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: 28 Mai 2007 10:07    Sujet du message: Répondre en citant

Bonjour, je ne l'ai jamais utilisé, mais je pense que la fonction Rule du package Calculus1 pourrait convenir pour votre problème.
J'ai trouvé sur maplesoft.com un maplet de dérivation pas à pas utilisant cette fonction:
http://www.maplesoft.com/applications/app_center_view.aspx?AID=1120
(s'inscrire d'abord en tant que membre)
A plus tard.
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: 28 Mai 2007 13:51    Sujet du message: Répondre en citant

Je définis F comme un opérateur linéaire par la fonction define, puis je donne les valeurs des F(D@@n)(y)) puis des F((D@@n)(y)*(D@@p)(y)).
On en déduira par linéarité les valeurs F(u,v) lorsque u,v sont des polynômes en D(y),....,(D@@k)(y)

Code:

> restart:
> define(F, 'linear');
> F(D(y)):=0:
> for n from 2 to 100 do
>     F((D@@n)(y)) := sum(binomial(n,k)*(D@@k)(y)*(D@@(n-k))(y),k=1..n-1) od:
> for n from 2 to 10 do
>   for p from 2 to 10 do
>      F((D@@n)(y)*(D@@p)(y)):=(D@@n)(y)*F((D@@p)(y))+(D@@p)(y)*F((D@@n)(y)) od od:
> F(-2*(D@@2)(y)+4*(D@@3)(y));

                            2             (2)
                     -4 D(y)  + 24 D(y) (D   )(y)

> F((D@@2)(y)*(D@@3)(y));

                    (2)    2            (3)         2
                6 (D   )(y)  D(y) + 2 (D   )(y) D(y)

> u:=-2*(D@@2)(y)+4*(D@@3)(y): v:=11*(D@@4)(y)-4*(D@@3)(y)+7*(D@@2)(y):
> F(u);

                            2             (2)
                     -4 D(y)  + 24 D(y) (D   )(y)

> F(v);

             (3)            (2)    2             (2)              2
   88 D(y) (D   )(y) + 66 (D   )(y)  - 24 D(y) (D   )(y) + 14 D(y)

> expand(F(expand(u*v)));   # calcul de F(u,v)

         (3)            (2)             (2)    3        (4)         2
  -368 (D   )(y) D(y) (D   )(y) - 132 (D   )(y)  - 44 (D   )(y) D(y)

                 (2)    2             (3)         2
         + 216 (D   )(y)  D(y) + 72 (D   )(y) D(y)

                (2)         2              (3)    2
         - 56 (D   )(y) D(y)  + 352 D(y) (D   )(y)

                 (3)       (2)    2         (4)            (2)
         + 264 (D   )(y) (D   )(y)  + 264 (D   )(y) D(y) (D   )(y)

> expand(u*F(v)+v*F(u));  # vérification

         (3)            (2)             (2)    3        (4)         2
  -368 (D   )(y) D(y) (D   )(y) - 132 (D   )(y)  - 44 (D   )(y) D(y)

                 (2)    2             (3)         2
         + 216 (D   )(y)  D(y) + 72 (D   )(y) D(y)

                (2)         2              (3)    2
         - 56 (D   )(y) D(y)  + 352 D(y) (D   )(y)

                 (3)       (2)    2         (4)            (2)
         + 264 (D   )(y) (D   )(y)  + 264 (D   )(y) D(y) (D   )(y)


C'est sans doute perfectible, mais c'est un premier pas.
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
Ksilver



Inscrit le: 27 Oct 2006
Messages: 5

MessagePosté le: 29 Mai 2007 11:02    Sujet du message: Répondre en citant

Merci beaucoup !!

je vais regarder ca des cette apres midi.
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 -> Programmation 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.