Re: [Maxima-lang-fr] problème avec charpoly
|
From: denis.thenail <denis.thenail@fr...> - 2009-12-18 03:40:30
|
Bonjour,
quel type de machine et quel OS travaillez-vous ?
Pour ce qui me concerne, c'est UBUNTU 9.10 sur un PC avec pour proc un
athlon 64 X2 4400+.
Merci d'avance
Stavros Macrakis a écrit :
> Salut,
>
> Je n'arrive toujours pas a reproduire le probleme, neanmoins je
> remarque (ce que j'aurais du remarquer avant...) que l'erreur a ete
> signalee par dimnary, donc c'est un probleme de formattage et non pas
> de calcul. Sans pouvoir resoudre ce probleme, je vous conseillerais
> simplement de ne pas imprimer l'expression %o6 (utiliser "$" au lieu
> de ";" a la fin de la commande %i6) ou bien de l'imprimer en mode
> display2d:false.
>
> -s
>
> 2009/12/16 denis.thenail <denis.thenail@...
> <mailto:denis.thenail@...>>
>
> Bonjour,
>
> la version de MAXIMA est 5.17.1
>
> j'ai construit la matrice en plusieurs étapes, c'est à dire
> plusieurs scripts
>
> voici celui qui à partir de la matrice tridiagonale construit la
> matrice à inverser :
> Ka5:matrix([X,-X,0,0,0,0]);
> Ka5:addrow(Ka5,[-X,2*X,-X,0,0,0]);
> Ka5:addrow(Ka5,[0,-X,2*X,-X,0,0]);
> Ka5:addrow(Ka5,[0,0,-X,2*X,-X,0]);
> Ka5:addrow(Ka5,[0,0,0,-X,2*X,-X]);
> Ka5:addrow(Ka5,[0,0,0,0,-X,X]);
> Ma5:matrix([2*Y,Y,0,0,0,0]);
> Ma5:addrow(Ma5,[Y,4*Y,Y,0,0,0]);
> Ma5:addrow(Ma5,[0,Y,4*Y,Y,0,0]);
> Ma5:addrow(Ma5,[0,0,Y,4*Y,Y, 0]);
> Ma5:addrow(Ma5,[0,0,0,Y,4*Y,Y]);
> Ma5:addrow(Ma5,[0,0,0,0,Y,2*Y]);
> Ia5:invert(Ma5);
> Aa5:Ia5 . Ka5;
>
> où Ka5 est la matrice de raideur et Ma5 une matrice de masse.
> A part ça, oui, je sais que je peux écrire la marice en une seule
> fois.
>
> voici le script complet
> Ka5:matrix([X,-X,0,0,0,0],[-X,2*X,-X,0,0,0],[0,-X,2*X,-X,0,0],[0,0,-X,2*X,-X,0],
> [0,0,0,-X,2*X,-X],[0,0,0,0,-X,X]);
> Ma5:matrix([2*Y,Y,0,0,0,0],[Y,4*Y,Y,0,0,0],[0,Y,4*Y,Y,0,0],[0,0,Y,4*Y,Y,
> 0],
> [0,0,0,Y,4*Y,Y],[0,0,0,0,Y,2*Y]);
> Ia5:invert(Ma5);
> Aa5:Ia5 . Ka5;
> charpoly(Aa5,l);
> factor(%);
> expand(%);
>
> et la réponse de MAXIMA :
> (%i1) batch("test.mac");
>
> batching #p/home/denis/MAXIMA/test.mac
> (%i2) Ka5 : matrix([X, - X, 0, 0, 0, 0], [- X, 2 X, - X, 0, 0, 0],
> [0, - X, 2 X, - X, 0, 0], [0, 0, - X, 2 X, - X, 0], [0, 0, 0, - X,
> 2 X, - X],
> [0, 0, 0, 0, - X, X])
> [ X - X 0 0 0 0 ]
> [ ]
> [ - X 2 X - X 0 0 0 ]
> [ ]
> [ 0 - X 2 X - X 0 0 ]
> (%o2) [ ]
> [ 0 0 - X 2 X - X 0 ]
> [ ]
> [ 0 0 0 - X 2 X - X ]
> [ ]
> [ 0 0 0 0 - X X ]
> (%i3) Ma5 : matrix([2 Y, Y, 0, 0, 0, 0], [Y, 4 Y, Y, 0, 0, 0],
> [0, Y, 4 Y, Y, 0, 0], [0, 0, Y, 4 Y, Y, 0], [0, 0, 0, Y, 4 Y, Y],
> [0, 0, 0, 0, Y, 2 Y])
> [ 2 Y Y 0 0 0 0 ]
> [ ]
> [ Y 4 Y Y 0 0 0 ]
> [ ]
> [ 0 Y 4 Y Y 0 0 ]
> (%o3) [ ]
> [ 0 0 Y 4 Y Y 0 ]
> [ ]
> [ 0 0 0 Y 4 Y Y ]
> [ ]
> [ 0 0 0 0 Y 2 Y ]
> (%i4) Ia5 : invert(Ma5)
> [ 362 97 26 7 2 1 ]
> [ ----- - ----- ----- - ----- ----- - ----- ]
> [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ]
> [ ]
> [ 97 194 52 14 4 2 ]
> [ - ----- ----- - ----- ----- - ----- ----- ]
> [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ]
> [ ]
> [ 26 52 182 49 14 7 ]
> [ ----- - ----- ----- - ----- ----- - ----- ]
> [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ]
> (%o4) [ ]
> [ 7 14 49 182 52 26 ]
> [ - ----- ----- - ----- ----- - ----- ----- ]
> [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ]
> [ ]
> [ 2 4 14 52 194 97 ]
> [ ----- - ----- ----- - ----- ----- - ----- ]
> [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ]
> [ ]
> [ 1 2 7 26 97 362 ]
> [ - ----- ----- - ----- ----- - ----- ----- ]
> [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ]
> (%i5) Aa5 : Ia5 . Ka5
> [ 153 X 194 X 52 X 14 X 4 X X ]
> [ ----- - ----- ----- - ----- ----- - ----- ]
> [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ]
> [ ]
> [ 97 X 179 X 104 X 28 X 8 X 2 X ]
> [ - ----- ----- - ----- ----- - ----- ----- ]
> [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ]
> [ ]
> [ 26 X 104 X 155 X 98 X 28 X 7 X ]
> [ ----- - ----- ----- - ----- ----- - ----- ]
> [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ]
> (%o5) [ ]
> [ 7 X 28 X 98 X 155 X 104 X 26 X ]
> [ - ----- ----- - ----- ----- - ----- ----- ]
> [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ]
> [ ]
> [ 2 X 8 X 28 X 104 X 179 X 97 X ]
> [ ----- - ----- ----- - ----- ----- - ----- ]
> [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ]
> [ ]
> [ X 4 X 14 X 52 X 194 X 153 X ]
> [ - ----- ----- - ----- ----- - ----- ----- ]
> [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ]
> (%i6) charpoly(Aa5, l)
> 153 X 179 X 155 X
> (%o6) (----- - l) ((----- - l) ((----- - l)
> 209 Y 209 Y 209 Y
> 2
> 155 X 153 X 179 X 18818 X
> ((----- - l) ((----- - l) (----- - l) - --------)
> 209 Y 209 Y 209 Y 2
> 43681 Y
> 179 X 153 X
> 2 52 X (----- - l) 2 104 X
> (----- - l)
> 20176 X 209 Y 5044 X 209 Y
> 26 X (-------- - ----------------) 104 X (-------- -
> -----------------)
> 2 209 Y 2 209 Y
> 43681 Y 43681 Y
> + ---------------------------------- +
> ------------------------------------)
> 209 Y 209 Y
> 179 X
> 14 X (----- - l) 2
> 155 X 209 Y 5432 X
> + (7 X (- (----- - l) (---------------- - --------)
> 209 Y 209 Y 2
> 43681 Y
> 179 X
> 2 52 X (----- - l)
> 20176 X 209 Y
> 98 X (-------- - ----------------)
> 2 209 Y
> 43681 Y
> - ----------------------------------))/(209 Y)
> 209 Y
> 153 X
> 28 X (----- - l) 2
> 155 X 209 Y 1358 X
> + (28 X (- (----- - l) (---------------- - --------)
> 209 Y 209 Y 2
> 43681 Y
> 153 X
> 2 104 X (----- - l)
> 5044 X 209 Y
> 98 X (-------- - -----------------)
> 2 209 Y
> 43681 Y
> - -----------------------------------))/(209 Y)
> 209 Y
> 2
> 153 X 179 X 18818 X
> 98 X ((----- - l) (----- - l) - --------)
> 209 Y 209 Y 2
> 43681 Y
> + (98 X (- -----------------------------------------
> 209 Y
> 179 X 153 X
> 14 X (----- - l) 2 28 X (----- - l)
> 2
> 209 Y 5432 X 209 Y 1358 X
> 26 X (---------------- - --------) 104 X (---------------- -
> --------)
> 209 Y 2 209 Y
> 2
> 43681 Y
> 43681 Y
> + ---------------------------------- +
> -----------------------------------))
> 209 Y 209 Y
> 179 X
> 2 52 X (----- - l)
> 20176 X 209 Y
> 28 X (-------- - ----------------)
> 2 209 Y
> 155 X 43681 Y
> /(209 Y)) + (2 X (- (----- - l) (----------------------------------
> 209 Y 209 Y
> 179 X
> 2 4 X (----- - l)
> 155 X 1552 X 209 Y
> - (----- - l) (-------- - ---------------))
> 209 Y 2 209 Y
> 43681 Y
> 179 X
> 14 X (----- - l) 2
> 155 X 209 Y 5432 X
> - (104 X (- (----- - l) (---------------- - --------)
> 209 Y 209 Y 2
> 43681 Y
> 179 X
> 2 52 X (----- - l)
> 20176 X 209 Y
> 98 X (-------- - ----------------)
> 2 209 Y
> 43681 Y
> - ----------------------------------))/(209 Y)
> 209 Y
> 179 X
> 14 X (----- - l) 2
> 209 Y 5432 X
> 28 X (---------------- - --------)
> 209 Y 2
> 43681 Y
> - (98 X (----------------------------------
> 209 Y
>
> Universal error handler called recursively (:ERROR NIL
> MAXIMA::DIMNARY "" "Couldn't protect")
> Universal error handler called recursively (:ERROR NIL
> MAXIMA::DIMNARY "" "Couldn't protect")
> Maxima encountered a Lisp error:
>
> Error in MAXIMA::DIMNARY [or a callee]: Caught fatal error [memory
> may be damaged]
>
> Automatically continuing.
> To reenable the Lisp debugger set *debugger-hook* to nil.
> (%i7) factor(%)
> Erreur de segmentation
>
> merci pour votre attention
>
>
> Stavros Macrakis a écrit :
>
> Je n'arrive pas a reproduire votre probleme. Pourriez-vous
> nous envoyer une copie complete de votre interaction avec
> Maxima, y compris la version?
>
> Par ailleurs, savez-vous que vous pouvez ecrire la matrice en
> une commande, a savoir:
>
> Ca5:matrix([153*Z/209,-194*Z/209,52*Z/209,-14*Z/209,4*Z/209,-X/209],
>
> [-97*Z/209,179*Z/209,-104*Z/209,28*Z/209,-8*Z/209,2*X/209],
> ... )
>
> Merci,
>
> -s
>
> 2009/12/15 denis.thenail <denis.thenail@...
> <mailto:denis.thenail@...> <mailto:denis.thenail@...
> <mailto:denis.thenail@...>>>
>
>
> Bonjour,
>
> je cherche le polynôme caractéristique d'une matrice (6,6).
> Cette
> matrice est pleine mais résulte d'un problème aux valeurs
> propres
> généralisé tridiagonal du type (K-omega^2M)(Phi)= (0)
> M définie postive, K semi-définie positive
> ramené sous la forme standard
> (M^{-1}K-omega^2I)(Phi)=(0)
> J'exécute avec Maxima :
>
> Ca5:matrix([153*Z/209,-194*Z/209,52*Z/209,-14*Z/209,4*Z/209,-X/209]);
>
> Ca5:addrow(Ca5,[-97*Z/209,179*Z/209,-104*Z/209,28*Z/209,-8*Z/209,2*X/209]);
>
> Ca5:addrow(Ca5,[26*Z/209,-104*Z/209,155*Z/209,-98*Z/209,28*Z/209,-7*Z/209]);
>
> Ca5:addrow(Ca5,[-7*Z/209,28*Z/209,-98*Z/209,155*Z/209,-104*Z/209,26*Z/209]);
>
> Ca5:addrow(Ca5,[2*Z/209,-8*Z/209,28*Z/209,-104*Z/209,179*Z/209,-97*Z/209]);
>
> Ca5:addrow(Ca5,[-Z/209,4*Z/209,-14*Z/209,52*Z/209,-194*Z/209,153*Z/209]);
> charpoly(Ca5,l);
> factor(%);
> expand(%);
>
> la commande charpoly entraîne :
> Universal error handler called recursively (:ERROR NIL
> MAXIMA::DIMNARY
> "" "Couldn't protect")
> Universal error handler called recursively (:ERROR NIL
> MAXIMA::DIMNARY
> "" "Couldn't protect")
> Maxima encountered a Lisp error:
>
> Error in MAXIMA::DIMNARY [or a callee]: Caught fatal error
> [memory may
> be damaged]
>
> Automatically continuing.
> To reenable the Lisp debugger set *debugger-hook* to nil.
>
> Merci de m'aider à améliorer ma pratique de MAXIMA,
>
> cordialement
>
> Dude007
>
>
>
>
> ------------------------------------------------------------------------------
> This SF.Net email is sponsored by the Verizon Developer
> Community
> Take advantage of Verizon's best-in-class app development
> support
> A streamlined, 14 day to market process makes app distribution
> fast and easy
> Join now and get one step closer to millions of Verizon
> customers
> http://p.sf.net/sfu/verizon-dev2dev
> _______________________________________________
> Maxima-lang-fr mailing list
> Maxima-lang-fr@...
> <mailto:Maxima-lang-fr@...>
> <mailto:Maxima-lang-fr@...
> <mailto:Maxima-lang-fr@...>>
>
> https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr
>
>
>
>
|