close all
clear all

%Základní promìnné 
a = 0.4;
r = 0.3;
m = 0.8;
l = 0.7;
t1 = 0.8;
t2 = 1.9;

% (Spuštìní simulace vytvoøené v prostøedí simulink (soubor simulace6_4.slx -
% slx je pøípona simulink, což je nástroj matlabu. Dostanete se tam, když
% zadáte do command windows pøíkaz simulink - tam mùžete napøíklad otevøít simulace 6_4 
% nebo vytvoøit svùj další)

sim('simulace6_4.slx'))

%Potøeboval jsem si nìjak vytvoøit graf pùlkružnice abych ho mohl vložit do
%animace

kruznicey = r*sin(0:0.01:pi()); %sinus s amplitudou r
kruznicex = r*cos(0:0.01:pi()) + m + r; %cosinus s amplitudou r a støedem posunutým o m + r

%Ze simulinku vyjdou promìnné xA, yA, xD, yD, vxD, vyD, vxA ve formátu
%timeseries, což si mùžete pøedstavit jako matici o dvou sloupcích. Jeden
%sloupec obsahuje èas a druhý hodnoty. Pøístup k èasùm nebo hodnotám
%získáte pøes promìnná.time nebo promìnná.data

cas = xA.time;


%Toto je animace, jednoduchý for cycle. Pro úkázku vždy vytáhne z timeseries
%jednu hodnotu - napø. Ax = xA.data(2) vytáhne druhou hodnotu ze sloupce
%data. 
for t = 1 : length(cas) %musím pøejít na poèty s celými èísly, abych mohl indexovat 
    %v datech (matice), tìžko bych tahal napøíklad xA.data(0.1) 
    Ax = xA.data(t);
    Ay = yA.data(t);
    Dx = xD.data(t);
    Dy = yD.data(t);
    
    %Nyní mám v pamìti x a y souøadnice bodù A,D pro urèité t
    
    %Pokud plotujete v [ ] tak výstupem je úseèka, všimnìte si, že v první
    %hranaté závorce jsou obì x souøadnice a v druhé y. Kruznicex a
    %kruznicey vykreslý do grafu pùlkruh podle tìch sin a cos nahoøe.
    %Nakonec se mají zvýraznit body A (souøadnice Ax,Ay) znaèkou + a stejnì
    %bod D.
    plot([Ax Dx], [Ay Dy], kruznicex, kruznicey, Ax, Ay, '-+', Dx, Dy, '-+')
    
    text(Ax-0.02,Ay+0.05,'A') %Pokud si chcete do grafu vložit písmeno nebo text... 
    text(Dx - 0.07,Dy +0.03,'D')%ty + a - jsou kvùli posunu písmena na èitelnìjší místo
    
    xlim([0 1.5]) %nastavení rozmezí os
    ylim([0 1.5])
   
    grid on %møížka - mám ji tam rád
    pbaspect([1 1 1]) %Mìøítko os, chci je 1:1:1 x:y:z .. i když tam z není, tak jistota je
    %jistota
    
    %Aby to po každém kroku udìlalo na chvilku pauzu
    pause(0.05);
end
