|
|

楼主 |
发表于 2017-12-16 17:04:00
|
显示全部楼层
是表达式哈哈 F,X0,W1,W2,W3……Wi 都是定值
以下是我写的代码,其中 i 都用n代替:
- import numpy as np
- import matplotlib.pyplot as plt
- import matplotlib.patches as mpatches
- n = int(input("请输入效数:"))
- x0 = int(input("请输入初始浓度:"))
- xn = int(input("请输入最终浓度:"))
- F = int(input("请输入每小时蒸发质量kg:"))
- def motion(F,x,w0,Wn):
- # (1)计算各效蒸发量和完成液浓度
- w = w0
- x_array = np.array([x])
- w_array = np.array([w])
-
- W = F*(1 - x0 / xn) #W为总蒸发量
-
- #通常各效蒸发量可按各效蒸发量相等的原则设定 W1=W2=……=Wn
- #并流加料时有自蒸发现象,涉及到比值,此处只列举了常用的二、三效并流加料
- if n == 2: #双效 比例为1:1
- W1 = W/2
- W2 = W/2
- elif n ==3: #三效 三者之比为1:1.1:1.2
- W1 = W / 3.3
- W2 = 1.1 * W1
- W3 = 1.2 * W1
-
- else: #W1,W2,W3……Wn都相等
- for i in range (1,n+1):
- W[i] = W / n
- i += 1
复制代码 |
|