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 

intégrer une matrice

 
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
nekcorp



Inscrit le: 08 Jan 2007
Messages: 5

MessagePosté le: 11 Avr 2007 19:03    Sujet du message: intégrer une matrice Répondre en citant

Bonjour voila je dois intégrer une matrice qui se résume simplement a intégrer chaqu'un des termes de la matrice. Cependant voila il n'i a pas de fonction toute faite pour intégrer une matrice et il faut pour chaque intégrale appeler le terme de la matrice et dire de l'intégrer, résultat je doit me taper 18 ligne de commande répétitive ce qui m'enchante pas vraiment. Je souhaiterais réslaiser le calcul a l'aide d'une boucle For faisant intervenir le i et j correspondant a la ligne et la colonne de la matrice. Voila ce que j'ai pour l'instant :
>q:=evalm(linalg[matrix](4,1,[psi1,psi2,psi3,psi4]));
> qT:=evalm(linalg[matrix](1,4,[psi1,psi2,psi3,psi4]));
> Phi:=evalm(q&*qT);
> L:=450;
> psi1:=Z->16/L^4*Z^4-32/L^3*Z^3+16/L^2*Z^2;
> psi2:=Z->16/L^4*Z^5-40/L^2*Z^4+32/L^2*Z^3-8/L*Z^2;
> psi3:=Z->-24/L^5*Z^5+52/L^4*Z^4-34/L^3*Z^3+7/L^2*Z^2;
> psi4:=Z->4/L^4*Z^5-8/L^3*Z^4+5/L^2*Z^3-1/L*Z^2;
> Phi[1,1]:=int(psi1(z)*psi1(z),z=0..168 .6);
> Phi[2,2]:=int(psi2(z)*psi2(z),z=0..168 .6);
> Phi[3,3]:=int(psi3(z)*psi3(z),z=0..168 .6);
> Phi[4,4]:=int(psi4(z)*psi4(z),z=0..168 .6);
> Phi[1,2]:=int(psi1(z)*psi2(z),z=0..168 .6);
> Phi[1,3]:=int(psi1(z)*psi3(z),z=0..168 .6);
> Phi[1,4]:=int(psi1(z)*psi4(z),z=0..168 .6);
> Phi[2,3]:=int(psi2(z)*psi3(z),z=0..168 .6);
> Phi[2,4]:=int(psi2(z)*psi4(z),z=0..168 .6);
> Phi[3,4]:=int(psi3(z)*psi4(z),z=0..168 .6);
> Phi[4,3]:=int(psi4(z)*psi3(z),z=0..168 .6);
> Phi[2,1]:=int(psi2(z)*psi1(z),z=0..168 .6);
> Phi[3,1]:=int(psi3(z)*psi1(z),z=0..168 .6);
> Phi[4,1]:=int(psi4(z)*psi1(z),z=0..168 .6);
> Phi[3,4]:=int(psi3(z)*psi4(z),z=0..168 .6);
> Phi[3,2]:=int(psi3(z)*psi2(z),z=0..168 .6);
> Phi[4,2]:=int(psi4(z)*psi2(z),z=0..168 .6);

Ma fonction psi est fonction de z.

En gros je souhaite obtenir un truc du genre :

>q:=evalm(linalg[matrix](4,1,[psi1,psi2,psi3,psi4]));
> qT:=evalm(linalg[matrix](1,4,[psi1,psi2,psi3,psi4]));
> Phi:=evalm(q&*qT);
> L:=450;

for i j from ..........
Phi[i,j]:=int(psi[i](z)*psi[j](z),z=0..168.6);

Mon problème est vraiment sérieux car je dois faire la même opération pour une matrice bien plus grande.

Cordialement.
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: 12 Avr 2007 8:47    Sujet du message: Répondre en citant

Bonjour, la procédure intMatrix(M,x,a,b) donne comme résultat la matrice obtenue en intégrant tous les coefficients de la matrice M, par rapport à la variable x, entre a et b.

Code:

> with(linalg):
 
> intMatrix:=proc(M::matrix, x::name, a, b)
> local i,j,IM,L;
>   IM:=NULL:
>   for i to rowdim(M) do
>      L:=NULL: 
>      for j to coldim(M) do
>         L:=L,simplify(int(M[i,j],x=a..b))
>      end do:
>      IM:=IM,[L]
>   end do:
>   evalm(matrix([IM]))
> end proc:
 
> psi1:=Z->16/L^4*Z^4-32/L^3*Z^3+16/L^2*Z^2;
> psi2:=Z->16/L^4*Z^5-40/L^2*Z^4+32/L^2*Z^3-8/L*Z^2;
> psi3:=Z->-24/L^5*Z^5+52/L^4*Z^4-34/L^3*Z^3+7/L^2*Z^2;
> psi4:=Z->4/L^4*Z^5-8/L^3*Z^4+5/L^2*Z^3-1/L*Z^2;
 
> q:=evalm(matrix(4,1,[psi1(Z),psi2(Z),psi3(Z),psi4(Z)]));
> qT:=evalm(matrix(1,4,[psi1(Z),psi2(Z),psi3(Z),psi4(Z)]));
> Phi:=evalm(q&*qT);
 
> intMatrix(Phi,Z,0,168.6);
 


Ceci devrait alléger quelque peu votre code.
A+
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
nekcorp



Inscrit le: 08 Jan 2007
Messages: 5

MessagePosté le: 12 Avr 2007 10:24    Sujet du message: Répondre en citant

Merci beaucoup pour ton aide.
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.