%% parametry loziska
n = 9;      %pocet kulicek
Bd = 6.35;  %prumer kulicek
Pd = 31.28;    %prumer roztecne kruznice
phi = 0; %stycny uhel

%% parametry motoru
RPM = 700;
fot = RPM/60; % mechanicka f
Rate = 2048; % poèet vzorkù
%% poruchove frekvence
ford = (n/2)*fot * (1-Bd/Pd*cosd((phi))); % frekvence defektu na vnejsi draze loziska

fird = (n/2)*fot * (1 + Bd/Pd*cosd(phi)); % frekvence defektu na vnitrni draze loziska

fbd = (n/2)*fot * (1 - (Bd/Pd)^2*(cosd(phi))^2);  % frekvence defektu kulicky

fc = (1/2)*fot*(1-Bd/Pd*cosd(phi));   % frekvence klece loziska
%%
q = plot(time, data);
xlabel('Èas [s]');
ylabel('Zrychlení [m/s^2]');
title('Výpadek fáze 50 Hz 705 RPM');
set(gca,'FontSize',24);
grid on;
hold on;


%% zpracovani dat
Fs = Rate;
len = size(data, 1);
f = (0:Fs/len:Fs-Fs/len);
Y = fft(data);


plot(f, abs(Y)/len*2);
xlim([0 fird*1.2]);
xlabel('Frekvence [Hz]');
ylabel('Amplituda');
title('FFT 30 s 50 Hz 705 ot/min');
grid on;
hold on;
set(gca,'FontSize',24);
plot([1 1] * fot, ylim, 'r');
plot([1 1] * fbd, ylim, 'r');
plot([1 1] * fc, ylim, 'r');
plot([1 1] * fird, ylim, 'r');
plot([1 1] * ford, ylim, 'r');


fbd = 50.13;
fot = 12.5;

lim = ylim;
text(fot, lim(2)*0.95, 'mechanická f', 'FontSize',20);
text(fbd, lim(2)*0.95, 'f kulièky', 'FontSize',20);
text(fc, lim(2)*0.95, 'f klece', 'FontSize',20);
text(fird, lim(2)*0.95, 'vnitøní dráha', 'FontSize',20);
text(ford, lim(2)*0.95, 'vnìjší dráha', 'FontSize',20);
set(gca,'FontSize',24);
