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 

Coefficients

 
Poster un nouveau sujet   Répondre au sujet    Apprendre Maple Index du Forum -> Séquences, listes, ensembles, tables ou tableaux...
Voir le sujet précédent :: Voir le sujet suivant  
Auteur Message
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 26 Mai 2005 0:38    Sujet du message: Coefficients Répondre en citant
Bonjour, devant créer une base orthonormée dans Rn avec un produit scalaire, j'avais pour idée de faire un petit programme pour m'entrainer. Il faudrait donc que je puisses creer des coefficients pour chaque vecteur, je pensais donc faire : e0:=1; e1:=X+a11; e2=X^2+a21X+a20 etc soit ei=Xî+somme ( aij*x^j,j variant de 0 a i-j) en maple ca donnerait :

Citation:
e[i]:=x^i+sum(a[i][j]*x^j,j=0..i-j);


mais voila lorsque je fait

Citation:
e[2];
il me repond e2 , sans m'écrire le polynome et ses coefficients,du coup le produit scalaire ne donne rien non plus ...
si quelqu'un pouvait m'aider, je ne sais pas du tout si la syntaxe est la bonne. Sad

Merci Smile
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 26 Mai 2005 16:23    Sujet du message: Répondre en citant
oui exact je m'en suis rendu compte pendant la nuit Smile
cela dit meme avec [i,j] j'ai toujours pareil;
le produit de e[1]*e[2]; donne e1e2 par exemple...
je peux lui demander la version developpée ?
le produit scalaire ne marche toujours pas
merci
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 26 Mai 2005 17:15    Sujet du message: Répondre en citant
ok super merci je pensais qu'en ecrivant juste e[i]=... ca marchait mais faut les calculer avant pour qu'il les mette en memoire comme definis, merci beaucoup Smile
ou alors : e:=unapply(x^i+sum(a[i,j]*x^j,j=0..i-1),i);

mais je sais pas si ca donnera pareil a la fin Smile merci en tout cas
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 26 Mai 2005 18:41    Sujet du message: Répondre en citant
super merci Smile
j'ai reussi a faire ma procédure elle marche nickel, enfin jusqu'a R5[X], apres a R6[X] c'est tres tres long, la procédure n'est pas du tout optimisée je suppose, avec toutes les boucles internes Smile

Citation:

base:=proc(ps,n)
local i,e,s,eq,k,h,l,eqns,m,v;
e[0]:=1;
s:=e[0];
eq:=NULL;
v:=e[0];

for i from 1 to n
do
e[i]:=x^i+sum(a[i,j]*x^j,j=0..i-1):
s:=s,e[i];
h[0,i]:=ps(e[0],e[i]):
eq:=eq,h[0,i];
od;

for k from 1 to n-1
do:

for l from k+1 to n
do
h[l,k]:=ps(e[l],e[k]):
eq:=eq,h[l,k]:
od:

od;

eqns:={eq};
solve(eqns);
assign(%);

for m from 2 to n
do
h:=[s][m]/sqrt(ps([s][m],[s][m]));
v:=v,h;
od;
[v];

end:


Si quelqu'un a des remarques pour l'optimiser, c'est ma seconde procédure ( la premiere etait une permettant de calculer le produit de matrices , l'equivalent de &* Smile ) Wink
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 27 Mai 2005 12:44    Sujet du message: Répondre en citant
un dernier petit probleme, lorsque je fais base(f,3) et que j'affiche [s] jai la base orthogonale de R3[X] , voulant la base orthonormale j'ai rajouté la derniere boucle afin de normer les vecteurs mais lorsque je fais base(f,3) il s'affiche la base de R2[X], j'ai du me tromper dans les indices mais je vois pas ou Rolling Eyes
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 27 Mai 2005 14:28    Sujet du message: Répondre en citant
super merci beacoup ca roule Smile
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 27 Mai 2005 14:37    Sujet du message: Répondre en citant
en fait non ca marche pas avec les modifications, il m'affiche des a10 a20 etc...
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 27 Mai 2005 14:39    Sujet du message: Répondre en citant
et la je comprends plus rien, meme avec mon ancienne procédure, lorsque e fait base(f,1) ca marche, base(f,2) idem mais si je fais restart et directement base(f,2) ca m'affiche les a0 a1 etc... alors qu'avant ca marchait
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
dureiken



Inscrit le: 26 Mai 2005
Messages: 9
MessagePosté le: 27 Mai 2005 16:48    Sujet du message: Répondre en citant
super merci beaucoup c'est impec la !
les voies du Maple sont impenetrables Smile
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 -> Séquences, listes, ensembles, tables ou tableaux... 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.