1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34
| plt.figure(1)
plt.subplot(211)
Qevap, =plt.plot(x_graph, Q[:,0], label='$Q_{evap}$')
Qabs, =plt.plot(x_graph, Q[:,1], label='$Q_{ads}$')
Qdes, =plt.plot(x_graph, Q[:,2], label='$Q_{des}$')
Qcond, =plt.plot(x_graph, Q[:,3], label='$Q_{cond}$')
plt.legend(handles=[Qevap, Qabs, Qdes, Qcond])
plt.ylabel(u'$Q (kW)$', fontsize=20, rotation=90)
plt.grid(True)
fig, ax1 = plt.subplots(212)
ax1.plot(x_graph, y1_graph, 'b-',x_graph, y2_graph, 'g-')
ax1.set_xlabel(u'$x_{NH3}$', fontsize=20, rotation=0)
ax1.set_ylabel(u'$COP$', fontsize=20, rotation=90)
for tl in ax1.get_yticklabels():
tl.set_color('b')
plt.grid(True)
ax2 = ax1.twinx()
ax2.plot(x_graph, eff, 'r-')
ax2.set_ylabel(u'$efficacite$', fontsize=20, rotation=90)
for tl in ax2.get_yticklabels():
tl.set_color('r')
#plt.title(':)')
plt.show() |
Partager