%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;
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
