% Program pro vypocet diagramu zatizeni
% Matice casu a vykonu
D=[2,600
   4,900
   7,1500
   4,1300
   3,1400
   3,1300
   1,1100]

% Soucet hodin intervalu davajici Dobu provozu
Tp=sum(D(:,1));
fprintf('Doba provozu Tp = %.1f hod\n',Tp)

% Vyrobena energie
W=D(:,1).'*D(:,2);
fprintf('Vyrobena energie W = %.1f kWh\n',W)

% Doba vyuziti maxima
T=W/max(D(:,2));
fprintf('Doba vyuziti maxima T = %.4f hod\n',T)

% Suma Pi^2 * Ti
D2=D(:,2).^2;
W2=D(:,1).'*D2;
fprintf('Suma Pi^2 * Ti = %.1f kW^2*hod\n',W2)

% Doba plnych ztrat
Tz=W2/max(D(:,2))^2;
fprintf('Doba plnych ztrat Tz = %.4f hod\n',Tz)

% Ekonomicka narocnost
En=Tz/T;
fprintf('Ekonomicka narocnost = %.6f []\n',En)

% Stredni vykon
Ps=W/Tp;
fprintf('Stredni vykon Ps = %.2f kW\n',Ps)
