clc
close all

%% ZADANÉ HODNOTY
%tíhové zrychlení [m/s^2]
g = 9.81
%moment setrvačnosti [kg*m^2]
I = 119.294e-006;
%hmotnost [kg]
m = 77.99e-003;
%vlastní frekvence [Hz]
frek_vl = 16;
%odporová zátěž [Ohm]
%Rz = 250;
%vzdálenost vazba - cívka [m]
l = 0.055;
%faktor kvality [-]
%Q = 200;

%% BUZENÍ

%počáteční výchylka [°]
q_st = 0;
%počáteční výchylka [rad]
q_rad = (pi*q_st)/180

%rameno setrvačnosti
i=sqrt(I/m)

%amplituda budících vibrací [m/sec^2]
amp_m = 12.6*g;%0.3*g;
%amplituda budících vibrací [°/sec^2]
amp_st = amp_m/i^2
%amplituda budících vibrací [rad/sec^2]
amp_rad = (amp_st*pi)/180

%frekvence budících vibrací [Hz]
frek_hz = frek_vl;
%frekvence budících vibrací [rad/sec]
frek_rad = 2*pi*frek_hz

%% NELINEÁRNÍ TUHOST

%hraniční výchylka pro linearitu tuhosti [°]
q_lin_st = 7;
%hraniční výchylka pro linearitu tuhosti [rad]
q_lin_rad = (q_lin_st*pi)/180;
%vychylka pro 10*K [°]
q_10k_st = 10;
%vychylka pro 10*K [rad]
q_10k_rad = (q_10k_st*pi)/180;

%% CÍVKA - známé hodnoty

%počet závitů cívky [-]
Nc = 1000;
%učinná délka závitu [m]
lc = 0.021;
% % %Magnetická indukce [T]
% B = 0.317;
%měrný odpor Cu [Ohm*m^2]
ro = 1.69e-008;
%poloměr záhybu cívky [m]
r_civ = 5e-003;
%průměr drátu [m]
d_dr = 1e-004;

%% VÝPOČTY CHARAKTERISTIK

%zobecněná hmotnost
m_zob=I;
%vlastní uhlová frekvence [rad/s]
w_vl = frek_vl*2*pi;
%zobecněná tuhost
K_zob = m_zob*(w_vl)^2;
%tuhost [Nm/rad]
K = q_lin_rad*K_zob

%% VÝPOČTY PARAMETRŮ CÍVKY

%obvod jednoho závitu [m]
o_zav = 2*lc+2*pi*r_civ;
%délka cívky [m]
l_civ = o_zav*Nc;
%plocha drátu [m^2]
S_dr = pi*(d_dr/2)^2;
%činný odpor cívky [Ohm]
Rc = ro*l_civ/S_dr

%% VÝPOČTY TLUMENÍ a ÚTLUMU

%kritické tlumení
b_k = 2*m_zob*w_vl
% %elektromagnetické tlumení
% b_e = ((2*B*lc*Nc*l)^2)/(Rc+Rz)
% %Pomerne elektromagnetické tlumeni
% b_p_em = b_e/b_k
%Pomerne mechanicke tlumeni
%b_p_mech = 1/(2*Q)

%% Parameters for AC/DC rectifier
blok_nap_diod=0.25; % typicke priepustne napatie Schottkyho diody
kapacita_C = 10e-6;% 180e-9;%220e-6; % [F]

%% A-F CHARAKTERISTIKA

% if q_st == 0
% 
% %spodní hraniční hodnota frekvence v Hz
% f_s_Hz = 14
% %horní hraniční hodnota frekvence v Hz
% f_h_Hz = 20
% 
% f_s_rad = 2*pi*f_s_Hz;
% f_h_rad = 2*pi*f_h_Hz;
% 
% 
% i =0;
% for frek_rad = f_s_rad:0.1:f_h_rad;
%     tspan = [0 10];
%     sim('model1s',tspan);
%     i=i+1;
%     hodnoty = ScopeData2.signals.values(2500:3000);
%     u_max (i)= max(ScopeData2.signals.values);
%     frek_hz (i) = frek_rad/(2*pi);
% end
% 
% figure
% plot(frek_hz, u_max, 'k')
% xlabel ('Frekvence [Hz]')
% ylabel ('Maximální napětí [V]')
% grid on
% 
% else
% end
% 
% pfer = u_max

%gh = max(ScopeData3.signals.values)
% 
% figure
% plot(a,b,'k','LineWidth',1)
% hold on
% plot(t,x,'r','LineWidth',1)
% grid on
% xlabel('Čas [s]')
% ylabel('Napětí [V]')
% legend('Upravená cívka','Původní cívka')
% axis([0 3.5 -4 4])


