notebook.community

Edit and run



In [ ]:

    
import matplotlib.pyplot as plt
x = []
y =[]
for i in range(0, 101):
    x.append(i / 100)
    y_i = - 0.4191 * (i/ 100) ** 4 + 1.5087 * (i/ 100) ** 3 - 2.3858 * (i/ 100) ** 2 + 2.2951 * (i/ 100) + 0.0005
    y.append(y_i)
plt.plot([0, 1], [0, 1])
plt.plot(x, y)
F = 1
# x_f = 0.555
for i in range(1, 10, 3):
    x_f = i / 10
    y_f = - 0.4191 * x_f ** 4 + 1.5087 * x_f ** 3 - 2.3858 * x_f ** 2 + 2.2951 * x_f + 0.0005
    x_1d = 0.999
    x_1w = 1 - x_1d
    D = F * (x_f - x_1w)/(x_1d - x_1w)
    f = F / D
    Rmin = (x_1d - y_f) / (y_f - x_f)
    R = Rmin * 1.3
    x_work = []
    y_work = []
    for j in range(1, 1000):
        x_work.append(j / 1000)
        if j / 1000 <= x_f:
            y_work_down = (R + f) / (R + 1) * j / 1000 - (1 - f) / (R + 1) * x_1w
            y_work.append(y_work_down)
        elif j / 1000 >= x_f:
            y_work_up = R / (R + 1) * j / 1000 + x_1d / (R + 1)
            y_work.append(y_work_up)
    plt.plot(x_work, y_work)
plt.axis([0, 1, 0, 1])
plt.show()



In [ ]:

    
import matplotlib.pyplot as plt
x = []
y =[]
for i in range(0, 101):
    x.append(i / 100)
    y_i = - 0.4191 * (i/ 100) ** 4 + 1.5087 * (i/ 100) ** 3 - 2.3858 * (i/ 100) ** 2 + 2.2951 * (i/ 100) + 0.0005
    y.append(y_i)
plt.plot([0, 1], [0, 1])
plt.plot(x, y)



def x_eq(y):
    global mix
    mixture = mix
    if mixture == 'bt':
        yx = 0.8551 * y ** 4 - 0.8917 * y ** 3 + 0.6754 * y ** 2 + 0.3567 * y + 0.0022
    elif mixture == 'te':
        yx = 0.6986 * y ** 4 - 1.067 * y ** 3 + 1.028 * y ** 2 + 0.3335 * y + 0.0031
    elif mixture == 'ea':
        yx = - 2.815 * y ** 5 + 7.403 * y ** 4 - 5.5997 * y ** 3 + 1.9874 * y ** 2 + 0.0262 * y + 0.0024
    return yx


def y_eq(x):
    global mix
    mixture = mix
    if mixture == 'bt':
        xy = - 0.4191 * x ** 4 + 1.5087 * x ** 3 - 2.3858 * x ** 2 + 2.2951 * x + 0.0005
    elif mixture == 'te':
        xy = - 0.7937 * x ** 4 + 2.0119 * x ** 3 - 2.2674 * x ** 2 + 2.0401 * x + 0.0059
    elif mixture == 'ea':
        xy = 4.0475 * x ** 5 - 12.916 * x ** 4 + 16.378 * x ** 3 - 10.649 * x ** 2 + 4.1377 * x + 0.0052
    return xy


def ch_comma_to_dot(string):
    if string.find(',') > -1:
        number = float(string[:string.find(',')] + '.' + string[string.find(',') + 1:])
    else:
        number = float(string)
    return number


def ch_d_to_c(number):
    string = str(round(number, 4))
    answer = string[:string.find('.')] + ',' + string[string.find('.') + 1:]
    return answer


#mix = input('смесь: ')
mix = 'bt'
x_1_temp = input('количество первого компонента в питании(доли): ')
x_1 = ch_comma_to_dot(x_1_temp)
#x_1d_temp = input('количество первого компонента в дистиллате(доли): ')
#x_1d = ch_comma_to_dot(x_1d_temp)
x_1d = 0.999
#F_temp = input('питание кмоль/ч: ')
#F = ch_comma_to_dot(F_temp)
F = 0.4
x_1w = 1 - x_1d
D = F * (x_1 - x_1w)/(x_1d - x_1w)
f = F / D

Rmin = (x_1d - y_eq(x_1)) / (y_eq(x_1) - x_1)
R = Rmin * 1.3
y = x_1d
x = x_1d
x_work = []
y_work = []
N = 0
x_pr = 0
x_last = 0
while x > x_1w:
    x_work.append(x)
    x = x_eq(y)
    x_last = x
    if round(x_last, 5) == round(x_pr, 5):
        x_work.pop()
        break
    if x > x_1:
        y_work_up = R / (R + 1) * x + x_1d / (R + 1)
        y = y_work_up
        y_work.append(y)
    elif x < x_1:
        y_work_down = (R + f) / (R + 1) * x - (1 - f) / (R + 1) * x_1w
        y = y_work_down
        y_work.append(y)
    x_pr = x_last
    N +=1
print(len(x_work), len(y_work))
plt.plot(x_work, y_work)
print(*(round(x, 4) for x in x_work))
print(*y_work)
print('Rmin =', ch_d_to_c(Rmin))
print('число тарело:', N)


plt.show()



