% Aufgabe 23-3, Aufgabenstellung b:

A = [5.6308  -4.6386   1.2155        0       0        0       0
   -4.6386  7.3013 -5.1020  1.3355       0        0       0
     1.2155  -5.1020   8.0215  -5.5991   1.4641        0       0
         0  1.3355 -5.5991  8.7937 -6.1318  1.6018       0
         0        0   1.4641  -6.1318   9.6203  -6.7016   1.7490
         0        0       0  1.6018 -6.7016  10.5040 -7.3103
         0        0       0        0   1.7490  -7.3103   9.3735] ;
 
B = [2   -1     0     0     0     0     0   ;
   -1     2   -1     0     0     0     0   ;
     0   -1     2   -1     0     0     0   ;
     0     0   -1     2   -1     0     0   ;
     0     0     0   -1     2   -1     0   ;
     0     0     0     0   -1     2   -1   ;
     0     0     0     0     0   -1     2   ] ;

[V D] = eig (A , B) ;
for i=1:7
   kappa(i) = D(i,i) ;       % Eigenwerte
end
Kappamin = min(kappa)        % Kleinster Eigenwert
piMalKappamin = pi*Kappamin

In das Command Window wird der kleinste Eigenwert (kappa) entsprechend der Aufgabenstellung ausgegeben, außerdem für einen Vergleich mit dem Ergebnis der Aufgabenstellung a) der mit π multiplizierte kleinste Eigenwert, aus dem sich der Näherungswert für die kritische Last nach

F_kr = 0,6976 E d ⁴ / l ²

berechnen lässt.

www.D@nkert.de

D

nkert.de