DMUG-Archiv 2009

Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

Re: (Numerisches?) Problem beim Bestimmen der Eigenvektoren

Hoi Andreas,

Es scheint tatsächlich ein Problem der Rechengenauigkeit zu sein. Durch
Verändern der Präzision erhält man das richtige Ergebnis (siehe Anhang).

In[20]:= Clear[w]
w = Eigensystem[SetPrecision[M, 4]][[2]]
Out[21]= {{0.05714, 0.05650, 0.05921, -0.9950},
        {0.5884, 0.5838, 0.5505, 0.09975},
        {0.7071, -0.7071, -0.00001116, -2.544*10^-10},
        {0.3917, 0.3917, -0.8326, -0.004860}}

die Angabe der Matrixelemente als Festkommazahlen mit vier Nachkommastellen reicht schlicht und einfach, mit anderen Worten

In[44]:= Clear[c1, c2, c3, c4, e1, e2, e3]
c1 = 34.8669`4;
c2 = 36.6194`4;
c3 = -3.7212`4;
c4 = -4.5547`4;
e1 = 1958.56`4;
e2 = 1950.11`4;
e3 = 2096.11`4;

und es funktioniert.

Gruss
Udo.


Verweise:
Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

DMUG DMUG-Archiv, http://www.mathematica.ch/archiv.html