# pohyb vysadkare s opozdenym otevrenim padaku
clear;
clc;
# vymazeme nastaveni vykreslovani
__gnuplot_raw__ "reset;";
# nastavuje opet vykresleni plnou carou
__gnuplot_set__ data style lines;
# obnoveni rozsahu os
__gnuplot_set__ xrange [:]
__gnuplot_set__ yrange [:]
###############################################
printf("Pohyb vysadkare s opozdenym otevrenim padaku\n");
% -----------------------------------------------
% definice poèáteèních podmínek
t0=0;
x0=0;
v0=0;
% -----------------------------------------------
% parametry prostredi
% -----------------------------------------------
g=9.81; 	% tíhové zrychlení
m=95;  	        % hmotnost parasutisty
rho=1.3;       	% hustota vzduchu 
C1=1.1;      	% koeficient odporu bez padaku
C2=1.33;       	% koeficient odporu s padakem
S1=1.1;         % prurez bez padaku
S2=45;          % prurez s padakem 
% -----------------------------------------------
% krok integrace a mez
dt = 0.01;
mez = 19;
% -----------------------------------------------
% vlastní numerická integrace - klasická Eulerova metoda
% -----------------------------------------------
cas_spusteni = clock();
t(1)=t0;
x(1)=x0;
v(1)=v0;
i=1;
% -----------------------------------------------
 while (t(i)<=mez)
   i=i+1;
   if (t(i-1)>10)
   a=g-C2*S2*rho*v(i-1)^2/2/m;
   else
   a=g-C1*S1*rho*v(i-1)^2/2/m;  
   endif  
   v(i)=v(i-1)+a*dt;
   x(i)=x(i-1)+v(i-1)*dt;
   t(i)=t(i-1)+dt;
 endwhile  
% -----------------------------------------------
celkovy_cas = etime(clock(), cas_spusteni);
disp(['Program bezel ',num2str(celkovy_cas),' sekund.']);
% -----------------------------------------------
% vykreslení zavislosti
clearplot
__gnuplot_set__ nokey
__gnuplot_set__ grid
__gnuplot_set__ tics out
__gnuplot_set__ mxtics
__gnuplot_set__ mytics
__gnuplot_set__ xlabel "t(s)"
__gnuplot_set__ ylabel "x(m)"
__gnuplot_set__ title "Parasutista s opozdenym otevrenim padaku"
% -----------------------------------------------
plot(t,x);
printf("Vykresleni zavislosti rychlosti na case\n");
printf("Stisknete libovolnou klavesu\n");
pause
__gnuplot_set__ xlabel "t(s)"
__gnuplot_set__ ylabel "v(m/s)"
plot(t,v); 
