clear all; clc; close all;



%TUHOSTI PRUŽIN ki[N/m]:

k2=216000;
k1=1274000;


%VZDÁLENOSTI[m]:
u1=4.5;
u2=4.5;
l=16.5;
h=2.5;


%---------------------------------------------------
% DEFINICE KONSTANT HMOTNOSTÍ[kg] a MOMENTÙ SETRVAÈNOSTI[kg.m^2]
%Kolo
ms=12350;
%  VZOREC: 1/12*ms*(b^2+h^2)   kde b je délka skøínì a h výška skøínì
Jsy=1/12*ms*((l^2)+(h^2));
m1=1100;
%--------------------------------------------------

%#####  2. VÝPOÈET VLASTNÍCH FREKVENCÍ

M = [ ms                          0                           0                      0 ;
       0                         Jsy                          0                      0 ;
       0                          0                           m1                     0 ;
       0                          0                           0                      m1;
] 

   
K = [ 8*k2                          4*k2*(u1-u2)                             -4*k2                       -4*k2;
     4*k2*(u1-u2)                 4*k2*(u1^2+u2^2)                           4*k2*u2                       -4*k2*u1 ;
      -4*k2                           4*k2*u2                              2*k1                      0 ;
      -4*k2                          -4*k2*u1                                0                      2*k1;
]   


%VLASTNÍ ÈÍSLA:


[U,V]=eig(K,M)

om1=sqrt(V(1,1))
om2=sqrt(V(2,2))
om3=sqrt(V(3,3))
om4=sqrt(V(4,4))


FREKVENCE=(1/(2*pi))*sqrt(V)

FREKVENCE1=FREKVENCE(1,1)
FREKVENCE2=FREKVENCE(2,2)
FREKVENCE3=FREKVENCE(3,3)
FREKVENCE4=FREKVENCE(4,4)

U1norm=U(:,1)/U(1,1)
U2norm=U(:,2)/U(1,2)
U3norm=U(:,3)/U(1,3)
U4norm=U(:,4)/U(1,4)






