function [] = urceni_stavu(tau,X,metoda)
% function [typ_poruchy, pomery, opravnenost] = urceni_stavu(tau,X,metoda)
% urceni_stavu(X(:,1),X(:,2),X(:,3),X(:,4),X(:,5),X(:,6),X(:,7),X(:,8),meto
% da)  volání   (t, X_CH4,X_C2H2,X_C2H4,X_H2,X_C2H6,X_CO,X_CO2,metoda)

%tau...čas měření, ve kterém chci provést analýzu
%typ poruchy...typ poruchy transformátoru
%pomery...pomery plynu pro mereni v case t (souřadnice bodu v obrazku)
%opravnenost...1 plyn prekrocil hodnotu z normy, 0 plyn neprekrocil hodnotu z normy  
%meze...meze z norem pro započetí analýzy pro plyny X_CH4,X_C2H2,X_C2H4,X_H2,X_C2H6,X_CO,X_CO2
clc
t=find(X(:,1)==tau);
meze=[30 2 60 50 20 400 3800];
plyny=[X(t,2) X(t,3) X(t,4) X(t,5) X(t,6) X(t,7) X(t,8)];
l=length(meze);
for i=1:l
    if plyny(i)>meze(i)
        opravnenost(i)=1;
    else opravnenost(i)=0;
    end
end
if metoda==1   
x=plyny(3)/plyny(5);
y=plyny(1)/plyny(4);
z=plyny(2)/plyny(3);
if x>=0 && x<0.2 && y>=0 && y<0.1 && z>=0 && z<0.01  
typ_poruchy='PD';
else typ_poruchy='ostatni';
end
if x>=0 && x<1 && y>=1 && y<inf && z>=0 && z<0.01  
typ_poruchy='T1';
else typ_poruchy='ostatni';
end
if x>=1 && x<4 && y>=1 && y<inf && z>=0 && z<0.1  
typ_poruchy='T2';
else typ_poruchy='ostatni';
end
if x>=4 && x<inf && y>=1 && y<inf && z>=0 && z<0.2  
typ_poruchy='T3';
else typ_poruchy='ostatni';
end
if x>=2 && x<inf && y>=0.1 && y<0.5 && z>=1 && z<2.5  
typ_poruchy='D1/D2';
else typ_poruchy='ostatni';
end
if x>=1 && x<inf && y>=0.1 && y<0.5 && z>=1 && z<inf  
typ_poruchy='D1';
else typ_poruchy='ostatni';
end
if x>=2 && x<inf && y>=0.1 && y<1 && z>=0.6 && z<2.5  
typ_poruchy='D2';
else typ_poruchy='ostatni';
end
end
if metoda==2
sumaplyny=plyny(1)+plyny(2)+plyny(3);
x=plyny(1)/sumaplyny;
y=plyny(3)/sumaplyny;
z=plyny(2)/sumaplyny;
if y>=0 && y<0.23 && z>=0.13 && z<1 
typ_poruchy='D1';
end
if y>=0.23 && y<0.4 && z>=0.13 && z<0.8 || y>=0.4 && y<0.71 && z>=0.29 && z<0.6
typ_poruchy='D2';
end
if y>=0.4 && y<0.85 && z>=0.15 && z<0.29 || y>=0.4 && y<0.5 && z>=0.13 && z<0.15 || y>=0 && y<0.5 && z>=0.02 && z<0.13
typ_poruchy='DT';
end
if y>=0.5 && y<1 && z>=0 && z<0.15 
typ_poruchy='T3';
end
if y>=0.2 && y<0.5 && z>=0 && z<0.04 
typ_poruchy='T2';
end
if y>=0.02 && y<0.2 && z>=0 && z<0.04 || y>=0 && y<0.02 && z>=0 && z<0.04&& x>=0.94 && x<0.98
typ_poruchy='T1';
end
if x>=0.98 && x<=1 
typ_poruchy='PD';
end
end
opravnenost
typ_poruchy
pomery=[x y z]
duval_obrazek(pomery);
end









