valeurs et vecteurs propres avec maxima

Tout ce qui concerne les différents logiciels mathématiques (Matlab, Maxima, Scilab, R, Maple, Mathematica, ...) et les calculatrices scientifiques.
Règles du forum
Merci d'éviter le style SMS dans vos messages et de penser à utiliser la fonction Recherche avant de poster un message. Pour joindre des fichiers à vos messages, consulter ce sujet.
> Penser à utiliser les balises Code pour poster du code.
> Préciser le nom du logiciel concerné dans le titre de votre message.

valeurs et vecteurs propres avec maxima

Messagepar lassaib » Samedi 12 Janvier 2013, 12:37

Bonjour, je veux calculer les valeurs et vecteurs propres d'une matrice ( a*v = e*s*v) où s est une matrice, e et v les valeurs et vecteurs propres respectifs. Mon code tourne mais ça ne termine pas les calculs sur la dernière commande.
Y a t il un grain de sable dans mon code ? Merci de votre aide.
Code: Tout sélectionner

/* [wxMaxima: comment start ]
Methode des valeurs propres
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
ze : 1.75$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
f : exp(-ze*s)*[1,s,u,s**2,t**2,u**2,s*u,u*t**2,u**3,(t*u)**2]$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Element de matrice de recouvrement
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
s1[n,m]:= f[n]*f[m]*u*(s**2-t**2)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
s[n,m]:=integrate(integrate(integrate(s1[n,m],t,0,u),u,0,s),s,0,inf)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Energie cinetique
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
k1[n,m]:=(diff(f[n],s)*diff(f[m],s)+diff(f[n],t)*diff(f[m],t)+
diff(f[n],u)*diff(f[m],u))*u*(s**2-t**2)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
k[n,m]:=integrate(integrate(integrate(k1[n,m],t,0,u),u,0,s),s,0,inf)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
p1[n,m]:=(diff(f[n],s)*diff(f[m],u)+diff(f[n],u)*diff(f[m],s))*s*(u**2-t**2) +
(diff(f[n],t)*diff(f[m],u)+diff(f[n],u)*diff(f[m],t))*t*(s**2-u**2)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
p[n,m]:=integrate(integrate(integrate(p1[n,m],t,0,u),u,0,s),s,0,inf)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Energie potentielle
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
z : 2$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
v1[n,m]:=f[n]*f[m]*(s**2-t**2-4*z*s*u)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
v[n,m]:=integrate(integrate(integrate(v1[n,m],t,0,u),u,0,s),s,0,inf)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Hamiltonien
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
h[n,m]:=k[n,m]+p[n,m]+v[n,m]$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Creation des matrices
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
hh : genmatrix(h,10,10);
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
ss : genmatrix(s,10,10);
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Matrice final
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
iss : invert(ss)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
aa : iss*hh$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Valeurs et vecteurs propres de aa vals : solve([determinant(aa)=0],[e]);
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
[vals,vecs] : eigenvectors(aa);
/* [wxMaxima: input   end   ] */


lassaib
Déca-utilisateur
 
Messages: 15
Inscription: Lundi 15 Février 2010, 13:55
Statut actuel: Post-bac | Doctorat

Publicité

Retourner vers Logiciels mathématiques

 


  • Articles en relation
    Réponses
    Vus
    Dernier message

Qui est en ligne

Utilisateurs parcourant ce forum: Aucun utilisateur enregistré et 2 invités