Energia_Sem.m

syms eta_ref beta T_cel T_re FR_Tau_a FR_UL T_amb Irrad T_fria T_gelada COP_chiller
syms eta_ther eta_el eta_1 eta_2 eta_frio eta_2chiller COP_e

eqn(1) = eta_el == eta_ref*(1-beta*(T_cel-25)); % T_ref=25;
eqn(2) = eta_ther==FR_Tau_a*(1-eta_el)-FR_UL*(T_cel-T_amb)/Irrad;
eqn(3) = eta_1==eta_el+eta_ther;
eqn(4) = eta_2==(1-(T_amb+273)/(T_cel+273))*eta_ther+eta_el;
eqn(5) = eta_frio==eta_ther*COP_chiller;
eqn(6) = eta_2chiller==((T_amb+273)/(T_gelada+273)-1)*eta_frio+eta_el;
% eqn(7) = eta_2==eta_ref;
eqn(7) = eta_2chiller==eta_ref;
eqn(8) = beta == 0.005; %0.002; % 
eqn(9) = T_fria==30;
eqn(10) = T_gelada==5;
eqn(11) = COP_chiller==0.55;
eqn(12) = T_amb==28;
eqn(13) = Irrad==500;
eqn(14) = FR_Tau_a==0.574/(1-0.123);
eqn(15) = FR_UL==4.85;
eqn(16) = eta_ref==0.2;
eqn(17) = COP_e == 3;
clear y x

for k = 1:5
 bet(k) = 1000+(k-1)/4*4000; % 200+(k-1)/4*800; %
 eqn(13) = Irrad==bet(k);
 for j=1:2
  eta_r(j) = 0.2; % 0.05+(j-1)/1*0.05; %
  eqn(16) = eta_ref==eta_r(j);
  for i=1:10
    x(i) = 30+(i-1)/9*90;
    if x(i) < 50
        eqn(11) = COP_chiller==0;
    else
        eqn(11) = COP_chiller==0.55;
    end
    eqn(7) = T_cel==x(i);
    S=vpasolve(eqn);
    if max(S.eta_2)>0 & isreal(max(S.eta_2))
        y(i,2*(k-1)+j)=max(S.eta_el)+max(S.eta_frio); % max(S.eta_1); %
    else
        y(i,2*(k-1)+j) = 0;
    end
  end
 end
end
y = double(y);
xx = 30:5:110;
yy = zeros(size(xx,2),size(y,2));
for i = 1:size(y,2)
     y1 = y(:,i);
     pp = csaps(x,y1);
     yy(:,i) = fnval(pp,xx)';
end
%Energia_Sem_graf(x,y*100)
Energia_Com_graf(x,y*100,'cool')