%programa preparacion datos toolbox ident respuesta Radiacion-Tsc
ws=input('Introduce workspace \n','s');
open(ws)
Tst=ans.Ts';
Irrt=ans.Irr';
VM1t=ans.VM1';
ident;
%Debemo importar variables, para ello presionamos en import, como variables
%de entrada introducimos [VM1t Irrt] y de salida Tst. presionamos OK
%A continuacion presionamos dentro del boton prepara modelo, remove means.
% Obtenemos el modelo parametrico, una vez obtenido presionamos dos veces
% sobre el y presionamos el boton present. Una vez hecho esto se cierre la
% ventana siguiente.
h=msgbox('Para continuar el programa, tras ejecutar ident, debe cerrar esta ventana');
waitfor(h);
sysident=ans;
sys=tf(sysident);
Ts_Irr=step(sys(:,2,:));
t=length(Ts_Irr);

for x=1:1:t;
    
    Ts(x)=Ts_Irr(x);
    Irr_s(x)=ones(1);
end

T_muestreo=1; %Tiempo de muestreo con el que se han sacado los datos con CubeDDE

V1=Irr_s(1);

x=1:1:t;

%for i=1:t
 %   if (VM1(i)==V1)
  %      nescalon=i;
  %end
    % VM1(i)=VM1(i)+30;
    %end

Tfinal=Ts(t);

Tinicial=Ts(1);

T5ciento=((Tfinal-Tinicial)*0.05)+Tinicial;

T63ciento=((Tfinal-Tinicial)*0.63)+Tinicial;


%plot(x,Ts(x),'',x,VM1(x),'')%,x,T5ciento,'',x,T63ciento);
%AXIS([0 1800 60 110]);legend('Tsc','VM1');

T5diferencia=abs(Ts(1)-T5ciento);
T5diferencia2=T5diferencia;
for i=2:1:t;
    T5diferencia=abs(Ts(i)-T5ciento);
    if (T5diferencia2>T5diferencia);
        T5diferencia2=T5diferencia;
        tiempo5=i*T_muestreo;
    end
end
    
T63diferencia=abs(Ts(1)-T63ciento);
T63diferencia2=T63diferencia;
for i=2:1:t;
    T63diferencia=abs(Ts(i)-T63ciento);
    if (T63diferencia2>T63diferencia)
        T63diferencia2=T63diferencia;
        tiempo63=i*T_muestreo;
    end
end

k=(Tfinal-Tinicial)/1;

tau=(tiempo63-tiempo5);

retraso=(tiempo5-(0*T_muestreo));

save 100_70simIrr_pru.mat;





