clc
clear all
syms r m k x1 x2
M=[20 0;0 10] %matriz de masa
K=[2 -1;-1 1]*1e5 %matriz de rigidez
[xx,omega]=eig(inv(M)*K); %valores propios
det=det(M*r^2 + K)
fr_natural=sqrt(omega)
r2=-omega
A=[M*r^2 + K];
fila1=A(1,:)*[x1 x2]'
syms r m k x1 x2 real
>> syms_real1
M =
20 0
0 10
K =
200000 -100000
-100000 100000
det =
200*r^4 + 4000000*r^2 + 10000000000
fr_natural =
130.6563 0
0 54.1196
r2 =
1.0e+04 *
-1.7071 0
0 -0.2929
fila1 =
x1*(20*r^2 + 200000) - 100000*x2