- MATLAB does numerical computation, so a round-off error compared to analytical computation is expected.
- The result of EIG doesn't give the eigenvalues in any particular order, so you may need to call sort on the eigenvalues and then reorder the eigenvectors accordingly.
eigen value what i am getting from matlab software that is not exact when i am compairing with hand calculated value and the mode shape i am getting is reverse
9 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
format long g
m1=50000;
m2=50000;
m3=50000;
mass=diag([m1 m2 m3]);
fprintf('The mass matrix \n');
disp(mass)
k1=50*10^6;
k2=50*10^6;
k3=50*10^6;
stiffness=[k1+k2 -k2 0;-k2 k2+k3 -k3;0 -k3 k3];
fprintf('The stiffness matrix \n');
disp(stiffness)
%Natural frequency(eigen value) and Mode shape(eigen vector) calculation FOR ANALYTICAL VALUE
[v,d]=eig(stiffness,mass);
fprintf('lambda value \n');
disp(d);
w=[sqrt(d)];
w1=w(1,1);
w2=w(2,2);
w3=w(3,3);
w=[w1,w2,w3];
fprintf('The natural frequencies of this structure are as follows(HZ) \n');
disp(w)
fprintf('The mode shape of the structure \n');
disp(v)
%Normalization of mode shape vectors
% for i=1:3
% v(:,i)=v(:,i)/v(3,i);
% end
% Individusal mode shapes vector
fprintf('individusal mode shapes \n');
v1=v(:,3)
v2=v(:,2)
v3=v(:,1)
0 Commenti
Risposte (1)
Christine Tobler
il 22 Apr 2022
The results look correct to me, can you say what you are expecting instead? Two points that are maybe relevant:
0 Commenti
Vedere anche
Categorie
Scopri di più su Linear Algebra in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!