In [52]:

    
eq_file = open('C:/Users/kamorozov/PycharmProjects/DDparser/equation.txt', 'r')
equations_yx = {}
equations_xy = {}
for temp_line in eq_file:
    if temp_line[0] == '#':
        next
    else:
        line = temp_line.split('\t')
        line[2] = line[2][:-1]
        eq = line[2].split(',')
        eq = '.'.join(eq)
        for i in range(2, 10):
            if eq.find('x' + str(i)) > -1:
                eq = eq[:eq.find('x' + str(i))] + ' * x ** ' + str(i) + eq[eq.find('x' + str(i)) + 2:]
        if eq.find('x +') > -1:
            eq = eq[:eq.find('x +')] + ' * x' + eq[eq.find('x +') + 1:]
        if line[1] == 'yx':
            while eq.find('x') > -1:
                eq = eq[:eq.find('x')] + 'y' + eq[eq.find('x') + 1:]
        line[2] = eq
        if line[1] == 'yx':
            equations_yx[line[0]] = line[2]
        else:
            equations_xy.update({line[0]: line[2]})
        print(line)
print()
print(equations_yx)
print()
print(equations_xy)
print(equations_xy.get('bt'))
x = 0.1
y = eval(equations_xy.get('bt'))
print(y)
print()
for mix in equations_xy:
    print(mix, end=' ')









    



['bt', 'yx', '0.8551 * y ** 4 - 0.8917 * y ** 3 + 0.6754 * y ** 2 + 0.3567 * y + 0.0022']
['te', 'yx', '0.6986 * y ** 4 - 1.067 * y ** 3 + 1.028 * y ** 2 + 0.3335 * y + 0.0031']
['ea', 'yx', '-6.2194 * y ** 6 + 16.525 * y ** 5 - 15.399 * y ** 4 + 6.9697 * y ** 3 - 1.1967 * y ** 2 + 0.3218y - 0.0004']
['C56', 'yx', '4.3539 * y ** 6 - 10.769 * y ** 5 + 10.58 * y ** 4 - 4.694 * y ** 3 + 1.2171 * y ** 2 + 0.3094 * y + 0.0013']
['C67', 'yx', '1.6943 * y ** 6 - 3.9056 * y ** 5 + 3.7748 * y ** 4 - 1.5298 * y ** 3 + 0.5488 * y ** 2 + 0.4165 * y + 0.0004']
['C78', 'yx', '0.5598 * y ** 5 - 0.7365 * y ** 4 + 0.6087 * y ** 3 + 0.046 * y ** 2 + 0.5223y - 0.0003']
['bt', 'xy', '-0.4191 * x ** 4 + 1.5087 * x ** 3 - 2.3858 * x ** 2 + 2.2951 * x + 0.0005']
['te', 'xy', '-0.7937 * x ** 4 + 2.0119 * x ** 3 - 2.2674 * x ** 2 + 2.0401 * x + 0.0059']
['ea', 'xy', '-5.0929 * x ** 6 + 19.326 * x ** 5 - 30.215 * x ** 4 + 25.512 * x ** 3 - 12.88 * x ** 2 + 4.3483 * x + 0.0011']
['C56', 'xy', '-0.5305 * x ** 4 + 1.8182 * x ** 3 - 2.784 * x ** 2 + 2.4936 * x + 0.0015']
['C67', 'xy', '-0.3479 * x ** 4 + 1.2762 * x ** 3 - 2.1562 * x ** 2 + 2.2266 * x + 0.0008']
['C78', 'xy', '0.4197 * x ** 3 - 1.3783 * x ** 2 + 1.9592 * x + 0.0006']

{'bt': '0.8551 * y ** 4 - 0.8917 * y ** 3 + 0.6754 * y ** 2 + 0.3567 * y + 0.0022', 'te': '0.6986 * y ** 4 - 1.067 * y ** 3 + 1.028 * y ** 2 + 0.3335 * y + 0.0031', 'ea': '-6.2194 * y ** 6 + 16.525 * y ** 5 - 15.399 * y ** 4 + 6.9697 * y ** 3 - 1.1967 * y ** 2 + 0.3218y - 0.0004', 'C56': '4.3539 * y ** 6 - 10.769 * y ** 5 + 10.58 * y ** 4 - 4.694 * y ** 3 + 1.2171 * y ** 2 + 0.3094 * y + 0.0013', 'C67': '1.6943 * y ** 6 - 3.9056 * y ** 5 + 3.7748 * y ** 4 - 1.5298 * y ** 3 + 0.5488 * y ** 2 + 0.4165 * y + 0.0004', 'C78': '0.5598 * y ** 5 - 0.7365 * y ** 4 + 0.6087 * y ** 3 + 0.046 * y ** 2 + 0.5223y - 0.0003'}

{'bt': '-0.4191 * x ** 4 + 1.5087 * x ** 3 - 2.3858 * x ** 2 + 2.2951 * x + 0.0005', 'te': '-0.7937 * x ** 4 + 2.0119 * x ** 3 - 2.2674 * x ** 2 + 2.0401 * x + 0.0059', 'ea': '-5.0929 * x ** 6 + 19.326 * x ** 5 - 30.215 * x ** 4 + 25.512 * x ** 3 - 12.88 * x ** 2 + 4.3483 * x + 0.0011', 'C56': '-0.5305 * x ** 4 + 1.8182 * x ** 3 - 2.784 * x ** 2 + 2.4936 * x + 0.0015', 'C67': '-0.3479 * x ** 4 + 1.2762 * x ** 3 - 2.1562 * x ** 2 + 2.2266 * x + 0.0008', 'C78': '0.4197 * x ** 3 - 1.3783 * x ** 2 + 1.9592 * x + 0.0006'}
-0.4191 * x ** 4 + 1.5087 * x ** 3 - 2.3858 * x ** 2 + 2.2951 * x + 0.0005
0.20761879

bt te ea C56 C67 C78



In [42]:

    
eq = '0,8551x4 - 0,8917x3 + 0,6754x2 + 0,3567x + 0,0022'
if eq.find('x' + str(3)) > -1:
    print(eq.find('x' + str(3)))
eq.count('x')



In [ ]: