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() |
