%% Dane
Pgen = 50000;
P0 = 70;
t0 = 500;
Pu = 6;
Pk = 0.06;
etha_w = 0.82;
etha_m = 0.98;
etha_g = 0.97;
%% Entlpie
i0 = XSteam('h_pT',P0, t0);
s0 = XSteam('s_pT', P0, t0);
iut = XSteam('h_ps', Pu, s0);
iu = i0-etha_w*(i0-iut);
su = XSteam('s_ph', Pu, iu);
ikt = XSteam('h_ps', Pk, su);
ik = iu-etha_w*(iu-ikt);
%% Wyznaczenie mocy pracy jałowej turbiny N0
N0 = (Pgen/(etha_m*etha_g))-Pgen;
%% Wykres Schneidera
mu = 0;
mk = 0;
P = [-5000:5:50000];
for k=1:5
if k==1
mu(k) = 0;
mk(k) = 0;
m1 = ((P+N0)/(i0-ik))+((iu-ik)/(i0-ik))*mu(k);
m2 = ((P+N0)/(i0-iu))-((iu-ik)/(i0-iu))*mk(k);
plot(P,m1,P,m2,'-r','LineWidth',2);
hold on;
else
mu(k) = mu(k-1)+10;
mk(k) = mk(k-1)+10;
end
Px(k) = 0.542*mu(k)-2.637;
P1 = [Px(k):10:50000]
m1 = ((P1+N0)/(i0-ik))+((iu-ik)/(i0-ik))*mu(k);
m2 = ((P1+N0)/(i0-iu))-((iu-ik)/(i0-iu))*mk(k);
plot(P1,m1,P1,m2,'-r');
hold on;


end
xlim([-5000 50000]);
ylim([0 70]);
grid on;
























%% Dane
Pgen = 50000;
P0 = 70;
t0 = 500;
Pu = 6;
Pk = 0.06;
etha_w = 0.82;
etha_m = 0.98;
etha_g = 0.97;
%% Entlpie
i0 = XSteam('h_pT',P0, t0);
s0 = XSteam('s_pT', P0, t0);
iut = XSteam('h_ps', Pu, s0);
iu = i0-etha_w*(i0-iut);
su = XSteam('s_ph', Pu, iu);
ikt = XSteam('h_ps', Pk, su);
ik = iu-etha_w*(iu-ikt);
%% Wyznaczenie mocy pracy jałowej turbiny N0
N0 = (Pgen/(etha_m*etha_g))-Pgen;
%% Wykres Schneidera
mu = 0;
mk = 0;
P = [-5000:5:50000];
for k=1:5
if k==1
mu(k) = 0;
mk(k) = 0;
m1 = ((P+N0)/(i0-ik))+((iu-ik)/(i0-ik))*mu(k);
m2 = ((P+N0)/(i0-iu))-((iu-ik)/(i0-iu))*mk(k);
plot(P,m1,P,m2,'-r','LineWidth',2);
hold on;
else
mu(k) = mu(k-1)+10;
mk(k) = mk(k-1)+10;
end
m1 = ((P+N0)/(i0-ik))+((iu-ik)/(i0-ik))*mu(k);
m2 = ((P+N0)/(i0-iu))-((iu-ik)/(i0-iu))*mk(k);
plot(P,m1,P,m2,'-r');
hold on;


end
xlim([-5000 50000]);
ylim([0 70]);
grid on;