% Skript aproximuje hodnoty soucinitele aerodynamickeho odporu vzduchu Cd
% polynomem vhodneho radu tak, aby polynom dobre aproximoval zadane body.
% Hodnoty soucinitele aerodynamickeho odporu vzduchu je treba ziskat z
% programu tatra.exe jako vysledne hodnoty pro parametr 1 nastaveny na 0.7,
% 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4.

% K ziskani hodnot Cd pro jednotlive hodnoty prvniho parametru vyzijte
% soubor tatra_Data/Nazev.tatra, kde Nazev je nazev vozu ve hre, napr.
% Red.tatra. Do souboru vepiste prislusnou hodnotu na pozici prvniho
% parametru, ve hre kliknete na tlacitko import a spuste hru. Vyslednou
% hodnotu Cd si opiste a opakujte proces pro vsechny testovane hodnoty
% parametru.

% x-ove souradnice
X = 0.7:0.1:1.4;
% namerene hodnoty Cd
Cd = [0.3075977 0.3044318 0.3012659 0.2981 0.2999 0.3017 0.3035 0.3053];

for i = 1:7
    % vypocet regresniho polynomu
    P = polyfit(X,Cd,i);

    % vycisleni regresniho polynomu v pozadovanych bodech
    t = 0.65:0.01:1.45;
    y = polyval(P,t);
    
    figure
    hold on
    box on

    xlabel('x')
    ylabel('y')

    % vykresleni regresniho polynomu
    plot(t,y,'LineWidth',2)
    % vykresleni namerenych bodu
    plot(X,Cd,'g+','LineWidth',2.5);
end