|
Apprendre Maple Site dédié au logiciel de calcul formel Maple
|
Voir le sujet précédent :: Voir le sujet suivant |
Auteur |
Message |
bendesarts
Inscrit le: 27 Fév 2009 Messages: 234
|
Posté le: 28 Fév 2009 12:20 Sujet du message: un peu d'automatique - partie réelle d'un complexe |
|
|
Bonjour à tous,
Je cherche à programmer sous maple les racines d'un polynome qui sont complexes purs. C'est pour un calcul de lieu d'Evans en automatique.
Voilà le calcul réalisé à la main
Voilà ce que j'ai fait sous maple mais qui ne fonctionne pas. J'aimerais réaliser sous maple le calcul que j'ai fait à la fin car je dois après utiliser ce type de calcul pour des fonctions plus complexes.
Je vous remercie d'avance pour votre aide
Bien cordialement,
Benjamin |
|
Revenir en haut de page |
|
|
ALS
Inscrit le: 11 Sep 2006 Messages: 647
|
Posté le: 01 Mar 2009 9:25 Sujet du message: Re: un peu d'automatique - partie réelle d'un complexe |
|
|
Bonjour, on doit aussi déclarer les autres variables comme réelles.
Voici la correction que je vous propose:
Code: |
> restart:
> assume(omega,real): assume(Km,real): assume(kp,real): assume(Ti,real): assume(Tn,real): assume(xi,real):
> eqcar:=Km*kp+Ti*s*(1+2*xi*Tn*s+Tn^2*s^2);
eqcar:=subs(s=I*omega,eqcar);
2 2
eqcar := Km~ kp~ + Ti~ s (1 + 2 xi~ Tn~ s + Tn~ s )
eqcar := Km~ kp~
2 2
+ Ti~ omega~ (1 + 2 I xi~ Tn~ omega~ - Tn~ omega~ ) I
> re:=Re(eqcar); im:=Im(eqcar);
2
re := Km~ kp~ - 2 Ti~ omega~ xi~ Tn~
2 2
im := Ti~ omega~ (1 - Tn~ omega~ )
> kp:=solve(re,kp);
2
2 Ti~ omega~ xi~ Tn~
kp := ---------------------
Km~
> solve(im,omega); omega:=%[2];
1 1
0, ---, - ---
Tn~ Tn~
1
omega := ---
Tn~
|
A plus tard sur ce forum. |
|
Revenir en haut de page |
|
|
bendesarts
Inscrit le: 27 Fév 2009 Messages: 234
|
Posté le: 01 Mar 2009 12:05 Sujet du message: |
|
|
Merci beaucoup pour votre aide
Par contre, j'ai encore un petit soucis pour finaliser le calcul ?
Je m'explique. La dernière étape est de remplacer omega que j'ai pris égal à 1/ Tn (solution physique) dans l'expression de kp
J'ai réaliser comme ceci :
subs(omega=1/Tn,kp);
Mais apparemment, j'obtiens ce résultat:
Autrement dit, maple n'a fait aucun remplacement.
Rappel :voilà où j'en suis dans les lignes de codes :
Benjamin |
|
Revenir en haut de page |
|
|
ALS
Inscrit le: 11 Sep 2006 Messages: 647
|
Posté le: 02 Mar 2009 8:24 Sujet du message: |
|
|
Bonjour, faites à la place
kp:=subs(omega=1/Tn,kp);
ALS |
|
Revenir en haut de page |
|
|
bendesarts
Inscrit le: 27 Fév 2009 Messages: 234
|
Posté le: 04 Mar 2009 14:50 Sujet du message: |
|
|
Ici, par contre, j'obtiens bien kp mais le w et le 1/Tn ne se simplifie pas dans l'expression de kp
Rappel des lignes de codes :
|
|
Revenir en haut de page |
|
|
ALS
Inscrit le: 11 Sep 2006 Messages: 647
|
Posté le: 05 Mar 2009 13:17 Sujet du message: |
|
|
Bonjour, j'ai légèrement modifié la fin:
Code: |
> restart:
> assume(omega,real): assume(Km,real): assume(kp,real): assume(Ti,real): assume(Tn,real): assume(xi,real):
> eqcar:=Km*kp+Ti*s*(1+2*xi*Tn*s+Tn^2*s^2);
> eqcar:=subs(s=I*omega,eqcar);
> re:=Re(eqcar); im:=Im(eqcar);
> kp:=solve(re,kp);
> solve(im,omega);
2 2
eqcar := Km~ kp~ + Ti~ s (1 + 2 xi~ Tn~ s + Tn~ s )
eqcar := Km~ kp~
2 2
+ Ti~ omega~ (1 + 2 I xi~ Tn~ omega~ - Tn~ omega~ ) I
2
re := Km~ kp~ - 2 Ti~ omega~ xi~ Tn~
2 2
im := Ti~ omega~ (1 - Tn~ omega~ )
2
2 Ti~ omega~ xi~ Tn~
kp := ---------------------
Km~
1 1
0, ---, - ---
Tn~ Tn~
> kp:=subs(omega=%[2],kp);
2 Ti~ xi~
kp := ---------
Tn~ Km~
> omega:=%%[2];
1
omega := ---
Tn~
|
Et voilà!! |
|
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.
|