#!/usr/bin/env python #coding=utf8 import numpy as np# from scipy.optimize import minimize# import matplotlib.pyplot as plt# import time# Po=1# Pzad=64# g=1.4# m=10# R=8.314# a=(g-1)/g T=293# a=0.286 " :" def E1(x1): return m*R*T*a*((x1/Po)**2-1) def E2(x1,x2): return m*R*T*a*((x2/x1)**2-1) def E3(x2,x3): return m*R*T*a*((x3/x2)**2-1) def E4(x3,x4): return m*R*T*a*((x4/x3)**2-1) def E5(x4): return m*R*T*a*((Pzad/x4)**2-1) start = time.time() def E(x1,x2,x3,x4):# return E1(x1)+E2(x1,x2)+E3(x2,x3)+E4(x3,x4)+E5(x4) def fun2(x): return E(*x) x0 =[4,15,24,40] # res = minimize(fun2, x0)# stop = time.time() print (" :",round(stop-start,3)) x1=round(res['x'][0],3)# I x2=round(res['x'][1],3)# II x3=round(res['x'][2],3)# III x4=round(res['x'][3],3)# IV x5=Pzad " :" x=[] x.append(x1) x.append(x2) x.append(x3) x.append(x4) x.append(x5) z1=E1(x1)# I z2=E2(x1,x2)+z1# I+II z3=E3(x2,x3)+z2# I+II +III z4=E4(x3,x4)+z3# I+II +III +IV z5=E5(x4)+z4# I+II +III +IV+V " :" y=[] y.append(z1) y.append(z2) y.append(z3) y.append(z4) y.append(z5) plt.title(' - %s %s \n\ %s %s '%(int(z5),len(y),Po,Pzad)) plt.ylabel(' ') plt.xlabel(' ') plt.plot(x[0],y[0],'o', label='I , - %s'%x1) plt.plot(x[1],y[1],'o', label='II , - %s'%x2) plt.plot(x[2],y[2],'o', label='III , - %s'%x3) plt.plot(x[3],y[3],'o', label='IV , - %s'%x4) plt.plot(x[4],y[4],'o', label='V , - %s'%x5) plt.plot(x,y,'r') plt.legend(loc='best') plt.grid(True) plt.show()
#!/usr/bin/env python #coding=utf8 import numpy as np# from scipy.optimize import minimize# import matplotlib.pyplot as plt# import time# Po=1# Pzad=64# g=1.4# m=10# R=8.314# a=(g-1)/g T=293# a=0.286 " :" def E1(x1): return m*R*T*a*((x1/Po)**2-1) def E2(x1,x2): return m*R*T*a*((x2/x1)**2-1) def E3(x2,x3): return m*R*T*a*((x3/x2)**2-1) def E4(x3,x4): return m*R*T*a*((x4/x3)**2-1) def E5(x4,x5): return m*R*T*a*((x5/x4)**2-1) def E6(x5): return m*R*T*a*((Pzad/x5)**2-1) start = time.time() def E(x1,x2,x3,x4,x5): return E1(x1)+E2(x1,x2)+E3(x2,x3)+E4(x3,x4)+E5(x4,x5)+E6(x5) def fun2(x): return E(*x) x0 =[4,15,20,30,40]# res = minimize(fun2, x0)# stop = time.time() print (" :",round(stop-start,3)) x1=round(res['x'][0],3)# I x2=round(res['x'][1],3)# II x3=round(res['x'][2],3)# III x4=round(res['x'][3],3)# IV x5=round(res['x'][4],3)# V x6=Pzad# VI " :" x=[] x.append(x1) x.append(x2) x.append(x3) x.append(x4) x.append(x5) x.append(x6) z1=E1(x1)# I z2=E2(x1,x2)+z1# I+II z3=E3(x2,x3)+z2# I+II+III z4=E4(x3,x4)+z3# I+II+III+IV z5=E5(x4,x5)+z4# I+II+III+IV+V z6=E6(x5)+z5# I+II+III+IV+V +VI " :" y=[] y.append(z1) y.append(z2) y.append(z3) y.append(z4) y.append(z5) y.append(z6) plt.title(' - %s %s \n\ %s %s '%(int(z6),len(y),Po,Pzad)) plt.ylabel(' ') plt.xlabel(' ') plt.plot(x[0],y[0],'o', label='I , - %s'%x1) plt.plot(x[1],y[1],'o', label='II , - %s'%x2) plt.plot(x[2],y[2],'o', label='III , - %s'%x3) plt.plot(x[3],y[3],'o', label='IV , - %s'%x4) plt.plot(x[4],y[4],'o', label='V , - %s'%x5) plt.plot(x[5],y[5],'o', label='VI , - %s'%x6) plt.plot(x,y,'r') plt.legend(loc='best') plt.grid(True) plt.show()
Source: https://habr.com/ru/post/347624/
All Articles