Buongiorno e grazie per l'ospitalità
il mio problema è che :
plottare eta solo se eta>0, non mi interessa
calcolare p02 con
- p2 =((1+((eta.*deltaideal)/(T01*cp))).^(k/(k-1))); se m>0
- p02=m^2 se m<0
il codice è questo:
Nmin = 20000; %Desired setpoint for engine revolution(rpm)
Nmax = 50000; %Maximum setpoint for engine revolution(rpm)
Nincr=15000; %Incremento numero di giri
N=Nmin;
while N<Nmax+1
U1 = (D1*pi.*N)/60;
U2 = (D2*pi.*N)/60;
alfa2b = atan((D1*tan(beta1b))/(sigma*D2));
deltaideal = sigma*(U2^2);
deltahii = 0.5.*((U1-((cot(beta1b).*m)/(ro1*Ai))).^2);
deltahid = 0.5.*(((sigma*D2*U1)/D1)-((m.*cot(alfa2b))/(ro1*Ad))).^2;
deltahfi = kf.*(m.^2);
deltahfd = kfd.*(m.^2);
deltan = deltanbf+deltanv+deltanc+deltand;
deltaloss = deltahii+deltahid+deltahfi+deltahfd;
eta = ((deltaideal)./(deltaideal+deltaloss))-deltan;
p2 =((1+((eta.*deltaideal)/(T01*cp))).^(k/(k-1)));
figure(1);
title('Rendimento')
xlabel('Mass Flow [kg/s]')
ylabel('eta')
grid on
plot(m,eta);
hold on
figure(2);
title('Pressione')
xlabel('Mass Flow [kg/s]')
ylabel('p02')
grid on
plot(m,p2);
hold on
N=N+Nincr
end