Matlab (calculate natural frequencies using eig function)

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 3K views
Junichirol
Messages
3
Reaction score
1

Homework Statement


Determine the corresponding natural frequencies
M = [1 0 0 0;0 1 0 0;0 0 1 0;0 0 0 1];
C = [2 0 0 0;0 0 0 0;0 0 0 0; 0 0 0 0];
K = [10 -5 0 0;-5 10 -5 0;0 -5 10 -5;0 0 -5 10];

Homework Equations


A = [zeros (4) eye(4); -M/K -M/C]
[V,D] = eig(A)

The Attempt at a Solution


just need help on A = [zeros (4) eye(4); -M/K -M/C].
-M/C =
NaN NaN NaN NaN
NaN NaN NaN NaN
NaN NaN NaN NaN
NaN NaN NaN -Inf

I also try with -inv(M*C)
-Inf -Inf -Inf -Inf
-Inf -Inf -Inf -Inf
-Inf -Inf -Inf -Inf
-Inf -Inf -Inf -Inf

This matrix caused the [V,D] = eig(A) cannot work.
Any step that i missed up?
 
Physics news on Phys.org
Junichirol said:

Homework Statement


Determine the corresponding natural frequencies
M = [1 0 0 0;0 1 0 0;0 0 1 0;0 0 0 1];
C = [2 0 0 0;0 0 0 0;0 0 0 0; 0 0 0 0];
K = [10 -5 0 0;-5 10 -5 0;0 -5 10 -5;0 0 -5 10];
Is that really the complete problem your were given?
 
DrClaude said:
Is that really the complete problem your were given?
For a system with the mass, damping, stiffness matrices given as below:
M = [1 0 0 0;0 1 0 0;0 0 1 0;0 0 0 1];
C = [2 0 0 0;0 0 0 0;0 0 0 0; 0 0 0 0];
K = [10 -5 0 0;-5 10 -5 0;0 -5 10 -5;0 0 -5 10];
Determine the corresponding natural frequencies by using Matlab eig function.